Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shallow water”

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.

At least 739 records · Page 41

The Role of Bedrock Circulation Depth and Porosity in Mountain Streamflow Response to Prolonged Drought

Quantitative understanding is lacking on how the depth of active groundwater circulation in bedrock affects mountain streamflow response to a multi-year drought. We use an integrated hydrological model to explore the sensitivity of a variety of streamflow metrics to bedrock circulation depth and porosity under a plausible extreme drought scenario lasting up to 5 years. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. With drought, a deeper circulation system with higher drainable porosity more effectively buffers minimum flow and significantly limits perennial stream loss in comparison to a shallow circulation system. Streamflow buffering is accomplished through extensive groundwater storage loss. However, deeper circulation systems experience prolonged recovery from drought in comparison to storage-limited shallow systems. Research highlights the importance of characterizing the deeper bedrock hydrogeology in mountainous watersheds to better understand and predict drought impacts on stream ecosystem health and water resource sustainability.

54 ENVIRONMENTAL SCIENCES↗

Do it yourself remote sensing: Generating an inexpensive, high tech, real science lake mapping project for the classroom

The utilization of modest equipment and software revealed bottom contours and water column conditions of a dynamic water body. Classroom discussions of field techniques and equipment capabilities followed by exercises with the data sets in cause-and-effect analysis all contributed to participatory education in the process of science. This project is presented as a case study of the value of engaging secondary and collegiate level students in planning, executing and appraising a real world investigation which they can directly relate to. A 1 km wide bay, experiencing marsh inflow, along an 8 km long lake situated 120 km north of Ottawa, Canada, on the glaciated Canadian Precambrian Shield was mapped in midsummer for submerged topography, bottom composition, temperature profile, turbudity, dissolved oxygen and biota distribution. Low level aerial photographs scanned into image processing software are permitting spatial classification of bottom variations in biology and geology. Instrumentation consisted of a portable sport fishing SONAR depth finder, an electronic lead line multiprobe with photocell, thermistor and dissolved oxygen sensors, a selective depth water sampler, portable pH meter, an underwater camera mounted on a home-made platform with a bottom-contact trigger and a disposable underwater camera for shallow survey work. Sampling transects were referenced using a Brunton hand transit triangulating several shore markers.

Metzger, Stephen M.↗

PBL Height from AIRS, GPS RO, and MERRA-2 Products in NASA GES DISC and Their 10 Year Seasonal Mean Intercomparison

Within the planetary boundary layer (PBL), surface forcing response, drag, turbulence, and vertical mixing are important processes and play a more critical role here than in the overlying “free atmosphere”. The PBL Height (PBLH) is an important parameter in climate models, weather forecasts, and air quality prediction. The NASA Goddard Earth Sciences Data and Information Services Center (GES DISC) provides data processing, archiving, and distribution services for numerous Earth science products. PBLH is a parameter in three products served by GES DISC, which are from the Atmospheric Infrared Sounder (AIRS), the Global Positioning System (GPS) radio occultation (RO) experiment, and the NASA reanalysis product Modern-Era Retrospective analysis for Research and Applications – 2 (MERRA-2). These products have different spatial and temporal resolutions and coverages, and their PBLH definitions are also different. To better serve the PBL research community, we have summarized the specifications of these products. A ten-year seasonal mean intercomparison is also conducted to provide further guidance to users. The intercomparison results show that MERRA-2 has a much shallower PBL than AIRS and GPS RO. An experimental study indicates the different PBLH definition in MERRA-2 caused smaller values of PBLH. The improvement of the water vapor retrieval in AIRS version 7 over version 6 results in the version 7 PBLH agreeing better with GPS RO and MERRA-2 than version 6, especially near the equator and low latitudes.

Feng Ding↗

Modeling Sea Ice Albedo and Transmittance Measurements With A Fully-Coupled Radiative Transfer Model

A rigorous treatment of the sea ice medium has been incorporated in the advanced Coupled Ocean-Atmosphere Radiative Transfer (COART) model. The inherent optical properties (IOPs) of brine pockets and air bubbles over the 0.25-4.0 μm spectral region are parameterized as a function of the vertical profile of the sea ice physical properties (temperature, salinity and density). We test the model performance using available albedo and transmittance measurements collected during the Impacts of Climate on the Ecosystems and Chemistry of the Arctic Pacific Environment (ICESCAPE) and the Surface Heat Budget of the Arctic Ocean (SHEBA) field campaigns. The observations are adequately simulated when at least three layers are used to represent bare (first-year and multi-year) ice, including a thin top layer characterized by low density and high scattering. Two layers can be sufficient to model isolated cases of multi-year ice, and apply well to ponded ice except for shallow ponds over thick ice. The albedo and transmittance of ponded ice in the visible are mainly determined by the optical properties of the ice underlying the water layer used to model the pond. Sensitivity results indicate that the air volume or ice density has the largest impact on the simulated fluxes. Possible contamination from light-absorbing impurities, such as black carbon or ice algae, is also implemented in the model and is able to effectively reduce the albedo and transmittance in the visible spectrum to further improve the model-observation agreement. The agreement between the observed and modeled spectra validates the parameterization of the sea ice IOPs, and endorses COART as an accurate tool for radiation studies in the cryosphere.

Coupled Ocean-9 Atmosphere Radiative Transfer (COA↗

Science Plan for the Deployment of the Third ARM Mobile Facility to the Southeastern United States at the Bankhead National Forest, Alabama (AMF3 BNF)

In 2018, the U.S. Department of Energy (DOE) held a workshop for the Atmospheric Radiation Measurement (ARM) (Mather and Voyles 2013) user facility to discuss critical climate challenges and locations where key ARM Mobile Facility (AMF) observational assets could impact Earth system modeling (ESM). As an outcome, the southeast United States (SE U.S.) was identified as a high-priority region to target climate-process studies that promote a deeper understanding of the climate system and bolster ARM interactions with the community to drive ESM advancement. The DOE ARM user facility is a globally recognized leader in deploying and operating strategically located observation sites around the world for studying the properties of aerosols and clouds and their interaction with radiation, precipitation, and the Earth’s surface. In partnering with the DOE Atmospheric System Research (ASR) program, ARM solicited a multi-agency Site Science Team approach to provide input and close interaction with ARM management towards a successful SE U.S. deployment of the ARM third Mobile Facility (AMF3) (Miller et al. 2016). These efforts included identifying key locations, science drivers and instruments, and measurement strategies to address the wider climate-process needs and ESM improvement. Community input served a vital role in establishing, refining, and informing the relevant drivers and decisions regarding this AMF3 deployment. The team has identified Northern Alabama (N. AL) as regionally representative to unlock the key opportunities that will improve our understanding and model representation of aerosol, cloud, and land surface processes and their couplings in the SE U.S. A defining aspect of the AMF3 deployment is its commitment to long-term (anticipated five-year) observations to mitigate potential seasonal-to-annual variability that often limits appropriate attribution of phenomena to local or larger-scale processes. The proposed location may leverage nearby surface networks and multi-agency and partner assets to enrich this multi-year deployment. One motivation is to understand the role of spatiotemporal variability (thermodynamic, land-surface) across aspects of the climate system, with our AMF3 team anticipating future demands on characterizing the relationships between local-to-regional cloud development and surface processes across a diverse patchwork of natural, managed, and urban landscapes as found throughout the N. AL regions. The main site targets an intact, representative, forested region – the Bankhead National Forest (BNF) – underscoring further team commitment to regionally important land atmosphere two-way interactive studies “from the canopy to the clouds”, with enhanced tower instrumentation augmenting traditional ARM capabilities adjacent to this site. Multiple supplemental sites will also be distributed across this region, prioritizing added needs for biodiversity. Anticipated high-priority cloud science themes will target N. AL as a regional SE U.S. hotbed for high-impact weather, convective cloud onset, and shallow to-deep cloud transitioning. Anticipated aerosol drivers will focus on chemical processes that control the evolution of organic aerosol, the seasonality and spatial distribution of water vapor and particle-phase water, and its role on aerosol optical properties. Anticipated land atmosphere drivers consider the two-way feedbacks between surface influence on aerosols, clouds, and precipitation properties and the associated radiative impacts on plant physiology and canopy-scale fluxes. Emphasis will include the study of the impact of surface processes on aerosols via precursor emission, and on clouds via moisture flux and thermal development.

Doppler lidar, aerosols, convection↗

Soil Water Retention and Hydraulic Conductivity Data and Model at Pump House in East River Watershed, Colorado 2019-2024

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors near Pump House at Mount Crested Butte in the East River Watershed. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format ER-X-Y, where ER refers to East River, X is the location identifier, and Y is the depth identifier at the same X (shallow Y=1). Specifically, ER-PHS, ER-LMC, ER-LMF, and ER-SMN are associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and ER-RBTn (upslope n=1) are sampling transects during the 2019 Rootball Campaign. The sample and location information can be found in metadata.csv. Sampling and Measurements Each sample falls into one of the three sampling methods – (1) intact cores, (2) repacked samples, or (3) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. Both intact cores and repacked samples were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores and in-situ soil sensors could be prioritized because these sampling methods are less destructive. While the repacked samples were packed to the target bulk density (estimated post-sampling, when sample volume was known), these samples had altered pore structures. Nevertheless, intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method also has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗

Soil Water Retention and Hydraulic Conductivity Data and Model at Snodgrass Mountain in East River Watershed, Colorado 2020-2025

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors at Snodgrass Mountain. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format SG-X-Y, where SG refers to Snodgrass Mountain, X is the location identifier, and Y is the depth identifier at the same X (shallow Y=1). Specifically, SG-EHS is associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and SG-ERTn (upslope n=1) are points along the Snodgrass electrical resistivity tomography transect not associated with the existing site names in the directory. The sample and location information can be found in metadata.csv. Sampling and Measurements Each sample falls into one of the three sampling methods – (1) intact cores, (2) repacked samples, or (3) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. Both intact cores and repacked samples were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores and in-situ soil sensors could be prioritized because these sampling methods are less destructive. While the repacked samples were packed to the target bulk density (estimated post-sampling, when sample volume was known), these samples had altered pore structures. Nevertheless, intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method also has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗

On the possibilty of clathrate hydrates on the Moon

One of the most important inferences of the Lunar Prospector mission data was the existence of subsurface water ice in the permanently shadowed craters near both lunar poles [Feldman et al., 1998]. We propose and substantiate an alternative explanation that hydrogen can exist in the shallow lunar subsurface in the form of clathrate hydrates: CH4 . 6H(2)o and/or CO2 . 6H(2)o.

Lunar subsurface calthrate hydrates↗

Hail Monitor Sensor

An inexpensive and simple hail monitor design has been developed that has a single piezoelectric ceramic disc and uses a metal plate as a sounding board. The structure is durable and able to withstand the launch environment. This design has several advantages over a multi-ceramic sensor, including reduced cost and complexity, increased durability, and improvement in impact response uniformity over the active surface. However, the most important characteristic of this design is the potential to use frequency discrimination between the spectrum created from raindrop impact and a hailstone impact. The sound of hail hitting a metal plate is distinctly different from the sound of rain hitting the same plate. This fortuitous behavior of the pyramid sensor may lead to a signal processing strategy, which is inherently more reliable than one depending on amplitude processing only. The initial concept has been im proved by forming a shallow pyramid structure so that hail is encouraged to bounce away from the sensor so as not to be counted more than once. The sloped surface also discourages water from collecting. Additionally, the final prototype version includes a mounting box for the piezo-ceramic, which is offset from the pyramid apex, thus helping to reduce non-uniform response (see Figure 2). The frequency spectra from a single raindrop impact and a single ice ball impact have been compared. The most notable feature of the frequency resonant peaks is the ratio of the 5.2 kHz to 3.1 kHz components. In the case of a raindrop, this ratio is very small. But in the case of an ice ball, the ratio is roughly one third. This frequency signature of ice balls should provide a robust method for discriminating raindrops from hailstones. Considering that hail size distributions (HSDs) and fall rates are roughly 1 percent that of rainfall, hailstone sizes range from a few tenths of a centimeter to several centimeters. There may be considerable size overlap between large rain and small hail. As hail occurs infrequently at KSC, the ideal HSD measurement sensor needs to have a collection area roughly 100 times greater than a raindrop-size distribution sensor or disdrometer. The sensitivity should be such that it can detect and count very small hail in the midst of intense rainfall consisting of large raindrop sizes. The dynamic range and durability should allow measurement of the largest hail sizes, and the operation and calibration strategy should consider the infrequent occurrence of hail fall over the KSC area.

Younquist, Robert↗

Potential Chemical Impacts of Subsurface CO 2 : An Integrated Experimental and Numerical Assessment for a Case Study of the Ogallala Aquifer

Leakage from geologic CO 2 sequestration (GCS) reservoirs to overlying underground sources of drinking water (USDW) is a tangible risk. This study is an integrated assessment that combines column experiments and reactive transport simulations of sediments sampled from the Ogallala aquifer above an active commercial-scale GCS site (the Farnsworth Unit in northern Texas). Experimental and simulation results suggest that carbonate mineral (calcite and dolomite) dissolution is the most significant reaction following CO 2 intrusion, and is also the dominant source of trace metal release. Cation exchange is another key mechanism controlling trace metal release by cation interference. Most of the trace metals, including Ba, Sr, As, Pb, and Zn, show a short-term release and quickly drop to the baseline values, suggesting low risk to the overlying USDW quality. Other trace metals, such as Mn and U, exhibit a tangible increase of their concentrations in the beginning, and drop to a higher level compared to the baseline, which may become a potential concern of long-term USDW quality change with CO2introduction. This study provides a comprehensive example of combining laboratory experiments and simulations for assessment of CO 2 -sediment interactions with combined release mechanisms in shallow groundwater aquifers. Data presented here provides useful insights for quantitative risk assessment and effective public education regarding CO 2 geological sequestration, and trace metal reactive transport studies in shallow groundwater aquifers.

54 ENVIRONMENTAL SCIENCES↗

An Unusual Inverted Saline Microbial Mat Community in an Interdune Sabkha in the Rub' Alkhali (the Empty Quarter), UAE: an Analog for Habitats on Present Mars

Salt flats (sabkha) are a recognized habitat for microbial life in desert environments and as analogs for habitats for life on Mars. Here we report on the physical setting and microbiology of interdune sabkhas among the large dunes in the Rub' al Khali (the Empty Quarter) in Liwa Oasis, United Arab Emirates. The salt flats, composed of gypsum and halite, between the dunes are moistened by relatively fresh ground water from below. The result is a salinity gradient that is inverted compared to most salt flat communities with the hypersaline layer at the top and freshwater layers below. We describe and characterize a rich photosynthetically-based microbial ecosystem that is protected from the arid outside environment below the translucent salt crust. Gases collected from sediments under shallow ponds in the sabkha contain methane in concentrations as high as 3400 ppm. The salt layer provides environmental protection to the habitat below and could preserve biomarkers and other evidence for life in the salt after it dries out. Chloride-filled depressions have been identified on Mars and although the surface flow of water is unlikely on Mars today, ground water is possible. Such a near surface system with modern groundwater flowing under ancient salt deposits could be present on Mars and could be accessed by surface rovers.

Hypersaline environments in deserts↗

Offshore Wind Composite Siting Exclusions Scenarios

This dataset consolidates criteria that restrict the siting of offshore wind energy installations, an update from Zuckerman et al. (2023). All scenarios exclude area from wind energy development due to legal or administrative restrictions (e.g., Marine Protected Areas), existing infrastructure, and technical constraints (e.g., 1,300 m bathymetry limit). A TIF data file and a PNG cartographic map of the data are provided for three scenarios ("Open", "Reference", and "Limited"), showing areas where offshore wind energy is restricted across the contiguous United States. The "Open" scenario only excludes areas with known preclusions for offshore wind deployment. The "Reference" scenario excludes all areas from the "Open" scenario, as well as additional best management practices such as known proposed protected areas, and state waters. The "Limited" scenario excludes all area from the "Reference" scenario. Additionally, it has more stringent siting criteria such as larger setbacks from shore and infrastructure, as well as shallower depth limits on floating technology. For further details and citation, please refer to the publication linked below: Zuckerman, Gabriel R., Anthony Lopez, Travis Williams, Rebecca Green, and Grant Buster. 2023. Impacts of Siting Considerations on Offshore Wind Technical Potential in the United States. Golden, CO: National Renewable Energy Laboratory. NREL/TP-6A20- 85088.

Array↗

Connecting the Radiative Influences of Aerosol upon the Mass Flux Profiles of Shallow Cumuli across the Southeast Atlantic Ocean Basin and its Boundaries (Final Report)

The Atlantic Ocean covers approximately 25% of Earth’s surface and the atmosphere above it is home to a complex array of clouds and aerosols that have important influences on regional and global weather and climate. These influences must be accurately depicted in short range, medium range, seasonal, and climate forecast models. Conditions over the tropical Atlantic are particularly complex due to continental scale plumes of dust from the Sahara Desert and smoke from agricultural burning in Africa that drift across the Atlantic Ocean basin toward the Americas. These plumes are often found meandering in the lower atmosphere above shallow tropical clouds that form above the ocean surface, presumably mingling with these clouds on occasion due to convective mixing processes. Elevated dust and smoke particles absorb incoming sunlight and substantially warm the marine atmosphere in the layer in which they are present. This warming may alter the thermal stability of the marine atmosphere and may throttle or enhance the development of clouds, change their internal structure and the rate at which they precipitate. Alternatively, it may isolate the lower atmosphere from drier layers above enabling water vapor to accumulate near the ocean surface potentially leading to the development of deeper convection. Our study was organized around the principal concept of determining the impact of African smoke and dust plumes upon cloud development at ASI and understanding how a popular shallow convection parameterization used in models responds to the presence of this aerosol. We analyzed observations collected during the US Department of Energy (DOE) Layered Atlantic Smoke Interactions with Clouds (LASIC) campaign using the Atmospheric Radiation Measurement Program’s Mobile Facility #1 (AMF-1), which was deployed to Ascension Island (ASI) in Southeastern Atlantic for a one-year period. Ascension Island is immediately downwind from an African source of these plumes, but far enough removed to enable the lower atmosphere to have reacted to their presence. A main goal of our study was to compute radiative heating rate profiles over Ascension Island and along the trajectory from the biomass burning regions along coastal Africa to Ascension Island. To compute the radiative heating rate due to aerosols and clouds, we employed observed profiles of temperature, humidity, and clouds from LASIC alongside aerosol optical properties from the Modern-era Retrospective analysis for Research and Applications Version 2 (MERRA-2), as input for the Rapid Radiation Transfer Model (RRTM). Radiative heating was also assessed across the southeast Atlantic Ocean using an ensemble of back trajectories from the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model. We were successful in this effort and the resulting publication is already being cited despite its relatively short lifetime in the literature. The second part of our study involved the process-level representation of the clouds observed over the Southeastern Atlantic in models. Our initial task, upon which the balance of this portion of the study depended, was to evaluate the representativeness of the clouds observed by AMF-1. What orographic influence did Ascension Island have on the measured cloud structure? We set out to answer this question by trying to separate observations that were clearly indicative of orographic forcing from those that were representative of open ocean. The siting of AMF-1, which was 340-m above sea-level on the slope of a steep escarpment, proved demonstrably problematic for cloud process measurements. Despite considerable effort using multiple approaches including artificial intelligence (AI), we were unable to successfully compensate for the orographic forcing present at the AMF-1 site on ASI and, hence, produce process-level summaries of the convective mass flux representative of open ocean over the Southeastern Atlantic. Even though the project has officially ended, we are making a final attempt using two new types of clustering algorithm (i.e., AI) on the recommendation of a former student, who is an expert in this area. Since we have only recently begun to test these new AI schemes, the recommendations from the second part of our project outlined below are based on our experience at the time of this report.

54 ENVIRONMENTAL SCIENCES↗

Crack Growth Rate and Fracture Toughness Tests on Irradiated Ex-Plant Materials

The performance of structural materials is critical for the safe and economic operation of light water reactors. Exposed to neutron irradiation during service, the reactor core internal materials can undergo significant microstructural and microchemical changes, leading to irradiation hardening and embrittlement. To ensure the structural integrity and functionality of nuclear reactor components during long-term operation, material degradation and damage mechanisms must be understood and managed adequately. In this work, irradiated materials harvested from the decommissioned Zorita reactor were studied for their cracking susceptibility and fracture resistance as a function of irradiation dose up to 47 displacement per atom (dpa). The material is a Type 304 stainless steel sectioned from the baffle plates of this pressurized water reactor with 38 years of service. Crack growth rate and fracture toughness J-resistance (J-R) curve tests were performed in low-corrosion-potential environments at ~315°C. All samples behaved similarly under cyclic loading, and no deteriorated corrosion-fatigue behavior was observed in the test environments. Under constant stress intensity factors, most samples did not show elevated crack growth rates, suggesting an adequate stress corrosion cracking resistance in the test environments. However, an unstable cracking behavior was observed in a 47-dpa sample, resulting in significantly higher crack growth rates than expected at high stress intensity factors. Crack instability was also observed in a 0.06-dpa sample but did not lead to a sustained high crack growth rate. The impact of neutron irradiation was more evident in the fracture toughness J-R curve tests. As the dose increased, the J-R curve declined considerably and became very shallow at high doses. A fully intergranular fracture morphology was also observed among the high-dose samples ruptured in an air atmosphere at room temperature. This brittle fracture mode in the absence of high temperature water environment confirmed a high degree of embrittlement of this material resulting from its service exposure to neutron irradiation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Daytime convective development over land: the role of surface forcing

Water availability at the Earth surface determines the portioning of the surface heat flux into its sensible and latent components, that is, the surface heat flux Bowen ratio. The two components affect differently the surface buoyancy flux and thus the development and growth of the convective boundary layer. As a result, the Bowen ratio has a strong impact on the daytime dry and moist convection development over land. Here, we use two canonical modeling test cases, one for the shallow convection and one for the shallow-to-deep convection transition, to document the impact of the surface Bowen ratio on daytime convection development. A simple approach is used where results from simulations featuring the original setup are contrasted with simulations where surface sensible heat flux takes on values of the latent heat flux and vice versa. Such a change illustrates the key impact of the surface water availability without changing the total surface heat flux. Because of the larger surface buoyancy flux, simulations with the reversed surface heat fluxes feature faster deepening of the convective boundary layer and wider clouds once moist convection develops. Mean cloud base width of cumulus clouds increases with the boundary layer depth. For shallow-to-deep convection transition, comparison of simulations with the observed large-scale flow that feature substantial lower-tropospheric shear with simulations with no mean horizontal flow suggests that the organization of convective circulations within the boundary layer strongly affects initial size and organization of shallow convection, as well as the subsequent width increase once convection deepens and eventually transitions to deep precipitating convection. A simple explanation is provided of why deeper well-mixed convective subcloud layer results in wider clouds. The key is the larger width of boundary layer coherent updraft structures when the mixed layer is deeper.

54 ENVIRONMENTAL SCIENCES↗

The FastEddy® Resident-GPU Accelerated Large-Eddy Simulation Framework: Moist Dynamics Extension, Validation and Sensitivities of Modeling Non-Precipitating Shallow Cumulus Clouds

Herein we describe the moist dynamics formulation implemented within the graphics processing unit-resident large-eddy simulation FastEddy® model, which includes a simple saturation adjustment scheme for condensation and evaporation processes. Two LES model intercomparison exercises for non-precipitating shallow cumulus clouds are simulated in order to validate this model extension, including a static forcing and a time-dependent forcing case. Overall, we find our dynamical, thermodynamical and microphysical quantities, along with turbulence variability and fluxes, to be commensurate with the corresponding model intercomparison results. In addition, sensitivities to specific model settings are investigated. Among these settings, it is shown that boundary layer and cloud layer structure and characteristics are sensitive to use of higher-order advection schemes impacting the vertical distribution of cloud content and associated turbulence statistics. Increasing the timescale of the saturation scheme leads to enhanced liquid water presence and decreases vertical velocity variance within the cloud deck. In some cases, these sensitivities agree with the model-to-model variability reported in the intercomparison exercises, highlighting the important role of specific model implementation choices in the context of shallow cumulus convection simulations. These analyses and findings also provide the basis for future extensions and applications of FastEddy® for modeling moist convection and precipitation scenarios.

54 ENVIRONMENTAL SCIENCES↗

Linking Stream Chemistry to Subsurface Redox Architecture

As drinking-water scarcity grows worldwide, we need to improve predictions of the quantity and quality of our water resources. An overarching problem for model improvement is that we do not know the geological structure of aquifers in sufficient detail. In this work, we demonstrate that mineral-water reactions imprint structure in the subsurface that impacts the flow and transport of some chemical species. Specifically, pyrite, a ubiquitous mineral, commonly oxidizes and depletes in the upper layers of the weathering profile in most humid watersheds, only remaining at depths of meters. We hypothesize that variations in concentrations (C) of pyrite-derived sulfate released into rivers as a function of discharge (q) reflect the rate-limiting step and depth of pyrite-oxidizing layers. We found that log C - log q behaviors thus differ in small and large watersheds in the Susquehanna River Basin as well as in selected watersheds in the Western United States. Although coal mining changes pyrite oxidation from closed to open system with respect to O 2 , patterns in stream chemistry as a function of discharge are consistent with deep and shallow pyrite oxidation zones in small and large watersheds respectively. Therefore, understanding the subsurface patterns of mineral reactions and how they affect the architecture of aquifers will elucidate patterns of changing river chemistry and our ability to manage water resources in the future under accelerated land use and climate change.

54 ENVIRONMENTAL SCIENCES↗

MODFLOW6 models used to evaluate potential stresses and hydrologic conditions driving water-level fluctuations in well ER-5-3-2, Frenchman Flat, Southern Nevada

The hydrograph for well ER-5-3-2 in Frenchman Flat, southern Nevada, has previously unexplained water-level fluctuations. Four, three-dimensional, groundwater models (MODFLOW 6) were developed to evaluate potential stresses and hydrologic conditions affecting the well ER-5-3-2 hydrograph. Four model scenarios were developed that simulated: (1) wellbore leakage without recharge, (2) wellbore leakage with recharge, (3) shallow (low transmissivity) and deep (high transmissivity) carbonate rocks, and (4) lateral heterogeneity of carbonate rocks. Input and output files for the four model scenarios are in the model and output directories, respectively. Hydraulic conductivity, specific storage, and wellbore-leakage rates (when simulated) were estimated with parameter estimation (PEST) by minimizing a weighted composite, sum-of-squares objective function. The objective function was informed by measurement and Tikhonov regularization observations. Measurement observations included drawdowns from the constant-rate aquifer test and water-level altitudes measured in well ER-5-3-2 from 2001-2021. Tikhonov regularization informed hydraulic conductivity and specific storage parameters that were insensitive to measurement observations, where homogeneity was the preferred relation. Batch files, executables, and MODFLOW 6, PEST, and post-processing utilities are in the ancillary directory. Supplementary data also are included in the ancillary directory, including site information, high-frequency water-level and aquifer-test data, transmissivity estimates, water-chemistry data, and water-temperature analyses. This USGS data release contains data, analyses, and model files for the simulations and analysis results described in U.S. Geological Survey Scientific Investigations Report (https://doi.org/10.3133/sir20225132).

54 ENVIRONMENTAL SCIENCES↗