An official website of the United States government
Abstract The land surface hydrology of the North American Great Lakes region regulates ecosystem water availability, lake levels, vegetation dynamics, and agricultural practices. In this study, we analyze the Great Lakes terrestrial water budget using the Noah‐MP land surface model to characterize the catchment hydrological regimes and identify the dominant quantities contributing to the variability in the land surface hydrology. We show that the Great Lakes domain is not hydrologically uniform and strong spatiotemporal differences exist in the regulators of the hydrological budget at daily, monthly, and annual timescales. Subseasonally, precipitation and soil moisture explain nearly all the terrestrial water budget variability in the southern basins, while the northern latitudes are snow‐dominated regimes. Seasonal assessments reveal greater differences among the basins. Precipitation, evaporation, and runoff are the dominant sources of variability at lower latitudes, while at higher latitudes, terrestrial water storage in the form of ground snowpack and soil moisture has the leading role. Differences in land cover categorizations, for example, croplands, forests, or urban zones, further induce spatial differences in the hydrological characteristics. This quantification of variability in the terrestrial water cycle embedded at different temporal scales is important to assess the impacts of changes in climate and land cover on catchment sensitivities across the diverse hydroclimate of the Great Lakes region.
more »
« less
Toward Understanding Parametric Controls on Runoff Sensitivity to Climate in the Community Land Model: A Case Study Over the Colorado River Headwaters
Abstract Crucial to the assessment of future water security is how the land model component of Earth System Models partition precipitation into evapotranspiration and runoff, and the sensitivity of this partitioning to climate. This sensitivity is not explicitly constrained in land models nor the model parameters important for this sensitivity identified. Here, we seek to understand parametric controls on runoff sensitivity to precipitation and temperature in a state‐of‐the‐science land model, the Community Land Model version 5 (CLM5). Process‐parameter interactions underlying these two climate sensitivities are investigated using the sophisticated variance‐based sensitivity analysis. This analysis focuses on three snow‐dominated basins in the Colorado River headwaters region, a prominent exemplar where land models display a wide disparity in runoff sensitivities. Runoff sensitivities are dominated by indirect or interaction effects between a few parameters of subsurface, snow, and plant processes. A focus on only one kind of parameters would therefore limit the ability to constrain the others. Surface runoff exhibits strong sensitivity to parameters of snow and subsurface processes. Constraining snow simulations would require explicit representation of the spatial variability across large elevation gradients. Subsurface runoff and soil evaporation exhibit very similar sensitivities. Model calibration against the subsurface runoff flux would therefore constrain soil evaporation. The push toward a mechanistic treatment of processes in CLM5 have dampened the sensitivity of parameters compared to earlier model versions. A focus on the sensitive parameters and processes identified here can help characterize and reduce uncertainty in water resource sensitivity to climate change.
more »
« less
- PAR ID:
- 10562868
- Publisher / Repository:
- Wiley
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 60
- Issue:
- 12
- ISSN:
- 0043-1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract. A common parameter in hydrological modeling frameworks is root zone water storage capacity (SR[L]), which mediates plant water availability during dry periods as well as the partitioning of rainfall between runoff and evapotranspiration. Recently, a simple flux-tracking-based approach was introduced to estimate the value of SR (Wang-Erlandsson et al., 2016). Here, we build upon this original method, which we argue may overestimate SR in snow-dominated catchments due to snow melt and evaporation processes. We propose a simple extension to the method presented by Wang-Erlandsson et al. (2016) and show that the approach provides a lower estimate of SR in snow-dominated watersheds. This SR dataset is available at a 1 km resolution for the continental USA, along with the full analysis code, on the Google Colab and Earth Engine platforms. We highlight differences between the original and new methods across the rain–snow transition in the Southern Sierra Nevada, California, USA. As climate warms and precipitation increasingly arrives as rain instead of snow, the subsurface may be an increasingly important reservoir for storing plant-available water between wet and dry seasons; therefore, improved estimates of SR will better clarify the future role of the subsurface as a storage reservoir that can sustain forests during seasonal dry periods and episodic drought.more » « less
-
null (Ed.)Abstract. Land models are essential tools for understanding and predicting terrestrial processes and climate–carbon feedbacks in the Earth system, but uncertainties in their future projections are poorly understood. Improvements in physical process realism and the representation of human influence arguably make models more comparable to reality but also increase the degrees of freedom in model configuration, leading to increased parametric uncertainty in projections. In this work we design and implement a machine learning approach to globally calibrate a subset of the parameters of the Community Land Model, version 5 (CLM5) to observations of carbon and water fluxes. We focus on parameters controlling biophysical features such as surface energy balance, hydrology, and carbon uptake. We first use parameter sensitivity simulations and a combination of objective metrics including ranked global mean sensitivity to multiple output variables and non-overlapping spatial pattern responses between parameters to narrow the parameter space and determine a subset of important CLM5 biophysical parameters for further analysis. Using a perturbed parameter ensemble, we then train a series of artificial feed-forward neural networks to emulate CLM5 output given parameter values as input. We use annual mean globally aggregated spatial variability in carbon and water fluxes as our emulation and calibration targets. Validation and out-of-sample tests are used to assess the predictive skill of the networks, and we utilize permutation feature importance and partial dependence methods to better interpret the results. The trained networks are then used to estimate global optimal parameter values with greater computational efficiency than achieved by hand tuning efforts and increased spatial scale relative to previous studies optimizing at a single site. By developing this methodology, our framework can help quantify the contribution of parameter uncertainty to overall uncertainty in land model projections.more » « less
-
Abstract. Land surface models (LSMs) are critical components of Earth system models (ESMs), enabling the simulation of energy and water fluxes that are essential for understanding climate systems. Soil hydraulic parameters, derived using pedotransfer functions (PTFs), are crucial for modeling soil–plant–water interactions; they introduce uncertainties in soil moisture simulations. However, a key knowledge gap exists in understanding how specific soil hydraulic properties contribute to these uncertainties and in identifying the regions most affected by them. This study conducts an intra-model sensitivity analysis within the Community Land Model version 5 (CLM5), examining how alternative soil parameter settings influence soil moisture variability across the contiguous United States (CONUS) using empirical orthogonal function (EOF) analysis. The EOF analysis revealed dominant spatial and temporal patterns of soil moisture across multiple experimental configurations, highlighting the impact of soil parameter variability on hydrological processes. The results showed significant discrepancies in soil moisture simulations, particularly in the central Great Plains, which may be attributed to the combination of arid climatic conditions and limitations in modeling saturated hydraulic conductivity and soil water retention curves. Seasonal soil moisture dynamics showed broad similarity to ERA5-Land patterns, with differences in magnitude and phase, indicating the importance of refined parameterization, particularly in the representation of infiltration and drainage processes. Comparisons with ERA5-Land, used here solely as a model-based reference for pattern consistency, revealed stronger similarity in regions with consistent climatic gradients but persistent differences in hydrologically complex areas, particularly in areas with arid climates, such as the Great Plains, where hydrological processes remain difficult to represent. Because CLM5 is forced by GSWP3, whereas ERA5-Land is an offline HTESSEL replay forced by ERA5, differences reflect both forcing and structural contrasts in addition to parameter effects. This research demonstrates the necessity to refine soil parameter representations, utilize high-resolution datasets, and consider climatic variability to inform the model development of LSMs. Importantly, these findings also pave the way for future efforts that incorporate dynamic soil properties into LSMs. This work illustrates the influence of soil properties on simulated variability. While the analysis documents their importance, a future direction will be to develop approaches that allow these properties to vary dynamically within LSMs. This study contributes to ongoing efforts toward more integrated modeling frameworks that capture soil–hydrology–climate interactions.more » « less
-
null (Ed.)Abstract. Plant activity in semi-arid ecosystems is largely controlled by pulses of precipitation, making them particularly vulnerable to increased aridity expected with climate change. Simple bucket-model hydrology schemes in land surface models (LSMs) have had limited ability in accurately capturing semi-arid water stores and fluxes. Recent, more complex, LSM hydrology models have not been widely evaluated against semi-arid ecosystem in situ data. We hypothesize that the failure of older LSM versions to represent evapotranspiration, ET, in arid lands is because simple bucket models do not capture realistic fluctuations in upper layer soil moisture. We therefore predict that including a discretized soil hydrology scheme based on a mechanistic description of moisture diffusion will result in an improvement in model ET when compared to data because the temporal variability of upper layer soil moisture content better corresponds to that of precipitation inputs. To test this prediction, we compared ORCHIDEE LSM simulations from (1) a simple conceptual 2-layer bucket scheme with fixed hydrological parameters; and (2) a 11-layer discretized mechanistic scheme of moisture diffusion in unsaturated soil based on Richards equations against daily and monthly soil moisture and ET observations, together with data-derived transpiration / evaporation, T / ET, ratios, from six semi-arid grass, shrub and forest sites in the southwestern USA. The 11-layer scheme also has modified calculations of surface runoff, bare soil evaporation, and water limitation to be compatible with the more complex hydrology configuration. To diagnose remaining discrepancies in the 11-layer model, we tested two further configurations: (i) the addition of a term that captures bare soil evaporation resistance to dry soil; and (ii) reduced bare soil fraction. We found that the more mechanistic 11-layer model results better representation of the daily and monthly ET observations. We show that is likely because of improved simulation of soil moisture in the upper layers of soil (top 5 cm). Some discrepancies between observed and modelled soil moisture and ET may allow us to prioritize future model development. Adding a soil resistance term generally decreased simulated E and increased soil moisture content, thus increasing T and T / ET ratios and reducing the negative T / ET model-data bias. By reducing the bare soil fraction in the model, we illustrated that modelled leaf T is too low at sparsely vegetated sites. We conclude that a discretized soil hydrology scheme and associated developments improves estimates of ET by allowing the model to more closely match the pulse precipitation dynamics of these semi-arid ecosystems; however, the partitioning of T from bare soil evaporation is not solved by this modification alone.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1029/2024WR037718
