Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial variability”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

From individual leaves to forest stands: importance of niche, distance decay, and stochasticity vary by ecosystem type and functional group for fungal community composition

Abstract Community assembly is influenced by environmental niche processes as well as stochastic processes that can be spatially dependent (e.g. dispersal limitation) or independent (e.g. priority effects). Here, we sampled senesced tree leaves as unit habitats to investigate fungal community assembly at two spatial scales: (i) small neighborhoods of overlapping leaves from differing tree species and (ii) forest stands of differing ecosystem types. Among forest stands, ecosystem type explained the most variation in community composition. Among adjacent leaves within stands, variability in fungal composition was surprisingly high. Leaf type was more important in stands with high soil fertility and dominated by differing tree mycorrhizal types (sugar maple vs. basswood or red oak), whereas distance decay was more important in oak-dominated forest stands with low soil fertility. Abundance of functional groups was explained by environmental factors, but predictors of taxonomic composition within differing functional groups were highly variable. These results suggest that fungal community assembly processes are clearest for functional group abundances and large spatial scales. Understanding fungal community assembly at smaller spatial scales will benefit from further study focusing on differences in drivers for different ecosystems and functional groups, as well as the importance of spatially independent factors such as priority effects.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Robust Anthropogenic Signal Identified in the Seasonal Cycle of Tropospheric Temperature

Previous work identified an anthropogenic fingerprint pattern in T AC (x, t), the amplitude of the seasonal cycle of mid- to upper-tropospheric temperature (TMT), but did not explicitly consider whether fingerprint identification in satellite T AC (x, t) data could have been influenced by real-world multidecadal internal variability (MIV). Here we address this question here using large ensembles (LEs) performed with five climate models. LEs provide many different sequences of internal variability noise superimposed on an underlying forced signal. Despite differences in historical external forcings, climate sensitivity, and MIV properties of the five models, their T AC (x, t) fingerprints are similar and statistically identifiable in 239 of the 240 LE realizations of historical climate change. Comparing simulated and observed variability spectra reveals that consistent fingerprint identification is unlikely to be biased by model underestimates of observed MIV. Even in the presence of large (factor of 3–4) intermodel and inter-realization differences in the amplitude of MIV, the anthropogenic fingerprints of seasonal cycle changes are robustly identifiable in models and satellite data. This is primarily due to the fact that the distinctive, global-scale fingerprint patterns are spatially dissimilar to the smaller-scale patterns of internal T AC (x, t) variability associated with the Atlantic multidecadal oscillation and El Niño–Southern Oscillation. The robustness of the seasonal cycle detection and attribution results shown here, taken together with the evidence from idealized aquaplanet simulations, suggest that basic physical processes are dictating a common pattern of forced T AC (x, t) changes in observations and in the five LEs. The key processes involved include GHG-induced expansion of the tropics, lapse-rate changes, land surface drying, and sea ice decrease.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The impact of multi-sensor land data assimilation on river discharge estimation

River discharge is one of the most critical renewable water resources. Accurately estimating river discharge with land surface models (LSMs) remains challenging due to the difficulty in estimating land water storages such as snow, soil moisture, and groundwater. While data assimilation (DA) ingesting optical, microwave, and gravity measurements from space can help constrain theses storage states, its impacts on runoff and eventually river discharge are not fully understood. In this study, by taking advantage of recently published land DA results that jointly assimilate eight different combinations of observations from the Moderate Resolution Imaging Spectroradiometer (MODIS), Gravity Recovery and Climate Experiment (GRACE), and Advanced Microwave Scanning Radiometer for EOS (AMSR-E), we quantify to what degree multi-sensor land DA improves the river discharge simulation skills over 40 global river basins, and investigate the complementary strengths of different satellite measurements on river discharge. To be more specific, river discharge is updated by feeding gridded runoff from the eight multi-sensor DA simulations into a vector-based river routing model named the Routing Application for Parallel computatIon of Discharge (RAPID). Our modeling results, including 7-year simulations at 177,458 river reaches globally, are used to study the seasonal to interannual variability of river discharge. It is found that assimilating GRACE has the greatest impact on global runoff patterns, leading to the most pronounced improvements in spatial river discharge in the middle and high latitudes with the R 2 increased by 0.16. The seasonal variation of spatial discharge is most skillful during the boreal summer. However, our evaluation also shows model and DA still struggle to generate reasonable variability and averaged discharge over permafrost regions. Finally, by assessing how different satellites add value to discharge forecasts, this study paves the way for more advanced multi-sensor satellite data assimilation to predict the terrestrial hydrological cycle.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Microstructural Topology as a Prescriptor for Quantum Coherence: Towards A Unified Framework for Decoherence in Superconducting Qubits

In superconducting quantum circuits, decoherence improvements are frequently obtained through process interventions that simultaneously modify surface chemistry, microstructural topology, and device geometry, leaving mechanistic attribution structurally underdetermined. Predictive materials engineering requires measurable structural statistics to be separated from geometry-dependent coupling coefficients into independently testable factors. We introduce the concept of classical and quantum microstructure. In that context, we formulate a channel-wise separable framework for decoherence in superconducting transmon qubits in which each loss channel is described by a reduced prescriptor. Here, a channel-specific microstructural state variable is determined independently of device geometry, and a geometry-dependent coupling functional is computable from field solutions without reference to surface chemistry. We derive this product form from a spatially resolved kernel representation and establish a perturbative separability criterion that defines the regime where independent variation of the variables is valid. The framework specifies five prescriptor classes for dominant loss pathways in transmon-class devices. Falsifiability is operationalized through a pre-committed 2x2 experimental protocol in which the variables must satisfy independent ratio checks within propagated uncertainty. A Minimum-Dataset Specification standardizes reporting for cross-laboratory inference. Part I establishes the conceptual and mathematical architecture; coordinated experimental validation is reserved for Part II.

Dravid, Vinayak P. [Northwestern U.]↗
DOE OSTI2025

Parallel-in-Time Solution of Hyperbolic PDE Systems via Characteristic-Variable Block Preconditioning

We consider the parallel-in-time solution of both linear and nonlinear hyperbolic partial differential equation (PDE) systems in one spatial dimension. In the nonlinear setting, the discretized equations are solved with a preconditioned residual iteration based on a global linearization. The linear(ized) equation systems are approximately solved parallel-in-time using a block preconditioner applied in the characteristic variables of the underlying linear(ized) hyperbolic PDE. This change of variables is motivated by the observation that intervariable coupling between characteristic variables is weak, at least locally where spatio-temporal variations in the eigenvectors of the associated flux Jacobian are sufficiently small, while that between the original variables is not. For an ℓ-dimensional system of PDEs, applying the preconditioner consists of solving a sequence of ℓ scalar linear(ized)-advection-like problems, each associated with a different characteristic wave-speed in the underlying linear(ized) PDE. Furthermore, we approximately solve these linear advection problems using multigrid reduction-in-time (MGRIT); however, any other suitable parallel-in-time method could be used. Numerical examples are shown for the (linear) acoustics equations in heterogeneous media and for the (nonlinear) shallow water equations and Euler equations of gas dynamics with shocks and rarefactions. For many test problems, the solver converges in just a handful of iterations and with mesh-independent convergence rates.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Large spatiotemporal variability in aerosol properties over central Argentina during the CACTI field campaign

Abstract. Few field campaigns with extensive aerosol measurements have been conducted over continental areas in the Southern Hemisphere. To address this data gap and better understand the interactions of convective clouds and the surrounding environment, extensive in situ and remote sensing measurements were collected during the Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign conducted between October 2018 and April 2019 over the Sierras de Córdoba range of central Argentina. This study describes measurements of aerosol number, size, composition, mixing state, and cloud condensation nuclei (CCN) collected on the ground and from a research aircraft during 7 weeks of the campaign. Large spatial and multiday variations in aerosol number, size, composition, and CCN were observed due to transport from upwind sources controlled by mesoscale to synoptic-scale meteorological conditions. Large vertical wind shears, back trajectories, single-particle measurements, and chemical transport model predictions indicate that different types of emissions and source regions, including biogenic emissions and biomass burning from the Amazon and anthropogenic emissions from Chile and eastern Argentina, contribute to aerosols observed during CACTI. Repeated aircraft measurements near the boundary layer top reveal strong spatial and temporal variations in CCN and demonstrate that understanding the complex co-variability of aerosol properties and clouds is critical to quantify the impact of aerosol–cloud interactions. In addition to quantifying aerosol properties in this data-sparse region, these measurements will be valuable to evaluate predictions over the midlatitudes of South America and improve parameterized aerosol processes in local, regional, and global models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A Process‐Model Perspective on Recent Changes in the Carbon Cycle of North America

Continental North America has been found to be a carbon (C) sink over recent decades by multiple studies employing a variety of estimation approaches. However, several key questions and uncertainties remain with these assessments. Here we used results from an ensemble of 19 state-of-the-art dynamic global vegetation models from the TRENDYv9 project to improve these estimates and study the drivers of its interannual variability. Our results show that North America has been a C sink with a magnitude of 0.37 ± 0.38 (mean and one standard deviation) PgC year —1 for the period 2000–2019 (0.31 and 0.44 PgC year —1 in each decade); split into 0.18 ± 0.12 PgC year —1 in Canada (0.15 and 0.20), 0.16 ± 0.17 in the United States (0.14 and 0.17), 0.02 ± 0.05 PgC year —1 in Mexico (0.02 and 0.02) and 0.01 ± 0.02 in Central America and the Caribbean (0.01 and 0.01). About 57% of the new C assimilated by terrestrial ecosystems is allocated into vegetation, 30% into soils, and 13% into litter. Losses of C due to fire account for 41% of the interannual variability of the mean net biome productivity for all North America in the model ensemble. Finally, we show that drought years (e.g., 2002) have the potential to shift the region to a small net C source in the simulations (—0.02 ± 0.46 PgC year —1 ). Our results highlight the importance of identifying the major drivers of the interannual variability of the continental-scale land C cycle along with the spatial distribution of local sink-source dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A framework for testing soil carbon dynamics post land-use transition in a multisector dynamics model

Soil carbon plays a crucial role in the global carbon cycle. Changes in land use can determine whether carbon is stored or is emitted into the atmosphere as carbon dioxide, which has broad implications for the human and Earth systems. These feedbacks to the carbon cycle and their socio-economic drivers are modelled by many global multisector dynamics models to project future possibilities for the human-Earth system. One notable model of this class is the Global Change Analysis Model (GCAM), which uses a simplified process to model soil organic carbon (SOC) content after land-use transition across 384 land units. While the current GCAM soil carbon framework is based on scientific principles, it has not been tested against experimental data. This work examines rates of SOC change from GCAM input data. Specifically, first order rate constants derived from model inputs were compared to values from two syntheses to assess GCAM’s accuracy. Welch’s t-tests and linear models were used to determine if rate constants were consistent across all tested geographical areas and land-use transition types. While we found that there was general agreement on the direction and magnitude (i.e., rate) of SOC change, the rate constant derived from GCAM and empirical values differed strongly in a subset of specific instances. These results indicate that GCAM’s current SOC dynamics during land use transition successfully capture broad patterns of change in this critical carbon pool, but should be interpreted with caution at finer spatial scales. One potential cause of these discrepancies is our highly aggregated variable, soil timescale, which could be made more granular to improve accuracy. When using economically rooted multisector dynamics models, such as GCAM, it is critical to understand such model limitations for representing specific Earth system processes.

carbon↗
DOE OSTI2022

Inverse Calculation of Burden Distribution Matrix Using B-spline Model Based PDF control in Blast Furnace Burden Charging Process

The inverse calculation of burden distribution matrix (BDM) is one of the most important challenges in the blast furnace operation in iron-making processes. In general, blast furnace consumes 65% of the total energy for the whole steel-making. Focusing on this practical challenge, this article proposes a new burden distribution spatial model in calculating burden charging process, and develops a B-spline approximation-based probability density function (PDF) control algorithm to assign the expected thickness distribution of burden layer and, thus, develops a new method for the required inverse calculation of BDM. First, a novel method for the thickness distribution of burden layer is given using B-spline model to produce an expected distribution shape subjected to a desired tracking within a specific spatial constraint. Then, according to the coexistence of continuous and bounded discrete variables in BDM, a novel hybrid optimization control method by combining integer programming and PDF tracking is further established for the effective inverse calculation of BDM. Finally, the proposed PDF-based iterative inverse calculation of BDM using B-spline models are tested using various data from industrial examples. Furthermore, the simulation results show that the proposed method is well suited to solve the BDM inverse calculation problem in practice.

42 ENGINEERING↗
DOE OSTI2023

Organic preservation of vase‐shaped microfossils from the late Tonian Chuar Group, Grand Canyon, Arizona, USA

Abstract Vase‐shaped microfossils (VSMs) are found globally in middle Neoproterozoic (800–730 Ma) marine strata and represent the earliest evidence for testate (shell‐forming) amoebozoans. VSM tests are hypothesized to have been originally organic in life but are most commonly preserved as secondary mineralized casts and molds. A few reports, however, suggest possible organic preservation. Here, we test the hypothesis that VSMs from shales of the lower Walcott Member of the Chuar Group, Grand Canyon, Arizona, contain original organic material, as reported by B. Bloeser in her pioneering studies of Chuar VSMs. We identified VSMs from two thin section samples of Walcott Member black shales in transmitted light microscopy and used scanning electron microscopy to image VSMs. Carbonaceous material is found within the internal cavity of all VSM tests from both samples and is interpreted as bitumen mobilized from Walcott shales likely during the Cretaceous. Energy dispersive X‐ray spectroscopy (EDS) and wavelength dispersive X‐ray spectroscopy (WDS) reveal that VSM test walls contain mostly carbon, iron, and sulfur, while silica is present only in the surrounding matrix. Raman spectroscopy was used to compare the thermal maturity of carbonaceous material within the samples and indicated the presence of pyrite and jarosite within fossil material. X‐ray absorption spectroscopy revealed the presence of reduced organic sulfur species within the carbonaceous test walls, the carbonaceous material found within test cavities, and in the sedimentary matrix, suggesting that organic matter sulfurization occurred within the Walcott shales. Our suite of spectroscopic analyses reveals that Walcott VSM test walls are organic and sometimes secondarily pyritized (with the pyrite variably oxidized to jarosite). Both preservation modes can occur at a millimeter spatial scale within sample material, and at times even within a single specimen. We propose that sulfurization within the Walcott Shales promoted organic preservation, and furthermore, the ratio of iron to labile VSM organic material controlled the extent of pyrite replacement. Based on our evidence, we conclude that the VSMs are preserved with original organic test material, and speculate that organic VSMs may often go unrecognized, given their light‐colored, translucent appearance in transmitted light.

58 GEOSCIENCES↗
DOE OSTI2021

Scaling plant responses to high temperature from cell to ecosystem

Plant survival and productivity decrease with exposure to less favourable environments. The negative effects of environmental factors such as water-deficit and salinity stress on vegetation are relatively well studied, primarily due to the large losses in crop yield and forest biomass caused by these stresses. Although the impact of combined heat and drought on plants has received considerable interest in recent years due to climate change, there has been less focus on ascertaining the impact of heat stress alone at a range of temporal and spatial scales. The steady increase in global mean temperatures and (more importantly) variability in extreme temperatures is translating into more frequent and extreme heat stress episodes, which in turn reduce plant productivity. Further, projected increases in both day and night temperatures and the predicted impact of these changes on agricultural and forest productivity have drawn attention to the need to better address the direct effects of heat on plants. In particular, there has been tremendous interest in understanding the short- and long-term effects of heat at the cellular and whole-plant levels, including work on crops of agricultural importance and species from less managed ecosystems. Hence, we made a concerted effort in this Special Issue to capture recent progress in our ability to quantify the impact of heat stress, and identify new tools and approaches that have been developed and implemented to capture a wide range of responses in plants. The series of reviews and research articles presented here captures our current understanding of plant physiological and molecular responses to heat stress, as well as covering phenotyping methods, genetic diversity and potential routes through which plants can acclimate and adapt to warming climates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Action Effectiveness Monitoring and Research of Dredged Material Placement at Woodland Islands (Final Report)

Benthic invertebrates play critical ecosystem roles including the breakdown of organic matter, sediment mixing, and nutrient cycling. In the lower Columbia River and estuary, benthic invertebrates provide a foraging resource for threatened and endangered juvenile salmon. Driven by the goal to create low velocity, shallow water, and riparian shrub habitats to benefit juvenile salmon, the USACE placed 237,000 CY of dredged material, resulting in the deposition of 13.5 acres of sand on the off-channel margins at Woodland Islands. The Pacific Northwest National Laboratory (PNNL) designed and implemented an action effectiveness research study to understand how dredged material placement at Woodland Islands affected sediment conditions and benthic invertebrate assemblages. A Before-After-Control-Impact, or BACI, study design was used to evaluate the response of environmental conditions and benthic invertebrates to dredged material placement. The design included one impact site and two control sites. Sampling occurred for two years prior to dredged material placement and two years after. Our study found that spatial variation was a significant factor for both environmental and biological response variables which suggests local conditions are important considerations for mechanisms affecting benthic assemblages. We found that off-channel habitats across all locations sampled were producing benthic invertebrates, many of which are common prey items for juvenile salmon and steelhead—e.g., insects, chironomids, crustaceans, and corophium. There was a significant BACI effect (i.e., an effect at the dredged material placement site, relative to conditions at the control sites) for concentration of carbon and ammonium in sediment, but not for phosphorous concentrations or for percent sand. The estimated abundance for the three invertebrate response variables—total abundance, total chironomid abundance, and total corophium abundance—was significantly lower at the impact site after dredged material placement, compared to the control sites. At Woodland Island, the estimated mean abundance for all invertebrates combined decreased 28% after dredged material placement. Estimated mean abundances of chironomid and corophium decreased by 8% and 88%, respectively. While invertebrate abundances were lower after dredged placement, the composition of benthic invertebrates was similar before and after placement suggesting that as the new habitat feature evolves recolonization will likely follow. These findings provide a foundation for understanding potential benefits and consequences of repurposing dredged material for habitat creation in the LCRE.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Hierarchical Fractional Advection-Dispersion Equation (FADE) to Quantify Anomalous Transport in River Corridor over a Broad Spectrum of Scales: Theory and Applications

Fractional calculus-based differential equations were found by previous studies to be promising tools in simulating local-scale anomalous diffusion for pollutants transport in natural geological media (geomedia), but efficient models are still needed for simulating anomalous transport over a broad spectrum of scales. This study proposed a hierarchical framework of fractional advection-dispersion equations (FADEs) for modeling pollutants moving in the river corridor at a full spectrum of scales. Applications showed that the fixed-index FADE could model bed sediment and manganese transport in streams at the geomorphologic unit scale, whereas the variable-index FADE well fitted bedload snapshots at the reach scale with spatially varying indices. Further analyses revealed that the selection of the FADEs depended on the scale, type of the geomedium (i.e., riverbed, aquifer, or soil), and the type of available observation dataset (i.e., the tracer snapshot or breakthrough curve (BTC)). When the pollutant BTC was used, a single-index FADE with scale-dependent parameters could fit the data by upscaling anomalous transport without mapping the sub-grid, intermediate multi-index anomalous diffusion. Pollutant transport in geomedia, therefore, may exhibit complex anomalous scaling in space (and/or time), and the identification of the FADE’s index for the reach-scale anomalous transport, which links the geomorphologic unit and watershed scales, is the core for reliable applications of fractional calculus in hydrology.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Inpainting radar missing data regions with deep learning

Abstract. Missing and low-quality data regions are a frequent problem for weather radars. They stem from a variety of sources: beam blockage, instrument failure, near-ground blind zones, and many others. Filling in missing data regions is often useful for estimating local atmospheric properties and the application of high-level data processing schemes without the need for preprocessing and error-handling steps – feature detection and tracking, for instance. Interpolation schemes are typically used for this task, though they tend to produce unrealistically spatially smoothed results that are not representative of the atmospheric turbulence and variability that are usually resolved by weather radars. Recently, generative adversarial networks (GANs) have achieved impressive results in the area of photo inpainting. Here, they are demonstrated as a tool for infilling radar missing data regions. These neural networks are capable of extending large-scale cloud and precipitation features that border missing data regions into the regions while hallucinating plausible small-scale variability. In other words, they can inpaint missing data with accurate large-scale features and plausible local small-scale features. This method is demonstrated on a scanning C-band and vertically pointing Ka-band radar that were deployed as part of the Cloud Aerosol and Complex Terrain Interactions (CACTI) field campaign. Three missing data scenarios are explored: infilling low-level blind zones and short outage periods for the Ka-band radar and infilling beam blockage areas for the C-band radar. Two deep-learning-based approaches are tested, a convolutional neural network (CNN) and a GAN that optimize pixel-level error or combined pixel-level error and adversarial loss respectively. Both deep-learning approaches significantly outperform traditional inpainting schemes under several pixel-level and perceptual quality metrics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Technical note: Using long short-term memory models to fill data gaps in hydrological monitoring networks

Abstract. Quantifying the spatiotemporal dynamics in subsurface hydrological flows over a long time window usually employs a network of monitoring wells. However, such observations are often spatially sparse with potential temporal gaps due to poor quality or instrument failure. In this study, we explore the ability of recurrent neural networks to fill gaps in a spatially distributed time-series dataset. We use a well network that monitors the dynamic and heterogeneous hydrologic exchanges between the Columbia River and its adjacent groundwater aquifer at the U.S. Department of Energy's Hanford site. This 10-year-long dataset contains hourly temperature, specific conductance, and groundwater table elevation measurements from 42 wells with gaps of various lengths. We employ a long short-term memory (LSTM) model to capture the temporal variations in the observed system behaviors needed for gap filling. The performance of the LSTM-based gap-filling method was evaluated against a traditional autoregressive integrated moving average (ARIMA) method in terms of error statistics and accuracy in capturing the temporal patterns of river corridor wells with various dynamics signatures. Our study demonstrates that the ARIMA models yield better average error statistics, although they tend to have larger errors during time windows with abrupt changes or high-frequency (daily and subdaily) variations. The LSTM-based models excel in capturing both high-frequency and low-frequency (monthly and seasonal) dynamics. However, the inclusion of high-frequency fluctuations may also lead to overly dynamic predictions in time windows that lack such fluctuations. The LSTM can take advantage of the spatial information from neighboring wells to improve the gap-filling accuracy, especially for long gaps in system states that vary at subdaily scales. While LSTM models require substantial training data and have limited extrapolation power beyond the conditions represented in the training data, they afford great flexibility to account for the spatial correlations, temporal correlations, and nonlinearity in data without a priori assumptions. Thus, LSTMs provide effective alternatives to fill in data gaps in spatially distributed time-series observations characterized by multiple dominant frequencies of variability, which are essential for advancing our understanding of dynamic complex systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

ReEDS Performance Improvement

The Regional Energy Deployment System (ReEDS) is an open-source, spatially explicit, long-term capacity expansion model for the bulk electric power system of the contiguous United States, encompassing multiple scenarios with technological and political assumptions (see https://github.com/NREL/ReEDS-2.0). With the increased needs for capabilities, higher temporal and spatial resolutions to model the evolution of the power system with modern technologies and low-carbon pathways, ReEDS' model solution times have increased significantly from 4-6 hours in 2018 to 18-48+ hours in 2023 . Also, the model size for commonly-run ReEDS scenarios reached 22 and 28 million equations and variables, respectively. These runtimes can be especially challenging under certain scenario settings (e.g., very high temporal or spatial resolution) or with limited computational power. In this presentation, we will discuss several methods we used to improve model runtime, including data preparation, model modification, and solver tuning. The implementation of these methods shrank the model size to 7.2 and 7.3 million equations and variables, respectively. Furthermore, this led to a 77% reduction in the model's run time for commonly-run ReEDS scenarios. We will discuss the process of identifying areas for solve time improvements and how the specific enhancements for the ReEDS model might be applied to other similar large-scale models.

ENERGY PLANNING, POLICY, AND ECONOMY,MATHEMATICS A↗
DOE OSTI2022

Uncertainty in Atlantic Multidecadal Oscillation derived from different observed datasets and their possible causes

As a leading mode of sea surface temperature (SST) variability over the North Atlantic in both observations and model simulations, the Atlantic Multidecadal Oscillation (AMO) can have a substantial influence on regional and global climate. By using Low-Frequency Component Analysis, we explore the uncertainties of the resulting AMO indices and the corresponding spatial patterns derived from three observational SST datasets. We found that the known coherent spatial pattern of the AMO at the basin scale over the North Atlantic appears in two out of the three datasets. Further analysis indicates that both the warming trend and the different techniques used to construct these observed gridded SSTs contribute to the AMO’s spatial coherence over the North Atlantic, especially during periods of sparse data sampling. The SST in the Extended Reconstructed SST dataset version 5 (ERSSTv5), changes from being systematically below the other datasets during the dense sampling periods on either side of the Second World War (WWII), to systematically above the other datasets during WWII, thereby introducing an artificial 10–20-year variability that affects the AMO’s spatial coherence. This coherence in the AMO’s spatial pattern is also affected by bias adjustment in ERSSTv5 at relative cool (i.e., non-summer) seasons, and by the heterogeneous North Atlantic warming pattern. The different AMO patterns can induce the different effects of wind, surface heat fluxes, and then drive ocean circulation and its heat transport convergence, especially for some seasons. For AMO indices, both the different detrending methods and different observational data result in uncertainty for the period 1935–1950. Such SST uncertainty is important to detect the relative role of the atmosphere and ocean in shaping the AMO.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Nitrogen Fertilization Restructured Spatial Patterns of Soil Organic Carbon and Total Nitrogen in Switchgrass and Gamagrass Croplands in Tennessee USA

Nitrogen (N) fertilizers can potentially alter spatial distribution of soil organic carbon (SOC) and total nitrogen (TN) concentrations in croplands such as switchgrass (SG: Panicum virgatum L.) and gamagrass (GG: Tripsacum dactyloides L.), but it remains unclear whether these effects are the same between crops and under different rates of fertilization. 13 C and 15 N are two important proxy measures of soil biogeochemistry, but they were rarely examined as to their spatial distributions in soil. Based on a three-year long fertilization experiment in Middle Tennessee, USA, the top mineral horizon soils (0–15 cm) were collected using a spatially explicit design within two 15-m 2 plots under three fertilization treatments in SG and GG croplands. A total of 288 samples were collected based on 12 plots and 24 samples in each plot. The fertilization treatments were no N input (NN), low N input (LN: 84 kg N ha -1 in urea) and high N input (HN: 168 kg N ha -1 in urea). The SOC, TN, SOC/TN (C: N), δ 13 C and δ 15 N were quantified and their within-plot variations and spatial distributions were achieved via descriptive and geostatistical methods. Results showed that SG generally displayed 10~120% higher plot-level variations in all variables than GG, and the plot-level variations were 20~77% higher in NN plots than LN and HN plots in SG but they were comparable in unfertilized and fertilized plots in GG. Relative to NN, LN and HN showed more significant surface trends and spatial structures in SOC and TN in both croplands, and the fertilization effect appeared more pronounced in SG. Spatial patterns in C: N, δ 13 C and δ 15 N were comparable among different fertilization treatments in both croplands. The descending within-plot variations were also identified among variables (SOC > TN > δ 15 N > C: N > δ 13 C). This study demonstrated that N fertilizations generally reduced the plot-level variance and simultaneously re-established spatial structures of SOC and TN in bioenergy croplands, which little varied with fertilization rate but was more responsive in switchgrass cropland.

59 BASIC BIOLOGICAL SCIENCES↗