Montane snowpack is a vital source of water supply in the Western United States. However, the future of snow in these regions in a changing climate is uncertain. Here, we use a large-ensemble approach to evaluate the consistency across 124 statistically downscaled snow water equivilent (SWE) projections between end-of-century (2076 – 2095) and early 21st century (2106 – 2035) periods. Comparisons were performed on dates corresponding with the end of winter (15 April) and spring snowmelt (15 May) in five western US montane domains. By benchmarking SWE climate change signals using the disparity between snow projections, we identified relationships between SWE projections that were repeatable across each domain, but shifted in elevation. In low to mid-elevations, 15 April average projected decreases to SWE of 48% or larger were greater than the disparity between models. Despite this, a significant portion of 15 April SWE volume (39 – 93%) existed in higher elevation regions where the disparities between snow projections exceeded the projected changes to SWE. Results also found that 15 April and 15 May projections were strongly correlated (r 0.82), suggesting that improvements to the spread and certainty of 15 April SWE projections would translate to improvements in later dates. The results of this study show that large-ensemble approaches can be used to measure coherence between snow projections and identify both 1) the highest-confidence changes to future snow water resources, and 2) the locations and periods where and when improvements to snow projections would most benefit future snow projections.

A document by Yifan Zhou. Click on the document to view its contents.

Current lightning predictions are uncertain because they either rely on empirical diagnostic relationships based on the present climate or use coarse-scale climate scenario simulations in which deep convection is parameterized. Previous studies demonstrated that simulations with convection-permitting resolutions (km-scale) improve lightning predictions compared to coarser-grid simulations using convection parameterization for different geographical locations but not over the boreal zone. In this study, lightning simulations with the NASA Unified-Weather Research and Forecasting (NU-WRF) model are evaluated over a 955x540 km2 domain including the Great Slave Lake in Canada for six lightning seasons. The simulations are performed at convection-parameterized (9 km) and convection-permitting (3 km) resolution using the Goddard 4ICE and the Thompson microphysics (MP) schemes. Four lightning indices are evaluated against observations from the Canadian Lightning Detection Network (CLDN), in terms of spatiotemporal frequency distribution, spatial pattern, daily climatology, and an event-based overall skill assessment. Concerning the model configuration, regardless of the spatial resolution, the Thompson scheme is superior to the Goddard 4ICE scheme in predicting the daily climatology but worse in predicting the spatial patterns of lightning occurrence. Several evaluation metrics indicate the benefit of working at a convection-permitting resolution. The relative performance of the different lightning indices depends on the evaluation criteria. Finally, this study demonstrates issues of the models to reproduce the observed spatial pattern of lightning well, which might be related to an insufficient representation of land surface heterogeneity in the study area.

Traditional watershed modeling often overlooks the role of vegetation dynamics. While past studies indicate possible improved hydrologic predictions by increasing the physical realism of vegetation dynamics in process-based models, there has been little quantitative evidence to support similar improvements in water quality predictions. To fill this knowledge-gap, we recently applied a modified Soil and Water Assessment Tool (SWAT) to quantify the extent of improvements that the assimilation of remotely sensed Leaf Area Index (LAI) would convey to streamflow, soil moisture, and nitrate load simulations across a 16,860 km2 agricultural watershed in the midwestern United States. We modified the SWAT source code to directly insert spatially distributed and temporally continuous LAI estimates from Moderate Resolution Imaging Spectroradiometer (MODIS). Compared to a “basic” traditional model with limited spatial information, our LAI assimilation model (i) significantly improved daily streamflow simulations during medium-to-low flow conditions, (ii) provided realistic spatial distributions of growing season soil moisture, and (iii) substantially reproduced the long-term observed variability of daily nitrate loads. Further analysis revealed that assimilation of MODIS LAI data corrected the model’s LAI overestimation tendency, which led to a proportionally increased rootzone soil moisture and decreased plant nitrogen uptake. With these new findings, our study confirms that assimilation of MODIS LAI data in watershed models can effectively improve both hydrology and water quality predictions.

The Snow Ensemble Uncertainty Project (SEUP) is an effort to establish a baseline characterization of snow water equivalent (SWE) uncertainty across North America with the goal of informing global snow observational needs. An ensemble-based modeling approach, encompassing a suite of current operational models, is used to assess the uncertainty in SWE and total snow storage (SWS) estimation during the 2009-2017 period. The highest modeled SWE uncertainty is observed in mountainous regions, likely due to the relatively deep snow, forcing uncertainties, and variability between the different models in resolving the snow processes over complex terrain. This highlights a need for high-resolution observations in mountains to capture the high spatial SWE variability. The greatest SWS is found in Tundra regions where, even though the spatiotemporal variability in modeled SWE is low, there is considerable uncertainty in the SWS estimates due to the large areal extent over which those estimates are spread. This highlights the need for high accuracy in snow estimations across the Tundra. In mid-latitude boreal forests, large uncertainties in both SWE and SWS indicate that vegetation-snow impacts are a critical area where focused improvements to modeled snow estimation efforts need to be made. Finally, the SEUP results indicate that SWE uncertainty is driving runoff uncertainty and measurements may be beneficial in reducing uncertainty in SWE and runoff, during the melt season at high latitudes (e.g., Tundra and Taiga regions) and in the Western mountain regions, whereas observations at (or near) peak SWE accumulation are more helpful over the mid-latitudes.

Irrigation is the largest human intervention in the water cycle that can modulate climate extremes, yet global irrigation water use (IWU) remains largely unknown. Microwave remote sensing offers a practical way to quantify IWU by monitoring changes in soil moisture caused by irrigation. This study evaluates the ability to quantify IWU by assimilating high-resolution (1km) SMAP-Sentinel 1 (SMAP-S1) remotely sensed soil moisture with a physically-based land surface model (LSM) using a particle batch smoother (PBS). A suite of synthetic experiments is devised to evaluate different error sources. Results from the synthetic experimentation show that unbiased simulations with known irrigation timing can produce an accurate irrigation estimate with a mean annual bias of 0.45% and the mean R2 of 96.5%, relative to observed IWU. Unknown irrigation timing can significantly deteriorate the model performance by increasing the mean annual bias to 23% and decreasing the mean R2 to 36%. In real-world experiments, the PBS data assimilation approach provides a mean bias of -18.6% when the timing of irrigation water use is known. This underestimation is possibly attributable to missing part of the irrigation signal. Yet, significantly higher irrigation was estimated over the irrigated pixels compared to the non-irrigated pixels, indicating that data assimilation can skillfully convey irrigation signals to LSMs. LSM calibration provides a 10% improvement to soil moistrue RMSE relative to the open-loop simulation. PBS data assimilation provides an additional 50% improvement to simulated soil moisture RMSE by correcting the model state and superimposing the optimal (unmodeled) irrigation on precipitation forcing.

It is known that representing wetland dynamics in land surface modeling improves models’ capacity to reproduce fluxes and land surface boundary conditions for atmospheric modeling in general circulation models. This study presents the development of the full coupling between the Noah-MP land surface model (LSM) and the HyMAP flood model in the NASA Land Information System and its application over the Inner Niger Delta (IND), a well-known hot-spot of strong land surface-atmosphere interactions in West Africa. Here, we define two experiments at 0.02º spatial resolution over the 2002-2018 period to quantify the impacts of the proposed developments on IND dynamics. One represents the one-way approach for simulating land surface and flooding processes (1-WAY), i.e., Noah-MP neglects surface water availability, and the proposed two-way coupling (2-WAY), where Noah-MP takes surface water availability into account in the vertical water and energy balance. Results show that accounting for two-way interactions between Noah-MP and HyMAP over IND improves all selected hydrological variables. Compared to 1-WAY, evapotranspiration derived from 2-WAY over flooding zones doubles, increased by 0.8mm/day, resulting in an additional water loss rate of ~18,900km3/year, ~40% drop of wetland extent during wet seasons and major improvement in water level variability at multiple locations. Significant soil moisture increase and surface temperature drop were also observed. Wetland outflows decreased by 35%, resulting in a substantial a Nash-Sutcliffe coefficient improvement, from -0.73 to 0.79. It is anticipated that future developments in global water monitoring and water‐related disaster warning systems will considerably benefit from these findings.

This first paper of the two-part series focuses on demonstrating the predictability of a hyper-resolution, offline terrestrial modeling system used for the High Mountain Asia (HMA) region. To this end, this study systematically evaluates four sets of model simulations at point scale, basin scale, and domain scale obtained from different spatial resolutions including 0.01 degree (∼ 1-km) and 0.25 degree (∼ 25-km). The assessment is conducted via comparisons against ground-based observations and satellite-derived reference products. The key variables of interest include surface net shortwave radiation, surface net longwave radiation, skin temperature, near-surface soil temperature, snow depth, snow water equivalent, and total runoff. In the evaluation against ground-based measurements, the superiority of the 0.01 degree estimates are mostly demonstrated across relatively complex terrain. Specifically, hyper-resolution modeling improves the skill in meteorological forcing estimates (except precipitation) by 9% relative to coarse-resolution estimates. The model forced by downscaled forcings in its entirety yields the highest predictability skill in model output states as well as precipitation, which improves the skill obtained by coarse-resolution estimates by 7%. These findings, on one hand, corroborate the importance of employing the hyper-resolution versus coarse-resolution modeling in areas characterized by complex terrain. On the other hand, by evaluating four sets of model simulations forced with different precipitation products, this study emphasizes the importance of accurate hyper-resolution precipitation products to drive model simulations.

Satellite observations of coastal Louisiana indicate an overall land loss over recent decades, which could be attributed to climate- and human-induced factors, including sea level rise (SLR). Climate-induced hydrological change (CHC) has impacted the way flood control structures are used, altering the spatiotemporal water distribution. Based on “what-if” scenarios, we determine relative impacts of SLR and CHC on increased flood risk over southern Louisiana and examine the role of water management via flood control structures in mitigating flood risk over the region. Our findings show that CHC has increased flood risk over the past 28 years. The number of affected people increases as extreme hydrological events become more exceptional. Water management reduces flood risk to urban areas and croplands, especially during exceptional hydrological events. For example, currently (i.e., 2016-2020 period), CHC-induced flooding puts an additional 73km2 of cropland under flood risk at least half of the time (median flood event) and 65km2 once a year (annual flood event), when compared to a past period (1993-1997). A ten- to twenty-fold increase relative to SLR-induced flooding. CHC also increases population vulnerability in southern Louisiana to flooding; additional 9900 residents currently live under flood risk at least half of the time, and that number increases to 27,400 for annual flood events. Residents vulnerable to SLR-induced flooding is lower (6000 and 3300 residents, respectively). Conclusions are that CHC is a major factor that should be accounted for flood resilience and that water management interventions can mitigate risks to human life and activities.