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.
A summary is given of the glacial and marine chronology of Mars. Hydrological models of oceans and ice sheets, the cratering record, hydrological cycling, and episodic glaciation are discussed. Evidence for a Noachian ocean is evaluated.
Ice velocities observed in 2005/06 at three GPS stations along the Sermeq Avannarleq flowline, West Greenland, are used to characterize an observed annual velocity cycle. We attempt to reproduce this annual ice velocity cycle using a 1-D ice-flow model with longitudinal stresses coupled to a 1-D hydrology model that governs an empirical basal sliding rule. Seasonal basal sliding velocity is parameterized as a perturbation of prescribed winter sliding velocity that is proportional to the rate of change of glacier water storage. The coupled model reproduces the broad features of the annual basal sliding cycle observed along this flowline, namely a summer speed-up event followed by a fall slowdown event. We also evaluate the hypothesis that the observed annual velocity cycle is due to the annual calving cycle at the terminus. We demonstrate that the ice acceleration due to a catastrophic calving event takes an order of magnitude longer to reach CU/ETH ('Swiss') Camp (46km upstream of the terminus) than is observed. The seasonal acceleration observed at Swiss Camp is therefore unlikely to be the result of velocity perturbations propagated upstream via longitudinal coupling. Instead we interpret this velocity cycle to reflect the local history of glacier water balance.
A snow hydrology has been implemented in an atmospheric general circulation model (GCM). The snow hydrology consists of parameterizations of snowfall and snow cover fraction, a prognostic calculation of snow temperature, and a model of the snow mass and hydrologic budgets. Previously, only snow albedo had been included by a specified snow line. A 3-year GCM simulation with this now more complete surface hydrology is compared to a previous GCM control run with the specified snow line, as well as with observations. In particular, the authors discuss comparisons of the atmospheric and surface hydrologic budgets and the surface energy budget for U.S. and Canadian areas. The new snow hydrology changes the annual cycle of the surface moisture and energy budgets in the model. There is a noticeable shift in the runoff maximum from winter in the control run to spring in the snow hydrology run. A substantial amount of GCM winter precipitation is now stored in the seasonal snowpack. Snow cover also acts as an important insulating layer between the atmosphere and the ground. Wintertime soil temperatures are much higher in the snow hydrology experiment than in the control experiment. Seasonal snow cover is important for dampening large fluctuations in GCM continental skin temperature during the Northern Hemisphere winter. Snow depths and snow extent show good agreement with observations over North America. The geographic distribution of maximum depths is not as well simulated by the model due, in part, to the coarse resolution of the model. The patterns of runoff are qualitatively and quantitatively similar to observed patterns of streamflow averaged over the continental United States. The seasonal cycles of precipitation and evaporation are also reasonably well simulated by the model, although their magnitudes are larger than is observed. This is due, in part, to a cold bias in this model, which results in a dry model atmosphere and enhances the hydrologic cycle everywhere.
Abstract. We introduce a hierarchical Bayesian model for the spatial distribution of rainfall corresponding to an extreme event of a specified duration that could be used with regional hydrologic models to perform a regional hydrologic risk analysis. An extreme event is defined if any gaging site in the watershed experiences an annual maximum rainfall event and the spatial field of rainfall at all sites corresponding to that occurrence is modeled. Applications to data from New York City demonstrate the effectiveness of the model for providing spatial scenarios that could be used for simulating loadings into the urban drainage system. Insights as to the homogeneity in spatial rainfall and its implications for modeling are provided by considering partial pooling in the hierarchical Bayesian framework.
Diverse vulnerabilities of Bangladesh's agricultural sector in 16 sub-regions are assessed using experiments designed to investigate climate impact factors in isolation and in combination. Climate information from a suite of global climate models (GCMs) is used to drive models assessing the agricultural impact of changes in temperature, precipitation, carbon dioxide concentrations, river floods, and sea level rise for the 2040-2069 period in comparison to a historical baseline. Using the multi-factor impacts analysis framework developed in Yu et al. (2010), this study provides new sub-regional vulnerability analyses and quantifies key uncertainties in climate and production. Rice (aman, boro, and aus seasons) and wheat production are simulated in each sub-region using the biophysical Crop Environment REsource Synthesis (CERES) models. These simulations are then combined with the MIKE BASIN hydrologic model for river floods in the Ganges-Brahmaputra-Meghna (GBM) Basins, and the MIKE21Two-Dimensional Estuary Model to determine coastal inundation under conditions of higher mean sea level. The impacts of each factor depend on GCM configurations, emissions pathways, sub-regions, and particular seasons and crops. Temperature increases generally reduce production across all scenarios. Precipitation changes can have either a positive or a negative impact, with a high degree of uncertainty across GCMs. Carbon dioxide impacts on crop production are positive and depend on the emissions pathway. Increasing river flood areas reduce production in affected sub-regions. Precipitation uncertainties from different GCMs and emissions scenarios are reduced when integrated across the large GBM Basins' hydrology. Agriculture in Southern Bangladesh is severely affected by sea level rise even when cyclonic surges are not fully considered, with impacts increasing under the higher emissions scenario.
Digital elevation models are now widely used to calculate terrain parameters to determine incoming solar and longwave radiation for use in surface climate models, interpretation of remote-sensing data, and parameters in hydrologic models. Because of the large number of points in an elevation grid, fast algorithms are useful to save computation time. A description is given of rapid methods for calculating slope and azimuth, solar illumination angle, horizons, and view factors for radiation from sky and terrain. Calculation time is reduced by fast algorithms and lookup tables.
Rising atmospheric carbon dioxide concentrations are expected to enhance photosynthesis and reduce crop water use. However, there is high uncertainty about the global implications of these effects for future crop production and agricultural water requirements under climate change. Here we combine results from networks of field experiments and global crop models to present a spatially explicit global perspective on crop water productivity (CWP, the ratio of crop yield to evapotranspiration) for wheat, maize, rice and soybean under elevated carbon dioxide and associated climate change projected for a high-end greenhouse gas emissions scenario. We find carbon dioxide effects increase global CWP by 10[0;47]%-27[7;37]% (median[interquartile range] across the model ensemble) by the 2080s depending on crop types, with particularly large increases in arid regions (by up to 48[25;56]% for rain fed wheat). If realized in the fields, the effects of elevated carbon dioxide could considerably mitigate global yield losses whilst reducing agricultural consumptive water use (4-17%). We identify regional disparities driven by differences in growing conditions across agro-ecosystems that could have implications for increasing food production without compromising water security. Finally, our results demonstrate the need to expand field experiments and encourage greater consistency in modeling the effects of rising carbon dioxide across crop and hydrological modeling communities.
High-resolution predictions of hydrometeorological variables are critical for supporting hydropower generation decisions and flood control at hydroelectric power plants. Traditional climate and hydrologic models rely on the numerical simulation of detailed physical processes. Therefore, running these simulations is time-, labor-, and computation-intensive. Improving the spatial and temporal resolution in these modeling outputs could lead to cubic increases in both the simulation time and computational demands, rendering high-resolution hydrometeorological predictions expensive and impractical. Many past studies apply the super resolution (SR) technique to downscale climate models using deep learners. However, deep learners are deemed “black-boxes,” as their derivation processes from low-resolution outputs to high-resolution outputs are often hidden. Their results are difficult for domain scientists to interpret and validate. Thus, there is a need for an exploratory machine learning approach that can partially integrate domain-specific theory and knowledge into the data-driven mapping process between simulation outputs of different spatial scales. The domain-specific theory and knowledge can be incorporated into the data model through an inductive approach in which process-related environmental variables are used and analyzed as key drivers (i.e., environmental surrogates) to reflect the complex physical processes. Many of these variables, such as land use land cover, soil types, topography, digital elevation, air temperature, and various watershed characteristics, can be directly measured through sensors or remote sensing techniques. Additionally, SR applications that can downscale hydrological and hydrodynamics models to efficiently produce high-resolution (1 m) flood depth grids are still rare. Since the flood depth grid can be used to support critical decisions for flood control operation at hydroelectric power plants, it is crucial to enable an SR-based capability for interpolating high-resolution flood inundation maps.
Effective flood simulation capabilities can tremendously support early warning and disaster prevention. To examine the applicability of a fully physics-based and high-performance flood simulation and forecasting modeling framework for a flood-prone region in Taiwan, we conduct a numerical experiment that couples the Weather Research and Forecasting (WRF) Model, WRF-Hydrological modeling system (WRF-Hydro), and the Two-Dimensional Runoff Inundation Toolkit for Operational Needs (TRITON) to perform integrated rainfall, streamflow, and flood simulations. Furthermore, we first use the coupled WRF and WRF-Hydro (WWH) to predict rainfall and streamflow and then drive TRITON with the predicted streamflow hydrographs to simulate flood depth and inundation area. With the refined spatial resolution and parameterization, this framework can better predict rainfall with reasonable spatial patterns. Although WWH could overestimate the amount of rainfall in some areas, the uncertain rainfall–streamflow predictions produce reasonable flood maps able to pinpoint regions at risk of flooding. In terms of model efficiency, the graphics processing unit–based computation can yield a speed-up factor as high as ∼13 compared to the central processing unit–based computation, promoting the efficacy of the coupled modeling framework in practical real-time flood forecasting.
High-resolution precipitation data, generated through dynamical downscaling (DD) or statistical downscaling (SD) of global climate model output, provide critical information for regional climate assessment and adaptation planning. Most downscaling development and validation have focused on accurate gridscale precipitation construction and ignored the spatial structure of precipitation across model grids and at the event scale. However, many applications, e.g., hydrologic modeling and the analysis using the downscaled precipitation, require a reasonable representation of the spatial structure of precipitation within watersheds. Therefore, a set of standard metrics to evaluate the representation of the spatial structure of individual storms across diverse downscaled precipitation products is desired. To address this need, we conducted an object-based evaluation of precipitation in decades-long DD and SD products over the contiguous United States (CONUS). Specifically, we evaluate their ability to reproduce various features of precipitation objects in the observations: total volume, precipitation area, peak intensity, and spatial structure. Multiple metrics (bias, Perkins score, and nonparametric statistical tests) are used to quantify model performance. Our evaluation reveals notable variations in performance among individual products across different climate zones and seasons, as well as between extreme and nonextreme events. In general, most DD products exhibit balanced performance across the four precipitation object features, while SD products vary more significantly in their performance across products. Based on this comprehensive evaluation, we provide guidance on choosing downscaled products for specific regions, seasons, and precipitation object features. These findings and recommendations can inform precipitation-relevant modeling and analysis over CONUS, guide future downscaling technique developments, and provide actionable information for climate impact assessment and adaptation.
High-resolution precipitation data, generated through dynamical downscaling (DD) or statistical downscaling (SD) of global climate model output, provide critical information for regional climate assessment and adaptation planning. Most downscaling development and validation have focused on accurate gridscale precipitation construction and ignored the spatial structure of precipitation across model grids and at the event scale. However, many applications, e.g., hydrologic modeling and the analysis using the downscaled precipitation, require a reasonable representation of the spatial structure of precipitation within watersheds. Therefore, a set of standard metrics to evaluate the representation of the spatial structure of individual storms across diverse downscaled precipitation products is desired. To address this need, we conducted an object-based evaluation of precipitation in decades-long DD and SD products over the contiguous United States (CONUS). Specifically, we evaluate their ability to reproduce various features of precipitation objects in the observations: total volume, precipitation area, peak intensity, and spatial structure. Multiple metrics (bias, Perkins score, and nonparametric statistical tests) are used to quantify model performance. Our evaluation reveals notable variations in performance among individual products across different climate zones and seasons, as well as between extreme and nonextreme events. In general, most DD products exhibit balanced performance across the four precipitation object features, while SD products vary more significantly in their performance across products. Based on this comprehensive evaluation, we provide guidance on choosing downscaled products for specific regions, seasons, and precipitation object features. These findings and recommendations can inform precipitation-relevant modeling and analysis over CONUS, guide future downscaling technique developments, and provide actionable information for climate impact assessment and adaptation.
It is well known that soil moisture is a characteristic of the land surface that strongly affects the partitioning of outgoing radiation into sensible and latent heat which significantly impacts both weather and climate. Detailed land surface schemes are now being coupled to mesoscale atmospheric models in order to represent the effect of soil moisture upon atmospheric simulations. However, there is little direct soil moisture data available to initialize these models on regional to continental scales. As a result, a Soil Hydrology Model (SHM) is currently being used to generate an indirect estimate of the soil moisture conditions over the continental United States at a grid resolution of 36 Km on a daily basis since 8 May 1995. The SHM is forced by analyses of atmospheric observations including precipitation and contains detailed information on slope soil and landcover characteristics.The purpose of this paper is to evaluate the utility of initializing a detailed coupled model with the soil moisture data produced by SHM.
Climate change (CC) is likely to affect the thousands of bar-built or barrier estuaries (here referred to as Small tidal inlets - STIs) around the world. Any such CC impacts on the stability of STIs, which governs the dynamics of STIs as well as that of the inlet-adjacent coastline, can result in significant socio-economic consequences due to the heavy human utilisation of these systems and their surrounds. This article demonstrates the application of a process based snap-shot modelling approach, using the coastal morphodynamic model Delft3D, to 3 case study sites representing the 3 main STI types; Permanently open, locationally stable inlets (Type 1), Permanently open, alongshore migrating inlets (Type 2) and Seasonally/Intermittently open, locationally stable inlets (Type 3). The 3 case study sites (Negombo lagoon - Type 1, Kalutara lagoon - Type 2, and Maha Oya river - Type 3) are all located along the southwest coast of Sri Lanka. After successful hydrodynamic and morphodynamic model validation at the 3 case study sites, CC impact assessment are undertaken for a high end greenhouse gas emission scenario. Future CC modified wave and riverflow conditions are derived from a regional scale application of spectral wave models (WaveWatch III and SWAN) and catchment scale applications of a hydrologic model (CLSM) respectively, both of which are forced with IPCC Global Climate Model output dynamically downscaled to approximately 50 km resolution over the study area with the stretched grid Conformal Cubic Atmospheric Model CCAM. Results show that while all 3 case study STIs will experience significant CC driven variations in their level of stability, none of them will change Type by the year 2100. Specifically, the level of stability of the Type 1 inlet will decrease from 'Good' to 'Fair to poor' by 2100, while the level of (locational) stability of the Type 2 inlet will also decrease with a doubling of the annual migration distance. Conversely, the stability of the Type 3 inlet will increase, with the time till inlet closure increasing by approximately 75%. The main contributor to the overall CC effect on the stability of all 3 STIs is CC driven variations in wave conditions and resulting changes in longshore sediment transport, not Sea level rise as commonly believed.
In ‘Satellite observations and modeling to understand the Lower Mekong River Basin streamflow variability’ [1] hydrological fluxes, meteorological variables, land cover land use maps, and soil characteristics and parameters data were compiled and processed for the Lower Mekong River Basin. In this work, daily streamflow time series data at nine gauges located at five different countries in the Mekong region (Thailand, Laos People׳s Democratic Republic (PDR), Myanmar, Cambodia, and Viet Nam) is presented. Satellite-based daily precipitation and air temperature (minimum & maximum) data is processed and provided over the entire basin as part of the dataset provided in this work. Moreover, land cover land use raster data that contains 18 classes that cover agriculture, urban, range and forests land cover land use classes for the basin is offered. In addition, a soil data that contains physical and chemical characteristics needed by physically based hydrological models to simulate the cycling of water and air is also provided.
This data package was generated to support the manuscript “Towards CONUS-Wide Machine Learning-Augmented Conceptually Interpretable Modeling of Catchment-Scale Precipitation-Storage-Runoff Dynamics.” It provides input files, model outputs, plotting data, scripts, notebooks, and documentation used to develop, evaluate, and reproduce Mass-Conserving Perceptron (MCP)-based hydrologic modeling experiments across 513 selected Catchment Attributes and Meteorology for Large-sample Studies in the United States (CAMELS-US) basins. The files are organized by modeling component and analysis purpose, including rainfall–runoff experiments, snow module experiments, coupled hydrologic-snow experiments, Long Short-Term Memory (LSTM) benchmark results, model skill metrics, initialization and epoch records, cell-state normalization files, Akaike Information Criterion (AIC)-based model comparison files, and data used to generate manuscript figures. Tabular files can be opened using standard spreadsheet software or Python/R data-analysis tools. Python scripts, Jupyter notebooks, and selected MATLAB scripts are included for model execution, postprocessing, plotting, and statistical analysis. Quality assurance and quality control were conducted through the source-data selection and modeling workflow. Meteorological forcing, streamflow, and static catchment attributes were derived from the CAMELS-US dataset, and snow water equivalent data were derived from the University of Arizona (UA) Snow Water Equivalent dataset. Selected basins and time periods were screened during the associated research workflow to avoid missing observations or poor-quality cases. Static geospatial features were processed primarily using Quantum Geographic Information System (QGIS) and Geospatial Data Abstraction Library (GDAL) workflows. Additional details are provided in the associated manuscript and documentation.
A new, revised three-dimensional (3-D) hydrostratigraphic framework model (HFM) for the Central and Western Pahute Mesa Corrective Action Units 101 and 102 was completed in 2019. The initial (Phase I) model was completed in 2002 and was documented in A Hydrostratigraphic Model and Alternatives for the Groundwater Flow and Contaminant Transport Model of Corrective Action Units 101 and 102: Central and Western Pahute Mesa, Nye County, Nevada (BN, 2002). Subsequent flow and transport modeling revealed the need to collect additional data to reduce uncertainties. The Phase II data collection and characterization effort included extensive well drilling and testing activities, and rebuilding of the Phase I 3-D HFM, resulting in the Phase II HFM (NSTec, 2014). Further drilling (Well ER-20-12), data collection, and interpretive work (e.g., Thirsty Canyon Lineament [Wurtz and Day, 2018]) were incorporated into the Pahute Mesa-Oasis Valley (PM-OV) HFM. The PM-OV HFM incorporates the area of interest from the Phase II HFM and adds additional areas, primarily to the north and west with minor additions to the south and east. The new, expanded area matches the PM-OV groundwater basin as defined in Fenelon et al. (2016) and includes Pahute Mesa, a former nuclear testing area at the Nevada National Security Site; and Oasis Valley, a groundwater discharge area downgradient from contaminant source areas on Pahute Mesa. The data and information necessary to build the PM-OV HFM are provided in Navarro (2019). The PM-OV HFM will be used by hydrologic modelers who are tasked with developing a model to determine how contaminants are transported by groundwater flow in an area of complex geology. The model area is large (more than 6,250 square kilometers) and geologically complex, including Paleozoic- to Mesozoic-age structures (e.g., thrust faults, Thirsty Canyon Lineament), at least seven Tertiary-age calderas, several intrusive bodies, and many relatively recent basin-and-range normal faults. Investigators from the U.S. Geological Survey; National Security Technologies, LLC; National Laboratories; and U.S. Department of Energy contractors have organized the more than 300 volcanic units and other rocks in the model area into 77 hydrostratigraphic units (HSUs). The volcanic rocks are subdivided into 42 aquifers, 11 confining units, and 9 composite units (containing both aquifer and confining unit rocks). The underlying pre-Tertiary rocks were divided into 5 HSUs, comprising 2 aquifers and 3 confining units. The model also includes 1 alluvial aquifer unit, 8 intrusive confining units, and 1 granite confining unit. The model depicts the thickness, extent, and geometric relationships of these HSUs (“layers” in the model) along with all the major structural features that control them, including calderas and faults. The most substantive differences between the Phase II HFM and the PM-OVHFM include (1) addition of the new data and resulting revised interpretations from the additional area encompassed by the PM-OV groundwater basin; (2) incorporation of 16 new faults; (3) 12 boreholes from the northern extension; (4) 2 additional caldera-collapse collars; and (5) correlation of the simulated faulting and hydrostratigraphy in the eastern part of the model domain with the delineations represented in the neighboring Rainier Mesa-Shoshone Mountain HFM-Prime and Yucca Flat-Climax Mine HFM. All PM-OV HFM delineations were based on published and unpublished data from geophysical investigations, surface geologic maps, and drill-hole records. Phase II approached model uncertainty using the model-building software to test model configurations and hydrologic consequences using a range of values for parameters such as fault offset amounts and HSU elevations. The PM-OV HFM expands on this approach to encompass the entire PM-OV groundwater basin.
Floods and associated landslides are one of the most widespread natural hazards on Earth, responsible for tens of thousands of deaths and billions of dollars in property damage every year. During 1993-2002, over 1000 of the more than 2,900 natural disasters reported were due to floods. These floods and associated landslides claimed over 90,000 lives, affected over 1.4 billion people and cost about $210 billion. The impact of these disasters is often felt most acutely in less developed regions. In many countries around the world, satellite-based precipitation estimation may be the best source of rainfall data due to lack of surface observing networks. Satellite observations can be of essential value in improving our understanding of the occurrence of hazardous events and possibly in lessening their impact on local economies and in reducing injuries, if they can be used to create reliable warning systems in cost-effective ways. This article addressed these opportunities and challenges by describing a combination of satellite-based real-time precipitation estimation with land surface characteristics as input, with empirical and numerical models to map potential of landslides and floods. In this article, a framework to detect floods and landslides related to heavy rain events in near-real-time is proposed. Key components of the framework are: a fine resolution precipitation acquisition system; a comprehensive land surface database; a hydrological modeling component; and landslide and debris flow model components. A key precipitation input dataset for the integrated applications is the NASA TRMM-based multi-satellite precipitation estimates. This dataset provides near real-time precipitation at a spatial-temporal resolution of 3 hours and 0.25deg x 0.25deg. By careful integration of remote sensing and in-situ observations, and assimilation of these observations into hydrological and landslide/debris flow models with surface topographic information, prediction of useful probabilistic maps of landslide and floods for emergency management in a timely manner is possible. Early results shows that the potential exists for successful application of satellite precipitation data in improving/developing global monitoring systems for flood/landslide disaster preparedness and management. The scientific and technological prototype can be first applied in a representative test-bed and then the information deliverables for the region can be tailored to the societal and economic needs of the represented affected countries.
This data package contains integrated hydrological model input and output files for Copper Creek, Colorado (24 km2), a tributary of the East River located in the headwaters of the Upper Colorado River Basin. The model code is the U.S. Geological Survey (USGS) Groundwater and Surface Water Flow (GSFLOW) model. The model contains a 100-m grid resolution and a daily timestep. The land surface model is dynamically linked to a three-dimensional groundwater flow model that allows for streamflow gaining and losing conditions. The groundwater model contains 12 model layers and extends 400 m below land surface. The original Copper Creek model was modified to contain geologic layers representing saprolite, shallow bedrock, and deep bedrock. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. For the shallow case, median flow depths occur in the shallow saprolite at depths <8 m, while the deep case promotes a median groundwater flow depth of 100 m. With this modeling framework we compare streamflow response to a plausible worst-case drought lasting up to five years. Streamflow metrics of analysis include average streamflow, fraction of stream network that is dry, no-flow duration, average groundwater flow to streams and time to recovery following the drought. Results and implications are presented in a paper submitted to Geophysical Research Letters titled, "The role of bedrock circulation depth and porosity in mountain streamflow response to prolonged drought" by Rosemary WH. Carroll, Andrew H. Manning and Kenneth H Williams. A Readme.txt file provides instructions on how to download all model files and execute each model scenario. In addition to the GSFLOW output/prms/copper_drought.csv file containing daily basin water stores and fluxes (refer to GSFLOW manual) and the output/prms/copper_drought_statvar.dat file with output defined in the gsflow3.control file (refer to GSFLOW Manual), output files also include spatially distributed daily values of total evapotranspiration, canopy evaporation, precipitation, snowfall, infiltration, snow water equivalent, potential evapotranspiration, recharge, sublimation, soil moisture, contributing interflow, water table elevations, changes in groundwater storage, groundwater evapotranspiration, interbasin groundwater flow (limited to the alluvium below the stream outlet), and surface-groundwater exchanges within the river system.