Search NASA⌕ Search

SEARCH · Search NASA

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Simulation of PV Variability as a Function of PV Generation and Plant Size

The deployment of photovoltaic (PV) systems continues to show significant expansion; however, this growth has brought added attention to issues around the variability of the solar resource. Both spatial and temporal variability exist. Temporal scales can range from the sub-second to multiyear, whereas spatial scales can range from a few meters to tens of kilometers. There are multiple methods described in the literature to quantify PV variability at various spatial and temporal scales. This study focuses on short-term temporal variability and uses similar approaches with the addition of PV plant size a parameter to quantify variability. The method employed here incorporates the normalization of clear-and cloudy-sky conditions and PV plant size to quantify nominal variability metrics. The distribution and fluctuations of these metrics provide relevant information that is useful for system operations. The National Solar Radiation Database (NSRDB) is used to simulate PV variability as a function of PV generation and plant size. Hypothetical but realistic system information at 33 locations is used to model PV generation by feeding NSRDB solar irradiance data to the National Renewable Energy Laboratory’s System Advisor Model (SAM). Over the selected region, it is found that the aggregated ramp rates for the 1-minute data are associated with standard deviations ranging from 0.002–0.055 on a daily basis; however, hourly intervals induce higher aggregated ramp rates than the other timescales. Even though minute-to-minute variations are significant for the 1-minute time-scale, the standard deviation aggregated into a daily metric is smaller because of the cancellation of values.

irradiance↗
DOE OSTI2021

Simulation of PV Variability as a Function of PV Generation and Plant Size: Preprint

The deployment of photovoltaic (PV) systems continues to show significant expansion; however, this growth has brought added attention to issues around the variability of the solar resource. Both spatial and temporal variability exist. Temporal scales can range from the sub-second to multiyear, whereas spatial scales can range from a few meters to tens of kilometers. There are multiple methods described in the literature to quantify PV variability at various spatial and temporal scales. This study focuses on short-term temporal variability and uses similar approaches with the addition of PV plant size a parameter to quantify variability. The method employed here incorporates the normalization of clear- and cloudy-sky conditions and PV plant size to quantify nominal variability metrics. The distribution and fluctuations of these metrics provide relevant information that is useful for system operations. The National Solar Radiation Database (NSRDB) is used to simulate PV variability as a function of PV generation and plant size. Hypothetical but realistic system information at 33 locations is used to model PV generation by feeding NSRDB solar irradiance data to the National Renewable Energy Laboratory’s System Advisor Model (SAM). Over the selected region, it is found that the aggregated ramp rates for the 1-minute data are associated with standard deviations ranging from 0.002–0.055 on a daily basis; however, hourly intervals induce higher aggregated ramp rates than the other timescales. Even though minute-to-minute variations are significant for the 1-minute timescale, the standard deviation aggregated into a daily metric is smaller because of the cancellation of values.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2021

An Automated Method for Mapping Giant Kelp Canopy Dynamics from UAV

Satellite and aerial imagery have been used extensively for mapping the abundance and distribution of giant kelp ( Macrocystis pyrifera ) in southern California. There is now great potential for using unoccupied aerial vehicles (UAVs) to map kelp canopy at very high resolutions. However, tides and currents have been shown to affect the amount of floating kelp canopy on the water surface, and the impacts of these processes on remotely sensed kelp estimates in this region have not been fully quantified. UAVs were used to map fine-scale changes in canopy area due to tidal height and current speed at kelp forests off the coast of Palos Verdes, CA and Santa Barbara, CA. An automated method for detecting kelp canopy was developed that was 67% accurate using red-green-blue (RGB) UAV imagery and 93% accurate using multispectral UAV imagery across a range of weather, ocean, and illumination conditions. Increases in tidal height of 1 m reduced the amount of floating kelp canopy by 15% in Santa Barbara and by over 30% in Palos Verdes. The effect of current speed on visible kelp canopy was inconclusive, but there was a trend towards lower canopy area with increased current speed. Therefore, while tidal height and current speed can introduce significant variability to estimates of kelp abundance, the magnitude of this variability is site specific. Still, UAVs are a valuable tool for mapping of kelp canopy and can provide greater spatial resolution and temporal coverage than is possible from many satellite sensors. This data can provide insight into the patterns and drivers of high frequency fluctuations in kelp abundance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Beryllium Isotopes in Rainfall for Stratospheric-Tropospheric Exchange Dynamics Experiment (BIRDIE) (Field Campaign Report)

The Beryllium Isotopes in Rainfall for Stratospheric-Tropospheric Exchange Dynamics Experiment or BIRDIE was a precipitation collection campaign designed to measure cosmogenic isotopes of beryllium in rainfall and other forms of precipitation across the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s Southern Great Plains (SGP) atmospheric observatory. Precipitation samples were collected over three weeks (May 21-June 10, 2022) at 11 SGP sites across a 120-km 2 domain in northern and central Oklahoma. The field campaign targeted summertime convection rain events. We aimed to capture any spatial and/or temporal variability in the concentrations of beryllium-10 (10Be) and beryllium-7 (7Be) across the ARM SGP site during these rainfall events. The ratio of 10Be to 7Be (10/7Be) has been shown to be a powerful tracer for transport processes in the stratosphere-troposphere exchange (STE) system (e.g., Koch and Rind 1998). Any observed variability is analyzed to see if it could be linked to STE processes, including overshooting of deep convection clouds that penetrate the lower stratosphere. Six wet weather system events were collected within our three-week span. Five out of six of these rainfall events are classified as moderately to strongly convective. Those events are studied in more detail. The five systems also follow each other over a short period of 11 days. In our rainfall samples we observed spatial and event variability in the concentrations of beryllium-10 and -7 that resulted in variability of 10/7Be ratios. Each weather event appeared to exhibit a unique 10/7Be pattern across the domain that was not directly linked to storm total precipitation amount. To better understand the causes for this spatial variability, cloud characteristics in the domain were looked at to assess cloud type, cloud height, and depth information. We also ran Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model back trajectories to identify possible air mass and moisture origins across our events. Detailed analysis is ongoing. The purpose of this technical report is to provide background on using cosmogenic isotopes as atmospheric tracers, describe the experimental design, and present the beryllium concentration and ratio data. Lastly, we note that BIRDIE coincided with two other deep convection field campaigns in the area: The National Aeronautics and Space Administration (NASA)’s Dynamics and Chemistry of the Summer Stratosphere (DCOTSS) experiment (https://dcotss.org/) centered in Salina, Kansas, and the Targeted Observation by Radars and (uncrewed aerial systems) UAS of Supercells (TORUS) experiment (https://catalog.eol.ucar.edu/torus_2022) across the U.S. Central Plains. Additionally, a Lawrence Livermore National Laboratory ( LLNL) team at DCOTSS attempted to measure 10Be/7Be ratios on aerosols collected in situ in the lower and mid-stratosphere over Kansas using high-altitude balloons. Precipitation samples were also collected on the ground in Salina, Kansas.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Beryllium Isotopes in Rainfall for Stratospheric-Tropospheric Exchange Dynamics Experiment (BIRDIE) (Final Report)

The Beryllium Isotopes in Rainfall for Stratospheric-Tropospheric Exchange Dynamics Experiment or BIRDIE was a precipitation collection campaign designed to measure cosmogenic isotopes of beryllium in rainfall and other forms of precipitation across the DOE Southern Great Plains (SGP) Atmospheric Radiation Measurement (ARM) site. Precipitation samples were collected over three weeks (May 21-June 10, 2022) at 11 SGP sites across a 120 km 2 domain in northern and central Oklahoma. The field campaign targeted summer-time convection rain events. We aimed to capture any spatial and/or temporal variability in the concentrations of beryllium-10 (10Be) and beryllium-7 (7Be) across the SGP ARM site during these rainfall events. The ratio of 10Be to 7Be (10/7Be) has been shown to be a powerful tracer for transport processes in the stratosphere-troposphere exchange (STE) system (e.g., Koch and Rind 1998). Any observed variability is analyzed to see if it could be linked to STE processes, including overshooting of deep convection clouds which penetrate the lower stratosphere. Six wet weather system events were collected within our three-week span. Five out of six of these rainfall events are classified as moderately to strongly convective. Those events are studied in more detail. The five systems also follow each other over a short time period of eleven days. In our rainfall samples we observed spatial and event variability in the concentrations of beryllium-10 and -7 which resulted in variability of 10/7 Be ratios. Each weather event appeared to exhibit a unique 10/7Be pattern across the domain that was not directly linked to storm total precipitation amount. To better understand the causes for this spatial variability, cloud characteristics in the domain were looked at to assess cloud type, cloud height and depth information. We also ran Hysplit back trajectories to identify possible air mass and moisture origins across our events. Detailed analysis is still on-going. The purpose of this technical report is to provide background on using cosmogenic isotopes as atmospheric tracers, describe the experimental design, and present the beryllium concentration and ratio data.

58 GEOSCIENCES↗
DOE OSTI2023

Machine Learning of Key Variables Impacting Extreme Precipitation in Various Regions of the Contiguous United States

Abstract Amplification in extreme precipitation intensity and frequency can cause severe flooding and impose significant social and economic consequences. Variations in extreme precipitation intensity, frequencies, and return periods can be attributed to many physical variables across spatial and temporal scales. Here we employ ensemble machine learning (ML) methods, namely random forest (RF), eXtreme Gradient Boosting (XGB), and artificial neural networks (ANN), to explore key contributing variables to monthly extreme precipitation intensity and frequency in six regions over the United States. We further establish emulators for return periods. Results show that the ML models for intensity perform better in regions with obvious seasonality (i.e., Northern Great Plains, Southern Great Plains, and West Coast) than the other three regions (Northeast, Southwest, and Rocky Mountains), while for frequency the models perform well for most regions. The Shapley additive explanation is used to help explain the relationships between extreme precipitation characteristics and identify top variables for RF and XGB. We find that latent heat flux, relative humidity, soil moisture, and large‐scale subsidence are key common variables across the regions for both monthly intensity and frequency, and their compound effects are non‐negligible. The developed ML models capture the probability and return period of extreme precipitation well for all regions and may be used for decision making (e.g., infrastructure planning and design).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

High resolution variability in wet deposition in the southeastern United States

Rainwater chemistry is determined by atmospheric pollutants and particles which vary spatially and temporally. Industrial and agricultural activities and meteorological events (e.g. sea breezes, severe weather, blowing dust) alter atmospheric particle and trace gas compositions. These gases and particles are scavenged by cloud and rain droplets that drive wet deposition. During an Intensive Operation Period (IOP) from April to October 2021, rainwater was collected at higher frequency intervals, usually daily, after precipitation events at three locations on the Savannah River Site (SRS). The farthest locations were separated by approximately 20 km. The mean concentration (μeq/L) of seven ions followed the Cl⁻ > SO 4 2− > Na⁺ > NO 3 ⁻ > K⁺ > Mg 2+ > Ca 2+ downward trend. Ion concentrations were compared to background ion concentrations from the National Atmospheric Deposition Program (NADP). The high frequency monthly averaged SRS data compared well with the monthly averaged NADP background but demonstrated extensive variability. In some months in 2021, the high frequency data compared better to the NADP site near the coast while in other months inland sites compared better. Strong spatial variability for ion concentrations was observed across SRS which was attributed to localized impacts in rainfall spatial variability. High frequency measurements allowed for comparison to regional weather patterns indicating influences from the Atlantic Ocean, Gulf of Mexico, and cities. This can account for spatial variability in the wet deposition flux. Sea breezes, Saharan dust, and anthropogenic sources were shown to impact wet deposition flux variability. Higher frequency precipitation chemistry sampling at numerous locations better captures ion concentration variability and improves measurement representativeness.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

The biogeography of soil and airborne fungi in the Southwestern USA in relation to climate and vegetation

To assess how fungal dispersal might respond to climate change, we examined how climate and geography influence the regional distribution of fungi in soil and air. Specifically, we hypothesized that neighboring fungal communities should be more similar than distant communities (i.e. spatially autocorrelated) and that fungal dispersal should be more limited in soil than in air. We collected soil and air samples from 60 sites across five states in the Southwestern USA. Then, we sequenced the ITS2 region to identify fungal taxa in each sample. Next, we used distance-based redundancy analysis to partition variation in fungal community composition between climate variables and spatial structure. Fungi were indeed spatially autocorrelated. Moreover, precipitation, maximum vapor pressure deficit, and soil moisture were significantly related to fungal community composition in soils. In comparison, only precipitation was significantly related to community composition in the air. After accounting for climate, the strength of spatial autocorrelation did not differ significantly in soilborne versus airborne fungi. Dispersal limitation was evident in soilborne fungi at short distances (<100 km) and was not observed at any distance in airborne fungi. Altogether, climate may influence which fungal taxa are present in soil and air, and fungi could feasibly wind disperse over regional scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A Probabilistic Model for Global EMIC Wave Activity Using Van Allen Probes Observations

Electromagnetic ion cyclotron (EMIC) waves play a key role in radiation belt dynamics through resonant interactions. However, their low occurrence probability, high variability, and spatial intermittency pose challenges for accurate modeling. In this study, we present a machine learning (ML)-based global EMIC wave model built on the entire data set from the Van Allen Probes mission. To capture the distinct statistical characteristics of wave occurrence and amplitude, the model is separated into two modules: an occurrence model trained using ML techniques, and a wave amplitude model sampled from observed probability distributions. The input parameters are limited to real-time or predictable variables to ensure practical applicability. Our model shows strong performance across the entire test set and demonstrates improved predictive capability over a baseline random occurrence model, particularly during quiet geomagnetic conditions. Evaluation during both quiet and active periods confirms the model's ability to represent the clustered and intermittent nature of EMIC wave activity. Furthermore, the model provides global estimates of wave power, enabling integration with radiation belt electron data and showing signatures consistent with wave-induced scattering. We found a good correlation between the global wave activity from the model and relativistic electron observation by Van Allen Probes, regardless of the availability of in situ wave observations. The modular structure of the model also allows for straightforward expansion for additional wave properties, such as wave frequency, which can be modeled independently. This flexible, event-sensitive approach offers a promising framework for data-driven radiation belt simulations and space weather applications.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Data from: Understanding the biogeochemical and spatial drivers of methane and carbon dioxide fluxes in a large temperate reservoir

This dataset contains spatially resolved measurements of CO₂ and CH₄ fluxes and associated environmental variables collected across 200 sites in Douglas Reservoir (Tennessee, USA) between July 29-August 2, 2024. Measurements include diffusive fluxes of CO₂ and CH₄, CH₄ ebullition, and biogeochemical and spatial variables such as dissolved oxygen, temperature, conductivity, chlorophyll-a, pH, water depth, and distance from the dam. Sampling was conducted using a spatially balanced design to capture longitudinal and depth-related gradients throughout the reservoir. The dataset is structured to support analyses of spatial variability, flux pathway comparisons, and modeling approaches (e.g., spatial statistics and machine learning) aimed at understanding controls on reservoir CO₂ and CH₄ fluxes and improving upscaling to whole-reservoir and regional estimates.

Neeper, Jamie [ORNL] (ORCID:0009000842342101)↗
DOE OSTI2023

Robust global detection of forced changes in mean and extreme precipitation despite observational disagreement on the magnitude of change

Detection and attribution (D&A) of forced precipitation change are challenging due to internal variability, limited spatial, and temporal coverage of observational records and model uncertainty. These factors result in a low signal-to-noise ratio of potential regional and even global trends. Here, we use a statistical method – ridge regression – to create physically interpretable fingerprints for the detection of forced changes in mean and extreme precipitation with a high signal-to-noise ratio. The fingerprints are constructed using Coupled Model Intercomparison Project phase 6 (CMIP6) multi-model output masked to match coverage of three gridded precipitation observational datasets – GHCNDEX, HadEX3, and GPCC – and are then applied to these observational datasets to assess the degree of forced change detectable in the real-world climate in the period 1951–2020. We show that the signature of forced change is detected in all three observational datasets for global metrics of mean and extreme precipitation. Forced changes are still detectable from changes in the spatial patterns of precipitation even if the global mean trend is removed from the data. This shows the detection of forced change in mean and extreme precipitation beyond a global mean trend is robust and increases confidence in the detection method's power as well as in climate models' ability to capture the relevant processes that contribute to large-scale patterns of change. We also find, however, that detectability depends on the observational dataset used. Not only coverage differences but also observational uncertainty contribute to dataset disagreement, exemplified by the times of emergence of forced change from internal variability ranging from 1998 to 2004 among datasets. Furthermore, different choices for the period over which the forced trend is computed result in different levels of agreement between observations and model projections. These sensitivities may explain apparent contradictions in recent studies on whether models under- or overestimate the observed forced increase in mean and extreme precipitation. Lastly, the detection fingerprints are found to rely primarily on the signal in the extratropical Northern Hemisphere, which is at least partly due to observational coverage but potentially also due to the presence of a more robust signal in the Northern Hemisphere in general.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

MSW Variability Mapping and Conversion to Biofuel

MSW (Municipal Solid Waste) is a form of biomass which consists of categorized components of waste/trash. The general categories are paper, yard trash, construction & debris, appliances, tires, glass, metals, aluminum & steel cans, plastics, organics, inorganics, and HHW (Household Hazardous Waste). This project focuses on the factors within a region or population that contribute to variability in the composition of MSW and in turn MSW’s convertibility to biofuel. A list of contributors was determined (Social Vulnerability Index, Access to Public Transportation, Racial Distribution, GDP, Personal Income) and then JMP was used to perform a Multivariate analysis to determine correlations and a Partial Least-Squares regression to determine Variable Importance Plots for each MSW category. In addition to data analysis, the convertibility of MSW to biofuel was studied via microwave pyrolysis system in order to separate and characterize the various gaseous and bio-oil products.

09 BIOMASS FUELS↗
DOE OSTI2020

EQSIM—A multidisciplinary framework for fault-to-structure earthquake simulations on exascale computers, part II: Regional simulations of building response

The existing observational database of the regional-scale distribution of strong ground motions and measured building response for major earthquakes continues to be quite sparse. As a result, details of the regional variability and spatial distribution of ground motions, and the corresponding distribution of risk to buildings and other infrastructure, are not comprehensively understood. Utilizing high-performance computing platforms, emerging high-resolution, physics-based ground motion simulations can now resolve frequencies of engineering interest and provide detailed synthetic ground motions at high spatial density. This provides an opportunity for new insight into the distribution of infrastructure seismic demands and risk. In the work presented herein, the EQSIM fault-to-structure computational framework described in a companion paper, McCallen et al., is employed to investigate the regional-scale response of buildings to large earthquakes. A representative M = 7.0 strike-slip event is used to explore the distribution and amplitude of building demand, and comparisons are made between building response computed with fault-to-structure simulations and building response computed with existing measured near-fault earthquake records. New information on the distribution and variability of building response from high-performance parallel simulations is described and analyzed, and favorable first comparisons between building response predicted with both fault-to-structure simulations and real ground motions records are presented.

58 GEOSCIENCES↗
DOE OSTI2024

Printed Potentiometric Ammonium Sensors for Agriculture Applications

Ammonium (NH 4 + ) concentration is critical to both nutrient availability and nitrogen (N) loss in soil ecosystems but can be highly variable across spatial and temporal scales. For this reason, effectively informing agricultural practices such as fertilizer management and understanding of mechanisms of soil N loss require sensor technologies to monitor ammonium concentrations in real time. Our work investigates the performance of fully printed ammonium ion-selective sensors used in diverse soil environments. Ammonium sensors consisting of a printed ammonium ion-selective electrode and a printed Ag/AgCl reference were fabricated and characterized in aqueous solutions and three different soil types (sand, peat, and clay) under the range of ion concentrations likely to be present in soil (0.01–100 mM). The response of ammonium sensors was further evaluated under variable gravimetric moisture content in the soil to reflect their reliability under field conditions. Ammonium sensors demonstrated a sensitivity of 53.6 ± 5.1 mV/decade when tested in aqueous solution, and a sensitivity of 55.7 ± 11 mV/dec, 57.5 ± 4.1 mV/dec, and 43.7 ± 4 mV/dec was measured in sand, clay, and peat soils, respectively.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Simulation of high-frequency dissolved oxygen dynamics in a shallow estuary, the Corsica River, Chesapeake Bay

Understanding shallow water biogeochemical dynamics is a challenge in coastal regions, due to the presence of highly variable land-water interface fluxes, tight coupling with sediment processes, tidal dynamics, and diurnal variability in biogeochemical processes. While the deployment of continuous monitoring devices has improved our understanding of high-frequency (12 - 24 hours) variability and spatial heterogeneity in shallow regions, mechanistic modeling of these dynamics has lagged behind conceptual and empirical models. The inherent complexity of shallow water systems is represented in the Corsica River estuary, a small basin within the Chesapeake Bay ecosystem, where abundant monitoring data have been collected from long-term monitoring stations, continuous monitoring sensors, synoptic sensor surveys, and measurements of sediment-water fluxes. A state-of-the-art modeling system, the Semi-implicit Cross-scale Hydroscience Integrated System Model (SCHISM), was applied to the Corsica domain with a high-resolution grid and nutrient loads from the most recent version of the Chesapeake Bay watershed model. The Corsica SCHISM model reproduced observed high-frequency variability in dissolved oxygen, as well as seasonal variability in chlorophyll-a and sediment-water fluxes. Time-series signal analyses using Empirical Model Decomposition and spectral analysis revealed that the diurnal and M2 tide frequencies are the dominant high-frequency modes and physical transport contributes a larger share to dissolved oxygen budgets than biogeochemical processes on an hourly time scale. Heterogeneity and patchiness in dissolved oxygen resulting from phytoplankton distributions and geometry-driven eddies amplify the physical transport effect, and on longer time scales oxygen is controlled more by photosynthesis and respiration. Our simulation demonstrates that interactions among physical and biological dynamics generate complex high-frequency variability in water quality and non-linear reposes to nutrient loading and environmental forcing in shallow water systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The ABCflux database: Arctic–boreal CO 2 flux observations and ancillary information aggregated to monthly time steps across terrestrial ecosystems

Past efforts to synthesize and quantify the magnitude and change in carbon dioxide (CO 2 ) fluxes in terrestrial ecosystems across the rapidly warming Arctic–boreal zone (ABZ) have provided valuable information but were limited in their geographical and temporal coverage. Furthermore, these efforts have been based on data aggregated over varying time periods, often with only minimal site ancillary data, thus limiting their potential to be used in large-scale carbon budget assessments. To bridge these gaps, we developed a standardized monthly database of Arctic–boreal CO 2 fluxes (ABCflux) that aggregates in situ measurements of terrestrial net ecosystem CO 2 exchange and its derived partitioned component fluxes: gross primary productivity and ecosystem respiration. The data span from 1989 to 2020 with over 70 supporting variables that describe key site conditions (e.g., vegetation and disturbance type), micrometeorological and environmental measurements (e.g., air and soil temperatures), and flux measurement techniques. Here, we describe these variables, the spatial and temporal distribution of observations, the main strengths and limitations of the database, and the potential research opportunities it enables. In total, ABCflux includes 244 sites and 6309 monthly observations; 136 sites and 2217 monthly observations represent tundra, and 108 sites and 4092 observations represent the boreal biome. The database includes fluxes estimated with chamber (19 % of the monthly observations), snow diffusion (3 %) and eddy covariance (78 %) techniques. The largest number of observations were collected during the climatological summer (June–August; 32 %), and fewer observations were available for autumn (September–October; 25 %), winter (December–February; 18 %), and spring (March–May; 25 %). ABCflux can be used in a wide array of empirical, remote sensing and modeling studies to improve understanding of the regional and temporal variability in CO 2 fluxes and to better estimate the terrestrial ABZ CO 2 budget.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Evaluating the Benefits of Bayesian Hierarchical Methods for Analyzing Heterogeneous Environmental Datasets: A Case Study of Marine Organic Carbon Fluxes

Large compilations of heterogeneous environmental observations are increasingly available as public databases, allowing researchers to test hypotheses across datasets. Statistical complexities arise when analyzing compiled data due to unbalanced spatial sampling, variable environmental context, mixed measurement techniques, and other reasons. Hierarchical Bayesian modeling is increasingly used in environmental science to describe these complexities, however few studies explicitly compare the utility of hierarchical Bayesian models to simpler and more commonly applied methods. Here we demonstrate the utility of the hierarchical Bayesian approach with application to a large compiled environmental dataset consisting of 5,741 marine vertical organic carbon flux observations from 407 sampling locations spanning eight biomes across the global ocean. We fit a global scale Bayesian hierarchical model that describes the vertical profile of organic carbon flux with depth. Profile parameters within a particular biome are assumed to share a common deviation from the global mean profile. Individual station-level parameters are then modeled as deviations from the common biome-level profile. The hierarchical approach is shown to have several benefits over simpler and more common data aggregation methods. First, the hierarchical approach avoids statistical complexities introduced due to unbalanced sampling and allows for flexible incorporation of spatial heterogeneitites in model parameters. Second, the hierarchical approach uses the whole dataset simultaneously to fit the model parameters which shares information across datasets and reduces the uncertainty up to 95% in individual profiles. Third, the Bayesian approach incorporates prior scientific information about model parameters; for example, the non-negativity of chemical concentrations or mass-balance, which we apply here. We explicitly quantify each of these properties in turn. We emphasize the generality of the hierarchical Bayesian approach for diverse environmental applications and its increasing feasibility for large datasets due to recent developments in Markov Chain Monte Carlo algorithms and easy-to-use high-level software implementations.

54 ENVIRONMENTAL SCIENCES↗