SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
An ensemble of land surface models and forcing data was developed to assess variability in SWE estimation over North America. In this study, the ensemble output was used to assess how SWE uncertainty impacts stream flow estimation. The analysis was conducted by major basins of North America over the 2009-2017 time period.
A new airborne synthetic aperture radar (SAR) system was recently developed for the estimation of snow water equivalent (SWE). The radar is part of the SWESARR (Snow Water Equivalent Synthetic Aperture Radar and Radiometer) instrument, an active passive microwave system specifically designed for the accurate estimation of SWE. The dual polarization (VV, VH) radar operates at three frequency bands (9.65 GHz, 13.6 GHz, and 17.25 GHz), with bandwidths of up to 200 MHz. The radar flew its first flight campaign in November 2019, along with SWESARR’s - already operational – radiometer. The radar collected comprehensive data sets over various terrains that show a successful system performance. The inst slated to participate in future SnowEx campaigns.
When compared with differences in snow accumulation predicted by widely used hydrological models, there is a much greater divergence among otherwise “good” models in their simulation of the snow ablation process. Here, we explore differences in the performance of the Variable Infiltration Capacity model (VIC), Noah land surface model with multiparameterization options (Noah-MP), the Catchment model, and the third-generation Simplified Simple Biosphere model (SiB3) in their ability to reproduce observed snow water equivalent (SWE) during the ablation season at 10 Snowpack Telemetry (SNOTEL) stations over 1992–2012. During the ablation period, net radiation generally has stronger correlations with observed melt rates than does air temperature. Average ablation rates tend to be higher (in both model predictions and observations) at stations with a large accumulation of SWE. The differences in the dates of last snow between models and observations range from several days to approximately a month (on average 5.1 days earlier than in observations). If the surface cover in the models is changed from observed vegetation to bare soil in all of the models, only the melt rate of the VIC model increases. The differences in responses of models to canopy removal are directly related to snowpack energy inputs, which are further affected by different algorithms for surface albedo and energy allocation across the models. We also find that the melt rates become higher in VIC and lower in Noah-MP if the shrub/grass present at the observation sites is switched to trees.
This study quantifies the level of observational accuracy required from a spaceborne light detection and ranging (LiDAR) snow depth retrieval mission for enabling beneficial impacts for snow estimation. The study is conducted over a region in Western Colorado using a suite of observing system simulation experiments(OSSEs).The Joint UK Land Environment Simulator, version 5.0 (JULES v5.0) is employed to simulate a suite of idealized LiDAR observations, considering a range of LiDAR snow depth retrieval errors, different hypothetical sensor swath widths, and the impact of cloud cover on observability. These simulated observations are then assimilated into the Noah land surface modelwith multi-parameterization options, version 3.6 (Noah-MP v3.6) model. This data assimilation setup is used to systematically evaluate the potential utility of LiDAR observations for improving modeled snow water equivalent (SWE)estimates and water budget variables such as runoff. Results from the OSSE runsshow that, in general, assimilation of synthetic LiDAR observations provide beneficial impacts when theLiDAR snow depth retrieval error standard deviation (σerror) is below 60 cm.Based on comparisons between the realistic (i.e., swath-limited and cloud-attenuated) case and the idealized (i.e., infinite swath width in the absence of cloud cover) case,this study concludes that observations with a conservative error standard deviation threshold of 40 cm (i.e., upper limit of the snow depth retrieval error that adds value to the SWE estimates viaassimilation) are needed for improving modeled snow estimates. More than a 33% reduction in SWEroot mean square errors and more than a 15% increase in correlation coefficientsare achieved when σ error ≤ 40 cmusing a 170-km sensor swath width in the presence of cloud attenuation effects. Further, the integrated hydrologic response, as represented by total (surface and subsurface) runoff estimates during the snow ablation season, are also enhanced when assimilating synthetic LiDAR snow depth retrievals with errors below this level.
he NASA LIS/WRF-Hydro system is a coupled modeling framework that combines the modeling and data assimilation (DA) capabilities of the NASA Land Information System (LIS) with the multi-scale surface hydrological modeling capabilities of the WRF-Hydro model, both of which are widely used in both operations and research. This coupled modeling framework builds on the linkage between land surface models (LSMs), which simulate surface boundary conditions in atmospheric models, and distributed hydrologic models, which simulate horizontal surface and sub-surface flow, adding new land DA capabilities. In the present study, we employ this modeling framework in the Tuolumne River basin in central California. We demonstrate the added value of the assimilation of NASA Airborne Snow Observatory (ASO) snow water equivalent (SWE) estimates in the Tuolumne basin. This analysis is performed in both LIS as an LSM column model and LIS/WRF-Hydro, with hydrologic routing. Results demonstrate that ASO DA in the basin reduced snow bias by as much as 30% from an open-loop (OL) simulation compared to three independent datasets. It also reduces downstream streamflow runoff biases by as much as 40%, and improves streamflow skill scores in both wet and dry years. Analysis of soil moisture and evapotranspiration (ET) also reveals the impacts of hydrologic routing from WRF-Hydro in the simulations, which would otherwise not be resolved in an LSM column model. By demonstrating the beneficial impact of SWE DA on the improving streamflow forecasts, the article outlines the importance of such observational inputs for reservoir operations and related water management applications.
Snow is a fundamental component of global and regional water budgets, particularly in mountainous areas and regions downstream that rely on snowmelt for water resources. Land surface models (LSMs) are commonly used to develop spatially distributed estimates of snow water equivalent (SWE) and runoff. However, LSMs are limited by uncertainties in model physics and parameters, among other factors. In this study, we describe the use of model calibration tools to improve snow simulations within the Noah-MP LSM as the first step in an observing system simulation experiment (OSSE). Noah-MP is calibrated against the University of Arizona (UA) SWE product over a western Colorado domain. With spatially varying calibrated parameters, we run calibrated and default Noah-MP simulations for water years 2010–20. By evaluating both simulations against the UA dataset, we show that calibration decreases domain averaged temporal RMSE and bias for snow depth from 0.15 to 0.13 m and from −0.036 to −0.0023 m, respectively, and improves the timing of snow ablation. Increased snow simulation performance also improves estimates of model-simulated runoff in four of six study basins, though only one has statistically significant improvement. Spatially distributed Noah-MP snow parameters perform better than default uniform values. We demonstrate that calibrating variables related to snow albedo calculations and rain–snow partitioning, among other processes, is a necessary step for creating a nature run that reasonably approximates true snow conditions for the OSSEs. Additionally, the inclusion of a snowfall scaling term can address biases in precipitation from meteorological forcing datasets, further improving the utility of LSMs for generating reliable spatiotemporal estimates of snow.
The NASA LIS/WRF-Hydro system is a coupled modeling framework that combines the modeling and data assimilation (DA) capabilities of the NASA Land Information System (LIS) with the multi-scale surface hydrological modeling capabilities of the WRF-Hydro model, both of which are widely used in both operations and research. This coupled modeling framework builds on the linkage between land surface models (LSMs), which simulate surface boundary conditions in atmospheric models, and distributed hydrologic models, which simulate horizontal surface and sub-surface flow, adding new land DA capabilities. In the present study, we employ this modeling framework in the Tuolumne River basin in central California. We demonstrate the added value of the assimilation of NASA Airborne Snow Observatory (ASO) snow water equivalent (SWE) estimates in the Tuolumne basin. This analysis is performed in both LIS as an LSM column model and LIS/WRF-Hydro, with hydrologic routing. Results demonstrate that ASO DA in the basin reduced snow bias by as much as 30% from an open-loop (OL) simulation compared to three independent datasets. It also reduces downstream streamflow runoff biases by as much as 40%, and improves streamflow skill scores in both wet and dry years. Analysis of soil moisture and evapotranspiration (ET) also reveals the impacts of hydrologic routing from WRF-Hydro in the simulations, which would otherwise not be resolved in an LSM column model. By demonstrating the beneficial impact of SWE DA on the improving streamflow forecasts, the article outlines the importance of such observational inputs for reservoir operations and related water management applications.
This second paper of the two-part series focuses on demonstrating the impact of assimilating satellite-based snow cover and freeze/thaw observations into the hyper-resolution, offline terrestrial modeling system used for the High Mountain Asia (HMA) region from 2003 to 2016. To this end, this study systematically evaluates a total of six sets of 0.01° (∼1 km) model simulations forced by different precipitation forcings, with and without the dual assimilation scheme enabled, at point-scale, basin-scale, and domain-scale. 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 (SWE), and total runoff. First, the point-scale assessment is mainly conducted via evaluating against ground-based measurements. In general, the assimilation enabled estimates are better than no-assimilation counterparts. Second, the basin-scale runoff assessment demonstrates that across three snow-dominated basins, the assimilation enabled experiment yields systematic improvements in all goodness-of-fit statistics through mitigating the negative effects brought by the fixed long-term precipitation correction factors. For example, when forced by the bias-corrected precipitation, the assimilation-enabled experiment improves the bias by 69%, the root-mean-squared error by 30%, and the unbiased root-mean-squared error by 18% (relative to the no-assimilation counterpart). Finally, the domainscale assessment is conducted via evaluating against satellite-based SWE and skin temperature products. Both sets of domain-scale analysis further corroborate the findings in the point-scale evaluations. Overall, this study suggests the benefits of the proposed multi-variate assimilation system in improving the cryospherichydrological process within a land surface model for use in HMA.
This second paper of the two-part series focuses on demonstrating the impact of assimilating satellite-based snow cover and freeze/thaw observations into the hyper-resolution, offline terrestrial modeling system used for the High Mountain Asia (HMA) region from 2003 to 2016. To this end, this study systematically evaluates a total of six sets of 0.01° (∼1 km) model simulations forced by different precipitation forcings, with and without the dual assimilation scheme enabled, at point-scale, basin-scale, and domain-scale. 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 (SWE), and total runoff. First, the point-scale assessment is mainly conducted via evaluating against ground-based measurements. In general, the assimilation enabled estimates are better than no-assimilation counterparts. Second, the basin-scale runoff assessment demonstrates that across three snow-dominated basins, the assimilation enabled experiment yields systematic improvements in all goodness-of-fit statistics through mitigating the negative effects brought by the fixed long-term precipitation correction factors. For example, when forced by the bias-corrected precipitation, the assimilation-enabled experiment improves the bias by 69%, the root-mean-squared error by 30%, and the unbiased root-mean-squared error by 18% (relative to the no-assimilation counterpart). Finally, the domain-scale assessment is conducted via evaluating against satellite-based SWE and skin temperature products. Both sets of domain-scale analysis further corroborate the findings in the point-scale evaluations. Overall, this study suggests the benefits of the proposed multi-variate assimilation system in improving the cryospheric-hydrological process within a land surface model for use in HMA.
Using information theory, our study quantifies the importance of selected indicators for the U.S. Drought Monitor (USDM) maps. We use the technique of mutual information (MI) to measure the importance of any indicator to the USDM, and because MI is derived solely from the data, our findings are independent of any model structure (conceptual, physically-based, or empirical). We also compare these MIs against the drought representation effectiveness ratings in the North America Drought Indices and Indicators Assessment (NADIIA) survey for Koeppen climate zones. This reveals: [1] agreement between some ratings and our MI values (high for example indicators like Standardized Precipitation-Evapotranspiration Index or SPEI); [2] some divergences (for example, soil moisture has high ratings but near-zero MIs for ESA-CCI soil moisture in the Western U.S., indicating the need of another remotely sensed soil moisture source); and [3] new insights into the importance of variables such as Snow Water Equivalent (SWE) that are not included in sources like NADIIA. Further analysis of the MI results yields findings related to: [1] hydrological mechanisms (summertime SWE domination during individual drought events through snowmelt into the water-scarce soil); [2] hydroclimatic types (the top pair of inputs in the Western and non-Western regions are SPEIs and soil moistures respectively); and [3] predictability (high for the California 2012-2017 event, with longer-timescale indicators dominating). Finally, the high MIs between multiple indicators jointly and the USDM indicate potentially high drought forecasting accuracies achievable using only model-based inputs, and the potential for global drought monitoring using only remotely sensed inputs, especially for locations having insufficient in situ observations.
Understanding the relationship between remotely sensed snow disappearance and seasonal water supply may become vital in coming years to supplement limited ground based, in situ measurements of snow in a changing climate. For the period 2001–2019, we investigated the relationship between satellite derived Day of Snow Disappearance (DSD)—the date at which snow has completely disappeared—and the seasonal water supply, i.e., the April—July total streamflow volume, for 15 snow dominated basins across the western U.S. A Monte Carlo framework was applied, using linear regression models to evaluate the predictive skill—defined here as a model’s ability to accurately predict seasonal flow volumes—of varied predictors, including DSD and in situ snow water equivalent (SWE), across a range of spring forecast dates. In all basins there is a statistically significant relationship between mean DSD and seasonal water supply (p ≤ 0.05), with mean DSD explaining roughly half of the variance. Satellite-based model skill improves later in the forecast season, surpassing the skill of in-situ-based (SWE) models in skill in 10 of the 15 basins by the latest forecast date. We found little to no correlation between model error and basin characteristics such as elevation and the ratio of snow water equivalent to total precipitation. Despite a relatively short data record, this exploratory analysis shows promise for improving seasonal water supply prediction, in particular for snow dominated basins lacking in situ observations.
Extremes events are occurring more frequently and are projected to increase in the future. This data set includes extreme event indices that were generated for the Colorado River basin using data from the Variable Infiltration Capacity hydrologic model. VIC simulations were run using six CMIP5 earth system models (ESMs), two different time periods, (historical, 1970-1999) and future (2070-2099), and for one Representative Concentration Pathway (RCP 8.5). ESMs available are MPI-ESM-LR, MIROC-ESM, IPSL-CM5A-LR, GFDL-ESM2M, and GFDL-ESM2G. The different indicators include temperature (min/max, freezing days, heating days), precipitation (dry days, max), runoff (min, max), soil moisture (min, max), evapotranspiration (max), snow water equivalent (SWE, max), wind speeds (max), and pressure (min). Detailed information on the indicators, file names, and how they were calculated is included in the README.txt document. Data sets were used to calculate results are in related references Bennett et al. 2020 and Talsma et al. 2021. Details on the VIC model configuration is described in related references Bennett et al. 2018 and Bennett et al. 2017.
Under a changing climate, increases in dry conditions and extreme heat events are projected to exacerbate wildfire risk in the northeastern U.S. In recent years, Vermont has observed higher annual temperatures, more frequent heatwaves, increased annual precipitation, extreme flood events, and decreased snowfall. The mechanisms through which environmental factors contribute to increased fire risk in humid environments, such as Vermont, are poorly understood. The team partnered with the National Weather Service, the Vermont Division of Forests, and the University of Vermont to investigate phenological trends and antecedent conditions influencing wildland fire risk. For the phenological analysis, from 2001 to 2023, vegetation data, phenological dates, and snow water equivalent (SWE) values were accessed from Landsat 5Thematic Mapper (TM), Landsat 7Enhanced Thematic Mapper Plus (ETM+), Landsat 8Operational Land Imager (OLI), Landsat 9OLI-2, the Moderate Resolution Imaging Spectroradiometer (MODIS), and the Snow Data Assimilation System (SNODAS), respectively. For the antecedent condition analysis, from 2008 to 2023, soil moisture, Environmental Stress Index, wind speed, minimum relative humidity, and daily precipitation data were obtained from the Global Land Data Assimilation System (GLDAS), SERVIR, gridMET Wind, gridMET Humidity, and NClimGrid, respectively. The study found that green-up dates over the study period remain relatively stable, while snowmelt dates appear increasingly variable. Minimum relative humidity was the most significant environmental variable correlated with wildfire risk in Vermont. Results from this study will inform the National Weather Service's preparation of fire forecasts before prescribed burns and support community outreach by the Vermont Agency of Natural Resources.
Estimation of riverbed profiles, also known as bathymetry, plays a vital role in many applications, such as safe and efficient inland navigation, prediction of bank erosion, land subsidence, and flood risk management. The high cost and complex logistics of direct bathymetry surveys, i.e, depth imaging, have encouraged the use of indirect measurements such as surface flow velocities. However, estimating high-resolution bathymetry from indirect measurements is an inverse problem that can be computationally challenging. Here, we propose a reduced-order model (ROM) based approach that utilizes a variational autoencoder (VAE), a type of deep neural network with a narrow layer in the middle, to compress bathymetry and flow velocity information and accelerate bathymetry inverse problems from flow velocity measurements. In our application, the shallow-water equations (SWE) with appropriate boundary conditions (BCs), e.g., the discharge and/or the free surface elevation, constitute the forward problem, to predict flow velocity. Then, ROMs of the SWEs are constructed on a nonlinear manifold of low dimensionality through a variational encoder and the bathymetry inversion problem is derived on the low-dimensional latent space in a Hierarchical Bayesian setting. Further, the reformulation allows variational inference with a small number (e.g., $\mathscr{O}$ (100) of ROM runs and efficient uncertainty quantification. We have tested our inversion approach on a one-mile reach of the Savannah River, GA, USA. Once the neural network is trained (offline stage), the proposed technique can perform the inversion operation orders of magnitude faster than traditional inversion methods that are commonly based on linear projections, such as principal component analysis (PCA), or the principal component geostatistical approach (PCGA). Furthermore, tests show that the algorithm can estimate the bathymetry with good accuracy even with sparse flow velocity measurements.
Storm-snow avalanches are challenging to forecast due to complex alpine terrain and during rapidly changing weather conditions. They can result in loss of lives and significant economic impact. We describe how a new device that continuously measures with high-frequency snowflake mass, size, density, and type, the Differential Emissivity Imaging Disdrometer (DEID), and show how the DEID can be used to aid avalanche forecasting when coupled with a storm-snow stability model. DEID measurements of snow accumulation, snow water equivalent (SWE), and snow density obtained during seventeen storms taken at the mid-Collins Snow-Study Plot at Alta Ski Area in Utah's Central Wasatch mountain range during winter 2020–2021 show excellent agreement with infrequent manual measurements. Additionally, two new variables, the Shape Density Index (SDI) and Complexity, are proposed and used to classify snowflake habit and estimate storm-snow shear strength. We illustrate how these DEID-derived data can be used to identify layers of concern in the storm snow such as density inversions, in real-time without digging snow pits. Furthermore, the DEID-data are used to run four variations of the SNOw Slope Stability model (SNOSS) for the storms investigated. The results are evaluated with data collected from tilt-board tests, infrasound measurements, and visual observations of avalanches. For a total fourteen storms analyzed, the DEID-driven SNOSS-modeled minimum stability index predicts the general stability of the storm-snow as indicated by observed avalanches, both natural and of unknown cause. Finally, the results provide a promising approach for nowcasting instabilities within storm-snow layers with a single instrument.
Abstract The Noah‐MP land surface model (LSM) relies on the Monin‐Obukhov (M‐O) Similarity Theory (MOST) to calculate land‐atmosphere exchanges of water, energy, and momentum fluxes. However, MOST flux‐profile relationships neglect canopy‐induced turbulence in the roughness sublayer (RSL) and parameterize within‐canopy turbulence in an ad hoc manner. We implement a new physics scheme (M‐O‐RSL) into Noah‐MP that explicitly parameterizes turbulence in RSL. We compare Noah‐MP simulations employing the M‐O‐RSL scheme (M‐O‐RSL simulations) and the default M‐O scheme (M‐O simulations) against observations obtained from 647 Snow Telemetry (SNOTEL) stations and two AmeriFlux stations in the western United States. M‐O‐RSL simulations of snow water equivalent (SWE) outperform M‐O simulations over 64% and 69% of SNOTEL sites in terms of root‐mean‐square‐error (RMSE) and correlation, respectively. The largest improvements in skill for M‐O‐RSL occur over closed shrubland sites, and the largest degradations in skill occur over deciduous broadleaf forest sites. Differences between M‐O and M‐O‐RSL simulated snowpack are primarily attributable to differences in aerodynamic conductance for heat underneath the canopy top, which modulates sensible heat flux. Differences between M‐O and M‐O‐RSL within‐canopy and below‐canopy sensible heat fluxes affect the amount of heat transported into snowpack and hence change snowmelt when temperatures are close to or above the melting point. The surface energy budget analysis over two AmeriFlux stations shows that differences between M‐O and M‐O‐RSL simulations can be smaller than other model biases (e.g., surface albedo). We intend for the M‐O‐RSL physics scheme to improve performance and uncertainty estimates in weather and hydrological applications that rely on Noah‐MP.
Localized projections of 21st‐century hydroclimate variables obtained from downscaling Global Climate Model (GCM) output are central to informing regional impact assessments and infrastructure planning. Regional GCM biases can be significant and, for dynamical downscaling, can be addressed either before (a priori) or after (a posteriori) downscaling. However, a priori bias correction (APBC) has generally unexplored effects on climate change signals. Here we analyze dynamically downscaled solutions of CMIP6 GCMs over the Western U.S., with and without APBC, and quantify APBC's impact on climate change signals relative to other irreducible uncertainty sources. For temperature and precipitation, the uncertainty introduced by APBC is negligible compared to that arising from GCM choice or internal variability. Furthermore, APBC greatly reduces regional models' unrealistically high snow‐water‐equivalent (SWE) biases that result directly from GCM errors. We leverage this finding to encourage the dynamical downscaling community to adopt APBC as a standard operating procedure.
Abstract This article presents a hydrological reconstruction of the Upper Colorado River Basin with an hourly temporal resolution, and 1-km spatial resolution from October 1982 to September 2019. The validated dataset includes a suite of hydrologic variables including streamflow, water table depth, snow water equivalent (SWE) and evapotranspiration (ET) simulated by an integrated hydrological model, ParFlow-CLM. The dataset was validated over the period with a combination of point observations and remotely sensed products. These datasets provide a long-term, natural-flow, simulation for one of the most over-allocated basins in the world.