Browse
Abstract
Portions of the northeastern United States (NE) have experienced drought every year since 2016. The U.S. Drought Monitor (USDM) has played an important role in drought characterization and management by providing weekly drought maps across the entire United States, including the NE. Unfortunately, the USDM lacks consistency between input variables leading to difficulties in defining boundaries between drought categories. This paper evaluates the National Water Model’s (NWM) ability to model streamflow and soil moisture, two important hydrological products that are frequently incorporated in drought indices. Using a 26-yr NWM retrospective simulation, comparisons were conducted between NWM output and observations of streamflow and soil moisture, as well as between drought categories derived from the NWM and observations and the USDM. Results indicate that NWM provides moderate predictions of streamflow at NE stations when comparing to historical observations, that NWM streamflow estimators are generally upwardly biased, and performance is worse at lower streamflow magnitudes. The NWM’s ability to predict soil moisture is worse than streamflow, with again a positive bias at most sites and strong variations in anomaly correlation across sites. When predicting drought categories, NWM streamflow is as strong a predictor of USDM drought categories as observed streamflow. Extending the NWM streamflow series using a maintenance of variance technique and only past records provides slight improvements over drought categories derived from the entire 26-yr retrospective simulation. Output from the NWM appears to have some skill in characterizing drought in the NE and provides a spatial resolution to improve the designation of drought boundaries.
Abstract
Portions of the northeastern United States (NE) have experienced drought every year since 2016. The U.S. Drought Monitor (USDM) has played an important role in drought characterization and management by providing weekly drought maps across the entire United States, including the NE. Unfortunately, the USDM lacks consistency between input variables leading to difficulties in defining boundaries between drought categories. This paper evaluates the National Water Model’s (NWM) ability to model streamflow and soil moisture, two important hydrological products that are frequently incorporated in drought indices. Using a 26-yr NWM retrospective simulation, comparisons were conducted between NWM output and observations of streamflow and soil moisture, as well as between drought categories derived from the NWM and observations and the USDM. Results indicate that NWM provides moderate predictions of streamflow at NE stations when comparing to historical observations, that NWM streamflow estimators are generally upwardly biased, and performance is worse at lower streamflow magnitudes. The NWM’s ability to predict soil moisture is worse than streamflow, with again a positive bias at most sites and strong variations in anomaly correlation across sites. When predicting drought categories, NWM streamflow is as strong a predictor of USDM drought categories as observed streamflow. Extending the NWM streamflow series using a maintenance of variance technique and only past records provides slight improvements over drought categories derived from the entire 26-yr retrospective simulation. Output from the NWM appears to have some skill in characterizing drought in the NE and provides a spatial resolution to improve the designation of drought boundaries.
Abstract
Using abundant rainfall gauge measurements and Global Precipitation Mission (GPM) data, spatial patterns of rainfall diurnal cycles and their seasonality over high mountain Asia (HMA) were examined. Spatial distributions of rainfall diurnal cycles over the HMA have a prominent seasonality regulated by circulations at different spatiotemporal scales, within which large regional contrasts are embedded. Rainfall diurnal variability is relatively weak in the premonsoon season, with larger amplitude over the western HMA, the southeastern HMA, as well as southern periphery regions, characterized by a dominant late afternoon to morning rainfall preference. The pattern of rainfall spatial distributions is closely related to the midlatitude westerlies. Both the mean rainfall and amplitudes of diurnal cycles become more pronounced with the advance of monsoon season but weaken during postmonsoon. The widespread late afternoon to night pattern over HMA migrating with seasonal atmospheric circulation is consistent with the lifetime of convective systems, which become active from the afternoon due to radiative heating and decay during the night. Stationary terrain-dependent night-to-morning rainfall patterns are visible in those east–west-orientated valleys over HMA and the Qaidam basin throughout the seasons. This salient geographical dependence is associated with local circulation produced by the strong differential thermal conditions over mountains and valleys, which can lift the warm moist air at the mouth of the valley and trigger nocturnal convection.
Significance Statement
The main purpose of this study is to explore how spatial patterns of rainfall diurnal cycles over high mountain Asia vary with the seasons. Our results show that the widespread late afternoon to night rainfall over high mountain Asia migrating with seasonal atmospheric circulation is consistent with the lifetime of convective systems. Stationary terrain-dependent night-to-morning rainfall patterns are visible in those east–west-orientated valleys over high mountain Asia and the Qaidam basin throughout the seasons. These results highlight the importance of large-scale atmospheric circulation and local circulation on precipitation, which is critical for water resources over high mountain Asia.
Abstract
Using abundant rainfall gauge measurements and Global Precipitation Mission (GPM) data, spatial patterns of rainfall diurnal cycles and their seasonality over high mountain Asia (HMA) were examined. Spatial distributions of rainfall diurnal cycles over the HMA have a prominent seasonality regulated by circulations at different spatiotemporal scales, within which large regional contrasts are embedded. Rainfall diurnal variability is relatively weak in the premonsoon season, with larger amplitude over the western HMA, the southeastern HMA, as well as southern periphery regions, characterized by a dominant late afternoon to morning rainfall preference. The pattern of rainfall spatial distributions is closely related to the midlatitude westerlies. Both the mean rainfall and amplitudes of diurnal cycles become more pronounced with the advance of monsoon season but weaken during postmonsoon. The widespread late afternoon to night pattern over HMA migrating with seasonal atmospheric circulation is consistent with the lifetime of convective systems, which become active from the afternoon due to radiative heating and decay during the night. Stationary terrain-dependent night-to-morning rainfall patterns are visible in those east–west-orientated valleys over HMA and the Qaidam basin throughout the seasons. This salient geographical dependence is associated with local circulation produced by the strong differential thermal conditions over mountains and valleys, which can lift the warm moist air at the mouth of the valley and trigger nocturnal convection.
Significance Statement
The main purpose of this study is to explore how spatial patterns of rainfall diurnal cycles over high mountain Asia vary with the seasons. Our results show that the widespread late afternoon to night rainfall over high mountain Asia migrating with seasonal atmospheric circulation is consistent with the lifetime of convective systems. Stationary terrain-dependent night-to-morning rainfall patterns are visible in those east–west-orientated valleys over high mountain Asia and the Qaidam basin throughout the seasons. These results highlight the importance of large-scale atmospheric circulation and local circulation on precipitation, which is critical for water resources over high mountain Asia.
Abstract
Although several flow routing (FR) algorithms are developed for hydrological modeling, it is still uncertain how the selection of algorithms may affect model results. This study aims to explore the similarity and dissimilarity in model results among different FR algorithms characterized by single flow direction (SD) and multiple flow direction (MD). The Coupled Hydro-Ecological Simulation System (CHESS) was incorporated with six different FR algorithms (D8, D∞, MD∞, MD8, MFD-md, and RMD∞) and then applied for modeling ecohydrological processes for a semiarid mountainous watershed in the western United States during 1991–2012. Comparisons were made between the model results at the catchment and the grid scale. After slightly adjusting one of the most sensitive soil parameters, all algorithms behave similarly in simulating stream hydrographs. When averaged for the watershed, the modeled ecohydrological variables mostly do not differ significantly (<5%) among the six FR algorithms. Nevertheless, the simulated ecohydrological variables are spatially more autocorrelated under the more dispersive MD algorithms. In addition, there exist significant (>5%) cell-level differences in modeled soil moisture among different FR algorithms, with propagated influences on the simulated evapotranspiration and vegetation growth variables. In hillslopes, the cell-level differences in model results tend to increase significantly as the flows move to the streams. Overall, this study proves that the watershed-level differences in model results among FR algorithms are low after model calibration, while significant differences still occur at the cell level. Thus, observational data are essential for testing which routing algorithm captures better the reality of local ecohydrological processes.
Significance Statement
The consideration of flow routing is essential for accurately simulating land surface ecohydrological processes. However, less is known about how the selection of flow routing algorithms may affect the model results. Based on model experiments, we found that the model results under different algorithms do not significantly differ from each other when averaged for the watershed. However, significant differences in model results exist at the individual cell level. These findings are useful for guiding future modeling-related research and also suggest the importance of field studies for testing which routing algorithm can better represent local ecohydrological processes.
Abstract
Although several flow routing (FR) algorithms are developed for hydrological modeling, it is still uncertain how the selection of algorithms may affect model results. This study aims to explore the similarity and dissimilarity in model results among different FR algorithms characterized by single flow direction (SD) and multiple flow direction (MD). The Coupled Hydro-Ecological Simulation System (CHESS) was incorporated with six different FR algorithms (D8, D∞, MD∞, MD8, MFD-md, and RMD∞) and then applied for modeling ecohydrological processes for a semiarid mountainous watershed in the western United States during 1991–2012. Comparisons were made between the model results at the catchment and the grid scale. After slightly adjusting one of the most sensitive soil parameters, all algorithms behave similarly in simulating stream hydrographs. When averaged for the watershed, the modeled ecohydrological variables mostly do not differ significantly (<5%) among the six FR algorithms. Nevertheless, the simulated ecohydrological variables are spatially more autocorrelated under the more dispersive MD algorithms. In addition, there exist significant (>5%) cell-level differences in modeled soil moisture among different FR algorithms, with propagated influences on the simulated evapotranspiration and vegetation growth variables. In hillslopes, the cell-level differences in model results tend to increase significantly as the flows move to the streams. Overall, this study proves that the watershed-level differences in model results among FR algorithms are low after model calibration, while significant differences still occur at the cell level. Thus, observational data are essential for testing which routing algorithm captures better the reality of local ecohydrological processes.
Significance Statement
The consideration of flow routing is essential for accurately simulating land surface ecohydrological processes. However, less is known about how the selection of flow routing algorithms may affect the model results. Based on model experiments, we found that the model results under different algorithms do not significantly differ from each other when averaged for the watershed. However, significant differences in model results exist at the individual cell level. These findings are useful for guiding future modeling-related research and also suggest the importance of field studies for testing which routing algorithm can better represent local ecohydrological processes.
Abstract
Most heavy precipitation events and extreme flooding over the U.S. Pacific coast can be linked to prevalent atmospheric river (AR) conditions. Thus, reliable quantitative precipitation estimation with a rich spatiotemporal resolution is vital for water management and early warning systems of flooding and landslides over these regions. At the same time, high-quality near-real-time measurements of AR precipitation remain challenging due to the complex topographic features of land surface and meteorological conditions of the region: specifically, orographic features occlude radar measurements while infrared-based algorithms face challenges, differentiating between both cold brightband (BB) precipitation and the warmer nonbrightband (NBB) precipitation. It should be noted that the latter precipitation is characterized by greater orographic enhancement. In this study, we evaluate the performance of a recently developed near-real-time satellite precipitation algorithm: Precipitation Estimation from Remotely Sensed Information Using Artificial Neural Networks (PERSIANN) Dynamic Infrared–Rain Rate-Now (PDIR-Now). This model is primarily dependent on infrared information from geostationary satellites as input; consequently, PDIR-Now has the advantage of short data latency, 15–60-min delay between observation to precipitation product delivery. The performance of PDIR-Now is analyzed with a focus on AR-related events for cases dominated by NBB and BB precipitation over the Russian River basin. In our investigations, we utilize S-band (3-GHz) precipitation profilers with Joss/Parsivel disdrometer measurements at the Middletown and Santa Rosa stations to classify BB and NBB precipitation events. In general, our analysis shows that PDIR-Now is more skillful in retrieving precipitation rates over both BB and NBB events across the topologically complex study area as compared to PERSIANN-Cloud Classification System (CCS). Also, we discuss the performance of well-known operational near-real-time precipitation products from 2017 to 2019. Conventional categorical and volumetric categorical indices, as well as continuous statistical metrics, are used to show the differences between various high-resolution precipitation products such as Multi-Radar Multi-Sensor (MRMS).
Abstract
Most heavy precipitation events and extreme flooding over the U.S. Pacific coast can be linked to prevalent atmospheric river (AR) conditions. Thus, reliable quantitative precipitation estimation with a rich spatiotemporal resolution is vital for water management and early warning systems of flooding and landslides over these regions. At the same time, high-quality near-real-time measurements of AR precipitation remain challenging due to the complex topographic features of land surface and meteorological conditions of the region: specifically, orographic features occlude radar measurements while infrared-based algorithms face challenges, differentiating between both cold brightband (BB) precipitation and the warmer nonbrightband (NBB) precipitation. It should be noted that the latter precipitation is characterized by greater orographic enhancement. In this study, we evaluate the performance of a recently developed near-real-time satellite precipitation algorithm: Precipitation Estimation from Remotely Sensed Information Using Artificial Neural Networks (PERSIANN) Dynamic Infrared–Rain Rate-Now (PDIR-Now). This model is primarily dependent on infrared information from geostationary satellites as input; consequently, PDIR-Now has the advantage of short data latency, 15–60-min delay between observation to precipitation product delivery. The performance of PDIR-Now is analyzed with a focus on AR-related events for cases dominated by NBB and BB precipitation over the Russian River basin. In our investigations, we utilize S-band (3-GHz) precipitation profilers with Joss/Parsivel disdrometer measurements at the Middletown and Santa Rosa stations to classify BB and NBB precipitation events. In general, our analysis shows that PDIR-Now is more skillful in retrieving precipitation rates over both BB and NBB events across the topologically complex study area as compared to PERSIANN-Cloud Classification System (CCS). Also, we discuss the performance of well-known operational near-real-time precipitation products from 2017 to 2019. Conventional categorical and volumetric categorical indices, as well as continuous statistical metrics, are used to show the differences between various high-resolution precipitation products such as Multi-Radar Multi-Sensor (MRMS).
Abstract
The land surface model is extensively used to simulate turbulence fluxes and hydrological and momentum variables at the land–atmosphere interface. In this study, the Community Land Model, version 5 (CLM5), driven by the 0.1° × 0.1° Chinese Meteorological Forcing Dataset (CMFD) and the field-surveyed soil parameters, is used to simulate land surface processes during 1979–2018. Various high-quality land surface datasets are adopted to assess the model simulations. In general, the CLM5 well captures the monthly variations of 0–10-cm soil moisture in subregions, particularly in the Tibetan Plateau, with an anomaly correlation coefficient between 0.56 and 0.88. However, the simulated soil moisture shows overall wet biases in the whole country, resulting from several reasons. The model simulation is skillful in replicating both the magnitude and spatial pattern when they are compared with the MODIS snow cover dataset. Compared with in situ measured soil temperature in multiple soil layers within 320-cm soil depth from 1980 to 2018, the simulations accurately capture spatial patterns, vertical profiles, and long-term warming trends. For land surface energy components, the simulations have a highly temporal correlation with the observation of Chinese Flux Observation and Research Network (ChinaFLUX) cropland and grassland sites, except for four forest sites, where biases exist in both atmospheric forcing variables and surface vegetation phenology in the model default input dataset. In summary, this study reveals the overall capability of CLM5 in reproducing land surface energy fluxes and hydrological variables over conterminous China, and the validation results may also provide some references for future model improvement and application.
Significance Statement
The offline Community Land Model, version 5 (CLM5), driven by a 0.1° × 0.1° (∼10 km) horizontal resolution atmospheric forcing dataset and a set of field-surveyed soil parameters, are used to simulate the land surface hydrological and heat fluxes in continental China for 1980–2018. The simulated hydrological variables and energy fluxes are validated with various sources of high-quality observation-based datasets. From our systematic evaluations, the current CLM5 high–resolution simulation accurately captures the spatial patterns and temporal variations in most of the water and energy balance components, although biases exist in some simulated variables. Overall, this study reveals the capability of the offline CLM5 simulation in conterminous China and provides the reference for future model improvement and application.
Abstract
The land surface model is extensively used to simulate turbulence fluxes and hydrological and momentum variables at the land–atmosphere interface. In this study, the Community Land Model, version 5 (CLM5), driven by the 0.1° × 0.1° Chinese Meteorological Forcing Dataset (CMFD) and the field-surveyed soil parameters, is used to simulate land surface processes during 1979–2018. Various high-quality land surface datasets are adopted to assess the model simulations. In general, the CLM5 well captures the monthly variations of 0–10-cm soil moisture in subregions, particularly in the Tibetan Plateau, with an anomaly correlation coefficient between 0.56 and 0.88. However, the simulated soil moisture shows overall wet biases in the whole country, resulting from several reasons. The model simulation is skillful in replicating both the magnitude and spatial pattern when they are compared with the MODIS snow cover dataset. Compared with in situ measured soil temperature in multiple soil layers within 320-cm soil depth from 1980 to 2018, the simulations accurately capture spatial patterns, vertical profiles, and long-term warming trends. For land surface energy components, the simulations have a highly temporal correlation with the observation of Chinese Flux Observation and Research Network (ChinaFLUX) cropland and grassland sites, except for four forest sites, where biases exist in both atmospheric forcing variables and surface vegetation phenology in the model default input dataset. In summary, this study reveals the overall capability of CLM5 in reproducing land surface energy fluxes and hydrological variables over conterminous China, and the validation results may also provide some references for future model improvement and application.
Significance Statement
The offline Community Land Model, version 5 (CLM5), driven by a 0.1° × 0.1° (∼10 km) horizontal resolution atmospheric forcing dataset and a set of field-surveyed soil parameters, are used to simulate the land surface hydrological and heat fluxes in continental China for 1980–2018. The simulated hydrological variables and energy fluxes are validated with various sources of high-quality observation-based datasets. From our systematic evaluations, the current CLM5 high–resolution simulation accurately captures the spatial patterns and temporal variations in most of the water and energy balance components, although biases exist in some simulated variables. Overall, this study reveals the capability of the offline CLM5 simulation in conterminous China and provides the reference for future model improvement and application.
Abstract
Root-zone soil moisture (RZSM) is an important variable in land–atmosphere interactions, notably affecting the global climate system. Contrary to satellite-based acquisition of surface soil moisture, RZSM is generally obtained from model-based simulations. In this study, in situ observations from the Naqu and Pali networks that represent different climatic conditions over the Tibetan Plateau (TP) and a triple collocation (TC) method are used to evaluate model-based RZSM products, including Global Land Evaporation Amsterdam Model (GLEAM) (versions 3.5a and 3.5b), Global Land Data Assimilation System (GLDAS) (versions 2.1 and 2.2), and the fifth-generation European Centre for Medium-Range Weather Forecasts reanalysis (ERA5). The evaluation results based on in situ observations indicate that all products tend to overestimate but could generally capture the temporal variation, and ERA5 exhibits the best performance with the highest R (0.875) and the lowest unbiased RMSE (ubRMSE; 0.015 m3 m−3) against in situ observations in the Naqu network. In the TC analysis, similar results are obtained: ERA5 has the best performance with the highest TC-derived R (0.785) over the entire TP, followed by GLEAM v3.5a (0.746) and GLDAS-2.1 (0.682). Meanwhile, GLEAM v3.5a and GLDAS-2.1 outperform GLEAM v3.5b and GLDAS-2.2 over the entire TP, respectively. Besides, possible error causes in evaluating these RZSM products are summarized, and the effectiveness of TC method is also evaluated with two dense networks, finding that TC method is reliable since TC-derived R is close to ground-derived R, with only 6.85% mean relative differences. These results using both in situ observations and TC method may provide a new perspective for the soil moisture product developers to further enhance the accuracy of model-based RZSM over the TP.
Significance Statement
The purpose of this study is to better understand the quality and applicability of GLEAM, GLDAS, and ERA5 RZSM products over the TP using both in situ observations and the triple collocation (TC) method, making it better applied to climate and hydrological research. This study provides four standard statistical metrics evaluation based on in situ observations, as well as the reliable metric, that is, correlation coefficient (R) derived from TC method, and highlights that TC-based evaluation could supplement the ground-based validation, especially over the data-scarce TP region.
Abstract
Root-zone soil moisture (RZSM) is an important variable in land–atmosphere interactions, notably affecting the global climate system. Contrary to satellite-based acquisition of surface soil moisture, RZSM is generally obtained from model-based simulations. In this study, in situ observations from the Naqu and Pali networks that represent different climatic conditions over the Tibetan Plateau (TP) and a triple collocation (TC) method are used to evaluate model-based RZSM products, including Global Land Evaporation Amsterdam Model (GLEAM) (versions 3.5a and 3.5b), Global Land Data Assimilation System (GLDAS) (versions 2.1 and 2.2), and the fifth-generation European Centre for Medium-Range Weather Forecasts reanalysis (ERA5). The evaluation results based on in situ observations indicate that all products tend to overestimate but could generally capture the temporal variation, and ERA5 exhibits the best performance with the highest R (0.875) and the lowest unbiased RMSE (ubRMSE; 0.015 m3 m−3) against in situ observations in the Naqu network. In the TC analysis, similar results are obtained: ERA5 has the best performance with the highest TC-derived R (0.785) over the entire TP, followed by GLEAM v3.5a (0.746) and GLDAS-2.1 (0.682). Meanwhile, GLEAM v3.5a and GLDAS-2.1 outperform GLEAM v3.5b and GLDAS-2.2 over the entire TP, respectively. Besides, possible error causes in evaluating these RZSM products are summarized, and the effectiveness of TC method is also evaluated with two dense networks, finding that TC method is reliable since TC-derived R is close to ground-derived R, with only 6.85% mean relative differences. These results using both in situ observations and TC method may provide a new perspective for the soil moisture product developers to further enhance the accuracy of model-based RZSM over the TP.
Significance Statement
The purpose of this study is to better understand the quality and applicability of GLEAM, GLDAS, and ERA5 RZSM products over the TP using both in situ observations and the triple collocation (TC) method, making it better applied to climate and hydrological research. This study provides four standard statistical metrics evaluation based on in situ observations, as well as the reliable metric, that is, correlation coefficient (R) derived from TC method, and highlights that TC-based evaluation could supplement the ground-based validation, especially over the data-scarce TP region.
Abstract
Time series of snowfall observations from over 500 stations in Oregon (OR) and Washington (WA) were generated for subregions of these states. Data problems encountered were as follows: 1) monthly totals in printed reports prior to 1940 that were not in the digital archive, 2) archived data listed as “missing” that were available, 3) digitized reports after 2010 eliminated good data, and 4) “zero” totals incorrectly entered in the official archive rather than “missing,” especially after 1980. Though addressing these was done, there is reduced confidence that some regional time series are representative of true long-term trends, especially for regions with few systematically reporting stations. For most regions characterized by consistent monitoring and with the most robust statistical reproducibility, we find no statistically significant trends in their periods of record (up to 131 years) for November–April seasonal totals through April 2020. This result includes the main snowfall regions of the Cascade Range. However, snowfall in some lower-elevation areas of OR and WA appear to have experienced declining trends, consistent with an increase in northeastern Pacific Ocean temperatures. Finally, previously constructed time series through April 2011 for regions in California are updated through April 2020 to include the recent, exceptionally low seasonal totals on the western slopes of the Sierra Nevada. This update indicates 2014/15 was the record lowest, 2013/14 was the 5th lowest, and 2012/13 was the 14th lowest of 142 years. Even so, the 1879–2020 linear trend in this key watershed region, though −2.6% decade−1, was not significantly different from zero due to high interannual variability and reconstruction uncertainty.
Abstract
Time series of snowfall observations from over 500 stations in Oregon (OR) and Washington (WA) were generated for subregions of these states. Data problems encountered were as follows: 1) monthly totals in printed reports prior to 1940 that were not in the digital archive, 2) archived data listed as “missing” that were available, 3) digitized reports after 2010 eliminated good data, and 4) “zero” totals incorrectly entered in the official archive rather than “missing,” especially after 1980. Though addressing these was done, there is reduced confidence that some regional time series are representative of true long-term trends, especially for regions with few systematically reporting stations. For most regions characterized by consistent monitoring and with the most robust statistical reproducibility, we find no statistically significant trends in their periods of record (up to 131 years) for November–April seasonal totals through April 2020. This result includes the main snowfall regions of the Cascade Range. However, snowfall in some lower-elevation areas of OR and WA appear to have experienced declining trends, consistent with an increase in northeastern Pacific Ocean temperatures. Finally, previously constructed time series through April 2011 for regions in California are updated through April 2020 to include the recent, exceptionally low seasonal totals on the western slopes of the Sierra Nevada. This update indicates 2014/15 was the record lowest, 2013/14 was the 5th lowest, and 2012/13 was the 14th lowest of 142 years. Even so, the 1879–2020 linear trend in this key watershed region, though −2.6% decade−1, was not significantly different from zero due to high interannual variability and reconstruction uncertainty.
Abstract
Dominant hydrological processes of a catchment could shift due to a changing climate. This climate-induced hydrological nonstationarity could affect the reliability of future runoff projection developed using a hydrological model calibrated for the historical period as the model or parameters may no longer be suitable under a different future hydroclimate. This paper explores whether competing parameterization approaches proposed to account for hydrological nonstationarity could improve the robustness of future runoff projection compared to the traditional approach where the model is calibrated targeting overall model performance over the entire historical period. The modeling experiments are carried out using climate and streamflow datasets from southeastern Australia, which has experienced a long drought and exhibited noticeable hydrological nonstationarity. The results show that robust multicriteria calibration based on the Pareto front can provide a more consistent model performance over contrasting hydroclimate conditions, but at a slight expense of increased bias over the entire historical period compared to the traditional approach. However, the robust calibration does not necessarily result in a more reliable projection of future runoff. This is because the systematic bias in any parameterization approach would propagate from the historical period to the future period and would largely be cancelled out when estimating the relative runoff change. Ensemble simulations combining results from different parameterization considerations could produce a more inclusive range of future runoff projection as it covers the uncertainties due to model parameterization.
Abstract
Dominant hydrological processes of a catchment could shift due to a changing climate. This climate-induced hydrological nonstationarity could affect the reliability of future runoff projection developed using a hydrological model calibrated for the historical period as the model or parameters may no longer be suitable under a different future hydroclimate. This paper explores whether competing parameterization approaches proposed to account for hydrological nonstationarity could improve the robustness of future runoff projection compared to the traditional approach where the model is calibrated targeting overall model performance over the entire historical period. The modeling experiments are carried out using climate and streamflow datasets from southeastern Australia, which has experienced a long drought and exhibited noticeable hydrological nonstationarity. The results show that robust multicriteria calibration based on the Pareto front can provide a more consistent model performance over contrasting hydroclimate conditions, but at a slight expense of increased bias over the entire historical period compared to the traditional approach. However, the robust calibration does not necessarily result in a more reliable projection of future runoff. This is because the systematic bias in any parameterization approach would propagate from the historical period to the future period and would largely be cancelled out when estimating the relative runoff change. Ensemble simulations combining results from different parameterization considerations could produce a more inclusive range of future runoff projection as it covers the uncertainties due to model parameterization.
Abstract
Long-term hydrological projections can vary substantially depending on the combination of meteorological forcing dataset, hydrologic model (HM), emissions scenario, and natural climate variability. Identifying dominant sources of model spread in an ensemble of hydrologic projections is critically important for developing reliable hydrological projections in support of flooding risk assessment and water resources management; however, it is not well understood due to the multifactor and multiscale complexities involved in the long-term hydrological projections. Therefore, a stepwise clustered Bayesian (SCB) ensemble method will be first developed to improve the performance of long-term hydrological projections. Meanwhile, a mixed-level factorial inference (MLFI) approach is employed to estimate multiple uncertainties in hydrological projections over the Jing River basin (JRB). MLFI is able to reveal the main and interactive effects of the anthropogenic emission and model choices on the SCB ensemble projections. The results suggest that the daily maximum temperature under RCP8.5 in the 2050s and 2080s is expected to respectively increase by 3.2° and 5.2°C, which are much higher than the increases under RCP4.5. The maximum increase of the RegCM driven by CanESM2 (CARM)-projected changes in streamflow for the 2050s and 2080s under RCP4.5 is 0.30 and 0.59 × 103 m s−3 in November, respectively. In addition, in a multimodel GCM–RCM–HM ensemble, hydroclimate is found to be most sensitive to the choice of GCM. Moreover, it is revealed that the percentage of contribution of anthropogenic emissions to the changes in monthly precipitation is relatively smaller, but it makes a more significant contribution to the total variance of changes in potential evapotranspiration and streamflow.
Significance Statement
Increasing concerns have been paid to climate change due to its aggravating impacts on the hydrologic regime, leading to water-related disasters. Such impacts can be investigated through long-term hydrological projection under climate change. However, it is not well understood what factor plays a dominant role in inducing extensive uncertainties associated with the long-term hydrological projections due to plausible meteorological forcings, multiple hydrologic models, and internal variability. The stepwise cluster Bayesian ensemble method and mixed-level factorial inference approach are employed to quantify the contribution of multiple uncertainty sources. We find that the total variance of changes in monthly precipitation, potential evapotranspiration, and streamflow can be mainly explained by the model choices. The identified dominant factor accounting for projection uncertainties is critically important for developing reliable hydrological projections in support of flooding risk assessment and water resources management. It is suggested that more reliable models should be taken into consideration in order to improve the projection robustness from a perspective of the Loess Plateau.
Abstract
Long-term hydrological projections can vary substantially depending on the combination of meteorological forcing dataset, hydrologic model (HM), emissions scenario, and natural climate variability. Identifying dominant sources of model spread in an ensemble of hydrologic projections is critically important for developing reliable hydrological projections in support of flooding risk assessment and water resources management; however, it is not well understood due to the multifactor and multiscale complexities involved in the long-term hydrological projections. Therefore, a stepwise clustered Bayesian (SCB) ensemble method will be first developed to improve the performance of long-term hydrological projections. Meanwhile, a mixed-level factorial inference (MLFI) approach is employed to estimate multiple uncertainties in hydrological projections over the Jing River basin (JRB). MLFI is able to reveal the main and interactive effects of the anthropogenic emission and model choices on the SCB ensemble projections. The results suggest that the daily maximum temperature under RCP8.5 in the 2050s and 2080s is expected to respectively increase by 3.2° and 5.2°C, which are much higher than the increases under RCP4.5. The maximum increase of the RegCM driven by CanESM2 (CARM)-projected changes in streamflow for the 2050s and 2080s under RCP4.5 is 0.30 and 0.59 × 103 m s−3 in November, respectively. In addition, in a multimodel GCM–RCM–HM ensemble, hydroclimate is found to be most sensitive to the choice of GCM. Moreover, it is revealed that the percentage of contribution of anthropogenic emissions to the changes in monthly precipitation is relatively smaller, but it makes a more significant contribution to the total variance of changes in potential evapotranspiration and streamflow.
Significance Statement
Increasing concerns have been paid to climate change due to its aggravating impacts on the hydrologic regime, leading to water-related disasters. Such impacts can be investigated through long-term hydrological projection under climate change. However, it is not well understood what factor plays a dominant role in inducing extensive uncertainties associated with the long-term hydrological projections due to plausible meteorological forcings, multiple hydrologic models, and internal variability. The stepwise cluster Bayesian ensemble method and mixed-level factorial inference approach are employed to quantify the contribution of multiple uncertainty sources. We find that the total variance of changes in monthly precipitation, potential evapotranspiration, and streamflow can be mainly explained by the model choices. The identified dominant factor accounting for projection uncertainties is critically important for developing reliable hydrological projections in support of flooding risk assessment and water resources management. It is suggested that more reliable models should be taken into consideration in order to improve the projection robustness from a perspective of the Loess Plateau.
Abstract
Fog plays a vital role in maintaining ecosystems in montane cloud forests. In these forests, a large amount of water on the surface of leaves and canopy (hereafter canopy water) evaporates during the morning. This biophysical process plays a critical factor in regulating afternoon fog formation. Recent studies have found that alterations in precipitation, temperature, humidity, and CO2 concentrations associated with future climate changes may affect terrestrial hydroclimatology, but the responses in cloud forests remain unclear. Utilizing numerical experiments with the Community Land Model, we explored changes in surface evaporative fluxes in Chi-Lan Mountain cloud forests in northeastern Taiwan under the RCP8.5 scenario with changes in the aforementioned various atmospheric variables. The results showed that increased rainfall intensity in climate change runs decreased the accumulation of canopy water, while larger water vapor concentrations led to more nighttime condensation on leaves. Elevated CO2 concentrations did not greatly impact canopy water amounts, but photosynthesis was enhanced, while transpiration was reduced and contributed to decreased latent heat fluxes, implying the importance of forest plant physiology in modulating land evaporative fluxes. Evapotranspiration decreased in Chi-Lan due to multiple combined factors, in contrast to the expected intensification in the global water cycle under global warming. The study, however, is restricted to an offline land surface model without land–atmosphere interactions and the interactions with adjacent grids, which deserves further analyses for the water cycle changes in the montane cloud forest regions.
Abstract
Fog plays a vital role in maintaining ecosystems in montane cloud forests. In these forests, a large amount of water on the surface of leaves and canopy (hereafter canopy water) evaporates during the morning. This biophysical process plays a critical factor in regulating afternoon fog formation. Recent studies have found that alterations in precipitation, temperature, humidity, and CO2 concentrations associated with future climate changes may affect terrestrial hydroclimatology, but the responses in cloud forests remain unclear. Utilizing numerical experiments with the Community Land Model, we explored changes in surface evaporative fluxes in Chi-Lan Mountain cloud forests in northeastern Taiwan under the RCP8.5 scenario with changes in the aforementioned various atmospheric variables. The results showed that increased rainfall intensity in climate change runs decreased the accumulation of canopy water, while larger water vapor concentrations led to more nighttime condensation on leaves. Elevated CO2 concentrations did not greatly impact canopy water amounts, but photosynthesis was enhanced, while transpiration was reduced and contributed to decreased latent heat fluxes, implying the importance of forest plant physiology in modulating land evaporative fluxes. Evapotranspiration decreased in Chi-Lan due to multiple combined factors, in contrast to the expected intensification in the global water cycle under global warming. The study, however, is restricted to an offline land surface model without land–atmosphere interactions and the interactions with adjacent grids, which deserves further analyses for the water cycle changes in the montane cloud forest regions.
