Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatiotemporal”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

DSCOVR/EPIC-derived global hourly and daily downward shortwave and photosynthetically active radiation data at 0.1° × 0.1° resolution

Downward shortwave radiation (SW) and photosynthetically active radiation (PAR) play crucial roles in Earth system dynamics. Spaceborne remote sensing techniques provide a unique means for mapping accurate spatiotemporally continuous SW–PAR, globally. However, any individual polar-orbiting or geostationary satellite cannot satisfy the desired high temporal resolution (sub-daily) and global coverage simultaneously, while integrating and fusing multisource data from complementary satellites/sensors is challenging because of co-registration, intercalibration, near real-time data delivery and the effects of discrepancies in orbital geometry. The Earth Polychromatic Imaging Camera (EPIC) on board the Deep Space Climate Observatory (DSCOVR), launched in February 2015, offers an unprecedented possibility to bridge the gap between high temporal resolution and global coverage and characterize the diurnal cycles of SW–PAR globally. In this study, we adopted a suite of well-validated data-driven machine-learning models to generate the first global land products of SW–PAR, from June 2015 to June 2019, based on DSCOVR/EPIC data. The derived products have high temporal resolution (hourly) and medium spatial resolution (0.1°×0.1°), and they include estimates of the direct and diffuse components of SW–PAR. We used independently widely distributed ground station data from the Baseline Surface Radiation Network (BSRN), the Surface Radiation Budget Network (SURFRAD), NOAA's Global Monitoring Division and the U.S. Department of Energy's Atmospheric System Research (ASR) program to evaluate the performance of our products, and we further analyzed and compared the spatiotemporal characteristics of the derived products with the benchmarking Clouds and the Earth's Radiant Energy System Synoptic (CERES) data. We found both the hourly and daily products to be consistent with ground-based observations (e.g., hourly and daily total SWs have low biases of -3.96 and -0.71 W m -2 and root-mean-square errors (RMSEs) of 103.50 and 35.40 W m -2 , respectively). The developed products capture the complex spatiotemporal patterns well and accurately track substantial diurnal, monthly, and seasonal variations in SW–PAR when compared to CERES data. They provide a reliable and valuable alternative for solar photovoltaic applications worldwide and can be used to improve our understanding of the diurnal and seasonal variabilities of the terrestrial water, carbon and energy fluxes at various spatial scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Using neural network ensembles to separate ocean biogeochemical and physical drivers of phytoplankton biogeography in Earth system models

Abstract. Earth system models (ESMs) are useful tools for predicting and understanding past and future aspects of the climate system. However, the biological and physical parameters used in ESMs can have wide variations in their estimates. Even small changes in these parameters can yield unexpected results without a clear explanation of how a particular outcome was reached. The standard method for estimating ESM sensitivity is to compare spatiotemporal distributions of variables from different runs of a single ESM. However, a potential pitfall of this method is that ESM output could match observational patterns because of compensating errors. For example, if a model predicts overly weak upwelling and low nutrient concentrations, it might compensate for this by allowing phytoplankton to have a high sensitivity to nutrients. Recently, we demonstrated that neural network ensembles (NNEs) are capable of extracting relationships between predictor and target variables within ocean biogeochemical models. Being able to view the relationships between variables, along with spatiotemporal distributions, allows for a more mechanistically based examination of ESM outputs. Here, we investigated whether we could apply NNEs to help us determine why different ESMs produce different spatiotemporal distributions of phytoplankton biomass. We tested this using three cases. The first and second case used different runs of the same ESM, except that the physical circulations differed between them in the first case, while the biological equations differed between them in the second. Our results indicated that the NNEs were capable of extracting the relationships between variables for different runs of a single ESM, allowing us to distinguish between differences due to changes in circulation (which do not change relationships) from changes in biogeochemical formulation (which do change relationships). In the third case, we applied NNEs to two different ESMs. The results of the third case highlighted the capability of NNEs to contrast the apparent relationships of different ESMs and some of the challenges it presents. Although applied specifically to the ocean components of an ESM, our study demonstrates that Earth system modelers can use NNEs to separate the contributions of different components of ESMs. Specifically, this allows modelers to compare the apparent relationships across different ESMs and observational datasets.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2006

First Directly Retrieved Global Distribution of Tropospheric Column Ozone from GOME: Comparison with the GEOS-CHEM Model

We present the first directly retrieved global distribution of tropospheric column ozone from Global Ozone Monitoring Experiment (GOME) ultraviolet measurements during December 1996 to November 1997. The retrievals clearly show signals due to convection, biomass burning, stratospheric influence, pollution, and transport. They are capable of capturing the spatiotemporal evolution of tropospheric column ozone in response to regional or short time-scale events such as the 1997-1998 El Nino event and a 10-20 DU change within a few days. The global distribution of tropospheric column ozone displays the well-known wave-1 pattern in the tropics, nearly zonal bands of enhanced tropospheric column ozone of 36-48 DU at 20degS-30degS during the austral spring and at 25degN-45degN during the boreal spring and summer, low tropospheric column ozone of <30 DU uniformly distributed south of 35 S during all seasons, and relatively high tropospheric column ozone of >33 DU at some northern high-latitudes during the spring. Simulation from a chemical transport model corroborates most of the above structures, with small biases of <+/-5 DU and consistent seasonal cycles in most regions, especially in the southern hemisphere. However, significant positive biases of 5-20 DU occur in some northern tropical and subtropical regions such as the Middle East during summer. Comparison of GOME with monthly-averaged Measurement of Ozone and Water Vapor by Airbus in-service Aircraft (MOZAIC) tropospheric column ozone for these regions usually shows good consistency within 1 a standard deviations and retrieval uncertainties. Some biases can be accounted for by inadequate sensitivity to lower tropospheric ozone, the different spatiotemporal sampling and the spatiotemporal variations in tropospheric column ozone.

Liu, Xiong↗
NASA NTRS2017

Scale Size-Dependent Characteristics of the Nightside Aurora

We have determined the spatiotemporal characteristics of the magnetosphere-ionosphere (M-I) coupling using auroral imaging. Observations at fixed positions for an extended period of time are provided by a ground-based all-sky imager measuring the 557.7 nanometer auroral emissions. We report on a single event of nightside aurora (at approximately 22 magnetic local time) preceding a substorm onset. To determine the spatiotemporal characteristics, we perform an innovative analysis of an all-sky imager movie (19 minutes duration, images at 3.31 hertz) that combines a two-dimensional spatial fast Fourier transform with a temporal correlation. We find a scale size-dependent variability where the largest scale sizes are stable on timescales of minutes while the small scale sizes are more variable. When comparing two smaller time intervals of different types of auroral displays, we find a variation in their characteristics. The characteristics averaged over the event are in remarkable agreement with the spatiotemporal characteristics of the nightside field-aligned currents during moderately disturbed times. Thus, two different electrodynamical parameters of the M-I coupling show similar behavior. This gives independent support to the claim of a system behavior that uses repeatable solutions to transfer energy and momentum from the magnetosphere to the ionosphere.

electrodynamical↗
NASA NTRS2021

Validation of IASI Satellite Ammonia Observations at the Pixel Scale Using In Situ Vertical Profiles

Satellite ammonia (NH 3 ) observations provide unprecedented insights into NH3 emissions, spatiotemporal variabilities and trends, but validation with in situ measurements remains lacking. Here, total columns from the Infrared Atmospheric Sounding Interferometer (IASI) were intercompared to boundary layer NH 3 profiles derived from aircraft- and surface-based measurements primarily in Colorado, USA, in the summer of 2014. IASI-NH 3 version 3 near real-time data set compared well to in situ derived columns (windows ±15 km around centroid, ±1 h around overpass time) with a correlation of 0.58, a slope of 0.78 ± 0.14 and an intercept of 2.1 × 10 15 ±1.5 × 10 15 molecules cm −2 . Agreement degrades at larger spatiotemporal windows, consistent with the short atmospheric lifetime of NH 3 . We also examined IASI version 3R data, which relies on temperature retrievals from the ERA Reanalysis, and a third product generated using aircraft-measured temperature profiles. The overall agreement improves slightly for both cases, and neither is biased within their combined measurement errors. Thus, spatiotemporal averaging of IASI over large windows can be used to reduce retrieval noise. Nonetheless, sampling artifacts of airborne NH 3 instruments result in significant uncertainties of the in situ-derived columns. For example, large validation differences exist between ascent and descent profiles, and the assumptions of the free tropospheric NH 3 profiles used above the aircraft ceiling significantly impact the validation. Because short-lived species like NH 3 largely reside within the boundary layer with complex vertical structures, more comprehensive validation is needed across a wide range of environments. More accurate and widespread in situ NH 3 data sets are therefore required for improved validations of satellite products.

Ammonia↗
NASA NTRSDate not supplied

How Well Satellite Remote Sensing Can Inform Surface Level Ozone Production Sensitivity to Precursor Trace Gas Emissions?

Surface-level ozone (O3) pollution mitigation strategies rely on a clear understanding of the atmospheric chemical processes involving O3, NOx (nitric oxide (NO) + nitrogen dioxide (NO2)) and volatile organic compounds (VOCs). A robust spatiotemporal classification of O3 photochemical regimes with relative chemical sensitivity of local O3 formation to emission reductions of NOx (NOx-limited regime) versus VOCs (radical-limited regime) is required to design anthropogenic emission reduction policies in order to improve surface air quality. Numerous studies found that the ratio of satellite-retrieved Vertical Column Densities (VCDs) of formaldehyde (HCHO) to NO2 can diagnose O3 production sensitivity. However, uncertainty still remains about the accuracy of such estimates using satellite-based retrievals. This study conducts an investigation to identify how well satellite data can inform surface-level O3 sensitivity by validating satellite-based estimates against aircraft-based retrievals taken during the Long Island Sound Tropospheric Ozone Study (LISTOS-2018) and Ozone Water-Land Environmental Transition Study (OWLETS-2) field campaigns. For this study, we use HCHO and NO2 retrievals from the Ozone Monitoring Instrument (OMI) onboard the Aura satellite, and high spatiotemporal resolution TROPOspheric Monitoring Instrument (TROPOMI) onboard the Sentinel-5 Precursor (S5P) satellite. The Community Multiscale Air Quality (CMAQ) modelling system is used to provide trace gas vertical profiles to recalculate air mass factor for deriving consistent satellite-based VCDs. To demonstrate the ability of satellites to diagnose O3 production regimes, we validate satellite-derived O3 sensitivity regimes, along with the OMI and TROPOMI products applied, using HCHO and NO2 retrievals from the aircraft-based instruments: Geostationary Trace gas and Aerosol Sensor Optimization (GeoTASO) spectrometer and Geostationary Coastal and Air Pollution Events Airborne Simulator (GCAS) instruments deployed during the LISTOS-2018 and OWLETS-2 campaigns. The results of this work will advance our fundamental knowledge of the spatiotemporal evolution of the complex O3-NOx-VOC chemistry, and will advance our understanding about the quality of remote-sensing, suborbital, and satellite data for accurately observing surface-level O3 production sensitivity.

Satellite↗
DOE OSTI2024

Temporal and Spatial Evolution of Non-Elastic Strain Accumulation in Stanstead Granite During Brittle Creep

Understanding the long-term behavior of brittle rocks requires fundamental consideration of time-dependent strain evolution and brittle creep processes. Previous studies have evaluated sub-critical crack growth during time-dependent deformation and damage evolution in brittle rocks; however, there is an incomplete knowledge of how damage evolves spatially and temporally within the body of intact rocks, where distributed regions of damage interact and coalesce during creep. This paper presents laboratory research focusing on evaluating brittle creep damage processes in Stanstead granite (SG) using 2-dimensional digital image correlation (2D-DIC). In the laboratory, the prismatic SG specimens were loaded beyond an estimated Crack Damage stress threshold (CD) level and then maintained a constant stress to initiate the creep process. DIC was used to characterize full-field spatiotemporal strain evolution, which was then interpreted in the context of local regions of “damage”, determined according to a strain-based criterion. Here, a method was proposed for identifying “existing” and “new” damage regions over specified intervals during the test, followed by spatial clustering of these regions to assess their spatiotemporal evolution. The clustering analysis results demonstrated the extension of existing damage regions was the main damage process during brittle creep, which is consistent with existing models of sub-critical crack growth. In addition, temporal analysis of tensile and shear strains on a point-by-point basis revealed both new damage formation and the strain concentration within existing damaged regions significantly contribute to overall specimen strain during primary creep. In contrast, during secondary creep, increases in specimen deformation are influenced by the accumulation of strains within already damaged regions.

58 GEOSCIENCES↗
DOE OSTI2020

Mechanoadaptive strain and functional osseointegration of dental implants in rats

Spatiotemporal implant-bone biomechanics and mechanoadaptive strains in peri-implant tissue are poorly understood. Physical and chemical characteristics of an implant-bone complex (IBC) were correlated in three-dimensional space (along the length and around a dental implant) to gather insights into time related integration of the implant with the cortical portion of a jaw bone in a rat. Rats (N = 9) were divided into three experimental groups with three rats per time point; 3-, 11-, and 24-day. All rats were fed crumbled hard pellets mixed with water (soft-food diet) for the first 3 days followed by a hard-food diet with intact hard-food pellets (groups of 11- and 24-day only). Biomechanics of the IBCs harvested from rats at each time point was evaluated by performing mechanical testing in situ in tandem with X-ray imaging. The effect of physical association (contact area) of a loaded implant with adapting peri-implant tissue, and resulting strain within was mapped by using digital volume correlation (DVC) technique. The IBC stiffness at respective time points was correlated with mechanical strain in peri-implant tissue. Results illustrated that IBC stiffness at 11-day was lower than that observed at 3-day. However, at 24-day, IBC stiffness recovered to that which was observed at 3-day. Furthermore, correlative microscopy and spectroscopy illustrated that the lower IBC stiffness was constituted by softer and less mineralized peri-implant tissue that contained varying expressions of osteoconductive elements. Lower IBC stiffness observed at 11-day was constituted by less mineralized peri-implant tissue with osteoconductive elements that included phosphorus (P) which was co-localized with higher expression of zinc (Zn), and lower expression of calcium (Ca). Higher IBC stiffness at 24-day was constituted by mineralized peri-implant tissue with higher expressions of osteoconductive elements including Ca and P, and lower expressions of Zn. These spatiotemporal correlative maps of peri-implant tissue architecture, heterogeneous distribution of mineral density, and elemental colocalization underscore mechanoadaptive physicochemical properties of peri-implant tissue that facilitate functional osseointegration of an implant. These results provided insights into 1) plausible "prescription" of mechanical loads as an osteoinductive "therapeutic dose" to encourage osteoconductive elements in the peri-implant tissue that would facilitate functional osseointegration of the implant; 2) a "critical temporal window" between 3 and 11 days, and perhaps it is this acute phase during which key candidate regenerative molecules can be harnessed to accelerate osseointegration of an implant under load.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Verifying infectious disease scenario planning for geographically diverse populations

In the face of the COVID-19 pandemic, the literature saw a spike in publications for epidemic models, and a renewed interest in capturing contact networks and geographic movement of populations. There remains a general lack of consensus in the modeling community around best practices for spatiotemporal epi-modeling, specifically as it pertains to the infection rate formulation and the underlying contact or mixing model. We mathematically verify several common modeling assumptions in the literature, to prove when certain choices can provide consistent results across different geographic resolutions, population densities and patterns, and mixing assumptions. The most common infection rate formulation, a computationally low cost per capita infection rate assumption, fails the consistency tests for heterogeneous populations and gravity-weighting assumptions. Future modeling efforts in spatiotemporal disease modeling should be wary of this limitation, particularly when working with more heterogeneous or sparse populations. Our results provide guidance for testing that a model preserves desirable properties even when model inputs mask potential problems due to symmetry or homogeneity. We also provide a recipe for performing this type of verification, strengthening decision support tools.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2019

Towards a Satellite Formaldehyde – in situ Hybrid Estimate for Organic Aerosol Abundance

Organic aerosol (OA) is one of the main components of the global particulate burden and intimately links natural and anthropogenic emissions with air quality and climate. It is challenging to accurately represent OA in global models. Direct quantification of global OA abundance is not possible with current remote sensing technology; however, it may be possible to exploit correlations of OA with remotely observable quantities to infer OA spatiotemporal distributions. In particular, formaldehyde (HCHO) and OA share common sources via both primary emissions and secondary production from oxidation of volatile organic compounds (VOCs). Here, we examine OA–HCHO correlations using data from summertime airborne campaigns investigating biogenic (NASA SEAC4RS and DC3), biomass burning (NASA SEAC4RS), and anthropogenic conditions (NOAA CalNex and NASA KORUS-AQ). In situ OA correlates well with HCHO (r=0.59–0.97), and the slope and intercept of this relationship depend on the chemical regime. For biogenic and anthropogenic regions, the OA–HCHO slopes are higher in low NOx conditions, because HCHO yields are lower and aerosol yields are likely higher. The OA–HCHO slope of wildfires is over 9 times higher than that for biogenic and anthropogenic sources. The OA–HCHO slope is higher for highly polluted anthropogenic sources (e.g., KORUS-AQ) than less polluted (e.g., CalNex) anthropogenic sources. Near-surface OAs over the continental US are estimated by combining the observed in situ relationships with HCHO column retrievals from NASA's Ozone Monitoring Instrument (OMI). HCHO vertical profiles used in OA estimates are from climatology a priori profiles in the OMI HCHO retrieval or output of specific period from a newer version of GEOS-Chem. Our OA estimates compare well with US EPA IMPROVE data obtained over summer months (e.g., slope =0.60–0.62, r=0.56 for August 2013), with correlation performance comparable to intensively validated GEOS-Chem (e.g., slope =0.57, r=0.56) with IMPROVE OA and superior to the satellite-derived total aerosol extinction (r=0.41) with IMPROVE OA. This indicates that OA estimates are not very sensitive to these HCHO vertical profiles and that a priori profiles from OMI HCHO retrieval have a similar performance to that of the newer model version in estimating OA. Improving the detection limit of satellite HCHO and expanding in situ airborne HCHO and OA coverage in future missions will improve the quality and spatiotemporal coverage of our OA estimates, potentially enabling constraints on global OA distribution.

anthropogenic emissions↗
DOE OSTI2022

Riverscape heterogeneity in estimated Chinook Salmon emergence phenology and implications for size and growth

Abstract Many salmonid‐bearing rivers exhibit thermal and hydrologic heterogeneity at multiple spatial and temporal scales, but how this translates into spatiotemporal patterns of fry emergence is poorly understood. Understanding this variability is important because emergence timing determines the biophysical conditions fish first experience (e.g., temperature, flow, and food supply), thereby influencing growth opportunities and survival during this critical life stage. We predicted spring Chinook Salmon ( Oncorhynchus tshawytscha ) emergence phenology across four northeastern Oregon subbasins over 5–9 years using empirical spawning and temperature data. We then related interannual emergence timing estimates to juvenile salmon size and growth rates at consistent sampling locations. There were clear longitudinal patterns of predicted emergence timing in each subbasin: The shape of these patterns was consistent among years, but not among subbasins. In two subbasins, emergence occurred progressively later with distance upstream, whereas in the other two subbasins emergence was earliest at upstream sites. Within each year, median emergence dates among sites within each subbasin ranged between 44 and 58 days. This spatial variation was comparable to interannual variation, with median emergence dates for a given location in each subbasin ranging between 47 and 74 days among years. Contrary to our expectations, juvenile salmon were not larger in years with earlier emergence, owing to slower estimated spring and summer growth rates compared to years with later emergence. Despite large interannual variation in estimated emergence dates, these results suggest that other factors (e.g., stream flow, temperature, and density‐dependence) were more important than growth duration in determining juvenile salmon growth rates and size among years. We demonstrated considerable spatial and interannual variation in emergence phenology within these subbasins. Understanding how this variation translates to spatiotemporal patterns of juvenile salmon habitat use, growth, and survival has important implications for guiding restoration efforts and understanding how climate change may impact these populations.

Kaylor, Matthew J.↗
DOE OSTI2022

Stream and ocean hydrodynamics mediate partial migration strategies in an amphidromous Hawaiian goby

Abstract Partial migration strategies, in which some individuals migrate but others do not, are widely observed in populations of migratory animals. Such patterns could arise via variation in migratory behaviors made by individual animals, via genetic variation in migratory predisposition, or simply by variation in migration opportunities mediated by environmental conditions. Here we use spatiotemporal variation in partial migration across populations of an amphidromous Hawaiian goby to test whether stream or ocean conditions favor completing its life cycle entirely within freshwater streams rather than undergoing an oceanic larval migration. Across 35 watersheds, microchemical analysis of otoliths revealed that most adult Awaous stamineus were freshwater residents (62% of n = 316 in 2009, 83% of n = 274 in 2011), but we found considerable variation among watersheds. We then tested the hypothesis that the prevalence of freshwater residency increases with the stability of stream flows and decreases with the availability of dispersal pathways arising from ocean hydrodynamics. We found that streams with low variation of daily discharge were home to a higher incidence of freshwater residents in each survey year. The magnitude of the shift in freshwater residency between survey years was positively associated with predicted interannual variability in the success of larval settlement in streams on each island based on passive drift in ocean currents. We built on these findings by developing a theoretical model of goby life history to further evaluate whether mediation of migration outcomes by stream and ocean hydrodynamics could be sufficient to explain the range of partial migration frequency observed across populations. The model illustrates that the proportion of larvae entering the ocean and differential survival of freshwater‐resident versus ocean‐going larvae are plausible mechanisms for range‐wide shifts in migration strategies. Thus, we propose that hydrologic variation in both ocean and stream environments contributes to spatiotemporal variation in the prevalence of migration phenotypes in A. stamineus . Our empirical and theoretical results suggest that the capacity for partial migration could enhance the persistence of metapopulations of diadromous fish when confronted with variable ocean and stream conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Characterizing climate pathways using feature importance on echo state networks

The 2022 National Defense Strategy of the United States listed climate change as a serious threat to national security. Climate intervention methods, such as stratospheric aerosol injection, have been proposed as mitigation strategies, but the downstream effects of such actions on a complex climate system are not well understood. The development of algorithmic techniques for quantifying relationships between source and impact variables related to a climate event (i.e., a climate pathway) would help inform policy decisions. Data-driven deep learning models have become powerful tools for modeling highly nonlinear relationships and may provide a route to characterize climate variable relationships. In this paper, we explore the use of an echo state network (ESN) for characterizing climate pathways. ESNs are a computationally efficient neural network variation designed for temporal data, and recent work proposes ESNs as a useful tool for forecasting spatiotemporal climate data. However, ESNs are noninterpretable black-box models along with other neural networks. The lack of model transparency poses a hurdle for understanding variable relationships. We address this issue by developing feature importance methods for ESNs in the context of spatiotemporal data to quantify variable relationships captured by the model. We conduct a simulation study to assess and compare the feature importance techniques, and we demonstrate the approach on reanalysis climate data. In the climate application, we consider a time period that includes the 1991 volcanic eruption of Mount Pinatubo. This event was a significant stratospheric aerosol injection, which acts as a proxy for an anthropogenic stratospheric aerosol injection. Furthermore, we are able to use the proposed approach to characterize relationships between pathway variables associated with this event that agree with relationships previously identified by climate scientists.

black-box models↗
DOE OSTI2020

Modeling regional precipitation over the Indus River basin of Pakistan using statistical downscaling

Complex processes govern spatiotemporal distribution of precipitation within the high-mountainous headwater regions (commonly known as the upper Indus basin (UIB)), of the Indus River basin of Pakistan. Reliable precipitation simulations particularly over the UIB present a major scientific challenge due to regional complexity and inadequate observational coverage. Here, we present a statistical downscaling approach to model observed precipitation of the entire Indus basin, with a focus on UIB within available data constraints. Taking advantage of recent high altitude (HA) observatories, we perform precipitation regionalization using K-means cluster analysis to demonstrate effectiveness of low-altitude stations to provide useful precipitation inferences over more uncertain and hydrologically important HA of the UIB. We further employ generalized linear models (GLM) with gamma and Tweedie distributions to identify major dynamic and thermodynamic drivers from a reanalysis dataset within a robust cross-validation framework that explain observed spatiotemporal precipitation patterns across the Indus basin. Final statistical models demonstrate higher predictability to resolve precipitation variability over wetter southern Himalayans and different lower Indus regions, by mainly using different dynamic predictors. The modeling framework also shows an adequate performance over more complex and uncertain trans-Himalayans and the northwestern regions of the UIB, particularly during the seasons dominated by the westerly circulations. However, the cryosphere-dominated trans-Himalayan regions, which largely govern the basin hydrology, require relatively complex models that contain dynamic and thermodynamic circulations. Furthermore, we also analyzed relevant atmospheric circulations during precipitation anomalies over the UIB, to evaluate physical consistency of the statistical models, as an additional measure of reliability. Overall, our results suggest that such circulation-based statistical downscaling has the potential to improve our understanding towards distinct features of the regional-scale precipitation across the upper and lower Indus basin. Additionally, such understanding should help to assess the response of this complex, data-scarce, and climate-sensitive river basin amid future climatic changes, to serve communal and scientific interests.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Proximity-based proteomics reveals the thylakoid lumen proteome in the cyanobacterium Synechococcus sp. PCC 7002

Abstract Cyanobacteria possess unique intracellular organization. Many proteomic studies have examined different features of cyanobacteria to learn about the intracellular structures and their respective functions. While these studies have made great progress in understanding cyanobacterial physiology, the conventional fractionation methods used to purify cellular structures have limitations; specifically, certain regions of cells cannot be purified with existing fractionation methods. Proximity-based proteomics techniques were developed to overcome the limitations of biochemical fractionation for proteomics. Proximity-based proteomics relies on spatiotemporal protein labeling followed by mass spectrometry of the labeled proteins to determine the proteome of the region of interest. We performed proximity-based proteomics in the cyanobacterium Synechococcus sp. PCC 7002 with the APEX2 enzyme, an engineered ascorbate peroxidase. We determined the proteome of the thylakoid lumen, a region of the cell that has remained challenging to study with existing methods, using a translational fusion between APEX2 and PsbU, a lumenal subunit of photosystem II. Our results demonstrate the power of APEX2 as a tool to study the cell biology of intracellular features and processes, including photosystem II assembly in cyanobacteria, with enhanced spatiotemporal resolution.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Geophysical methods reveal the soil architecture and subsurface stratigraphic heterogeneities across land-lake interfaces along Lake Erie

Abstract Purpose The land-lake interface is a unique zone where terrestrial and aquatic ecosystems meet, forming part of the Earth’s most geochemically and biologically active zones. The unique characteristics of this interface are yet to be properly understood due to the inherently high spatiotemporal variability of subsurface properties, which are difficult to capture with the traditional soil sampling methods. Geophysical methods offer non-invasive techniques to capture variabilities in soil properties at a high resolution across various spatiotemporal scales. Methods We combined electromagnetic induction (EMI), electrical resistivity tomography (ERT), and ground penetrating radar (GPR) with data from soil cores and in situ sensors to investigate hydrostratigraphic heterogeneities across land-lake interfaces along the western basin of Lake Erie. Results EMI revealed high spatial heterogeneities in ECa distribution across the land-lake interfaces, with higher values in the wetland and transition zones compared to the upland zone. Soil ECa maps matched soil maps from a public database with the hydric soil units delineated as high conductivity zones (ECa > 40 mS/m). ERT and GPR showed vertical variation in soil properties with clear stratigraphic boundaries, and correlation of ERT profiles with lithologs from piezometers revealed the stratigraphic units of silt–clay and till sequence down to 3.5 m depth which are consistent with the surficial geology of the study area. Conclusions These results validate the use of multiple geophysical methods for extrapolating soil properties and mapping stratigraphic structures at land-lake interfaces, thereby providing the missing information required to improve the earth system model (ESM) of coastal interfaces.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Scaling kinetic Monte-Carlo simulations of grain growth with combined convolutional and graph neural networks

Graph neural networks (GNN) have emerged as a promising machine learning method for microstructure simulations such as grain growth. However, accurate modeling of realistic grain boundary networks requires large simulation cells, which GNN has difficulty scaling up to. To alleviate the computational costs and memory footprint of GNN, we suggest a hybrid architecture combining a convolutional neural network (CNN) based bijective autoencoder to compress the spatial dimensions, and a GNN that evolves the microstructure in the latent space of reduced spatial sizes. Our results demonstrate that the new design significantly reduces computational costs with using fewer message passing layer (from 12 down to 3) compared with GNN alone. The reduction in computational cost becomes more pronounced as the spatial size increases, indicating strong computational scalability. For the largest mesh evaluated (160 3 ), our method reduces memory usage and runtime in inference by 117× and 115×, respectively, compared with GNN-only baseline. More importantly, it shows higher accuracy and stronger spatiotemporal capability than the GNN-only baseline, especially in long-term testing. Such combination of scalability and accuracy is essential for simulating realistic material microstructures over extended time scales. The improvements can be attributed to the bijective autoencoder’s ability to compress information losslessly from spatial domain into a high dimensional feature space, thereby producing more expressive latent features for the GNN to learn from, while also contributing its own spatiotemporal modeling capability. Training data are generated from stochastic grain growth simulations, providing realistic variability for learning robust microstructure evolution. Comprehensive system validation confirms that the model is accurate, robust, and scalable.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Can upscaling ground nadir SIF to eddy covariance footprint improve the relationship between SIF and GPP in croplands?

Ground solar-induced chlorophyll fluorescence (SIF) is important for the mechanistic understanding of the dynamics of vegetation gross primary production (GPP) at fine spatiotemporal scales. However, eddy covariance (EC) observations generally cover larger footprint areas than ground SIF observations (a bare fiber with nadir), and this footprint mismatch between nadir SIF and GPP could complicate the canopy SIF-GPP relationships. Here, we upscaled nadir SIF observations to EC footprint and investigated the change in SIF-GPP relationships after the upscaling in cropland. We included 13 site-years data in our study, with seven site-years corn, four site-years soybeans, and two site-years miscanthus, all located in the US Corn Belt. All sites’ crop nadir SIF observations collected from the automated FluoSpec2 system (a hemispheric-nadir system) were upscaled to the GPP footprint-based SIF using vegetation indices (VIs) calculated from high spatiotemporal satellite reflectance data. We found that SIF-GPP relationships were not substantially changed after upscaling nadir SIF to GPP footprint at our crop sites planted with corn, soybean, and miscanthus, with R 2 change after the upscaling ranging from -0.007 to 0.051 and root mean square error (RMSE) difference from -0.658 to 0.095 umol m –2 s –1 relative to original nadir SIF-GPP relationships across all the site-years. The variation of the SIF-GPP relationship within each species across different site-years was similar between the original nadir SIF and upscaled SIF. Different VIs, EC footprint models, and satellite data led to marginal differences in the SIF-GPP relationships when upscaling nadir SIF to EC footprint. Furthermore, our study provided a methodological framework to correct this spatial mismatch between ground nadir SIF and GPP observations for croplands and potentially for other ecosystems. Our results also demonstrated that the spatial mismatch between ground nadir SIF and GPP might not significantly affect the SIF-GPP relationship in cropland that are largely homogeneous.

54 ENVIRONMENTAL SCIENCES↗