Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Climate Zones”

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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Project Title: Adaptive Vertical Grid Enhancement for E3SM (Final Technical Report)

During the first phase of the project (2018/05-2021/05), we proved the hypothesis that, without tuning, high vertical resolution improves the representation of unresolved, parameterized low altitude clouds such as stratocumulus clouds in global climate models. We evaluated the effectiveness of the newly developed computational framework, Framework for Improvement by Vertical Enhancement (FIVE) by implementing it into E3SM. FIVE aims to improve representation of the target cloud, low level cloud for this project, with use of a locally high vertical resolution for computation of selected processes. FIVE was assessed via direct comparison among E3SM, E3SM-FIVE, and satellite data as well as a newly developed ensemble artificial neural network methodology. We demonstrated that E3SM-FIVE remarkably improves marine stratocumulus representation. To advance the FIVE methodology, we took a first step toward the adaptive vertical grid in the FIVE method by modifying the E3SM single column model. For future application of FIVE to the global storm resolving model with a few km horizontal grid resolution, we developed the vertical-meridional two-dimensional (2D) Hadley circulation model, which allows one to run simulations with much higher vertical resolution and significantly longer duration than current three-dimensional global storm resolving model. During the second phase of the project (2021/05-2023/11; period for this technical report), we sought a way to appropriately transform FIVE to a regionalized version of FIVE, which allows one to select locations with/without FIVE. In parallel, we demonstrated that E3SM-FIVE with a 25-km mesh significantly reduced the bias associated with coastal stratocumulus. Furthermore, we showed that the same improvement was able to be obtained with the use of both the Regional Refined Model (RRM) targeted to one of the stratocumulus regions with 25-km resolution and FIVE. We explored a new implementation strategy of FIVE with our 2D Hadley circulation model. In parallel we applied the 2D Hadley circulation model to study the double intertropical convergence zone bias with storm resolving model resolution and aerosol-cloud interactions within the Hadley circulation with near large eddy simulation resolution. Last, we have been investigating Arctic mixed-phase stratocumulus clouds with large eddy simulations in order to improve their representation in global models. During the course of the entire project, 6 journal articles were published, and 3 manuscripts are currently in preparation. We gave a total of 30 presentations at various domestic and international conferences. In addition we gave 6 seminar talks at universities and a national laboratory. E3SM- FIVE is available as a branch of E3SM.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Stable Methane Isotopologues From Northern Lakes Suggest That Ebullition Is Dominated by Sub‐Lake Scale Processes

Abstract Stable isotopes have emerged as popular study targets when investigating emission of methane (CH 4 ) from lakes. Yet little is known on how isotopic patterns conform to variations in emission magnitudes—a highly relevant question. Here, we present a large multiyear data set on stable isotopes of CH 4 ebullition (bubbling) from three small adjacent subarctic lakes. The δ 13 C‐CH 4 and δD‐CH 4 range from −78.4‰ to −53.1‰ and from −369.8‰ to −218.8‰, respectively, and vary greatly among the lakes. The signatures suggest dominant hydrogenotrophic methanogenesis, particularly in the deep zones, but there are also signals of seemingly acetoclastic production in some high fluxing shallow areas, possibly fueled by in situ vegetation, but in‐sediment anaerobic CH 4 oxidation cannot be ruled out as an alternative cause. The observed patterns, however, are not consistent across the lakes. Neither do they correspond to the spatiotemporal variations in the measured bubble CH 4 fluxes. Patterns of acetoclastic and hydrogenotrophic production plus oxidation demonstrate that gains and losses of sediment CH 4 are dominated by sub‐lake scale processes. The δD‐CH 4 in the bubbles was significantly different depending on measurement month, likely due to evaporation effects. On a larger scale, our isotopic data, combined with those from other lakes, show a significant difference in bubble δD‐CH 4 between postglacial and thermokarst lakes, an important result for emission inventories. Although this characteristic theoretically assists in source partitioning studies, most hypothetical future shifts in δD‐CH 4 due to high‐latitude lake area or production pathway are too small to lead to atmospheric changes detectable with current technology.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2007

Microbial Extremophiles in Aspect of Limits of Life

During Earth's evolution accompanied by geophysical and climatic changes a number of ecosystems have been formed. These ecosystems differ by the broad variety of physicochemical and biological factors composing our environment. Traditionally, pH and salinity are considered as geochemical extremes, as opposed to the temperature, pressure and radiation that are referred to as physical extremes (Van den Burg, 2003). Life inhabits all possible places on Earth interacting with the environment and within itself (cross species relations). In nature it is very rare when an ecotope is inhabited by a single species. As a rule, most ecosystems contain the functionally related and evolutionarily adjusted communities (consortia and populations). In contrast to the multicellular structure of eukaryotes (tissues, organs, systems of organs, whole organism), the highest organized form of prokaryotic life in nature is the benthic colonization in biofilms and microbial mats. In these complex structures all microbial cells of different species are distributed in space and time according to their functions and to physicochemical gradients that allow more effective system support, self-protection, and energy distribution. In vitro, of course, the most primitive organized structure for bacterial and archaeal cultures is the colony, the size, shape, color, consistency, and other characteristics of which could carry varies specifics on species or subspecies levels. In table 1 all known types of microbial communities are shown (Pikuta et a]., 2005). In deep underground (lithospheric) and deep-sea ecosystems an additional factor - pressure, and irradiation - could also be included in the list of microbial communities. Currently the beststudied ecosystems are: human body (due to the medical importance), and fresh water and marine ecosystems (due to the reason of an environmental safety). For a long time, extremophiles were terra incognita, since the environments with aggressive parameters (compared to the human body temperature, pH, mineralization, and pressure) were considered a priori as a dead zone.

Pikuta, Elena V.↗
NASA NTRSDate not supplied

The Influence of Sea ice on Arctic Cloud Properties: What Can We Learn by Applying a In Situ Observational Strategy to Satellite Data?

Sea ice is declining because of anthropogenic climate change. This change alters many aspects of the Arctic climate system, including the way that the surface and atmosphere interact. Atmosphere-surface coupling processes represent an important cloud feedback mechanism that can alter the Arctic surface energy budget. For the Arctic, it has been hypothesized that a reduction in sea ice cover could lead to an increase in clouds. If this process were to occur as originally hypothesized, sea ice loss in all seasons would lead to an increase in clouds. Recent observational studies find a cloud response to sea ice loss in non-summer months and no cloud response in summer months. However, previous studies rely on inter-annual variability and reanalysis to control for the influence of meteorology, reducing the confidence in the resulting conclusions. We adopt a phenomenological, event-based approach that does not need to use meteorological reanalysis. The approach analyzes cloud properties derived from CALIPSO-CloudSat over sea ice and adjacent ice-free footprints by compositing individual satellite ground tracks that cross the Arctic sea ice edge. The underlying assumption, which we verify, is that footprints that are close to each other in space and time experience similar large-scale meteorological conditions. Our results show larger cloud fraction and more cloud liquid water over ice-free than over sea ice footprints and provide additional evidence for a seasonal dependence of cloud-sea ice coupling that is in line with previous work. We find a different result where the maximum cloud property differences between sea ice and ice-free ocean occurs in spring, not fall as earlier studies suggest. We argue that these cloud differences between sea ice and ice-free ocean are primarily caused by the influence of the surface type on the thermodynamic stability of the lower troposphere and not principally from an increase in surface evaporation. In addition, we explore the sensitivity of these results to marginal ice zone width, season, Atlantic vs. Pacific sector, and the cloud property dependence on the distance from the sea ice edge. Overall, our results provide further evidence that cloud-sea ice coupling processes are not offsetting the observed surface energy budget perturbation due to sea ice loss in summer; thus, the cloud response to declining sea ice appears to contribute to amplified Arctic warming.

Patrick C Taylor↗
NASA NTRS1994

Size distribution and scattering phase function of aerosol particles retrieved from sky brightness measurements

Ground-based measurements of the solar transmission and sky radiance in a horizontal plane through the Sun are taken in several geographical regions and aerosol types: dust in a desert transition zone in Israel, sulfate particles in Eastern and Western Europe, tropical aerosol in Brazil, and mixed continental/maritime aerosol in California. Stratospheric aerosol was introduced after the eruption of Mount Pinatubo in June 1991. Therefore measurements taken before the eruption are used to analyze the properties of tropospheric aerosol; measurements from 1992 are also used to detect the particle size and concentration of stratospheric aerosol. The measurements are used to retrieve the size distribution and the scattering phase function at large scattering angles of the undisturbed aerosol particles. The retrieved properties represent an average on the entire atmospheric column. A comparison between the retrieved phase function for a scattering angle of 120 deg, with phase function predicted from the retrieved size distribution, is used to test the assumption of particle homogeneity and sphericity in radiative transfer models (Mie theory). The effect was found to be small (20% +/- 15%). For the stratospheric aerosol (sulfates), as expected, the phase function was very well predicted using the Mie theory. A model with a power law distribution, based on the spectral dependence of the optical thickness, alpha, cannot estimate accurately the phase function (up to 50% error for lambda = 0.87 microns). Before the Pinatubo eruption the ratio between the volumes of sulfate and coarse particles was very well correlated with alpha. The Pinatubo stratospheric aerosol destroyed this correlation. The aerosol optical properties are compared with analysis of the size, shape, and composition of the individual particles by electron microscopy of in situ samples. The measured volume size distribution before the injection of stratospheric aerosol consistently show two modes, sulfate particles with r(sub m) less than 0.2 microns and coarse paritcles with r(sub m) greater than 0.7 microns. The 'window' in the tropospheric aerosol in this radius range was used to observe a stable stratospheric aerosol in 1992, with r(sub m) approximately 0.5 microns. A combination of such optical thickness and sky measurements can be used to assess the direct forcing and the climatic impact of aerosol. Systematic inversion for the key aerosol types (sulfates, smoke, dust, and maritime aerosol) of the size distribution and phase function can give the relationship between the aerosol physical and optical properties that can be used to compute the radiative forcing. This forcing can be validated in dedicated field experiments.

Kaufman, Y. J.↗
NASA NTRS1997

Cabauw Experimental Results from the Project for Intercomparison of Land-Surface Parameterization Schemes

In the Project for Intercomparison of Land-Surface Parameterization Schemes phase 2a experiment, meteorological data for the year 1987 from Cabauw, the Netherlands, were used as inputs to 23 land-surface flux schemes designed for use in climate and weather models. Schemes were evaluated by comparing their outputs with long-term measurements of surface sensible heat fluxes into the atmosphere and the ground, and of upward longwave radiation and total net radiative fluxes, and also comparing them with latent heat fluxes derived from a surface energy balance. Tuning of schemes by use of the observed flux data was not permitted. On an annual basis, the predicted surface radiative temperature exhibits a range of 2 K across schemes, consistent with the range of about 10 W/m in predicted surface net radiation. Most modeled values of monthly net radiation differ from the observations by less than the estimated maximum monthly observational error (+/- 10 W/sq m). However, modeled radiative surface temperature appears to have a systematic positive bias in most schemes; this might be explained by an error in assumed emissivity and by models' neglect of canopy thermal heterogeneity. Annual means of sensible and latent heat fluxes, into which net radiation is partitioned, have ranges across schemes of 30 W/sq m and 25 W/sq m, respectively. Annual totals of evapotranspiration and runoff, into which the precipitation is partitioned, both have ranges of 315 mm. These ranges in annual heat and water fluxes were approximately halved upon exclusion of the three schemes that have no stomatal resistance under non-water-stressed conditions. Many schemes tend to underestimate latent heat flux and overestimate sensible heat flux in summer, with a reverse tendency in winter. For six schemes, root-mean-square deviations of predictions from monthly observations are less than the estimated upper bounds on observation errors (5 W/m for sensible beat flux and 10 W/m for latent heat flux). Actual runoff at the site is believed to be dominated by vertical drainage to ground water, but several schemes produced significant amounts of runoff as overland flow or interflow. There is a range across schemes of 184 mm (40% of total pore volume) in the simulated annual mean root-zone soil moisture. Unfortunately, no measurements of soil moisture were available for model evaluation. A theoretical analysis suggested that differences in boundary conditions used in various schemes are not sufficient to explain the large variance in soil moisture. However, many of the extreme values of soil moisture could be explained in terms of the particulars of experimental setup or excessive evapotranspiration.

Chen, Tian Hong↗
DOE OSTI2021

The Relationship Between Topography, Bedrock Weathering, and Water Storage Across a Sequence of Ridges and Valleys

Abstract Bedrock weathering regulates nutrient mobilization, water storage, and soil production. Relative to the mobile soil layer, little is known about the relationship between topography and bedrock weathering. Here, we identify a common pattern of weathering and water storage across a sequence of three ridges and valleys in the sedimentary Great Valley Sequence in Northern California that share a tectonic and climate history. Deep drilling, downhole logging, and characterization of chemistry and porosity reveal two weathering fronts. The shallower front is ∼7 m deep at the ridge of all three hillslopes, and marks the extent of pervasive fracturing and oxidation of pyrite and organic carbon. A deeper weathering front marks the extent of open fractures and discoloration. This front is 11 m deep under two ridges of similar ridge‐valley spacing, but 17.5 m deep under a ridge with nearly twice the ridge‐valley spacing. Hence, at ridge tops, the fraction of the hillslope relief that is weathered scales with hillslope length. In all three hillslopes, below this second weathering front, closed fractures and unweathered bedrock extend about one‐half the hilltop elevation above the adjacent channels. Neutron probe surveys reveal that seasonally dynamic moisture is stored to approximately the same depth as the shallow weathering front. Under the channels that bound our study hillslopes, the two weathering fronts coincide and occur within centimeters of the ground surface. Our findings provide evidence for feedbacks between erosion and weathering in mountainous landscapes that result in systematic subsurface structuring and water routing.

58 GEOSCIENCES↗
NASA NTRS2011

Groundwater and Terrestrial Water Storage

Most people think of groundwater as a resource, but it is also a useful indicator of climate variability and human impacts on the environment. Groundwater storage varies slowly relative to other non-frozen components of the water cycle, encapsulating long period variations and trends in surface meteorology. On seasonal to interannual timescales, groundwater is as dynamic as soil moisture, and it has been shown that groundwater storage changes have contributed to sea level variations. Groundwater monitoring well measurements are too sporadic and poorly assembled outside of the United States and a few other nations to permit direct global assessment of groundwater variability. However, observational estimates of terrestrial water storage (TWS) variations from the GRACE satellites largely represent groundwater storage variations on an interannual basis, save for high latitude/altitude (dominated by snow and ice) and wet tropical (surface water) regions. A figure maps changes in mean annual TWS from 2009 to 2010, based on GRACE, reflecting hydroclimatic conditions in 2010. Severe droughts impacted Russia and the Amazon, and drier than normal weather also affected the Indochinese peninsula, parts of central and southern Africa, and western Australia. Groundwater depletion continued in northern India, while heavy rains in California helped to replenish aquifers that have been depleted by drought and withdrawals for irrigation, though they are still below normal levels. Droughts in northern Argentina and western China similarly abated. Wet weather raised aquifer levels broadly across western Europe. Rains in eastern Australia caused flooding to the north and helped to mitigate a decade long drought in the south. Significant reductions in TWS seen in the coast of Alaska and the Patagonian Andes represent ongoing glacier melt, not groundwater depletion. Figures plot time series of zonal mean and global GRACE derived non-seasonal TWS anomalies (deviation from the mean of each month of the year) excluding Greenland and Antarctica. The two figures show that 2010 was the driest year since 2003. The drought in the Amazon was largely responsible, but an excess of water in 2009 seems to have buffered that drought to some extent. The drying trend in the 25-55 deg S zone is a combination of Patagonian glacier melt and drought in parts of Australia.

Rodell, Matthew↗
DOE OSTI2021

Temporal, Spatial, and Temperature Controls on Organic Carbon Mineralization and Methanogenesis in Arctic High-Centered Polygon Soils

Warming temperatures in continuous permafrost zones of the Arctic will alter both hydrological and geochemical soil conditions, which are strongly linked with heterotrophic microbial carbon (C) cycling. Heterogeneous permafrost landscapes are often dominated by polygonal features formed by expanding ice wedges: water accumulates in low centered polygons (LCPs), and water drains outward to surrounding troughs in high centered polygons (HCPs). These geospatial differences in hydrology cause gradients in biogeochemistry, soil C storage potential, and thermal properties. Presently, data quantifying carbon dioxide (CO 2 ) and methane (CH 4 ) release from HCP soils are needed to support modeling and evaluation of warming-induced CO 2 and CH 4 fluxes from tundra soils. This study quantifies the distribution of microbial CO 2 and CH 4 release in HCPs over a range of temperatures and draws comparisons to previous LCP studies. Arctic tundra soils were initially characterized for geochemical and hydraulic properties. Laboratory incubations at −2, +4, and +8°C were used to quantify temporal trends in CO 2 and CH 4 production from homogenized active layer organic and mineral soils in HCP centers and troughs, and methanogen abundance was estimated from mcrA gene measurements. Results showed that soil water availability, organic C, and redox conditions influence temporal dynamics and magnitude of gas production from HCP active layer soils during warming. At early incubation times (2–9 days), higher CO 2 emissions were observed from HCP trough soils than from HCP center soils, but increased CO 2 production occurred in center soils at later times (>20 days). HCP center soils did not support methanogenesis, but CH 4 -producing trough soils did indicate methanogen presence. Consistent with previous LCP studies, HCP organic soils showed increased CO 2 and CH 4 production with elevated water content, but HCP trough mineral soils produced more CH 4 than LCP mineral soils. HCP mineral soils also released substantial CO 2 but did not show a strong trend in CO 2 and CH 4 release with water content. Knowledge of temporal and spatial variability in microbial C mineralization rates of Arctic soils in response to warming are key to constraining uncertainties in predictive climate models.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Quantifying Subsurface Biogeochemical Variability in a High Altitude Watershed During Winter Isolation

Shallow subsurface microbial and geochemical processes in watersheds are dynamic and responsive to seasonal and long-term environmental change. At high elevation or latitude, such changes may occur during winter months when normal sampling is impossible, impractical, unsafe, or may actually alter the processes being studied. Yet accurate modeling of biogeochemical parameters in these environments requires sampling that can capture events and persistent cold climate trends. Gaps in current models exist where data have not been collected during extended snow- or ice-covered periods. The goals of our exploratory study were to examine the biogeochemical processes occurring in a well-studied high altitude watershed during winter when the system is largely inaccessible. We tested the hypothesis that during snow cover in the East River (ER) watershed, episodic excursions of microbial community structure and biogeochemical processes and concentrations fluctuated from values extrapolated from pre- and post-snow time periods. Through our research we identified previously unresolved subsurface-surface transport phenomena and a heretofore under-appreciated methane cycle occurring during wintertime months, thus providing valuable insight into watershed processes. Our research showed that autonomous sampling techniques, first used in marine systems, is extremely versatile and can be adapted for use in harsh continental settings to examine temporal changes in freshwater chemistry and microbial community assemblage. Further, we found that by sampling during periods when systems are hard to access, will lend insight and better accountability for global budgets for methane, an important greenhouse gas. By strengthening our knowledge of wintertime biogeochemical relationships in high altitude watersheds our research leads to more complete descriptions of cryptic processes and ultimately will inform and improve the reactive transport models that simulate processes in these systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Performance assessment of a real water source heat pump within a hardware-in-the-loop (HIL) testing environment

Over the last decade, the global fight against climate change through electrification has led to an increase in research on building heating, ventilation, and air conditioning (HVAC) systems that utilize intelligent control algorithms to provide demand-side grid service while also maintaining the thermal comfort of building occupants. As the pivotal point between building electricity consumption and indoor thermal comfort, high-efficiency electrical vapor-compression heat pumps are at the center of these emerging studies, and various grid-interactive and occupant-comfort control algorithms have been developed for them. The impact of these algorithms on heat pump operation and performance when subjected to different weather conditions, building loads, and grid requests calls for investigation and verification via experimental testing with actual heat pumps integrated with real-time building and grid responses. This study introduces a Water-Source Heat Pump (WSHP) Hardware-in-The-Loop (HIL) Test Facility that is the first of its kind. This testbed utilizes a 2-ton variable speed water-to-air heat pump that is capable of interacting with a virtual environment currently comprised of an EnergyPlus (E+) building simulation, an agent-based occupant behavioral model, and a single U-tube ground-loop heat exchanger (GLHE) model. Detailed descriptions of the testbed’s physical design and operation, virtual environment, as well as their mutual communication is provided. An uncertainty analysis is also performed under manufacturer specified heating and cooling design conditions. This analysis shows that the total load across the WSHP’s demand side heat exchanger, i.e., the sum of its latent and sensible components, can be measured with a relative uncertainty of ± 10.4% and ± 3.6% in cooling and heating mode respectively. The WSHP’s coefficient of performance (COP) can be measured with relative uncertainties of ± 10.4% in cooling mode, and ± 3.7% in heating mode. A preliminary 24-h experimental demonstration is then performed utilizing the DOE prototype small commercial office building model in E+. The simulation takes place in Atlanta, GA on the date of 08/26/15 from 12:00 AM to 11:59 PM using TMY3 weather data. Here, the results from this demonstration show that over the course of this experiment the simulated outputs of zone dry-bulb temperature, zone humidity ratio, and WSHP inlet water temperature can be tracked by testbed emulators up to a root mean squared error (RMSE) of ± 0.27 °C, ± 0.376 g/kg, and ± 0.85 °C respectively. The WSHP’s dynamic behavioral characteristics and performance are also captured, and correspond well with the authors’ previous understanding of heat pump efficiency as a function of evaporator and condenser fluid inlet conditions respectively.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

CERF: IM3 Projected Western US Power Plant Locations

Overview The Capacity Expansion Regional Feasibility (CERF) model is an open-source geospatial python package that provides new power plant locations at a 1km resolution. The model ingests U.S. state or regional-scale electricity system capacity expansion plans, such as those produced by the Global Change Analysis Model (GCAM-USA), and identifies feasible, site-specific locations for individual new power plants (renewable and non-renewable). CERF combines high-resolution geospatial suitability analyses with an economic algorithm that selects individual plant siting locations based on grid interconnection costs and the locational marginal value of new generation. The model incorporates a wide range of dynamic constraints and opportunities, such as protected lands, population density, existing infrastructure, and water availability. This dataset provides CERF power plant siting results for IM3 Phase 2 simulations across eight different scenarios for the Western US through 2055. The scenarios include combinations of two Shared Socioeconomic Pathways (SSP3 and SSP5) with four high-resolution climate projections specific to the United States (see, https://tgw-data.msdlive.org/). These climate projections include "hotter" and "cooler" variants for two Representative Concentration Pathways (RCP4.5 and RCP8.5). The resulting eight simulations are: rcp45cooler_ssp3 rcp45cooler_ssp5 rcp45hotter_ssp3 rcp45hotter_ssp5 rcp85cooler_ssp3 rcp85cooler_ssp5 rcp85hotter_ssp3 rcp85hotter_ssp5 CERF siting results in this dataset correspond to capacity expansion plans in the GCAM-USA IM3 Phase 2 simulation data and are available for each of the above scenarios. Data Details Temporal Range: 2015-2055 in 5-year timesteps. Note that 2015 is the experiment base year and 2020 and beyond represent model simulation years. Spatial Range: Plant locations are provided for the eleven states in the Western US including Arizona, California, Colorado, Idaho, Montana, New Mexico, Nevada, Oregon, Utah, Washington, and Wyoming. Spatial Resolution: 1 km-squared, provided in x and y coordinates Geospatial Projection: Albers Equal Area Conic (ESRI:102003) File Type: csv The dataset contains subdirectories for each of the eight scenarios described in the overview. Each scenario folder contains two subfolders with the following information: 1. Power Plant Data This directory contains a single .csv file of power plant locations for both pre-existing (non-CERF sited plants in operation in 2015) and new (CERF-sited) power plants across the temporal range along with additional CERF model output parameters for CERF-sited plants. Plant with a siting year earlier than 2020 correspond to facilities that are operational leading into the first timestep CERF simulation. For a more detailed description of CERF model output parameters, see the CERF model documentation. Note that the cerf_plant_id parameter is unique within each scenario file but not across scenario files. Parameter Descriptions scenario - Name of scenario cerf_plant_id - Unique siting identifier cerf_sited - If True, indicates that plant was sited by CERF model. If False, indicates pre-existing facility region_name - Name of region (state) tech_id - Technology ID tech_name - Full generation technology name inclusive of cooling type (if applicable) and additional characteristics tech_simple - Simplified generation technology type unit_size_mw - Power plant unit size (MW) xcoord - X coordinate in the default CRS (meters) ycoord - Y coordinate in the default CRS (meters) index - Index position in the flattend 2D array buffer_in_km - Exclusion buffer around site (km) sited_year - Year of siting retirement_year - Year of retirement lmp_zone - Locational marginal price (LMP) zone ID locational_marginal_price_usd_per_mwh - Locational marginal price ($/MWh) generation_mwh_per_year - Generation output (MWh/yr) operating_cost_usd_per_year - Cost of plant operations ($/yr) net_operational_value - Net operational value based on LMP and and operating costs ($/yr) interconnection_cost - Cost of interconnection for transmission & gas pipeline (if applicable) net_locational_cost -- Difference of interconnection cost and operating value ($/yr) capacity_factor_fraction - Capacity factor (fraction) carbon_capture_rate_fraction - Carbon capture rate (fraction) fuel_co2_content_tons_per_btu - Fuel CO2 content (tons/Btu) fuel_price_usd_per_mmbtu - Fuel price ($/MMBtu) fuel_price_esc_rate_fraction - Fuel price escalation rate (fraction) heat_rate_btu_per_kWh - Heat rate (Btu/kWh) lifetime_yrs - Technology lifetime for annuity (years) operational_life_yrs - Operational lifetime for retirement (years) variable_om_usd_per_mwh - Variable operation and maintenance costs of yearly capacity use ($/MWh) variable_om_esc_rate_fraction - Variable operation and maintenance costs escalation rate (fraction) carbon_tax_usd_per_ton - Carbon tax ($/ton) carbon_tax_esc_rate_fraction - Carbon tax escalation rate (fraction) 2. Storage Data This directory contains information on new and pre-existing energy storage facilities operational in each timestep along with various storage operational parameters. The 2015 timestep provides pre-existing energy storage data and corresponds with facilities that are operational leading into the first model simulation timestep. Note that coordinates in the storage files correspond to the interconnection point on the grid (substation location), not individual energy storage locations. Energy storage is added in a cumulative process at each given interconnection point. That is, each individual file provides the total operational storage capacity interconnected to the specified substation for the given timestep, inclusive of previously installed storage at that location and new storage installed in that timestep at that location. Parameters scenario - Name of scenario timestep - Simulation timestep name - Unique storage identifier s_typ - Type of energy storage technology (battery or pumped storage hydro) s_node - Node ID of interconnecting substation xcoord - X coordinate in the default CRS (meters) ycoord - Y coordinate in the default CRS (meters) charge_rate - Maximum charge rate (power capacity) of storage system (MW) discharge_rate - Maximum discharge rate (power capacity) of storage system (MW) duration - Duration of storage system (hours) max_SoC - Allowed maximum state of charge (energy capacity) of storage system (MWh) min_SoC -Allowed minimum state of charge (energy capacity) of storage system (MWh) charge_eff - Efficiency of charge (fraction between 0 and 1) discharge_eff - Efficiency of discharge (fraction between 0 and 1) Acknowledgment IM3 is a multi-institutional effort led by Pacific Northwest National Laboratory and supported by the U.S. Department of Energy's Office of Science as part of research in MultiSector Dynamics, Earth and Environmental Systems Modeling Program.

CERF↗
DOE OSTI2025

A Scale‐Adaptive Urban Hydrologic Framework: Incorporating Network‐Level Storm Drainage Pipes Representation

Abstract Below‐ground urban stormwater networks (BUSNs) significantly influence urban flood dynamics, yet their representation at the watershed or larger scales remains challenging. We introduce a scalable urban hydrologic framework that centers on a novel network‐level BUSN representation, balancing the needs for physical basis, parameter parsimony, and computational efficiency. Our framework conceptualizes an urban watershed into four interacting zones: hillslopes (natural), storm‐sewersheds (urban), a sub‐network channel (tributaries), and a main channel. We develop an innovative Graph Theory‐based algorithm to derive network‐level BUSN parameters from publicly available datasets, enabling efficient, scalable parameterization. We demonstrate this framework's applicability at nine representative watersheds in the Houston metropolitan region, USA, with urban imperviousness ranging from 0% to 64% and drainage areas ranging from 24 to 302 . Our model achieves satisfying computational efficiency, completing hourly time step simulations for 18 years in less than 5 sec per watershed on a standard PC. Validation against observed daily streamflow confirms that the model can capture small‐to‐large flood peaks and seasonal and annual water balance over these watersheds. Comparisons with the National Water Model show better performance in predicting flood peaks and overall water balance, underscoring the promises of our new framework for urban hydrologic modeling at large scales. Furthermore, analysis reveals nonlinear relationships between BUSNs' designed capacities and flood reduction effects. Our approach bridges the gap between detailed hydraulic and large‐scale hydrologic models, providing a valuable tool for urban flood prediction and management across broader spatial and temporal scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Forced and Unforced Decadal Behavior of the Interhemispheric SST Contrast during the Instrumental Period (1881–2012): Contextualizing the Late 1960s–Early 1970s Shift

The sea surface temperature (SST) contrast between the Northern Hemisphere (NH) and Southern Hemisphere (SH) influences the location of the intertropical convergence zone (ITCZ) and the intensity of the monsoon systems. This study examines the contributions of external forcing and unforced internal variability to the interhemispheric SST contrast in HadSST3 and ERSSTv5 observations, and 10 models from phase 5 of the Coupled Model Intercomparison Project (CMIP5) from 1881 to 2012. Using multimodel mean fingerprints, a significant influence of anthropogenic, but not natural, forcing is detected in the interhemispheric SST contrast, with the observed response larger than that of the model mean in ERSSTv5. The forced response consists of asymmetric NH–SH SST cooling from the mid-twentieth century to around 1980, followed by opposite NH–SH SST warming. The remaining best-estimate residual or unforced component is marked by NH–SH SST maxima in the 1930s and mid-1960s, and a rapid NH–SH SST decrease around 1970. Examination of decadal shifts in the observed interhemispheric SST contrast highlights the shift around 1970 as the most prominent from 1881 to 2012. Both NH and SH SST variability contributed to the shift, which appears not to be attributable to external forcings. Most models examined fail to capture such large-magnitude shifts in their control simulations, although some models with high interhemispheric SST variability are able to produce them. Large-magnitude shifts produced by the control simulations feature disparate spatial SST patterns, some of which are consistent with changes typically associated with the Atlantic meridional overturning circulation (AMOC).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Thermodynamics of Minerals Stable near Earth's Surface

Project “Thermodynamics of Minerals Stable Near the Earth’s Surface” uses specialized calorimetric techniques to obtain new thermochemical data for solid phases of interest in the Earth’s “critical zone”, the near surface regions of the planet most directly affecting and affected by human activity. Major activity was in three areas: manganese oxide materials, aluminum clusters which are precursors to aluminum oxyhydroxide precipitates, and zeolites. All these classes of minerals are potential carriers of heavy metals and radionuclides in the environment, both within their crystal structures and as adsorbed layers or colloids. Their thermodynamic stability, as well as the sequence of stability of competing phases, governs the aqueous solubility of contaminants (and nutrients) in soil, groundwater, and sediments. The thermodynamics of these phases is strongly affected by particle size and crystallinity, since they often occur as very fine particles which do not coarsen in the near-surface environment. Many of these materials exhibit extensive polymorphism, multicomponent solid solution, and variable hydration states. Thus, one must define and systematize important thermodynamic variables, since one cannot make measurements on each and every possible phase, particle size, and composition. A major finding of the work on manganese oxides was labile water-induced oxidation-reduction reactions. These will be studied further by a combination of thermodynamic and synchrotron-based spectroscopic methods, the latter in collaboration with Peggy O’Day at UC Merced. The adsorption of other species, including simple organics, on such redox-active surfaces will also be investigated by a combination of thermochemical and spectroscopic studies. Studies of the energetics of water adsorption on mineral surfaces will be expanded to emphasize the adsorption of CO2 and small organic molecules, relevant to soil, groundwater, and the subsurface environment. Gas adsorption calorimetry, immersion calorimetry, and solution calorimetry will be combined with appropriate analytical and spectroscopic techniques to obtain a comprehensive picture of the interactions and their energetics. The evolution of aqueous aluminum clusters to form aluminum oxyhydroxide precipitates will be studied in detail by a combination of calorimetric and structural methods applied to aqueous samples aged at room temperature. The energetics of new classes of clusters will be explored. The project addresses several common themes relevant to the broad mission of the Department of Energy. Thermochemical data for minerals stable near the Earth’s surface, though still incomplete, are essential for understanding both natural and anthropogenic processes that affect our very existence on this planet. Such data contribute to the fundamental underpinning for mitigation and management strategies for the extraction and use of energy. Fundamental thermodynamic and mechanistic understanding is required to model processes ranging from heavy metal pollution to nuclear waste disposal to climate change.

58 GEOSCIENCES↗
DOE OSTI2024

Aircraft Measurements from a U.S. Western Wildfire Demonstrating Day and Night Differences in the Chemical Composition and Optical Properties of Biomass Burning Aerosols and Their Atmospheric Evolution

The composition and transformations of biomass burning aerosols (BBA) have been measured onboard the NOAA Twin Otter research aircraft during the Fire Influence on Regional to Global Environments and Air Quality field study. Here, we analyze real-time aerosol mass spectrometry measurements across three flights during the afternoon, late afternoon, and night of August 28, 2019, for one midsized wildfire. Analysis of several metrics showed that the aerosol composition and optical properties varied depending on the burning conditions at the fire zone and the time of day the BBA was emitted, with substantial variations in the available sunlight. The total aerosol mass loadings were dominated by organic components with a much smaller contribution from inorganic species. A gradual buildup of organic material was observed during the afternoon as the plume aged, indicating the condensation of photochemically formed low-volatility oxidized organic compounds. Highly hygroscopic ammonium nitrate was the main inorganic component, suggesting potential water content in BBA particles and the likelihood of their aqueous-phase reactivity. Depletions of particle-phase NO 3 – and Cl – relative to carbon monoxide were observed in the late afternoon and nighttime plumes, respectively, aligning with known gas-particle partitioning thermodynamics and the heterogeneous chemistry of dissolved nitrate and chloride. The wavelength-dependent light absorption by aerosol species was higher for the plume sampled at night and showed no significant changes with plume age, despite observed trends in composition and mass downwind. These differences in particle composition and optical properties demonstrate that the processes involved in BBA aging are not uniform for the same wildfire over the course of the day and depend highly on when the BBA was emitted, as well as the burning phase at the emissions source.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Identifying recharge sources and their impacts on a North Central New Mexico shallow aquifer using unsupervised machine learning

In this article, shallow aquifers are important but highly variable resources in arid to semi-arid regions. Limited shallow aquifer volume results in high sensitivity to recharge fluctuations, which can impact the local fauna and flora, and transport of contaminants in the aquifer or vadose zone. Aquifer response to external forcing (e.g., precipitation) is usually solved by estimating aquifer parameters and running physics-based models to match known fluctuations of hydraulic head. However, this technique is time and computationally expensive. Furthermore, high aquifer complexity decreases precision in physics-based models. Alternatively supervised machine learning is used to predict aquifer dynamics. However, these techniques rely on input data and struggle to interpret aquifer response for missing sources (i.e., snowpack data). To counter these problems, we propose an unsupervised machine learning technique (NMFk) to estimate the impact of different sources on aquifer recharge. NMFk is used to understand the influence of external forcing on shallow aquifer recharge in the Pajarito Plateau (Los Alamos, NM, USA). The results show how NMFk can be used to reduce the data dimension in a complex field dataset to three recharge signals that cause fluctuations within the field data. Here, the source signals are interpreted as rainfall, snowmelt, and a delayed aquifer response to the previous two signals. These results evidence how heterogeneous aquifers delimited by canyons incised into the Pajarito Plateau respond in similar ways to the source signals identified by NMFk. Furthermore, results show the importance of the local geology where faults act as sinks, and anthropogenic disturbances can facilitate infiltration amplifying the interpreted signal.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2004

Annual Cycles of Multiyear Sea Ice Coverage of the Arctic Ocean: 1999-2003

For the years 1999-2003, we estimate the time-varying perennial ice zone (PIZ) coverage and construct the annual cycles of multiyear (MY, including second year) ice coverage of the Arctic Ocean using QuikSCAT backscatter, MY fractions from RADARSAT, and the record of ice export from satellite passive microwave observations. An area balance approach extends the winter MY coverage from QuikSCAT to the remainder of the year. From these estimates, the coverage of MY ice at the beginning of each year is 3774 x 10(exp 3) sq km (2000), 3896 x 10(exp 3) sq km (2001), 4475 x 10(exp 3) sq km (2002), and 4122 x 10(exp 3) sq km (2003). Uncertainties in coverage are approx.150 x 10(exp 3) sq km. In the mean, on 1 January, MY ice covers approx.60% of the Arctic Ocean. Ice export reduces this coverage to approx.55% by 1 May. From the multiple annual cycles, the area of first-year (FY) ice that survives the intervening summers are 1192 x 10(exp 3) sq km (2000), 1509 x 10(exp 3) sq km (2001), and 582 x 10(exp 3) sq km (2002). In order for the MY coverage to remain constant from year to year, these replenishment areas must balance the overall area export and melt during the summer. The effect of the record minimum in Arctic sea ice area during the summer of 2002 is seen in the lowest area of surviving FY ice of the three summers. In addition to the spatial coverage, the location of the PIZ is important. One consequence of the unusual location of the PIZ at the end of the summer of 2002 is the preconditioning for enhanced export of MY ice into the Barents and Kara seas. Differences between the minimums in summer sea ice coverage from our estimates and passive microwave observations are discussed.

Arctic Ocean↗