Search

Northern Mexico is home to more than 32 million people and is of significant agricultural and economic importance for the country. The region includes three distinct hydroclimatic regions, all of which regularly experience severe dryness and flooding and are highly susceptible to future changes in precipitation. To date, little work has been done to characterize future trends in either mean or extreme precipitation over northern Mexico. To fill this gap, we investigate projected precipitation trends over the region in the NA-CORDEX ensemble of dynamically downscaled simulations. We first verify that these simulations accurately reproduce observed precipitation over northern Mexico, as derived from the Multi-Source Weighted-Ensemble Precipitation (MSWEP) product, demonstrating that the NA-CORDEX ensemble is appropriate for studying precipitation trends over the region. By the end of the century, simulations forced with a high-emissions scenario project that both mean and extreme precipitation will decrease to the west and increase to the east of the Sierra Madre highlands, decreasing the zonal gradient in precipitation. We also find that the North American monsoon, which is responsible for a substantial fraction of the precipitation over the region, is likely to start later and last approximately three weeks longer. The frequency of extreme precipitation events is expected to double throughout the region, exacerbating the flood risk for vulnerable communities in northern Mexico. Collectively, these results suggest that the extreme precipitation-related dangers that the region faces, such as flooding, will increase significantly by the end of the century, with implications for the agricultural sector, economy, and infrastructure. Significance Statement: Northern Mexico regularly experiences severe flooding and its important agricultural sector can be heavily impacted by variations in precipitation. Using high-resolution climate model simulations that have been tested against observations, we find that these hydroclimate extremes are likely to be exacerbated in a warming climate; the dry (wet) season is projected to receive significantly less (more) precipitation (approximately ±10% by the end of the century). Simulations suggest that some of the changes in precipitation over the region can be related to the North American monsoon, with the monsoon starting later in the year and lasting several weeks longer. Our results also suggest that the frequency of extreme precipitation will increase, although this increase is smaller than that projected for other regions, with the strongest storms becoming 20% more frequent per degree of warming. These results suggest that this region may experience significant changes to its hydroclimate through the end of the century that will require significant resilience planning. [ABSTRACT FROM AUTHOR]

Macroscale hydrological/land surface models are important tools for assessing historical and predicting future characteristics of extreme hydrological events, yet quantitative understandings of how these large-scale models perform in simulating extreme hydrological characteristics remain limited. Here we evaluate simulated high and low flows from 23 macroscale models within three modeling experiments (i.e., 14 climate models from CMIP6, 6 global hydrological models from ISIMIP2a, and 3 land surface models from GLDAS) against observation in 633 unimpaired catchments globally over 1971–2010. Our findings reveal limitations in simulating extreme flow characteristics by these models. Specifically, we find that (i) most models overestimate high-flow magnitudes (bias range: from +15% to +70%) and underestimate low-flow magnitudes (bias range: from −80% to −20%); (ii) interannual variability in high and low flows is reasonably reproduced by ISIMIP2a and GLDAS models but poorly reproduced by CMIP6 models; (iii) no model consistently replicates the observed trend direction in high and low flows in over two-thirds of the catchments, and most models overestimate high-flow trends and underestimate low-flow trends; and (iv) CMIP6 and GLDAS models show timing biases, with early high flows and late low flows, while ISIMIP2a models exhibit the opposite pattern. Furthermore, all models performed better in more humid environments and noncold regions, with model structure and parameterization contributing more to uncertainties than climatic forcings. Overall, our results demonstrate that extreme flow characteristics simulated from current state-of-the-art macroscale models still contain large uncertainties and provide important guidance regarding the robustness of assessing extreme hydrometeorological events based on these modeling outputs. Significance Statement: Macroscale hydrological and land surface models represent crucial tools for assessing historical trends and making predictions about future hydrological changes. Nevertheless, our current understanding of the quantitative performance of these large-scale models in simulating extreme hydrological characteristics remains limited. Here, we evaluate simulated high and low flows from 23 state-of-the-art macroscale models against observation in 633 unimpaired catchments globally over 1971–2010. Our results reveal important limitations in the extreme flow characteristics simulated from these models and provide important guidance regarding the robustness of assessing extreme hydrometeorological events based on these modeling outputs. The model evaluation performed herein serves as a pivotal, offering valuable insights to inform the development of the next generation of macroscale hydrological and land surface models. [ABSTRACT FROM AUTHOR]

Statistical processing of numerical model output has been a part of both weather forecasting and climate applications for decades. Statistical techniques are used to correct systematic biases in atmospheric model outputs and to represent local effects that are unresolved by the model, referred to as downscaling. Many downscaling techniques have been developed, and it has been difficult to systematically explore the implications of the individual decisions made in the development of downscaling methods. Here we describe a unified framework that enables the user to evaluate multiple decisions made in the methods used to statistically postprocess output from weather and climate models. The Ensemble Generalized Analog Regression Downscaling (En-GARD) method enables the user to select any number of input variables, predictors, mathematical transformations, and combinations for use in parametric or nonparametric downscaling approaches. En-GARD enables explicitly predicting both the probability of event occurrence and the event magnitude. Outputs from En-GARD include errors in model fit, enabling the production of an ensemble of projections through sampling of the probability distributions of each climate variable. We apply En-GARD to regional climate model simulations to evaluate the relative importance of different downscaling method choices on simulations of the current and future climate. We show that choice of predictor variables is the most important decision affecting downscaled future climate outputs, while having little impact on the fidelity of downscaled outcomes for current climate. We also show that weak statistical relationships prevent such approaches from predicting large changes in extreme events on a daily time scale. [ABSTRACT FROM AUTHOR]

Flood-related risks to people and property are expected to increase in the future due to environmental and demographic changes. It is important to quantify and effectively communicate flood hazards and exposure to inform the design and implementation of flood risk management strategies. Here we develop an integrated modeling framework to assess projected changes in regional riverine flood inundation risks. The framework samples climate model outputs to force a hydrologic model and generate streamflow projections. Together with a statistical and hydraulic model, we use the projected streamflow to map the uncertainty of flood inundation projections for extreme flood events. We implement the framework for rivers across the state of Pennsylvania, United States. Our projections suggest that flood hazards and exposure across Pennsylvania are overall increasing with future climate change. Specific regions, including the main stem Susquehanna River, lower portion of the Allegheny basin, and central portion of Delaware River basin, demonstrate higher flood inundation risks. In our analysis, the climate uncertainty dominates the overall uncertainty surrounding the flood inundation projection chain. The combined hydrologic and hydraulic uncertainties can account for as much as 37% of the total uncertainty. We discuss how this framework can provide regional and dynamic flood-risk assessments and help to inform the design of risk-management strategies. [ABSTRACT FROM AUTHOR]

Traditional multimodel methods for estimating future changes in precipitation intensity, duration, and frequency (IDF) curves rely on mean or median of models' IDF estimates. Such multimodel estimates are impaired by large estimation uncertainty, shadowing their efficacy in planning efforts. Here, assuming that each climate model is one representation of the underlying data generating process, i.e., the Earth system, we propose a novel extension of current methods through pooling model data: (i) evaluate performance of climate models in simulating the spatial and temporal variability of the observed annual maximum precipitation (AMP), (ii) bias-correct and pool historical and future AMP data of reasonably performing models, and (iii) compute IDF estimates in a nonstationary framework from pooled historical and future model data. Pooling enhances fitting of the extreme value distribution to the data and assumes that data from reasonably performing models represent samples from the "true" underlying data generating distribution. Through Monte Carlo simulations with synthetic data, we show that return periods derived from pooled data have smaller biases and lesser uncertainty than those derived from ensembles of individual model data. We apply this method to NA-CORDEX models to estimate changes in 24-h precipitation intensity–frequency (PIF) estimates over the Susquehanna watershed and Florida peninsula. Our approach identifies significant future changes at more stations compared to median-based PIF estimates. The analysis suggests that almost all stations over the Susquehanna and at least two-thirds of the stations over the Florida peninsula will observe significant increases in 24-h precipitation for 2–100-yr return periods. [ABSTRACT FROM AUTHOR]

Surface temperature is one of the key parameters for estimating regional evapotranspiration (ET) based on the Surface Energy Balance System (SEBS) model using remote sensing data. However, continuous daily remote sensing surface temperature data are often not available due to the weather and environmental conditions. This paper proposed a scheme to obtain reliable ET that estimating ET using WRF-simulated surface skin temperature (TSK) and then modifying the deviation using the normalized difference vegetation index (NDVI). This study aims to explore whether the model data can be a viable option when the continuous-time-series remote sensing surface temperature is missing for estimating ET. Comparison results show that the correlation between WRF TSK and the measured temperature of the 2-cm soil (Ts) is better than MODIS land surface temperature (LST) in the study area, while the correlation between MODIS LST and the measured surface radiation temperature (IRT) is better than WRF TSK. The MODIS LST is significantly higher than Ts, and the WRF TSK is closer to Ts than MODIS LST. However, the ET calculated using WRF TSK was not good, exhibiting relatively high ET in the whole area and a poor correlation with the measurements, whereby R2, RMSE, and the percent bias (PBIAS) were equal to 0.1256, 5.2783 mm, and −202.17%, respectively. According to the principle of land surface process simulation in WRF, this paper proposes using NDVI to modify ET calculated using TSK. The comparison between the modified ET and the measurements exhibited a relatively good correlation, with R2 = 0.7532, RMSE = 1.0993 mm, and PBIAS = −17.9%. Therefore, the model surface temperature data can be used to estimate continuous-time-series regional ET when NDVI is used to modify the deviation, which indicates the surface temperature data simulated by the WRF Model can become the optional data for estimating ET and compensate for the shortcoming of poor time continuity of remote sensing data, further expanding the application prospects of meteorological model data in the remote sensing field. [ABSTRACT FROM AUTHOR]

How the globally uniform component of sea surface temperature (SST) warming influences rainfall in the African Sahel remains insufficiently studied, despite mean SST warming being among the most robustly simulated and theoretically grounded features of anthropogenic climate change. A prior study using the NOAA Geophysical Fluid Dynamics Laboratory (GFDL) AM2.1 atmospheric general circulation model (AGCM) demonstrated that uniform SST warming strengthens the prevailing northerly advection of dry Saharan air into the Sahel. The present study uses uniform SST warming simulations performed with 7 GFDL and 10 CMIP5 AGCMs to assess the robustness of this drying mechanism across models and uses observations to assess the physical credibility of the severe drying response in AM2.1. In all 17 AGCMs, mean SST warming enhances the free-tropospheric meridional moisture gradient spanning the Sahel and with it the Saharan dry-air advection. Energetically, this is partially balanced by anomalous subsidence, yielding decreased precipitation in 14 of the 17 models. Anomalous subsidence and precipitation are tightly linked across the GFDL models but not the CMIP5 models, precluding the use of this relationship as the start of a causal chain ending in an emergent observational constraint. For AM2.1, cloud–rainfall covariances generate radiative feedbacks on drying through the subsidence mechanism and through surface hydrology that are excessive compared to observations at the interannual time scale. These feedbacks also act in the equilibrium response to uniform warming, calling into question the Sahel's severe drying response to warming in all coupled models using AM2.1. [ABSTRACT FROM AUTHOR]

Global climate models play an important role in quantifying past and projecting future changes in drought. Previous studies have pointed to shortcomings in these models for simulating droughts, but systematic evaluation of their level of agreement has been limited. Here, historical simulations (1950-2004) for 20 models from the latest Coupled Model Intercomparison Project (CMIP5) were analyzed for a variety of drought metrics and thresholds using a standardized drought index. Model agreement was investigated for different types of drought (precipitation, runoff, and soil moisture) and how this varied with drought severity and duration. At the global scale, climate models were shown to agree well on most precipitation drought metrics, but systematically underestimated precipitation drought intensity compared to observations. Conversely, simulated runoff and soil moisture droughts varied significantly across models, particularly for intensity. Differences in precipitation simulations were found to explain model differences in runoff and soil moisture drought metrics over some regions, but predominantly with respect to drought intensity. This suggests it is insufficient to evaluate models for precipitation droughts to increase confidence in model performance for other types of drought. This study shows large but metric-dependent discrepancies in CMIP5 for modeling different types of droughts that relate strongly to the component models (i.e., atmospheric or land surface scheme) used in the coupled modeling systems. Our results point to a need to consider multiple models in drought impact studies to account for high model uncertainties. [ABSTRACT FROM AUTHOR]

Climate models project significant twenty-first-century declines in water availability over the American West from anthropogenic warming. However, the physical mechanisms underpinning this response are poorly characterized, as are the uncertainties from vegetation's modulation of evaporative losses. To understand the drivers and uncertainties of future hydroclimate in the American West, a 35-member single model ensemble is used to examine the response of summer soil moisture and runoff to anthropogenic forcing. Widespread dry season soil moisture declines occur across the region despite increases in total water-year precipitation and ubiquitous increases in plant water-use efficiency. These modeled soil moisture declines are initially forced by significant snowpack losses that directly diminish summer soil water, even in regions where water-year precipitation increases. When snowpack priming is coupled with a warming- and CO2-induced shift in phenology and increased primary production, widespread increases in leaf area further reduces summer soil moisture and runoff by outpacing decreased stomatal conductance from high CO2. The net effects lead to the co-occurrence of both a 'greener' and 'drier' future across the western United States. Because simulated vegetation exerts a large influence on predicted changes in water availability in the American West, these findings highlight the importance of reducing the substantial uncertainties in the ecological processes increasingly incorporated into numerical Earth system models. [ABSTRACT FROM AUTHOR]

The article discusses the Weather Research and Forecasting (WRF) Model which is an atmospheric model designed for research and numerical weather prediction supported by the National Center for Atmospheric Research (NCAR). It states that WRF provides specialty capabilities for a range of Earth system prediction applications including wildland fires, hydrology and hurricanes.

Most of the currently employed procedures for bias correction and statistical downscaling primarily consider a univariate approach by developing a statistical relationship between large-scale precipitation/temperature with the local-scale precipitation/temperature, ignoring the interdependency between the two variables. In this study, a multivariate approach, asynchronous canonical correlation analysis (ACCA), is proposed and applied to global climate model (GCM) historic simulations and hindcasts from phase 5 of the Coupled Model Intercomparison Project (CMIP5) to downscale monthly precipitation and temperature over the conterminous United States. ACCA is first applied to the CNRM-CM5 GCM historical simulations for the period 1950-99 and compared with the bias-corrected dataset based on quantile mapping from the Bureau of Reclamation. ACCA is also applied to CNRM-CM5 hindcasts and compared with univariate asynchronous regression (ASR), which applies regular regression to sorted GCM and observed variables. ACCA performs better than ASR and quantile mapping in preserving the cross correlation at grid points where the observed cross correlations are significant while reducing fractional biases in mean and standard deviation. Results also show that preservation of cross correlation increases the bias in standard deviation slightly, but estimates observed precipitation and temperature with increased likelihood, particularly for months exhibiting significant cross correlation. ACCA also better estimates the joint likelihood of observed precipitation and temperature under hindcasts since hindcasts estimate the observed variability in precipitation better. Implications of preserving cross correlations across climate variables for projecting runoff and other land surface fluxes are also discussed. [ABSTRACT FROM AUTHOR]

Global climate models (GCMs) are usually used for future climate projections. Model output from GCMs needs to be downscaled and stochastic weather generators such as Climate Generator (CLIGEN) are tools to downscale GCM output and to produce synthetic weather sequences that are statistically similar to the observed weather data. Two methods of adjusting CLIGEN parameters were developed to reproduce precipitation sequences for southeastern Australia (SEA), where significant changes in annual precipitation had occurred, and for southwestern Western Australia (SWWA), where the precipitation has shown a significant decreasing trend since the 1920s. The adjustment methods have been validated using observed precipitation data for these regions. However, CLIGEN outputs ultimately will be used as input to other simulation models. The objective of this research was to further validate the methods of CLIGEN parameter adjustment using conceptual hydrological models to simulate streamflow and to compare the streamflow using observed and CLIGEN-generated precipitation data. Six precipitation sites from SEA and SWWA were selected and synthetic time series of daily precipitation were generated for these sites. Conceptual hydrological models, namely, the Australian Water Balance Model and SimHyd, were used for flow simulation and were calibrated using recorded daily streamflow data from six gauging stations in SEA and SWWA. Both monthly and annual streamflow show statistically similar patterns using observed and CLIGEN-generated precipitation data. The adjustment methods for CLIGEN parameters are further validated and can be used to reproduce the significant changes, both abrupt and gradually decreasing, in streamflow for these two climatically contrasting regions of Australia. [ABSTRACT FROM AUTHOR]

Gridded climate datasets are produced in many parts of the world by applying various interpolation methods to weather observations, to which are sometimes added secondary information (in addition to geographic location) such as topography and radar or atmospheric model outputs. For a region of interest, the choice of a dataset for a given study can be a significant challenge given the lack of information on the similarities and differences that exist between datasets, or about the benefits that one dataset may present relative to another. This study aims to provide information on the spatial and temporal differences between gridded precipitation datasets and their implication for hydrological modeling. Three gridded datasets for the province of Quebec are considered: the Natural Resources Canada (NRCan) dataset, the Canadian Precipitation Analysis (CaPA) dataset, and the dataset from the Ministère du Développement Durable, de I'Environnement et de la Lutte contre les Changements Climatiques du Québec (MDDELCC). Using statistical metrics and diagrams, these precipitation datasets are compared with each other. Hydrological responses of 181 Quebec watersheds with respect to each gridded precipitation dataset are also analyzed using the hydrological model HSAMI. The results indicate strong similarities in the southern parts and disparities in the central and northern parts of the province of Quebec. Analysis of hydrological simulations indicates that the CaPA dataset offers the best results, particularly for watersheds located in the central and northern parts of the province. MDDELCC shows the best performance in watersheds located on the south shore of the St. Lawrence River and comes out as the overall second-best option. [ABSTRACT FROM AUTHOR]

A simple conceptual model of surface specific humidity change over land is described, based on the effect of increased moisture advection from the oceans in response to sea surface temperature (SST) warming. In this model, future q over land is determined by scaling the present-day pattern of land q by the fractional increase in the oceanic moisture source. Simple model estimates agree well with climate model projections of future (mean spatial correlation coefficient 0.87), so over both land and ocean can be viewed primarily as a thermodynamic process controlled by SST warming. Precipitation change is also affected by , and the new simple model can be included in a decomposition of tropical precipitation change, where it provides increased physical understanding of the processes that drive over land. Confidence in the thermodynamic part of extreme precipitation change over land is increased by this improved understanding, and this should scale approximately with Clausius-Clapeyron oceanic q increases under SST warming. Residuals of actual climate model from simple model estimates are often associated with regions of large circulation change, and can be thought of as the 'dynamical' part of specific humidity change. The simple model is used to explore intermodel uncertainty in , and there are substantial contributions to uncertainty from both the thermodynamic (simple model) and dynamical (residual) terms. The largest cause of intermodel uncertainty within the thermodynamic term is uncertainty in the magnitude of global mean SST warming. [ABSTRACT FROM AUTHOR]

In the last decade operational probabilistic ensemble flood forecasts have become common in supporting decision-making processes leading to risk reduction. Ensemble forecasts can assess uncertainty, but they are limited to the uncertainty in a specific modeling system. Many of the current operational flood prediction systems use a multimodel approach to better represent the uncertainty arising from insufficient model structure. This study presents a multimodel approach to building a global flood prediction system using multiple atmospheric reanalysis datasets for river initial conditions and multiple TIGGE forcing inputs to the ECMWF land surface model. A sensitivity study is carried out to clarify the effect of using archive ensemble meteorological predictions and uncoupled land surface models. The probabilistic discharge forecasts derived from the different atmospheric models are compared with those from the multimodel combination. The potential for further improving forecast skill by bias correction and Bayesian model averaging is examined. The results show that the impact of the different TIGGE input variables in the HTESSEL/Catchment-Based Macroscale Floodplain model (CaMa-Flood) setup is rather limited other than for precipitation. This provides a sufficient basis for evaluation of the multimodel discharge predictions. The results also highlight that the three applied reanalysis datasets have different error characteristics that allow for large potential gains with a multimodel combination. It is shown that large improvements to the forecast performance for all models can be achieved through appropriate statistical postprocessing (bias and spread correction). A simple multimodel combination generally improves the forecasts, while a more advanced combination using Bayesian model averaging provides further benefits. [ABSTRACT FROM AUTHOR]

Univariate bias correction algorithms, such as quantile mapping, are used to address systematic biases in climate model output. Intervariable dependence structure (e.g., between different quantities like temperature and precipitation or between sites) is typically ignored, which can have an impact on subsequent calculations that depend on multiple climate variables. A novel multivariate bias correction (MBC) algorithm is introduced as a multidimensional analog of univariate quantile mapping. Two variants are presented. MBCp and MBCr respectively correct Pearson correlation and Spearman rank correlation dependence structure, with marginal distributions in both constrained to match observed distributions via quantile mapping. MBC is demonstrated on two case studies: 1) bivariate bias correction of monthly temperature and precipitation output from a large ensemble of climate models and 2) multivariate correction of vertical humidity and wind profiles, including subsequent calculation of vertically integrated water vapor transport and detection of atmospheric rivers. The energy distance is recommended as an omnibus measure of performance for model selection. As expected, substantial improvements in performance relative to quantile mapping are found in each case. For reference, characteristics of the MBC algorithm are compared against existing bivariate and multivariate bias correction techniques. MBC performs competitively and fills a role as a flexible, general purpose multivariate bias correction algorithm. [ABSTRACT FROM AUTHOR]

Future changes in floods and water availability across China under representative concentration pathway 2.6 (RCP2.6) and RCP8.5 are studied by analyzing discharge simulations from the Inter-Sectoral Impact Model Intercomparison Project (ISI-MIP) with the consideration of uncertainties among global climate models (GCMs) and hydrologic models. Floods and water availability derived from ISI-MIP simulations are compared against observations. The uncertainties among models are quantified by model agreement. Only model agreement >50% is considered to generate reliable projections of floods and water availability and their relationships with climate change. The results show five major points. First, ISI-MIP simulations have acceptable ability in modeling floods and water availability. The spatial patterns of changes in floods and water availability highly depend on the outputs of GCMs. Uncertainties from GCMs/hydrologic models predominate the uncertainties in the wet/dry areas in eastern/northwestern China. Second, the magnitudes of floods throughout China increase during 2070-99 under RCP8.5 relative to those with the same return periods during 1971-2000. The increase rates of larger floods are higher than those of the smaller ones. Third, water availability decreases/increases in southern/northern China under RCP8.5, but changes negligibly under RCP2.6. Fourth, more severe floods in the future are driven by more intense precipitation extremes over China. The negligible change in mean precipitation and the increase in actual evapotranspiration reduce the water availability in southern China. Fifth, model agreements are higher in simulated floods than water availability because increasing precipitation extremes are more consistent among different GCM outputs compared to mean precipitation. [ABSTRACT FROM AUTHOR]

Projections of regional climate, net basin supply (NBS), and water levels are developed for the mid- and late twenty-first century across the Laurentian Great Lakes basin. Two state-of-the-art global climate models (GCMs) are dynamically downscaled using a regional climate model (RCM) interactively coupled to a one-dimensional lake model, and then a hydrologic routing model is forced with time series of perturbed NBS. The dynamical downscaling and coupling with a lake model to represent the Great Lakes create added value beyond the parent GCM in terms of simulated seasonal cycles of temperature, precipitation, and surface fluxes. However, limitations related to this rudimentary treatment of the Great Lakes result in warm summer biases in lake temperatures, excessive ice cover, and an abnormally early peak in lake evaporation. While the downscaling of both GCMs led to consistent projections of increases in annual air temperature, precipitation, and all NBS components (overlake precipitation, basinwide runoff, and lake evaporation), the resulting projected water level trends are opposite in sign. Clearly, it is not sufficient to correctly simulate the signs of the projected change in each NBS component; one must also account for their relative magnitudes. The potential risk of more frequent episodes of lake levels below the low water datum, a critical shipping threshold, is explored. [ABSTRACT FROM AUTHOR]

Projections of regional climate, net basin supply (NBS), and water levels are developed for the mid- and late twenty-first century across the Laurentian Great Lakes basin. Two state-of-the-art global climate models (GCMs) are dynamically downscaled using a regional climate model (RCM) interactively coupled to a one-dimensional lake model, and then a hydrologic routing model is forced with time series of perturbed NBS. The dynamical downscaling and coupling with a lake model to represent the Great Lakes create added value beyond the parent GCM in terms of simulated seasonal cycles of temperature, precipitation, and surface fluxes. However, limitations related to this rudimentary treatment of the Great Lakes result in warm summer biases in lake temperatures, excessive ice cover, and an abnormally early peak in lake evaporation. While the downscaling of both GCMs led to consistent projections of increases in annual air temperature, precipitation, and all NBS components (overlake precipitation, basinwide runoff, and lake evaporation), the resulting projected water level trends are opposite in sign. Clearly, it is not sufficient to correctly simulate the signs of the projected change in each NBS component; one must also account for their relative magnitudes. The potential risk of more frequent episodes of lake levels below the low water datum, a critical shipping threshold, is explored. [ABSTRACT FROM AUTHOR]

The bias correction of climate model outputs is based on the main assumption of the time invariance of the bias, in which the statistical relationship between observations and climate model outputs in the historical period stays constant in the future period. The present study aims to assess statistical bias correction under nonstationary bias conditions and its implications on the simulated streamflow over two snowmelt-driven Canadian catchments. A pseudoreality approach is employed in order to derive a proxy of future observations. In this approach, CRCM-ECHAM5 ensemble simulations are used as pseudoreality observations to correct for bias in the CRCM-CGCM3 ensemble simulations in the reference (1971-2000) period. The climate model simulations are then bias corrected in the future (2041-70) period and compared with the future pseudoreality observations. This process demonstrates that biases (precipitation and temperature) remain after the bias correction. In a second step, the uncorrected and bias-corrected CRCM-CGCM3 simulations are used to drive the Soil and Water Assessment Tool (SWAT) hydrological model in both periods. The bias correction decreases the error on mean monthly streamflow over the reference period; such findings are more mixed over the future period. The results of various hydrological indicators show that the climate change signal on streamflow obtained with uncorrected and bias-corrected simulations is similar in both its magnitude and its direction for the mean monthly streamflow only. Regarding the indicators of extreme hydrological events, more mixed results are found with site dependence. All in all, bias correction under nonstationary bias is an additional source of uncertainty that cannot be neglected in hydrological climate change impact studies. [ABSTRACT FROM AUTHOR]

The aridity of a terrestrial climate is often quantified using the dimensionless ratio of annual precipitation ( P) to annual potential evapotranspiration (PET). In this study, the climatological patterns and greenhouse warming responses of terrestrial P, Penman-Monteith PET, and are compared among 16 modern global climate models. The large-scale climatological values and implied biome types often disagree widely among models, with large systematic differences from observational estimates. In addition, the PET climatologies often differ by several tens of percent when computed using monthly versus 3-hourly inputs. With greenhouse warming, land P does not systematically increase or decrease, except at high latitudes. Therefore, because of moderate, ubiquitous PET increases, decreases (drying) are much more widespread than increases (wetting) in the tropics, subtropics, and midlatitudes in most models, confirming and expanding on earlier findings. The PET increases are also somewhat sensitive to the time resolution of the inputs, although not as systematically as for the PET climatologies. The changes in the balance between P and PET are also quantified using an alternative aridity index, the ratio , which has a one-to-one but nonlinear correspondence with . It is argued that the magnitudes of changes are more uniformly relevant than the magnitudes of changes, which tend to be much higher in wetter regions. The ratio and its changes are also found to be excellent statistical predictors of the land surface evaporative fraction and its changes. [ABSTRACT FROM AUTHOR]

This study compared two types of approaches to downscale seasonal precipitation ( P) and 2-m air temperature (T2M) forecasts from the North American Multimodel Ensemble (NMME) over the states of Alabama, Georgia, and Florida in the southeastern United States (SEUS). Each NMME model forecast was evaluated. Two multimodel ensemble (MME) schemes were tested by assigning equal weight to all forecast members (SuperEns) or by assigning equal weights to each model's ensemble mean (MeanEns). One type of downscaling approach used was a model output statistics (MOS) method, which was based on direct spatial disaggregation and bias correction of the NMME P and T2M forecasts using the quantile mapping technique [spatial disaggregation with bias correction (SDBC)]. The other type of approach used was a perfect prognosis (PP) approach using nonparametric locally weighted polynomial regression (LWPR) models, which used the NMME forecasts of Niño-3.4 sea surface temperatures (SSTs) to predict local-scale P and T2M. Both SDBC and LWPR downscaled P showed skill in winter but no skill or limited skill in summer at all lead times for all NMME models. The SDBC downscaled T2M were skillful only for the Climate Forecast System, version 2 (CFSv2), model even at far lead times, whereas the LWPR downscaled T2M showed limited skill or no skill for all NMME models. In many cases, the LWPR method showed significantly higher skill than the SDBC. After bias correction, the SuperEns mostly showed higher skill than the MeanEns and most of the single models, but its skill did not outperform the best single model. [ABSTRACT FROM AUTHOR]

Extreme value theory is used as a diagnostic for two high-resolution (12-km parameterized convection and 1.5-km explicit convection) Met Office regional climate model (RCM) simulations. On subdaily time scales, the 12-km simulation has weaker June-August (JJA) short-return-period return levels than the 1.5-km RCM, yet the 12-km RCM has overly large high return levels. Comparisons with observations indicate that the 1.5-km RCM is more successful than the 12-km RCM in representing (multi)hourly JJA very extreme events. As accumulation periods increase toward daily time scales, the erroneous 12-km precipitation extremes become more comparable with the observations and the 1.5-km RCM. The 12-km RCM fails to capture the observed low sensitivity of the growth rate to accumulation period changes, which is successfully captured by the 1.5-km RCM. Both simulations have comparable December-February (DJF) extremes, but the DJF extremes are generally weaker than in JJA at daily or shorter time scales. Case studies indicate that 'gridpoint storms' are one of the causes of unrealistic very extreme events in the 12-km RCM. Caution is needed in interpreting the realism of 12-km RCM JJA extremes, including short-return-period events, which have return values closer to observations. There is clear evidence that the 1.5-km RCM has a higher degree of realism than the 12-km RCM in the simulation of JJA extremes. [ABSTRACT FROM AUTHOR]

Seasonal predictions of precipitation and surface air temperature from the Climate Forecast System, version 2 (CFSv2), are evaluated against gridded daily observations from 1982 to 2007 over 17 hydroclimatic regions in China. The seasonal predictive skill is quantified with skill scores including correlation coefficient, RMSE, and mean bias for spatially averaged seasonal precipitation and temperature forecasts for each region. The evaluation focuses on identifying regions and seasons where significant skill exists, thus potentially contributing to skill in hydrological prediction. The authors find that the predictive skill of CFSv2 precipitation and temperature forecasts has a stronger dependence on seasons and regions than on lead times. Both temperature and precipitation forecasts show higher skill from late summer [July-September (JAS)] to late autumn [October-December (OND)] and from winter [December-February (DJF)] to spring [March-May (MAM)]. The skill of CFSv2 precipitation forecasts is low during summer [June-August (JJA)] and winter (DJF) over all of China because of low potential predictability of the East Asian summer monsoon and the East Asian winter monsoon for China. As expected, temperature predictive skill is much higher than precipitation predictive skill in all regions. As observed precipitation shows significant correlation with the Oceanic Niño index over western, southwestern, and central China, the authors found that CFSv2 precipitation forecasts generally show similar correlation pattern, suggesting that CFSv2 precipitation forecasts can capture ENSO signals. This evaluation suggests that using CFSv2 forecasts for seasonal hydrological prediction over China is promising and challenging. [ABSTRACT FROM AUTHOR]

In this study, the authors contrast four century-long meteorological datasets comprising of two sets of observations [Climate Research Unit (CRU) and Parameter–Elevation Regressions on Independent Slopes Model (PRISM)] and two atmospheric reanalyses [Twentieth Century Reanalysis (20CR) and Florida Climate Institute–Florida State University Land–Atmosphere Regional Reanalysis version 1.0 (FLAReS1.0)] to diagnose the El Niño–Southern Oscillation (ENSO) forced variations on the streamflow in 28 watersheds spread across the southeastern United States (SEUS). The datasets are used to force three different lumped (calibrated) hydrological models with precipitation from these four sources of century-long datasets separately to obtain the median prediction from 1800 (=3 models × 600 simulations per model per watershed per season) multimodel estimates of seasonal mean streamflow across the 28 watersheds in the SEUS for each winter season from 1906 to 2005. The authors then compare and contrast themean streamflow and its variability estimates fromall three of the century-long climate forcings. The multimodel strategy of simulating the seasonal mean streamflow is to reduce the hydrological model uncertainty. The authors focus on the boreal winter season when ENSO influence on the SEUS climate variations is well known. The authors find that the atmospheric reanalysis over the SEUS is able to reasonably capture the ENSO teleconnections as depicted in the CRU and PRISM precipitation datasets. Even the observed decadal modulation of this teleconnection by Atlantic multidecadal oscillation (AMO) is broadly captured. The streamflow in the 28 watersheds also show similar consistency across the four datasets in that the positive correlations of the boreal winter Niño-3.4 SST anomalies with corresponding anomalies of streamflow, the associated shift in the probability density function of the streamflow with the change in phase of ENSO, and the decadal modulation of the ENSO teleconnection by the AMO are sustained in the streamflow simulations forced by all four climate datasets (CRU, PRISM, 20CR, and FLAReS1.0). However, the ENSO signal in the streamflow is consistently much stronger in the southern watersheds (over Florida) of the SEUS across all four climate datasets. During the negative phase of the AMO, however, there is a clear shift of the ENSO teleconnections with streamflow, with winter streamflows in northern watersheds (over the Carolinas) exhibiting much stronger correlations with the ENSO Niño-3.4 index relative to the southern watersheds of the SEUS. This study clearly indicates that the proposed methodology using FLAReS1.0 serves as a viable alternative to reconstruct twentieth-century SEUS seasonal winter hydrology that captures the interannual variations of ENSO and associated decadal variations forced by the AMO. However, it is found that the FLAReS1.0 forced streamflow is far from adequate in simulating the streamflow dynamics of the watershed over the SEUS at a daily time scale. [ABSTRACT FROM AUTHOR]

A high-resolution climate model (4-km horizontal grid spacing) is used to examine the following question: How will long-term changes in climate impact the partitioning of annual precipitation between evapotranspiration and runoff in the Colorado Headwaters? This question is examined using a climate sensitivity approach in which eight years of current climate is compared to a future climate created by modifying the current climate signal with perturbation from the NCAR Community Climate System Model, version 3 (CCSM3), model forced by the A1B scenario for greenhouse gases out to 2050. The current climate period is shown to agree well with Snowpack Telemetry (SNOTEL) surface observations of precipitation ( P) and snowpack, as well as streamflow and AmeriFlux evapotranspiration (ET) observations. The results show that the annual evaporative fraction (ET/ P) for the Colorado Headwaters is 0.81 for the current climate and 0.83 for the future climate, indicating increasing aridity in the future despite a positive increase of precipitation. Runoff decreased by an average of 6%, reflecting the increased aridity. Precipitation increased in the future winter by 12%, but decreased in the summer as a result of increased low-level inhibition to convection. The fraction of precipitation that fell as snow decreased from 0.83 in the current climate to 0.74 in the future. Future snowpack did not change significantly until January. From January to March the snowpack increased above ~3000 m MSL and decreased below that level. Snowpack decreased at all elevations in the future from April to July. The peak snowpack and runoff over the headwaters occurred 2-3 weeks earlier in the future simulation, in agreement with previous studies. [ABSTRACT FROM AUTHOR]

A flash flood and landslide in the Leh region of the Indus Valley in the Indian state of Jammu and Kashmir on 5-6 August 2010 resulted in hundreds of deaths and great property damage. Observations have led to the hypothesis that the storm, which formed over the Tibetan Plateau, was steered over the steep edge of the plateau by 500-hPa winds and then energized by the ingestion of lower-level moist air, which was approaching from the Arabian Sea and Bay of Bengal and rose up the Himalayan barrier. A coupled land surface and atmospheric model simulation validates this hypothesized storm scenario, with the model storm taking the form of a traveling mesoscale squall line with a leading convective line, trailing stratiform region, and midlevel inflow jet. In this region, the development of a mesoscale storm over high terrain is highly unusual, especially one in the form of a propagating squall line system. This unusual storm occurrence and behavior could serve as a warning sign in flash flood prediction. The coupled atmosphere and land surface model showed that the excessive runoff leading to the flood and landslide were favored by the occurrence of this unusual meteorological event coinciding temporally and spatially with favorable hydrologic conditions. Additionally, the model simulations showed that previous rainstorms had moistened the soil during the entire season and especially over the few days leading up to the Leh flood, so the normally arid mountainsides were likely not able to rapidly absorb the additional rainfall of the sudden 5 August squall line. [ABSTRACT FROM AUTHOR]

In many hydrological studies, the main limiting factor in model performance is the low meteorological data quality. In some cases, no meteorological records even exist. Installing weather stations becomes a necessity in these areas when water resource management becomes an issue. The objective of this study is to propose a new experimental and exploratory method for determining the optimal density of a weather station network when being used for long-term hydrological modeling. Data from the Canadian Regional Climate Model at 15-km resolution (CRCM15) were used to create a virtual network of stations with long and complete series of meteorological data over the Toulnustouc River basin in central Québec, Canada. The weather stations to be fed to HSAMI, Hydro-Québec's lumped rainfall-runoff hydrological model, were selected in order to minimize the number of stations while maintaining the best hydrological performance possible using a multi-objective optimization algorithm. It was shown that the number of stations making up the network on the Toulnustouc River basin should be at least two but not higher than four. If the stations are positioned optimally, there is little to no gain to be made with a denser network. The optimization algorithm clearly identified that combinations of two or three stations can result in better hydrological performance than if a high-density network was fed to the model. Thus, the major conclusion of this study is that if weather stations are positioned at optimal locations, a very few number of them are required to model runoff with as good as or better performance than when a high-density network is used. [ABSTRACT FROM AUTHOR]

Projections of storm-track changes over the continental United States and southern Canada made by 23 models from phase 5 of the Coupled Model Intercomparison Project (CMIP5) have been compared to changes projected by 11 models from phase 3 of CMIP (CMIP3). Overall, under representative concentration pathway 8.5 (RCP8.5) forcing, CMIP5 models project much more significant decreases in North American storm-track activity than CMIP3 models under the Special Report on Emission Scenarios (SRES) A2 scenario, with the largest decrease in summer and the smallest decrease in spring. The decrease is found both in temporal variance and cyclone statistics, with the frequency of strong cyclones projected to decrease by 15.9%, 6.6%, 32.6%, and 16.9% for winter, spring, summer, and fall, respectively. There is a strong consensus among the 23 models regarding the sign of the projected change, with less than 20% of the models projecting changes in the opposite sign in any of the storm-track parameters examined. Nevertheless, there are also significant model-to-model differences in the magnitude of the projected changes. Projected changes in mean flow baroclinicity have also been examined. Model-to-model differences in the projected storm-track change are found to correlate significantly with model-to-model differences in the projected change in a locally defined mean available potential energy (MAPE) across the ensemble of 34 CMIP5 and CMIP3 models, suggesting that the differences in the projected change in local MAPE can partly account for not only the model-to-model differences but also the differences between CMIP5 and CMIP3 projections. Examination of projected precipitation change suggests that models projecting larger decrease in North American storm-track activity also project a farther northward intrusion of the decrease in subtropical precipitation. [ABSTRACT FROM AUTHOR]

In this study, the authors applied version 4 of the Community Land Model (CLM4) integrated with an uncertainty quantification (UQ) framework to 20 selected watersheds from the Model Parameter Estimation Experiment (MOPEX) spanning a wide range of climate and site conditions to investigate the sensitivity of runoff simulations to major hydrologic parameters and to assess the fidelity of CLM4, as the land component of the Community Earth System Model (CESM), in capturing realistic hydrological responses. They found that for runoff simulations, the most significant parameters are those related to the subsurface runoff parameterizations. Soil texture-related parameters and surface runoff parameters are of secondary significance. Moreover, climate and soil conditions play important roles in the parameter sensitivity. In general, water-limited hydrologic regime and finer soil texture result in stronger sensitivity of output variables, such as runoff and its surface and subsurface components, to the input parameters in CLM4. This study evaluated the parameter identifiability of hydrological parameters from streamflow observations at selected MOPEX basins and demonstrated the feasibility of parameter inversion/calibration for CLM4 to improve runoff simulations. The results suggest that in order to calibrate CLM4 hydrologic parameters, model reduction is needed to include only the identifiable parameters in the unknowns. With the reduced parameter set dimensionality, the inverse problem is less ill posed. [ABSTRACT FROM AUTHOR]

The impact of climate change on three small- to medium-sized river catchments (Ammer, Mulde, and Ruhr) in Germany is investigated for the near future (2021-50) following the Intergovernmental Panel on Climate Change (IPCC) Special Report on Emissions Scenarios (SRES) A1B scenario. A 10-member en-semble of hydrological model (HM) simulations, based on two high-resolution regional climate models (RCMs) driven by two global climate models (GCMs), with three realizations of ECHAM5 (E5) and one realization of the Canadian Centre for Climate Modelling and Analysis version 3 (CCCma3; C3) is estab-lished. All GCM simulations are downscaled by the RCM Community Land Model (CLM), and one re-alization of E5 is downscaled also with the RCM Weather Research and Forecasting Model (WRF). This concerted 7-km, high-resolution RCM ensemble provides a sound basis for runoff simulations of small catchments and is currently unique for Germany. The hydrology for each catchment is simulated in an overlapping scheme, with two of the three HMs used in the project. The resulting ensemble hence contains for each chain link (GCM-realization-RCM-HM) at least two members and allows the investigation of quali-tative and limited quantitative indications of the existence and uncertainty range of the change signal. The ensemble spread in the climate change signal is large and varies with catchment and season, and the results show that most of the uncertainty of the change signal arises from the natural variability in winter and from the RCMs in summer. [ABSTRACT FROM AUTHOR]

Permafrost is a characteristic aspect of the terrestrial Arctic and the fate of near-surface permafrost over the next century is likely to exert strong controls on Arctic hydrology and biogeochemistry. Using output from the fifth phase of the Coupled Model Intercomparison Project (CMIP5), the authors assess its ability to simulate present-day and future permafrost. Permafrost extent diagnosed directly from each climate model's soil temperature is a function of the modeled surface climate as well as the ability of the land surface model to represent permafrost physics. For each CMIP5 model these two effects are separated by using indirect estimators of permafrost driven by climatic indices and compared to permafrost extent directly diagnosed via soil temperatures. Several robust conclusions can be drawn from this analysis. Significant air temperature and snow depth biases exist in some model's climates, which degrade both directly and indirectly diagnosed permafrost conditions. The range of directly calculated present-day (1986-2005) permafrost area is extremely large (~4-25 × 106 km2). Several land models contain structural weaknesses that limit their skill in simulating cold region subsurface processes. The sensitivity of future permafrost extent to temperature change over the present-day observed permafrost region averages (1.67 ± 0.7) × 106 km2 °C−1 but is a function of the spatial and temporal distribution of climate change. Because of sizable differences in future climates for the representative concentration pathway (RCP) emission scenarios, a wide variety of future permafrost states is predicted by 2100. Conservatively, the models suggest that for RCP4.5, permafrost will retreat from the present-day discontinuous zone. Under RCP8.5, sustainable permafrost will be most probable only in the Canadian Archipelago, Russian Arctic coast, and east Siberian uplands. [ABSTRACT FROM AUTHOR]

Irrigation is an important human activity that may impact local and regional climate, but current climate model simulations and data assimilation systems generally do not explicitly include it. The European Centre for Medium-Range Weather Forecasts (ECMWF) Interim Re-Analysis (ERA-Interim) shows more irrigation signal in surface evapotranspiration (ET) than the Modern-Era Retrospective Analysis for Research and Applications (MERRA) because ERA-Interim adjusts soil moisture according to the observed surface temperature and humidity while MERRA has no explicit consideration of irrigation at the surface. But, when compared with the results from a hydrological model with detailed considerations of agriculture, the ET from both reanalyses show large deficiencies in capturing the impact of irrigation. Here, a back-trajectory method is used to estimate the contribution of irrigation to precipitation over local and surrounding regions, using MERRA with observation-based corrections and added irrigation-caused ET increase from the hydrological model. Results show substantial contributions of irrigation to precipitation over heavily irrigated regions in Asia, but the precipitation increase is much less than the ET increase over most areas, indicating that irrigation could lead to water deficits over these regions. For the same increase in ET, precipitation increases are larger over wetter areas where convection is more easily triggered, but the percentage increase in precipitation is similar for different areas. There are substantial regional differences in the patterns of irrigation impact, but, for all the studied regions, the highest percentage contribution to precipitation is over local land. [ABSTRACT FROM AUTHOR]

Precipitation simulated by regional climate models (RCMs) is generally biased with respect to observations, especially at the local scale of a few tens of kilometers. This study investigates how well two different RCMs are able to reproduce the spatial correlation patterns of observed summer precipitation for the central United States. On local scales, gridded precipitation observations and simulated precipitation are compared for the period of the 1987 First International Satellite Land Surface Climatological Project Field Experiment (FIFE) campaign. The results show that spatial correlation length scales on the order of 130 km are found in both observed data and RCM simulations. When simulations and observations are aggregated to different grid sizes, the pattern correlation significantly decreases when the aggregation length is less than roughly 100 km. Furthermore, the intermodel standard deviation between simulations with different domains or resolutions increases for aggregation lengths below ~130 km. Below this length scale there is a high level of randomness in the precise location of precipitation events. Conversely, spatial correlation values increase above this length scale, reflecting larger predictive certainty of the RCMs at larger scales. The findings on aggregated grid scales are shown to be largely independent of the underlying RCMs grid resolutions but not of the overall size of RCM domain. With regard to hydrological modeling applications, these findings indicate that precipitation extracted from the present RCM simulations at a catchment scale below the intermodel standard deviation length cannot be expected to accurately match observations. [ABSTRACT FROM AUTHOR]

Recent advances in global flow routing schemes have shown the importance of using high-resolution topography for representing floodplain inundation dynamics more reliably. This study presents and evaluates the Hydrological Modeling and Analysis Platform (HyMAP), which is a global flow routing scheme specifically designed to bridge the gap between current state-of-the-art global flow routing schemes by combining their main features and introducing new features to better capture floodplain dynamics. The ultimate goals of HyMAP are to provide the scientific community with a novel scheme suited to the assimilation of satellite altimetry data for global water discharge forecasts and a model that can be potentially coupled with atmospheric models. In this first model evaluation, HyMAP is coupled with the Interactions between Soil-Biosphere-Atmosphere (ISBA) land surface model in order to simulate the surface water dynamics in the Amazon basin. The model is evaluated over the 1986-2006 period against an unprecedented source of information, including in situ and satellite-based datasets of water discharge and level, flow velocity, and floodplain extent. Results show that the model can satisfactorily simulate the large-scale features of the water surface dynamics of the Amazon River basin. Among all stream gauges considered, 23% have Nash-Sutcliffe coefficients (NS) higher than 0.50 and 68% above zero. About 28% of the stations have volume errors lower than 15%. Simulated discharges at Óbidos had NS = 0.89. Time series of simulated floodplains at the basin scale agrees well with satellite-based estimates, with a relative error of 7% and correlation of 0.89. These results indicate nonnegligible improvements in comparison to previous studies for the same region. [ABSTRACT FROM AUTHOR]

The Colorado River is the primary water source for much of the rapidly growing southwestern United States. Recent studies have projected reductions in Colorado River flows from less than 10% to almost 50% by midcentury because of climate change-a range that has clouded potential management responses. These differences in projections are attributable to variations in climate model projections but also to differing land surface model (LSM) sensitivities. This second contribution to uncertainty-specifically, variations in LSM runoff change with respect to precipitation (elasticities) and temperature (sensitivities)-are evaluated here through comparisons of multidecadal simulations from five commonly used LSMs (Catchment, Community Land Model, Noah, Sacramento Soil Moisture Accounting model, and Variable Infiltration Capacity model) all applied over the Colorado River basin at 1/8° latitude by longitude spatial resolution. The annual elasticity of modeled runoff (fractional change in annual runoff divided by fractional change in annual precipitation) at Lees Ferry ranges from two to six for the different LSMs. Elasticities generally are higher in lower precipitation and/or runoff regimes; hence, the highest values are for models biased low in runoff production, and the range of elasticities is reduced to two to three when adjusted to current runoff climatology. Annual temperature sensitivities (percent change in annual runoff per degree change in annual temperature) range from declines of 2% to as much as 9% per degree Celsius increase at Lees Ferry. For some LSMs, small areas, primarily at midelevation, have increasing runoff with increasing temperature; however, on a spatial basis, most sensitivities are negative. [ABSTRACT FROM AUTHOR]

This study assesses the performance of the regional climate model (RCM) simulations from the North American Regional Climate Change Assessment Program (NARCCAP) for the Upper Colorado River basin (UCRB), U.S. Rocky Mountains. The UCRB is a major contributor to the Colorado River's runoff. Its significant snow-dominated hydrological regime makes it highly sensitive to climatic changes, and future water shortage in this region is likely. The RCMs are evaluated with a clear RCM output user's perspective and a main focus on snow. Snow water equivalent (SWE) and snow duration, as well as air temperature and precipitation from five RCMs, are compared with snowpack telemetry (SNOTEL) observations, with National Centers for Environmental Prediction (NCEP)-National Center for Atmospheric Research (NCAR) Reanalysis II (R2), which provides the boundary conditions for the RCM simulations, and with North American Regional Reanalysis (NARR). Overall, most RCMs were able to significantly improve on the results from the NCEP-NCAR reanalysis. However, in comparison with spatially aggregated point observations and NARR, the RCMs are generally too dry, too warm, simulate too little SWE, and have a too-short snow cover duration with a too-late start and a too-early end of a significant snow cover. The intermodel biases found are partly associated with inadequately resolved topography (at the spatial resolution of the RCMs), imperfect observational data, different forcing techniques (spectral nudging versus no nudging), and the different land surface schemes (LSS). Attributing the found biases to specific features of the RCMs remains difficult or even impossible without detailed knowledge of the physical and technical specification of the models. [ABSTRACT FROM AUTHOR]

Three different complexity snow schemes implemented in the ECMWF land surface scheme Hydrology Tiled ECMWF Scheme of Surface Exchanges over Land (HTESSEL) are evaluated within the EC-EARTH climate model. The snow schemes are (i) the original HTESSEL single-bulk-layer snow scheme, (ii) a new snow scheme in operations at ECMWF since September 2009, and (iii) a multilayer version of the previous. In offline site simulations, the multilayer scheme outperforms the single-layer schemes in deep snowpack conditions through its ability to simulate sporadic melting events thanks to the lower thermal inertial of the uppermost layer. Coupled atmosphere-land/snow simulations performed by the EC-EARTH climate model are validated against remote sensed snow cover and surface albedo. The original snow scheme has a systematic early melting linked to an underestimation of surface albedo during spring that was partially reduced with the new snow schemes. A key process to improve the realism of the near-surface atmospheric temperature and at the same time the soil freezing is the thermal insulation of the snowpack (tightly coupled with the accuracy of snow mass and density simulations). The multilayer snow scheme outperforms the single-layer schemes in open deep snowpack (such as prairies or tundra in northern latitudes) and is instead comparable in shallow snowpack conditions. However, the representation of orography in current climate models implies limitations for accurately simulating the snowpack, particularly over complex terrain regions such as the Rockies and the Himalayas. [ABSTRACT FROM AUTHOR]

Anthropogenic activities have been significantly perturbing global freshwater flows and groundwater reserves. Despite numerous advances in the development of land surface models (LSMs) and global terrestrial hydrological models (GHMs), relatively few studies have attempted to simulate the impacts of anthropogenic activities on the terrestrial water cycle using the framework of LSMs. From the comparison of simulated terrestrial water storage with the Gravity Recovery and Climate Experiment (GRACE) satellite observations it is found that a process-based LSM, the Minimal Advanced Treatments of Surface Interaction and Runoff (MATSIRO), outperforms the bucket-model-based GHM called H08 in simulating hydrologic variables, particularly in water-limited regions. Therefore, the water regulation modules of H08 are incorporated into MATSIRO. Further, a new irrigation scheme based on the soil moisture deficit is developed. Incorporation of anthropogenic water regulation modules significantly improves river discharge simulation in the heavily regulated global river basins. Simulated irrigation water withdrawal for the year 2000 (2462 km3 yr−1) agrees well with the estimates provided by the Food and Agriculture Organization (FAO). Results indicate that irrigation changes surface energy balance, causing a maximum increase of ~50 W m−2 in latent heat flux averaged over June-August. Moreover, unsustainable anthropogenic water use in 2000 is estimated to be ~450 km3 yr−1, which corresponds well with documented records of groundwater overdraft, representing an encouraging improvement over the previous modeling studies. Globally, unsustainable water use accounts for ~40%% of blue water used for irrigation. The representation of anthropogenic activities in MATSIRO makes the model a suitable tool for assessing potential anthropogenic impacts on global water resources and hydrology. [ABSTRACT FROM AUTHOR]

The representation of hydrological processes in land surface schemes (LSSs) has recently been improved. In this study, the usability of GCM runoff for river discharge modeling is evaluated by validating the mean, timing, and amplitude of the modeled annual discharge cycles against observations. River discharge was calculated for six large rivers using runoff, precipitation, and actual evaporation from the GCMs ECHAM5 and Hadley Centre Global Environmental Model version 2 (HadGEM2). Four methods were applied: 1) accumulation of GCM runoff, 2) routing of GCM runoff, 3) routing of GCM runoff combined with temporal storage of subsurface runoff, and 4) offline hydrological modeling with the global distributed hydrological model PCRaster Global Water Balance (PCR-GLOBWB) using meteorological data from the GCMs as forcing. The quality of discharge generated by all four methods is highly influenced by the quality of the GCM data. In small catchments, the methods that include runoff routing perform equally well, although offline modeling with PRC-GLOBWB outperforms the other methods for ECHAM5 data. For larger catchments, routing introduces realistic travel times, decreased day-to-day variability, and it reduces extremes. Complexity of the LSS of both GCMs is comparable to the complexity of the hydrological model. However, in HadGEM2 the absence of subgrid variability for saturated hydraulic conductivity results in a large subsurface runoff flux and a low seasonal variability in the annual discharge cycle. The analysis of these two GCMs shows that when LSSs are tuned to reproduce realistic water partitioning at the grid scale and a routing scheme is also included, discharge variability and change derived from GCM runoff could be as useful as changes derived from runoff obtained from offline simulations using large-scale hydrological models. [ABSTRACT FROM AUTHOR]

Satellite and gridded meteorological data can be used to estimate evaporation ( E) from land surfaces using simple diagnostic models. Two satellite datasets indicate a positive trend (first time derivative) in global available energy from 1983 to 2006, suggesting that positive trends in evaporation may occur in 'wet' regions where energy supply limits evaporation. However, decadal trends in evaporation estimated from water balances of 110 wet catchments do not match trends in evaporation estimated using three alternative methods: 1) , a model-tree ensemble approach that uses statistical relationships between E measured across the global network of flux stations, meteorological drivers, and remotely sensed fraction of absorbed photosynthetically active radiation; 2) , a Budyko-style hydrometeorological model; and 3) , the Penman-Monteith energy-balance equation coupled with a simple biophysical model for surface conductance. Key model inputs for the estimation of and are remotely sensed radiation and gridded meteorological fields and it is concluded that these data are, as yet, not sufficiently accurate to explain trends in E for wet regions. This provides a significant challenge for satellite-based energy-balance methods. Trends in for 87 'dry' catchments are strongly correlated to trends in precipitation ( R2 == 0.85). These trends were best captured by , which explicitly includes precipitation and available energy as model inputs. [ABSTRACT FROM AUTHOR]

The ability of two versions of the Noah land surface model (LSM) to simulate the water cycle of the Little Washita River experimental watershed is evaluated. One version that uses the standard hydrological parameterizations of Noah 2.7 (STD) is compared another version that replaces STD''s subsurface hydrology with a simple aquifer model and topography-related surface and subsurface runoff parameterizations (GW). Simulations on a distributed grid at fine resolution are compared to the long-term distribution of observed daily-mean runoff, the spatial statistics of observed soil moisture, and locally observed latent heat flux. The evaluation targets the typical behavior of ensembles of models that use realistic, near-optimal sets of parameters important to runoff. STD and GW overestimate the ratio of runoff to evapotranspiration. In the subset of STD and GW runs that best reproduce the timing and the volume of streamflow, the surface-to-subsurface runoff ratio is overestimated and simulated streamflow is much flashier than observations. Both models'' soil columns wet and dry too quickly, implying that there are structural shortcomings in the formulation of STD that cannot be overcome by adding GW''s increased complexity to the model. In its current formulation, GW extremely underestimates baseflow''s contribution to total runoff and requires a shallow water table to function realistically. In the catchment (depth to water table >10 m), GW functions as a simple bucket model. Because model parameters are likely scale and site dependent, the need for even ''physically based'' models to be extensively calibrated for all domains on which they are applied is underscored. [ABSTRACT FROM AUTHOR]

In earth system models, the partitioning of precipitation among the variations of continental water storage, evapotranspiration, and freshwater runoff to the ocean has a major influence on the terrestrial water and energy budgets and thereby on simulated climate on a wide range of scales. The evaluation of continental hydrology is therefore a crucial task that requires offline simulations driven by realistic atmospheric forcing to avoid the systematic biases commonly found in global atmospheric models. Generally, this evaluation is done mainly by comparison with in situ river discharge data, which does not guarantee that the spatiotemporal distribution of water storage and evapotranspiration is correctly simulated. In this context, the Interactions between Soil, Biosphere, and Atmosphere–Total Runoff Integrating Pathways (ISBA-TRIP) continental hydrological system of the Centre National de Recherches Météorologiques is evaluated by using the additional constraint of terrestrial water storage (TWS) variations derived from three independent gravity field retrievals (datasets) from the Gravity Recovery and Climate Experiment (GRACE). On the one hand, the results show that, in general, ISBA-TRIP captures the seasonal and the interannual variability in both TWS and discharges. GRACE provides an additional constraint on the simulated hydrology and consolidates the former evaluation only based on river discharge observations. On the other hand, results indicate that river storage variations represent a significant contribution to GRACE measurements. While this remark highlights the need to improve the TRIP river routing model for a more useful comparison with GRACE [Decharme et al. (Part II of the present study)], it also suggests that low-resolution gravimetry products do not necessarily represent a strong additional constraint for model evaluation, especially in downstream areas of large river basins where long-term discharge data are available. [ABSTRACT FROM AUTHOR]

The present study assesses the potential of the U.S. Climate Variability and Predictability (CLIVAR) Drought Working Group (DWG) models in simulating interannual precipitation variability over North America, especially the Great Plains. It also provides targets for the idealized DWG model experiments investigating drought origin. The century-long Atmospheric Model Intercomparison Project (AMIP) simulations produced by version 3.5 of NCAR’s Community Atmosphere Model (CAM3.5), the Lamont-Doherty Earth Observatory’s Community Climate Model (CCM3), and NASA’s Seasonal-to-Interannual Prediction Project (NSIPP-1) atmospheric models are analyzed; CCM3 and NSIPP-1 models have 16- and 14-ensemble simulations, respectively, while CAM3.5 only has 1. The standard deviation of summer precipitation is different in AMIP simulations. The maximum over the central United States seen in observations is placed farther to the west in simulations. Over the central plains the models exhibit modest skill in simulating low-frequency precipitation variability, a Palmer drought severity index proxy. The presence of a linear trend increases correlations in the period 1950–99 when compared with those for the whole century. The SST links of the Great Plains drought index have features in common with observations over both the Pacific and Atlantic Oceans. Interestingly, summer-to-fall precipitation regressions of the warm Trend, cold Pacific, and warm Atlantic modes of annual mean SST variability (used in forcing the DWG idealized model experiments) tend to dry the southwestern, midwestern, and southeastern regions of the United States in the observations and, to a lesser extent, in the simulations. The similarity of the idealized SST-forced droughts in DWG modeling experiments with AMIP precipitation regressions of the corresponding SST principal components, evident especially in the case of the cold Pacific pattern, suggests that the routinely conducted AMIP simulations could have served as an effective proxy for the more elaborated suite of DWG modeling experiments. [ABSTRACT FROM AUTHOR]

The variance budget of land surface hydrological quantities is analyzed in the second Atmospheric Model Intercomparison Project (AMIP2) simulation made with the Canadian Centre for Climate Modelling and Analysis (CCCma) third-generation general circulation model (AGCM3). The land surface parameterization in this model is the comparatively sophisticated Canadian Land Surface Scheme (CLASS). Second-order statistics, namely variances and covariances, are evaluated, and simulated variances are compared with observationally based estimates. The soil moisture variance is related to second-order statistics of surface hydrological quantities. The persistence time scale of soil moisture anomalies is also evaluated. Model values of precipitation and evapotranspiration variability compare reasonably well with observationally based and reanalysis estimates. Soil moisture variability is compared with that simulated by the Variable Infiltration Capacity-2 Layer (VIC-2L) hydrological model driven with observed meteorological data. An equation is developed linking the variances and covariances of precipitation, evapotranspiration, and runoff to soil moisture variance via a transfer function. The transfer function is connected to soil moisture persistence in terms of lagged autocorrelation. Soil moisture persistence time scales are shorter in the Tropics and longer at high latitudes as is consistent with the relationship between soil moisture persistence and the latitudinal structure of potential evaporation found in earlier studies. In the Tropics, although the persistence of soil moisture anomalies is short and values of the transfer function small, high values of soil moisture variance are obtained because of high precipitation variability. At high latitudes, by contrast, high soil moisture variability is obtained despite modest precipitation variability since the persistence time scale of soil moisture anomalies is long. Model evapotranspiration estimates show little variability and soil moisture variability is dominated by precipitation and runoff, which account for about 90% of the soil moisture variance over land surface areas. [ABSTRACT FROM AUTHOR]

The seasonal and annual climatological behavior of selected components of the hydrological cycle are presented from coupled and uncoupled configurations of the atmospheric component of the Community Climate System Model (CCSM) Community Atmosphere Model version 3 (CAM3). The formulations of processes that play a role in the hydrological cycle are significantly more complex when compared with earlier versions of the atmospheric model. Major features of the simulated hydrological cycle are compared against available observational data, and the strengths and weaknesses are discussed in the context of specified sea surface temperature and fully coupled model simulations. The magnitude of the CAM3 hydrological cycle is weaker than in earlier versions of the model, and is more consistent with observational estimates. Major features of the exchange of water with the surface, and the vertically integrated storage of water in the atmosphere, are generally well captured on seasonal and longer time scales. The water cycle response to ENSO events is also very realistic. The simulation, however, continues to exhibit a number of long-standing biases, such as a tendency to produce double ITCZ-like structures in the deep Tropics, and to overestimate precipitation rates poleward of the extratropical storm tracks. The lower-tropospheric dry bias, associated with the parameterized treatment of convection, also remains a simulation deficiency. Several of these biases are exacerbated when the atmosphere is coupled to fully interactive surface models, although the larger-scale behavior of the hydrological cycle remains nearly identical to simulations with prescribed distributions of sea surface temperature and sea ice. [ABSTRACT FROM AUTHOR]

A Global Land Data Assimilation System (GLDAS) has been developed. Its purpose is to ingest satellite- and ground-based observational data products, using advanced land surface modeling and data assimilation techniques, in order to generate optimal fields of land surface states and fluxes. GLDAS is unique in that it is an uncoupled land surface modeling system that drives multiple models, integrates a huge quantity of observation-based data, runs globally at high resolution (0.25°), and produces results in near–real time (typically within 48 h of the present). GLDAS is also a test bed for innovative modeling and assimilation capabilities. A vegetation-based “tiling” approach is used to simulate subgrid-scale variability, with a 1-km global vegetation dataset as its basis. Soil and elevation parameters are based on high-resolution global datasets. Observation-based precipitation and downward radiation and output fields from the best available global coupled atmospheric data assimilation systems are employed as forcing data. The high-quality, global land surface fields provided by GLDAS will be used to initialize weather and climate prediction models and will promote various hydrometeorological studies and applications. The ongoing GLDAS archive (started in 2001) of modeled and observed, global, surface meteorological data, parameter maps, and output is publicly available. [ABSTRACT FROM AUTHOR]