Search NASA⌕ Search

SEARCH · Search NASA

Results for “surface modeling”

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 343 records · Page 19

Evaluation of Rooftop Solar Potential in Chernihiv and Lviv, Ukraine, and Efficacy of High-Resolution 3D Data Digital Twins [Slides]

NREL evaluated rooftop solar photovoltaic (PV) siting opportunities in the cities of Chernihiv and Lviv, Ukraine, leveraging very-high-resolution 3D elevation data to calculate technical potential. The study assessed total rooftop solar capacity and annual energy production. In Chernihiv and Lviv, 116,503 buildings were analyzed for their rooftop solar potential. The buildings in the study areas include a mixture of residential, commercial, and industrial buildings, characterized by diverse roof shapes and sizes. The total estimated rooftop solar capacity is 332 MWDC in Chernihiv and 873 MWDC in Lviv with annual energy production up to 376.2 GWhDC in Chernihiv and 995.5 GWhDC in Lviv. This accounting provides a clear estimate of the potential rooftop solar installations that could be realized under optimal conditions, considering both technical constraints and the geographic distribution of available rooftop space. Related to Russia's invasion of Ukraine, NREL estimates loss from buildings damaged or destroyed of 2,754 buildings, 20.15 MWDC of capacity, and 22,869 MWhDC of annual energy production in Chernihiv as well as a loss of 1,316 buildings, 34.49 MWDC of capacity, and 39,369 MWhDC of annual energy production in Lviv. The study also assessed the feasibility of adapting this methodology on a national scale using either simulated digital surface models (DSMs) or a digital twin approach. While the very-high-resolution DSM provided more precise results, the simulated DSM demonstrated reasonable accuracy for broader applications in modeling aggregated distributed solar supply. The feasibility of using a digital twin for Ukraine's national rooftop solar potential is considered promising, with certain limitations in areas with highly variable building stock and heavy war damage.

14 SOLAR ENERGY↗

Satellite-based Investigation of Power-Line Vegetation Encroachment in the US (SILVANUS)

Rapid wide-area assessments of vegetation encroachment on transmission and distribution line rights-of-way (ROW) is a highly desirable capability for understanding risks to the power grid during severe weather and wildfire events. Conventional assessments are time-consuming and expensive due to the need for in-situ inspections and the use of aerial assets. Performing conventional assessments on a wide area would require immense resources and time that might not be available within the horizon of an expected adverse event. Developing a capability to accurately assess vegetation encroachment into ROWs will enable faster analysis of potential grid vulnerabilities in NAERM. This project sought to develop a prototype capability for rapid ROW vegetation encroachment assessments by using Puerto Rico as a test case. Puerto Rico is a heavily forested island territory frequently impacted by tropical cyclones that threaten the electric grid by downing trees across transmission and distribution lines. Multispectral satellite imagery (MSI) enable very high resolution (i.e., 0.5 - 2 meter) assessment of vegetation conditions at scale and with revisit times appropriate for regular monitoring (e.g., weekly to quarterly, depending on cloud cover) of the entire island. Synthetic aperture radar (SAR) data from satellites was also investigated as solution to the cloud-cover issue as they are active sensors that emit and receive a microwave signal rather than relying on solar illumination, and are therefore unaffected by cloud cover and can collect data during day or night. Finally, MSI-derived digital surface models (DSMs) map the height of objects relative to sea level, and were assessed for their ability to estimate the height of forest canopies relative to coincident transmission lines. The results of the mapping investigation are described in this report.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Evaluation of Multi-Fidelity Soil Moisture Products Across the Continental United States

We have aggregated the most recent soil moisture datasets from a diverse range of sources, encompassing the Continental United States (CONUS). These sources encompass gridded data from remote sensing products, reanalysis products, machine learning-based projects, and land surface modeling products. Additionally, we have obtained and processed in-situ soil moisture observations from the International Soil Moisture Network. The collected datasets exhibit variations in both temporal and spatial resolutions. Among the 20 datasets, six are available at a spatial resolution of 0.25 degrees, while three are at a coarser spatial resolution of 25 km. To minimize spatial interpolation, we conducted data uncertainty evaluations at the 0.25-degree spatial resolution. For our data evaluations, we maintained a monthly temporal resolution, which effectively captures soil moisture seasonality and interannual variability. Our data processing strategy preserves the raw data and interpolated data at their original temporal resolutions. Datasets with higher temporal resolutions, including daily, three-hourly, and hourly datasets, are set aside for subsequent analyses. These analyses will delve into topics such as soil moisture changes and recovery during extreme weather events. Furthermore, we have processed auxiliary data to enhance our evaluation, leveraging tools such as Google Earth Engine. This includes incorporating topography data, land use land cover data, Köppen-Geiger climate classification, and more to provide a comprehensive assessment from multiple sources.

Li, Lingcheng↗

Crystal Structure Prediction of Binary Alloys via Deep Potential

Predicting crystal structure has been a challenging problem in physics and materials science for a long time. A reliable energy calculation engine combined with an efficient global search algorithm, such as particle swarm optimization algorithm or genetic algorithm, is needed to conduct crystal structure prediction. In recent years, machine learning-based interatomic potential energy surface models have been proposed, potentially allowing us to perform crystal structure prediction for systems with the accuracy of density functional theory (DFT) and the speed of empirical force fields. In this paper, we employ a previously developed Deep Potential model to predict the intermetallic compound of the aluminum–magnesium system, and find six meta-stable phases with negative or nearly zero formation energy. In particular, Mg 12 Al 8 shows excellent ductility and Mg 5 Al 27 has a high Young's modulus. Based on our benchmark results, we propose a relatively robust structure screening criterion that selects potentially stable structures from the Deep Potential-based convex hull and performs DFT refinement. By using this criterion, the computational cost needed to construct the convex hull with ab initio accuracy can be dramatically reduced.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Pore-Scale Dynamics of Liquid CO2–Water Displacement in 2D Axisymmetric Porous Micromodels Under Strong Drainage and Weak Imbibition Conditions: High-Speed μPIV Measurements

Resolving pore-scale transient flow dynamics is crucial to understanding the physics underlying multiphase flow in porous media and informing large-scale predictive models. Surface properties of the porous matrix play an important role in controlling such physics, yet interfacial mechanisms remain poorly understood, in part due to a lack of direct observations. This study reports on an experimental investigation of the pore-scale flow dynamics of liquid CO 2 and water in two-dimensional (2D) circular porous micromodels with different surface characteristics employing high-speed microscopic particle image velocimetry (μPIV). The design of the micromodel minimized side boundary effects due to the limited size of the domain. The high-speed μPIV technique resolved the spatial and temporal dynamics of multiphase flow of CO 2 and water under reservoir-relevant conditions, for both drainage and imbibition scenarios. When CO 2 displaced water in a hydrophilic micromodel (i.e., drainage), unstable capillary fingering occurred and the pore flow was dominated by successive pore-scale burst events (i.e., Haines jumps). When the same experiment was repeated in a nearly neutral wetting micromodel (i.e., weak imbibition), flow instability and fluctuations were virtually eliminated, leading to a more compact displacement pattern. Energy balance analysis indicates that the conversion efficiency between surface energy and external work is less than 30%, and that kinetic energy is a disproportionately smaller contributor to the energy budget. This is true even during a Haines jump event, which induces velocities typically two orders of magnitude higher than the bulk velocity. These novel measurements further enabled direct observations of the meniscus displacement, revealing a significant alteration of the pore filling mechanisms during drainage and imbibition. While the former typically featured burst events, which often occur only at one of the several throats connecting a pore, the latter is typically dominated by a cooperative filling mechanism involving simultaneous invasion of a pore from multiple throats. This cooperative filling mechanism leads to merging of two interfaces and releases surface energy, causing instantaneous high-speed events that are similar, yet fundamentally different from, burst events. Finally, pore-scale velocity fields were statistically analyzed to provide a quantitative measure of the role of capillary effects in these pore flows.

Li, Yaofa↗

Numerical Optimization of Spray-Guided Spark Assistance for Cold Idle Operation in a Heavy-Duty Gasoline Compression Ignition Engine

This article describes the results of a response surface model (RSM)-based numerical optimization campaign for spray-guided spark assistance at cold operations in a heavy-duty gasoline compression ignition (GCI) engine. On the basis of an earlier work on spark-assisted GCI cold combustion, a space-filling design of experiments (DoE) method was first undertaken to investigate a multitude of hardware design variables and engine operating parameters. The main design variables included the number of injector nozzles, fuel split quantities and injection timings, and spark timing. The objective variables were engine combustion efficiency (ŋ c ), maximum pressure rise rate (MPRR), and engine-out nitrogen oxide (NO x ) emissions. A total of 150 design candidates were automatically generated using the Sobol sequence method provided by the commercial software package, CAESES. Then, closed-cycle computational fluid dynamic (CFD) spark-assisted GCI simulations under cold idling operations were performed. The outcomes from the CFD-DoE design campaign were utilized to construct high-fidelity RSMs that allowed for further design optimization of the spark plug- and fuel injector-related design variables, along with fuel injection strategy parameters. A merit function with respect to objective variables was formulated with an appropriate weight assignment on each objective variable. Finally, the best design candidate was identified from the RSM-based optimization process and further validated in the CFD analysis. The best design candidate showed the potential to significantly improve combustion efficiency (ŋ c > 90%) over the baseline at cold idle while satisfying MPRR and NOx emissions constraints (MPRR < 5 bar/CAD and NO x < 4.5 g/kWh).

33 ADVANCED PROPULSION SYSTEMS↗

BICEP/Keck. XVII. Line-of-sight Distortion Analysis: Estimates of Gravitational Lensing, Anisotropic Cosmic Birefringence, Patchy Reionization, and Systematic Errors

We present estimates of line-of-sight distortion fields derived from the 95 and 150 GHz data taken by BICEP2, BICEP3, and the Keck Array up to the 2018 observing season, leading to cosmological constraints and a study of instrumental and astrophysical systematics. Cosmological constraints are derived from three of the distortion fields concerning gravitational lensing from large-scale structure, polarization rotation from magnetic fields or an axion-like field, and the screening effect of patchy reionization. We measure an amplitude of the lensing power spectrum ${A}_{L}^{\phi \phi }=0.95\pm 0.20$. We constrain polarization rotation, expressed as the coupling constant of a Chern–Simons electromagnetic term g aγ ≤ 2.6 × 10 -2 /H I , where H I is the inflationary Hubble parameter, and an amplitude of primordial magnetic fields smoothed over 1 Mpc B 1Mpc ≤ 6.6 nG at 95 GHz. We constrain the rms of optical depth fluctuations in a simple "crinkly surface" model of patchy reionization, finding A τ < 0.19 (2σ) for the coherence scale of L c = 100. We show that all of the distortion fields of the 95 and 150 GHz polarization maps are consistent with simulations including lensed ΛCDM, dust, and noise, with no evidence for instrumental systematics. In some cases, the EB and TB quadratic estimators presented here are more sensitive than our previous map-based null tests at identifying and rejecting spurious B-modes that might arise from instrumental effects. Finally, we verify that the standard deprojection filtering in the BICEP/Keck data processing is effective at removing temperature to polarization leakage.

79 ASTRONOMY AND ASTROPHYSICS↗

Machine-learning-based investigation of the variables affecting summertime lightning occurrence over the Southern Great Plains

Lightning is affected by many factors, many of which are not routinely measured, well understood, or accounted for in physical models. Several commonly used machine learning (ML) models have been applied to analyze the relationship between Atmospheric Radiation Measurement (ARM) data and lightning data from the Earth Networks Total Lightning Network (ENTLN) in order to identify important variables affecting lightning occurrence in the vicinity of the Southern Great Plains (SGP) ARM site during the summer months (June, July, August and September) of 2012 to 2020. Testing various ML models, we found that the random forest model is the best predictor among common classifiers. When convective clouds were detected, it predicts lightning occurrence with an accuracy of 76.9 % and an area under the curve (AUC) of 0.850. Using this model, we further ranked the variables in terms of their effectiveness in nowcasting lightning and identified geometric cloud thickness, rain rate and convective available potential energy (CAPE) as the most effective predictors. The contrast in meteorological variables between no-lightning and frequent-lightning periods was examined for hours with CAPE values conducive to thunderstorm formation. Besides the variables considered for the ML models, surface variables and mid-altitude variables (e.g., equivalent potential temperature and minimum equivalent potential temperature, respectively) have statistically significant contrasts between no-lightning and frequent-lightning hours. For example, the minimum equivalent potential temperature from 700 to 500 hPa is significantly lower during frequent-lightning hours compared with no-lightning hours. Finally, a notable positive relationship between the intracloud (IC) flash fraction and the square root of CAPE ($\sqrt{CAPE}$) was found, suggesting that stronger updrafts increase the height of the electrification zone, resulting in fewer flashes reaching the surface and consequently a greater IC flash fraction.

54 ENVIRONMENTAL SCIENCES↗

Quantifying the role of moss in terrestrial ecosystem carbon dynamics in northern high latitudes

Mosses are ubiquitous in northern terrestrial ecosystems, and play an important role in regional carbon, water and energy cycling. Current global land surface models that do not consider mosses may bias the quantification of regional carbon dynamics. Here we incorporate mosses as a new plant functional type into the process-based Terrestrial Ecosystem Model (TEM 5.0), to develop a new model (TEM_Moss). The new model explicitly quantifies the interactions between vascular plants and mosses and their competition for energy, water, and nutrients. Compared to the estimates using TEM 5.0, the new model estimates that the regional terrestrial soils currently store 132.7 Pg more C and will store 157.5 and 179.1 Pg more C under the RCP8.5 and RCP2.6 scenarios, respectively, by the end of the 21st century. Ensemble regional simulations forced with different parameters for the 21st century with TEM_Moss predict that the region will accumulate 161.1±142.1 Pg C under the RCP2.6 scenario and 186.7±166.1 Pg C under the RCP8.5 scenario over the century. Our study highlights the necessity of coupling moss into Earth system models to adequately quantify terrestrial carbon–climate feedbacks in the Arctic.

54 ENVIRONMENTAL SCIENCES↗

X-BASE: the first terrestrial carbon and water flux products from an extended data-driven scaling framework, FLUXCOM-X

Mapping in situ eddy covariance measurements of terrestrial land–atmosphere fluxes to the globe is a key method for diagnosing the Earth system from a data-driven perspective. We describe the first global products (called X-BASE) from a newly implemented upscaling framework, FLUXCOM-X, representing an advancement from the previous generation of FLUXCOM products in terms of flexibility and technical capabilities. The X-BASE products are comprised of estimates of CO 2 net ecosystem exchange (NEE), gross primary productivity (GPP), evapotranspiration (ET), and for the first time a novel, fully data-driven global transpiration product (ETT), at high spatial (0.05°) and temporal (hourly) resolution. X-BASE estimates the global NEE at −5.75 ± 0.33 Pg C yr −1 for the period 2001–2020, showing a much higher consistency with independent atmospheric carbon cycle constraints compared to the previous versions of FLUXCOM. The improvement of global NEE was likely only possible thanks to the international effort to increase the precision and consistency of eddy covariance collection and processing pipelines, as well as to the extension of the measurements to more site years resulting in a wider coverage of bioclimatic conditions. However, X-BASE global net ecosystem exchange shows a very low interannual variability, which is common to state-of-the-art data-driven flux products and remains a scientific challenge. With 125 ± 2.1 Pg C yr −1 for the same period, X-BASE GPP is slightly higher than previous FLUXCOM estimates, mostly in temperate and boreal areas. X-BASE evapotranspiration amounts to 74.7×10 3 ± 0.9×10 3 km 3 globally for the years 2001–2020 but exceeds precipitation in many dry areas, likely indicating overestimation in these regions. On average 57 % of evapotranspiration is estimated to be transpiration, in good agreement with isotope-based approaches, but higher than estimates from many land surface models. Despite considerable improvements to the previous upscaling products, many further opportunities for development exist. Pathways of exploration include methodological choices in the selection and processing of eddy covariance and satellite observations, their ingestion into the framework, and the configuration of machine learning methods. For this, the new FLUXCOM-X framework was specifically designed to have the necessary flexibility to experiment, diagnose, and converge to more accurate global flux estimates.

Nelson, Jacob A.↗

Development of observation-based global multilayer soil moisture products for 1970 to 2016

Abstract. Soil moisture (SM) datasets are critical to understanding the global water, energy, and biogeochemical cycles and benefit extensive societal applications. However, individual sources of SM data (e.g., in situ and satellite observations, reanalysis, offline land surface model simulations, Earth system model – ESM – simulations) have source-specific limitations and biases related to the spatiotemporal continuity, resolutions, and modeling and retrieval assumptions. Here, we developed seven global, gap-free, long-term (1970–2016), multilayer (0–10, 10–30, 30–50, and 50–100 cm) SM products at monthly 0.5∘ resolution (available at https://doi.org/10.6084/m9.figshare.13661312.v1; Wang and Mao, 2021) by synthesizing a wide range of SM datasets using three statistical methods (unweighted averaging, optimal linear combination, and emergent constraint). The merged products outperformed their source datasets when evaluated with in situ observations (mean bias from −0.044 to 0.033 m3 m−3, root mean square errors from 0.076 to 0.104 m3 m−3, Pearson correlations from 0.35 to 0.67) and multiple gridded datasets that did not enter merging because of insufficient spatial, temporal, or soil layer coverage. Three of the new SM products, which were produced by applying any of the three merging methods to the source datasets excluding the ESMs, had lower bias and root mean square errors and higher correlations than the ESM-dependent merged products. The ESM-independent products also showed a better ability to capture historical large-scale drought events than the ESM-dependent products. The merged products generally showed reasonable temporal homogeneity and physically plausible global sensitivities to observed meteorological factors, except that the ESM-dependent products underestimated the low-frequency temporal variability in SM and overestimated the high-frequency variability for the 50–100 cm depth. Based on these evaluation results, the three ESM-independent products were finally recommended for future applications because of their better performances than the ESM-dependent ones. Despite uncertainties in the raw SM datasets and fusion methods, these hybrid products create added value over existing SM datasets because of the performance improvement and harmonized spatial, temporal, and vertical coverages, and they provide a new foundation for scientific investigation and resource management.

54 ENVIRONMENTAL SCIENCES↗

Implementation and evaluation of the unified stomatal optimization approach in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES)

Abstract. Stomata play a central role in regulating the exchange of carbon dioxide and water vapor between ecosystems and the atmosphere. Their function is represented in land surface models (LSMs) by conductance models. The Functionally Assembled Terrestrial Ecosystem Simulator (FATES) is a dynamic vegetation demography model that can simulate both detailed plant demographic and physiological dynamics. To evaluate the effect of stomatal conductance model formulation on forest water and carbon fluxes in FATES, we implemented an optimality-based stomatal conductance model – the Medlyn (MED) model – that simulates the relationship between photosynthesis (A) and stomatal conductance to water vapor (gsw) as an alternative to the FATES default Ball–Woodrow–Berry (BWB) model. To evaluate how the behavior of FATES is affected by stomatal model choice, we conducted a model sensitivity analysis to explore the response of gsw to climate forcing, including atmospheric CO2 concentration, air temperature, radiation, and vapor pressure deficit in the air (VPDa). We found that modeled gsw values varied greatly between the BWB and MED formulations due to the different default stomatal slope parameters (g1). After harmonizing g1 and holding the stomatal intercept parameter (g0) constant for both model formulations, we found that the divergence in modeled gsw was limited to conditions when the VPDa exceeded 1.5 kPa. We then evaluated model simulation results against measurements from a wet evergreen forest in Panama. Results showed that both the MED and BWB model formulations were able to capture the magnitude and diurnal changes of measured gsw and A but underestimated both by about 30 % when the soil was predicted to be very dry. Comparison of modeled soil water content from FATES to a reanalysis product showed that FATES captured soil drying well, but translation of drying soil to modeled physiology reduced the models' ability to match observations. Our study suggests that the parameterization of stomatal conductance models and current model response to drought are the critical areas for improving model simulation of CO2 and water fluxes in tropical forests.

54 ENVIRONMENTAL SCIENCES↗

Algorithmically detected rain-on-snow flood events in different climate datasets: a case study of the Susquehanna River basin

Abstract. Rain-on-snow (RoS) events in regions of ephemeral snowpack – such as the northeastern United States – can be key drivers of cool-season flooding. We describe an automated algorithm for detecting basin-scale RoS events in gridded climate data by generating an area-averaged time series and then searching for periods of concurrent precipitation, surface runoff, and snowmelt exceeding predefined thresholds. When evaluated using historical data over the Susquehanna River basin (SRB), the technique credibly finds RoS events in published literature and flags events that are followed by anomalously high streamflow as measured by gauge data along the river. When comparing four different datasets representing the same 21-year period, we find large differences in RoS event magnitude and frequency, primarily driven by differences in estimated surface runoff and snowmelt. Using dataset-specific thresholds improves agreement between datasets but does not account for all discrepancies. We show that factors such as meteorological forcing and coupling frequency, as well as choice of land surface model, play roles in how data products capture these compound extremes and suggest care is to be taken when climate datasets are used by stakeholders for operational decision-making.

54 ENVIRONMENTAL SCIENCES↗

CLM5 CAMELS Basins Ensemble

Land surface models such as Community Land Model Version 5 (CLM5) are essential tools for simulating the behaviors of the terrestrial system. Despite the extensive application of CLM5, limited attention has been paid to the underlying uncertainties associated with its hydrologic parameters and the implications that these uncertainties have on water resources applications. To address this long-standing issue, we conduct a comprehensive hydrologic parameter uncertainty characterization (UC) of CLM5 over the hydroclimatic gradients of the Conterminous United States using five meteorological datasets. Key datasets produced from the UC experiment include a benchmark dataset of CLM5 default hydrological performance, parameter sensitivity identified for 28 hydrological metrics, and large ensemble outputs for hydrological predictions. The presented datasets can assist CLM5 calibration and to support broad applications such as evaluating vulnerabilities to droughts and floods. The dataset can be used to identify under what hydroclimate conditions parametric uncertainties demonstrate substantial effects on hydrological predictions and clarify where further investigations are needed to understand how land runoff uncertainties interact with other Earth system processes. Please refer to the included README.pdf for a description of the included files. Note that raw CLM5 model outputs for each forcing dataset are hosted on a Globus endpoint: https://app.globus.org/file-manager?destination_id=d22ef858-27b0-11ed-a910-fd3165076336.

Yan, Hongxiang↗

CLM5 CAMELS Basins Ensemble

Land surface models such as Community Land Model Version 5 (CLM5) are essential tools for simulating the behaviors of the terrestrial system. Despite the extensive application of CLM5, limited attention has been paid to the underlying uncertainties associated with its hydrologic parameters and the implications that these uncertainties have on water resources applications. To address this long-standing issue, we conduct a comprehensive hydrologic parameter uncertainty characterization (UC) of CLM5 over the hydroclimatic gradients of the Conterminous United States using five meteorological datasets. Key datasets produced from the UC experiment include a benchmark dataset of CLM5 default hydrological performance, parameter sensitivity identified for 28 hydrological metrics, and large ensemble outputs for hydrological predictions. The presented datasets can assist CLM5 calibration and to support broad applications such as evaluating vulnerabilities to droughts and floods. The dataset can be used to identify under what hydroclimate conditions parametric uncertainties demonstrate substantial effects on hydrological predictions and clarify where further investigations are needed to understand how land runoff uncertainties interact with other Earth system processes. Please refer to the README.pdf for a description of the included files. Note that raw CLM5 model outputs for each forcing dataset are hosted on a Globus endpoint: https://app.globus.org/file-manager?destination_id=d22ef858-27b0-11ed-a910-fd3165076336.

Yan, Hongxiang↗

Responsible Adoption of Artificial Intelligence (AI) in Electric Grid Operations

The future of the grid will be powered by AI—or undermined by it. Artificial intelligence is rapidly reshaping grid operations, improving fault detection, forecasting accuracy, and real-time optimization. As AI systems move closer to operational decision loops, however, they introduce new consequence pathways: expanded attack surfaces, model integrity risks, regulatory exposure, and human-automation challenges. This talk presents a consequence-driven framework for deploying AI responsibly in the electric grid. Attendees will gain practical strategies to strengthen resilience, boost reliability, and deploy AI securely — ensuring the grid of the future is not only smarter but safer.

25 - ENERGY STORAGE↗

A Multicomponent Reactive Transport Model for Integrated Surface‐Subsurface Hydrology Problems

Abstract Despite the widespread use of integrated hydrology models in a variety of applications, consideration of multicomponent reactive transport is still not common. The implementation of these processes requires coupling transport at the surface‐subsurface interface and efficient solution of the non‐linear geochemical model that is consistent with the integrated hydrology solution. The Advanced Terrestrial Simulator provides a flexible multiphysics framework that facilitated this process. In this work, the integrated reactive transport process kernel (PK) was weakly coupled to the integrated hydrology PK. In turn, integrated transport and reactions were coupled using an operator splitting approach. This splitting enabled an explicit solution of the integrated transport problem, including a novel algorithm to calculate exchange fluxes across the surface‐subsurface interface and a point‐by‐point solution of the geochemical problem. Geochemical capabilities were added using well‐established external codes, but rather than using a custom interface to each, a generic interface was used that clearly specifies the variables and operations used by the chemistry PK. The implementation is demonstrated with two example simulations: transport of a tracer in a soil column as it saturates over time and water ponds on the surface and reactive transport in a hillslope driven by successive wet‐dry cycles that result in infiltration, runoff and exfiltration processes.

54 ENVIRONMENTAL SCIENCES↗

The Modeled Seasonal Cycles of Surface N 2 O Fluxes and Atmospheric N 2 O

Nitrous oxide (N 2 O) is a greenhouse gas and stratospheric ozone-depleting substance with large and growing anthropogenic emissions. Previous studies identified the influx of N 2 O-depleted air from the stratosphere to partly cause the seasonality in tropospheric N 2 O (aN 2 O), but other contributions remain unclear. Here, we combine surface fluxes from eight land and four ocean models from phase 2 of the Nitrogen/N 2 O Model Intercomparison Project with tropospheric transport modeling to simulate aN 2 O at eight remote air sampling sites for modern and pre-industrial periods. Models show general agreement on the seasonal phasing of zonal-average N 2 O fluxes for most sites, but seasonal peak-to-peak amplitudes differ several-fold across models. The modeled seasonal amplitude of surface aN 2 O ranges from 0.25 to 0.80 ppb (interquartile ranges 21%–52% of median) for land, 0.14–0.25 ppb (17%–68%) for ocean, and 0.28–0.77 ppb (23%–52%) for combined flux contributions. The observed seasonal amplitude ranges from 0.34 to 1.08 ppb for these sites. The stratospheric contributions to aN 2 O, inferred by the difference between the surface-troposphere model and observations, show 16%–126% larger amplitudes and minima delayed by ~1 month compared to Northern Hemisphere site observations. Land fluxes and their seasonal amplitude have increased since the pre-industrial era and are projected to grow further under anthropogenic activities. Our results demonstrate the increasing importance of land fluxes for aN 2 O seasonality. Considering the large model spread, in situ aN 2 O observations and atmospheric transport-chemistry models will provide opportunities for constraining terrestrial and oceanic biosphere models, critical for projecting carbon-nitrogen cycles under ongoing global warming.

54 ENVIRONMENTAL SCIENCES↗