Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mathematics and Computing, Environmental sciences”

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 325 records · Page 18

Approximate inverse-based block preconditioners in poroelasticity

We focus on the fully implicit solution of the linear systems arising from a three-field mixed finite element approximation of Biot’s poroleasticity equations. The objective is to develop algebraic block preconditioners for the efficient solution of such systems by Krylov subspace methods. In this work, we investigate the use of approximate inverse-based techniques to decouple the native system of equations and obtain explicit sparse approximations of the Schur complements related to the physics-based partitioning of the unknowns by field type. Here, the proposed methods are tested in various numerical experiments including real-world applications dealing with petroleum and geotechnical engineering.

54 ENVIRONMENTAL SCIENCES↗

Conservative Numerical Schemes with Optimal Dispersive Wave Relations: Part II. Numerical Evaluations

A new energy and enstrophy conserving scheme (EEC) for the shallow water equations is proposed and evaluated using a suite of test cases over the global spherical or bounded domain. The evaluation is organized around a set of pre-defined properties: accuracy of individual operators, accuracy of the whole scheme, conservation of key quantities, control of the divergence variable, representation of the energy and enstrophy spectra, and simulation of nonlinear dynamics. The results confirm that the scheme is between the first and second order accurate, and conserves the total energy and potential enstrophy up to the time truncation errors. Here, the scheme is capable of producing more physically realistic energy and enstrophy spectra, indicating that it can help prevent the unphysical energy cascade towards the finest resolvable scales. With an optimal representation of the dispersive wave relations, the scheme is able to keep the flow close to being non-divergent, and maintain the geostrophically balanced structures with large-scale geophysical flows over long-term simulations.

54 ENVIRONMENTAL SCIENCES↗

NRAP-open-IAM: A flexible open-source integrated-assessment-model for geologic carbon storage risk assessment and management

Large-scale implementation of geologic carbon storage (GCS) to help reduce atmospheric greenhouse gas emissions requires stakeholder confidence that injected CO2 will remain contained and that potential subsurface environmental risks are acceptably small and manageable. The U.S. Department of Energy’s National Risk Assessment Partnership (NRAP) has developed an open-source integrated assessment model (NRAP-Open-IAM) to help address questions about a potential GCS site’s ability to effectively contain injected CO 2 and protect groundwater and other overlying environmentally sensitive receptors. NRAP-Open-IAM allows a user to: (1) incorporate relevant site geologic and injection scenario data; (2) characterize important site features and events;(3) couple fast prediction models of various system components of the engineered geologic system; and (4) execute stochastic, dynamic simulation of whole GCS system performance, leakage risk assessment, and uncertainty quantification. NRAP-Open-IAM is available on GitLab (https://gitlab.com/NRAP/OpenIAM), and is accompanied by multiple application examples and detailed user and developer guides.

54 ENVIRONMENTAL SCIENCES↗

Toward more mechanistic representations of biogeochemical processes in river networks: Implementation and demonstration of a multiscale model

A recently introduced multiscale model for representing the effects of small-scale hyporheic-zone biogeochemical processes is extended from the reach scale to river network scales. The model uses advection-dispersion-reaction equations for the channel network and one-dimensional advection-reaction subgrid models for the hyporheic zone. We summarize implementation in the integrated surface/subsurface hydrology modeling system Amanzi-ATS. The extension uses topologically defined meshes to represent stream/river networks and associates a hyporheic subgrid model with each channel grid cell. Biogeochemical reaction modeling capability residing in community software is accessed through an application programming interface. The implementation is verified against independent numerical solutions on a single reach. Mesh convergence studies show that commonly used semi-distributed representations can introduce significant spatial discretization error. Denitrification of farm runoff in a subbasin of the Portage River Basin in Ohio, USA is used to demonstrate the general-purpose reactive transport capability.

54 ENVIRONMENTAL SCIENCES↗

Machine learning to discover mineral trapping signatures due to CO 2 injection

Mineral trapping is pursued as a geological CO 2 sequestration (GCS) mechanism because it permanently stores CO 2 in solid phases or minerals. However, CO 2 mineral-trapping mechanisms are poorly understood due to (1) lack of sufficient field and laboratory data characterizing these complex processes, and (2) challenges to develop site-specific reactive-transport models coupling fluid flow and geochemical reactions occurring at various temporal (from milliseconds to years) and spatial (from pore (millimeters) to field (kilometers)) scales. Reactive transport with additional complexities such as heterogeneity can make the simulation outputs even more difficult to interpret because of complex nonlinearity and multi-scale interdependencies. Furthermore, the values of model outputs such as concentrations can vary by several orders of magnitude, making it harder to correlate and characterize the impact of the variables via traditional data interpretation techniques such as exploratory data analyses. Recently, machine learning (ML) has shown promise in feature discovery and in highlighting hidden mechanisms that cannot be obtained by existing data-analytics and statistical methods. In this study, we applied an unsupervised ML approach, non-negative matrix factorization with custom -means clustering (NMF) to the data generated by reactive-transport simulations of GCS. The reactive-transport data consisted of 19 attributes, including four physio-chemical variables (pH, porosity, aqueous CO 2 , and sequestered CO 2 ), six chemical species (K + , Na + , HCO, Ca 2+ , Mg 2+ , Fe 2+ ), and four carbonate minerals (calcite, dolomite, siderite, and ankerite), a feldspar mineral (albite), and four clay minerals (illite, clinochlore, kaolinite, and smectite) over a period of 200 years of simulation time. Furthermore, the simulation data used was for Morrow B sandstone at the Farnsworth hydrocarbon unit in Texas. Data are sampled at two locations within the model domain: (1) at the injection well and (2) 200 m west of the injection well. The injection was performed for a period of 10 years. Using NMF, we estimated the temporal interdependencies among the 19 attributes over a span of 200 years. We found that NMF was able to identify four reaction stages and their dominant attributes; these cannot be directly discerned through traditional visualization (e.g., line plots, Pareto analysis, Glyph-based visualization methods) or exploratory data analysis tools of the simulation data. The four stages were: reactions in the injection phase followed by short-, mid-, and long-term reactions. The NMF analysis also revealed that 10 among the 19 attributes are dominant. These dominant attributes for mineral trapping include calcite, dolomite at injection well, siderite at 200 m away from the injection well, clinochlore, kaolinite, Na + , K + , Ca 2+ , Mg 2+ , pH, and aqeuous CO 2 . Finally, at late times (65–200 years), our results showed that calcite plays a major role in mineral trapping with insignificant contribution from siderite, ankerite, and clay minerals. These findings make the proposed unsupervised ML-model attractive for reactive-transport sensing towards real-time GCS monitoring.

54 ENVIRONMENTAL SCIENCES↗

Baroclinic Residual Circulation and Mass Transport Due to Internal Tides

Baroclinic (BC) tidal residual circulation due to internal tides is investigated around islands over a shallow ridge using a numerical ocean model. Internal tides enhance vertical mixing over shallow slopes, leading to horizontal density gradients that drive BC residual circulation along the main thermocline. For a strongly stratified summer case, the vertical diffusivity estimated by the Mellor and Yamada turbulence closure model exceeds 1 × 10 -2 m 2 s -1 , and the velocity of BC residual circulations reaches 0.2 m s -1 . The magnitude of BC residual circulation is larger than that of barotropic residual tidal circulation, implying that BC residual circulation due to internal tides plays an important role in coastal ocean circulation. Furthermore, BC residual flow accounts for an equal percentage (5%) of the total tidal kinetic energy as the barotropic residual flow under summer stratification conditions. Results from a coupled sediment resuspension and transport model show the growth of intermediate nepheloid layers formed by strong bottom shear stress and BC residual flow. The residual component contributes largely (23% of the total in summer) to the total suspended sediment flux. Seasonal variability is explored, with weaker winter stratification leading to reduced mixing, and thus weaker BC residual circulation and sediment flux. Lastly, the magnitude of the BC residual circulation is also shown to be proportional to the square of the tidal amplitude.

54 ENVIRONMENTAL SCIENCES↗

Improving Seasonal Forecast Using Probabilistic Deep Learning

The path toward realizing the potential of seasonal forecasting and its socioeconomic benefits relies on improving general circulation model (GCM) based dynamical forecast systems. To improve dynamical seasonal forecasts, it is crucial to set up forecast benchmarks, and clarify forecast limitations posed by model initialization errors, formulation deficiencies, and internal climate variability. With huge costs in generating large forecast ensembles, and limited observations for forecast verification, the seasonal forecast benchmarking and diagnosing task proves challenging. Here, we develop a probabilistic deep learning-based statistical forecast methodology, drawing on a wealth of climate simulations to enhance seasonal forecast capability and forecast diagnosis. By explicitly modeling the internal climate variability and GCM formulation differences, the proposed Conditional Generative Forecasting (CGF) methodology enables bypassing crucial barriers in dynamical forecast, and offers a top-down viewpoint to examine how complicated GCMs encode the seasonal predictability information. We apply the CGF methodology for global seasonal forecast of precipitation and 2 m air temperature, based on a unique data set consisting 52,201 years of climate simulation. Results show that the CGF methodology can faithfully represent the seasonal predictability information encoded in GCMs. We successfully apply this learned relationship in real-world seasonal forecast, achieving competitive performance compared to dynamical forecasts. Using this CGF as benchmark, we reveal the impact of insufficient forecast spread sampling that limits the skill of the considered dynamical forecast system. Finally, we introduce different strategies for composing ensembles using the CGF methodology, highlighting the potential for leveraging the strengths of multiple GCMs to achieve advantgeous seasonal forecast.

54 ENVIRONMENTAL SCIENCES↗

Simulating the Transport and Rupture of Pollen in the Atmosphere

Pollen, one type of primary biological aerosol particle (PBAP), is emitted from the terrestrial biosphere and can undergo physical changes in the atmosphere via particle rupture. To examine the fate of pollen and its atmospheric processing, a pollen emission and transport scheme is coupled to the Weather Research and Forecasting Model with Chemistry (WRF-Chem). We simulate the emission of pollen and its impacts on the cloud properties and precipitation in the Southern Great Plains from 12 to 19 April 2013, a period with both high pollen emissions and convective activity. We conduct a suite of ensemble runs that simulate primary pollen and three different pollen rupture mechanisms that generate subpollen particles, including (a) high humidity-induced surface rupture, (b) high humidity-induced in-atmosphere plus surface rupture, and (c) lightning-induced rupture, where in-cloud and cloud-to-ground lightning strikes trigger pollen rupture events. When relative humidity is high (>80%), coarse primary pollen (~1 μg m –3 ) is converted into fine subpollen particles (~1.2e –4 μg m –3 ), which produces 80% more subpollen particles than lightning-induced rupture. The in-atmosphere humidity-driven rupture predominantly produces subpollen particles, which is further enhanced during a frontal thunderstorm. During strong convection, vertical updrafts lift primary pollen and subpollen particles (~0.5e –4 μg m –3 ) to the upper troposphere (~12 km) and laterally transports the ruptured pollen in the anvil top outflow. In regions of high pollen and strong convection, ruptured pollen can influence warm cloud formation by decreasing low cloud (<4 km) cloud water mixing ratios and increasing ice phase hydrometeors aloft (>10 km)

54 ENVIRONMENTAL SCIENCES↗

The Influence of Climate Variability and Future Climate Change on Atlantic Hurricane Season Length

Abstract Atlantic hurricane season length is important for emergency management preparation, motivating the need to understand its variability and change. We investigated the influence of ocean variability on Atlantic hurricane season length in observations and a future climate simulated by the Energy Exascale Earth System Model (E3SM). We found that multiple factors influence hurricane season length, through their influence on season start and end. Warm western subtropical Atlantic sea‐surface temperature anomalies (SSTAs) during boreal spring (before the official hurricane season start) drive early starts to the hurricane season, and vice versa for cool SSTAs. Meanwhile, La Niña in autumn (before the official hurricane season end) drives late ends to the hurricane season, and vice versa for El Niño. E3SM projects a 27‐day increase in future Atlantic hurricane season length given La Niña and warm northern tropical Atlantic SSTAs. This research documents sources of predictability for Atlantic hurricane season length.

54 ENVIRONMENTAL SCIENCES↗

Effects of Pollen on Hydrometeors and Precipitation in a Convective System

Anemophilous (wind-driven) pollen is one type of primary biological aerosol particle, which can rupture under high humidity conditions and form smaller sub-pollen particles (SPPs). Both pollen and SPPs can reach the upper troposphere under convective conditions, acting as cloud condensation nuclei (CCN) and ice nucleating particles (INPs), thus influencing cloud formation and precipitation. However, the impacts of these biological aerosols on cold cloud formation and local climate remain unclear as there are large uncertainties on their emission flux and ice nucleating abilities. Here, we incorporate pollen emission and rupture processes in the Weather Research and Forecasting Model with Chemistry (WRF-Chem) simulations and update the Morrison microphysics scheme within WRF-Chem using aerosol-aware INP parameterizations to account for pollen in addition to other anthropogenic and biogenic aerosol. INP parameterizations for pollen and SPP are derived from laboratory experiments. When including pollen rupture rates as observed in a series of chamber studies, SPP concentrations increase, leading to an increase of cloud ice and water by up to 50% and potentially extending the duration of the convective system. Among all simulated hydrometeors, graupel and raindrops exhibit the largest enhancements from the inclusion of SPPs, with intensifying precipitation at the backside of the convective system and a greater spatial extent. Sensitivity simulations indicate that SPPs have a greater effect on cloud microphysical processes than whole pollen grains, and further observational evidence is needed to constrain these processes.

54 ENVIRONMENTAL SCIENCES↗

Global energy spectrum of the general oceanic circulation

Abstract Advent of satellite altimetry brought into focus the pervasiveness of mesoscale eddies $${{{{{{{\bf{{{{{{{{\mathcal{O}}}}}}}}}}}}}}}}({100})$$ O ( 100 ) km in size, which are the ocean’s analogue of weather systems and are often regarded as the spectral peak of kinetic energy (KE). Yet, understanding of the ocean’s spatial scales has been derived mostly from Fourier analysis in small "representative” regions that cannot capture the vast dynamic range at planetary scales. Here, we use a coarse-graining method to analyze scales much larger than what had been possible before. Spectra spanning over three decades of length-scales reveal the Antarctic Circumpolar Current as the spectral peak of the global extra-tropical circulation, at ≈ 10 4 km, and a previously unobserved power-law scaling over scales larger than 10 3 km. A smaller spectral peak exists at ≈ 300 km associated with mesoscales, which, due to their wider spread in wavenumber space, account for more than 50% of resolved surface KE globally. Seasonal cycles of length-scales exhibit a characteristic lag-time of ≈ 40 days per octave of length-scales such that in both hemispheres, KE at 10 2 km peaks in spring while KE at 10 3 km peaks in late summer. These results provide a new window for understanding the multiscale oceanic circulation within Earth’s climate system, including the largest planetary scales.

54 ENVIRONMENTAL SCIENCES↗

Differential credibility assessment for statistical downscaling

Climate science is increasingly using (i) ensembles of climate projections from multiple models derived using different assumptions and/or scenarios and (ii) process-oriented diagnostics of model fidelity. Efforts to assign differential credibility to projections and/or models are also rapidly advancing. A framework to quantify and depict the credibility of statistically downscaled model output is presented and demonstrated. Here, the approach employs transfer functions in the form of robust and resilient generalized linear models applied to downscale daily minimum and maximum temperature anomalies at 10 locations using predictors drawn from ERA-Interim reanalysis and two global climate models (GCM; GFDL-ESM2M and MPI-ESM-LR). The downscaled time series are used to derive several impact relevant CLIMDEX temperature indices that are assigned credibility based on (1) the reproduction of relevant large-scale predictors by the GCMs (i.e. fraction of regression beta-weights derived from predictors that are well-reproduced) and (2) the degree of variance in the observations reproduced in the downscaled series following application of a new variance inflation technique. Credibility of the downscaled predictands varies across locations, between the two GCM and is generally higher for minimum temperature than maximum temperature. The differential credibility assessment framework demonstrated here is easy to use and flexible. It can be applied as is to inform decision makers regarding projection confidence, and/or extended to include other components of the transfer functions, and/or used to weight members of a statistically downscaled ensemble.

54 ENVIRONMENTAL SCIENCES↗

CMIP6-based Multi-model Hydroclimate Projection over the Conterminous US, Version 1.1

This dataset presents a suite of high-resolution downscaled hydro-climate projections over the conterminous United States (CONUS) based on multiple selected Global Climate Models (GCMs) from the Coupled Models Intercomparison Project phase 6 (CMIP6). The CMIP6 GCMs are downscaled using either statistical (DBCCA) or dynamical (RegCM) downscaling approaches based on two meteorological reference datasets (Daymet and Livneh). Subsequently, the downscaled precipitation, temperature, and wind speed are employed to drive two calibrated hydrologic models (VIC and PRMS), enabling the simulation of projected future hydrologic responses across the CONUS. Each ensemble member covers the 1980-2019 baseline and 2020-2059 near-future periods under the high-end (SSP585) emission scenario. Moreover, utilizing only DBCCA and Daymet, the projections are further extended to the 2060-2099 far-future period and across three additional emission scenarios (SSP370, SSP245, and SSP126). This dataset is formulated to support the SECURE Water Act Section 9505 Assessment for the US Department of Energy (DOE) Water Power Technologies Office (WPTO). For further details on this dataset, please refer to Kao et al. (2022) and Rastogi et al. (2022).

13 HYDRO ENERGY↗

Gridded Sub-daily Climate Forcings for North America Based on Daymet and GSWP3 (Daymet-GSWP3)

To support high spatial and temporal resolution land surface modeling, this dataset provides 3-hourly time step historic weather forcing at 1-km spatial resolution for the entire North America. The latest Daymet V4 data provides gridded historic daily weather observations at 1-km spatial resolution from 1980 to 2014. Using sub-daily temporal information from the Global Soil Wetness Project Phase 3 (GSWP3), Daymet was further temporally downscaled to 3-hourly time steps and provided in the format required for land surface model simulations. The process of temporal downscaling preserves the relative magnitude in each sub-daily time step from GSWP3 while maintaining the total and average values from Daymet for each day. This results in a blended 1980-2014 Daymet-GSWP3 dataset. Available variables include surface air temperature, precipitation, specific humidity, shortwave and longwave radiation, wind speed, and pressure. These data can be used as a high-resolution meteorological forcing dataset to support high-resolution land surface modeling where accurate meteorological forcing datasets built from historic observations and/or reanalysis datasets are desirable.

54 ENVIRONMENTAL SCIENCES↗

Gridded Sub-daily Climate Forcings for North America Based on Daymet and GSWP3 (Daymet-GSWP3)

To support high spatial and temporal resolution land surface modeling, this dataset provides 3-hourly time step historic weather forcing at 1-km spatial resolution for the entire North America. The latest Daymet V4 data provides gridded historic daily weather observations at 1-km spatial resolution from 1980 to 2014. Using sub-daily temporal information from the Global Soil Wetness Project Phase 3 (GSWP3), Daymet was further temporally downscaled to 3-hourly time steps and provided in the format required for land surface model simulations. The process of temporal downscaling preserves the relative magnitude in each sub-daily time step from GSWP3 while maintaining the total and average values from Daymet for each day. This results in a blended 1980-2014 Daymet-GSWP3 dataset. Available variables include surface air temperature, precipitation, specific humidity, shortwave and longwave radiation, wind speed, and pressure. These data can be used as a high-resolution meteorological forcing dataset to support high-resolution land surface modeling where accurate meteorological forcing datasets built from historic observations and/or reanalysis datasets are desirable.

54 ENVIRONMENTAL SCIENCES↗

High-Resolution Fire Weather Index Data for the Conterminous US (1980–2099), Version 1

This dataset presents a suite of high-resolution fire weather index datasets calculated from observation (gridMet, Livneh, Daymet V4), reanalysis (AgERA5), downscaled hydro-climate projections over the conterminous United States (CONUS) based on multiple selected Global Climate Models (GCMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6). Aside from the daily FWI datasets, we also include a set of FWI extreme indicators at annual, seasonal, and monthly scales, including 1) fwixx: maximum FWI; 2) fwisa: mean FWI. For annual fwisa, it refers to the season with the maximum seasonal average; 3) fwils: Length of fire season over a specified period, where fire season is defined as the days exceeding the median value of the normalized FWI during the reference period (1980-1984); 4) fwixd: Number of extreme fire weather days over a specified period, where extreme day is defined as the day with FWI > the 95th percentile of the FWI during the reference period (1980-1984). All FWI datasets cover 1980-2020 baseline and the model simulated products including the downscaled products additionally include 2021-2099 near-future periods under the high-end (SSP585) emission scenario.

54 ENVIRONMENTAL SCIENCES↗

CONUS-wide Projected Flood Frequency and Uncertainty Estimates, Version 1.0

This dataset presents a large-ensemble of CONUS-wide projected flood frequency and uncertainty estimates across ~2.7 million NHDPlusV2 river reaches over the CONUS. The framework producing this dataset leverages a multi-model, uncertainty-aware modeling framework that allows evaluating shifts in flood frequences at the stream reach level across the CONUS. CONUS-wide ensemble streamflow projections generated from hydrologic simulations driven by downscaled and bias-corrected Coupled Model Intercomparison Project Phase 6 (CMIP6) outputs are used to derive these flood frequency and uncertainty estimates over the period 1980 - 2099. A spatially consistent regional L-moment algorithm is applied across clusters defined by the US Hydrologic Unit Code Subregions (HUC4s and HUC8s) and NHDPlusV2 stream orders to estimate flood frequencies. The dataset also includes at-site based flood estimates that allow for the comparison between local and regional approach-based estimates, assess projected changes, and characterize their uncertainties. For more reliable estimation of rare flood frequencies such as 500 and 1000-year return periods, super-ensemble based estimates are also included in the dataset. This dataset is derived to support the "Impact-Informed Dam Safety Risk Assessment for Securing Hydropower Assests" project for the US Department of Energy (DOE) Hydropower and Hydrokinetic Office (H2O). For further details, refer to Kao et al. (2022), Ghimire et al. (2023), Ghimire et al. (2025), and Hosking and Wallis (1997).

Ghimire, Ganesh [ORNL] (ORCID:0000000242843941)↗