Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Global Sensitivity Analysis”

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 775 records · Page 43

Data for: Patel et al. Carbon flux estimates are sensitive to data source: A comparison of field and lab temperature sensitivity data

This dataset contains data and code used for the paper "Carbon flux estimates are sensitive to data source: A comparison of field and lab temperature sensitivity data" [DOI COMING SOON]A large literature exists on mechanisms driving soil production of the greenhouse gases CO2 and CH4. Measurements of these gases’ fluxes are often performed using closed-chamber incubations in the laboratory or in situ, i.e., in the field. Although it is common knowledge that measurements obtained through field studies vs. laboratory incubations can diverge because of the vastly different conditions of these environments, few studies have systematically examined these patterns. It is crucial to understand the magnitude and reasons for any differences, as these data are used to parametrize and benchmark ecosystem- to global-scale models, which are then susceptible to the biases of the source data. Here, we specifically examine how greenhouse gas measurements may be influenced by whether the measurement/incubation was conducted in the field vs. laboratory, focusing on CO2 and CH4 measurements. We use Q10 of greenhouse gas flux (temperature sensitivity) for our analyses, because of the ubiquity of this metric in biological and Earth system sciences and its importance to many modeling frameworks. We predicted that laboratory measurements would be less variable, but also less representative of true field conditions. However, there was greater variability in the Q10 values calculated from lab-based measurements of CO2 fluxes, because lab experiments explore extremes rarely seen in situ, and reflect the physical and chemical disturbances occurring during sampling, transport, and incubation. Overall, respiration Q10 values were significantly greater in laboratory incubations (mean = 4.19) than field measurements (mean = 3.05), with strong influences of incubation temperature and climate region/biome. However, this was in part because field measurements typically represent total respiration (Rs), whereas lab incubations typically represent heterotrophic respiration (Rh), making direct comparisons difficult to interpret. Focusing only on Rh-derived Q10, these values showed almost identical distributions across laboratory (n = 1110) and field (n = 581) experiments, providing strong support for using the former as an experimental proxy for the latter, although we caution that geographic biases in the extant data make this conclusion tentative. Due to a smaller sample size of CH4 Q10 data, we were unable to perform a comparable robust analysis, but we expect similar interactions with soil temperature, moisture, and environmental/climatic variables. Our results here suggest the need for more concerted efforts to document and standardize these data, including sample and site metadata. This dataset contains a compressed (.zip) archive of the data and R scripts used for this manuscript. The dataset includes files in .csv format, which can be accessed and processed using MS Excel or R. This archive can also be accessed on GitHub at https://github.com/kaizadp/field_lab_q10 (DOI: 10.5281/zenodo.7106554).

54 ENVIRONMENTAL SCIENCES↗

New Mastcam Multispectral Rock Classes in Sulfate-Bearing Strata, Gale Crater, Mars

Orbital data of Gale Crater, Mars has identified a transition in the stratigraphy of Mt. Sharp indicating an environmental change from a wetter one that accommodated clay mineral formation to a drier environment that led to the precipitation of sulfates. Mt. Sharp is not the only location on Mars for which this transition has been observed which suggests a global environmental change around the Noachian-Hesperian transition. Also visible from orbit is a dark-toned marker band that lies within the sulfate-bearing strata hypothesized to be a volcanic ash deposit. The Mars Science Laboratory (MSL) Curiosity rover is currently exploring this region of Mt. Sharp. In-situ analyses by the rover show changes in morphology and chemistry between the sulfate-bearing layers below and above the marker band and the marker band itself. The MSL instrument payload includes the Mast Camera (Mastcam), which collects multispectral images to provide context for other instruments and broad mineralogic interpretations. In addition to mafic Fe-bearing minerals, Mastcam is sensitive to oxides and some hydrated minerals, important indicators of alteration. In our previous work, we compiled a database of Mastcam spectra through Curiosity’s exploration of Vera Rubin ridge (sols 0-2302), from which 9 spectral classes of rocks were identified (Fig. 1D). Each class has diagnostic spectral features that reflect a common mineralogic interpretation. More recently, as Curiosity traversed through the clay-rich region, Glen Torridon (GT), there was some variation in the rock spectral classes potentially due to nontronite and other phyllosilicates, but not enough to warrant new spectral classes (Fig. 1B). We hypothesize that new classes will appear in accordance with sulfate-bearing strata and the marker band, specifically appearing as a ‘downturn’ in the last two wavelength filters that is inherent, but not unique, to polyhydrated Mg-sulfates. Here, we expand upon our analysis of Mastcam multispectral observations through Curiosity’s second encounter with the marker band (up to sol 3672), comparing the spectral diversity within the sulfate-bearing unit and the marker band to the spectral classes encountered previously in the traverse. These analyses will provide a basis for comparison and aid in selecting targets for in-situ investigations as Curiosity continues its ascent of Mt. Sharp.

A M Eng↗

ICESat Observations of Arctic Sea Ice: A First Look

Analysis of near-coincident ICESat and RADARSAT imagery shows that the retrieved elevations from the laser altimeter are sensitive to new openings (containing thin ice or open water) in the sea ice cover as well as to surface relief of old and first-year ice. The precision of the elevation estimates, measured over relatively flat sea ice, is approx. 2 cm. Using the thickness of thin-ice in recent openings to estimate sea level references, we obtain the sea-ice freeboard along the altimeter tracks. This step is necessitated by the large uncertainties in the sea surface topography compared to that required for accurate determination of freeboard. Unknown snow depth introduces the largest uncertainty in the conversion of freeboard to ice thickness. Surface roughness is also derived, for the first time, from the variability of successive elevation estimates along the altimeter track. Overall, these ICESat measurements provide an unprecedented view of the Arctic Ocean ice cover at length scales at and above the spatial dimension of the altimeter footprint of approx. 70 m.

Arctic Ocean↗

An Updated Life Cycle Assessment of Utility-Scale Solar Photovoltaic Systems Installed in the United States

Given the high deployment targets for solar photovoltaics (PV) needed to meet U.S. decarbonization goals, and the limited carbon budget remaining to limit global temperature rise, accurate accounting of the energy-use and greenhouse-gas emissions over the life-cycle of PV systems is needed. In the United States, most PV systems are large utility-scale systems which use single-axis trackers and central inverters, which are not commonly examined in existing life-cycle assessment (LCA) literature. In this study, we present a cradle-to-grave LCA of a typical silicon U.S. utility PV (UPV) installation which is consistent with the utility system features documented in the annual NREL PV system cost benchmark reports. We analyze and present results for four main metrics: cumulative energy demand (CED), greenhouse gas (GHG) emissions, energy payback time (EPBT), and carbon payback time (CPBT). We consider six primary manufacturing options: three based on an imported PV module supply chain (comparing low-carbon imports, high-carbon imports, and average imports), and three based on a potential domestic PV module supply chain (comparing low-carbon U.S. regions, high-carbon U.S. regions, and average U.S. regions). These manufacturing options were then paired with installation locations to create six main cases: low-carbon options were installed in Phoenix (Arizona), high-carbon options were installed in Seattle (Washington), and average options were installed in Fredonia (Kansas). These locations were selected to represent a range of irradiance and grid mixes in the United States, in order to illustrate the likely range of EPBTs and CPBTs possible across the United States. For all six cases, a sensitivity analysis for end-of-life (EOL) handling was explored to capture current and future management options: landfilling, partial recycling, and high-quality recycling. For the purposes of this report, the benchmark system was defined to use an average imported supply chain with partial recycling, installed in Fredonia (Kansas). CED results show ratios at or below 0.1 MJ oil-eq /MJgenerated which demonstrates efficient use of primary energy resources (below a 1:1 ratio), and represents a slight improvement over previous results in literature. GHG emissions per kWh range from 10-36 g CO 2 e, which are consistent with or lower than previous results published by NREL and IEA-PVPS. We use a graphical approach for calculating EPBT and CPBT in this report, which improves upon methods typically used in literature by accounting for non-linearity and avoiding data quality issues associated with long-term projections. EPBT was determined to vary from 0.5 to 1.2 years, with a benchmark EPBT of 0.6 years; CPBT was shown to vary from 0.8 to 20 years, with benchmark CPBT of 2.1 years, which is lower than other estimates from recent literature (typically >2 years).

14 SOLAR ENERGY↗

Radiative Susceptibility of Cloudy Atmospheres to Droplet Number Perturbations: 1. Theoretical Analysis and Examples from MODIS

Theoretical and satellite-based assessments of the sensitivity of broadband shortwave radiative fluxes in cloudy atmospheres to small perturbations in the cloud droplet number concentration (N) of liquid water clouds under constant water conditions are performed. Two approaches to study this sensitivity are adopted: absolute increases in N, for which the radiative response is referred to as absolute cloud susceptibility, and relative increases in N or relative cloud susceptibility. Estimating the former is more challenging as it requires an assumed value for either cloud liquid water content or geometrical thickness; both susceptibilities require an assumed relationship between the droplet volume and effective radius. Expanding upon previous susceptibility studies, present radiative calculations include the effect of AN perturbations on droplet asymmetry parameter and single-scattering albedo, in addition to extinction. Absolute cloud susceptibility has a strong nonlinear dependence on the droplet effective radius as expected, while relative cloud susceptibility is primarily dependent on optical thickness. Molecular absorption and reflecting surfaces both reduce the relative contribution of the cloud to the top-of-atmosphere (TOA) flux and therefore also reduce the TOA albedo susceptibility. Transmittance susceptibilities are negative with absolute values similar to albedo susceptibility, while atmospheric absorptance susceptibilities are about an order of magnitude smaller than albedo susceptibilities and can be either positive or negative. Observation-based susceptibility calculations are derived from MODIS pixel-level retrievals of liquid water cloud optical thickness, effective radius, and cloud top temperature; two data granule examples are shown. Susceptibility quantifies the aerosol indirect effect sensitivity in a way that can be easily computed from model fields. As such, susceptibilities derived from MODIS observations provide a higher-order test of model cloud properties used for indirect effect studies. MODIS-derived global distributions of cloud susceptibility and radiative forcing calculations are presented in a companion paper.

Platnick, Steven↗

FY22 Status Report on the ART-GCR CMVB and CNWG International Collaborations

This work presents the numerical model for the high temperature test reactor (HTTR) loss of forced cooling (LOFC) experiment with the INL codes Griffin, BISON, and RELAP-7, based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Promising results were obtained, with the overall behavior of the reactor successfully captured. Changes in the heat transfer coupling, as compared to the fiscal year (FY)-21 model, enabled drastic improvement of the steady-state solution, both in terms of computational time and global energy discrepancies. The former was reduced by a factor of roughly 60, whereas the latter decreased from 11% to less than 2%. Furthermore, the discrepancy in the steady-state multiplication factor was improved from +2,300 and +2,900 pcm (for the 30 and 9 MW cases) to -700 and +1,200 pcm, respectively, and now falls well within the large measurement uncertainties stemming from graphite impurities. Validation of the Monte Carlo model used to generate cross sections was also performed against available measurements. Though significant, the discrepancies remain acceptable overall in light of the large uncertainty. Numerous improvements are still needed to better compare with the experiment involving the 9 MW case and to instill greater confidence in the model’s ability to accurately predict the 30 MW behavior. Specifically, the power levels predicted by the 9 MW transient simulation following re-criticality remain low, pointing to an underestimation of the passive cooling of the core. A key aspect of future work will be to better understand the flow pattern during the LOFC event, particularly to determine if natural or forced convection is occurring inside the reactor pressure vessel (RPV). More generally, additional validation data would be immensely useful for further enhancing the numerical model and better matching the experiments. In addition, a more sophisticated thermal-hydraulics model that simulates all the channels as a single system model should be considered to take into account the rest of the primary loop. Finally, even if the results are in better agreement with the experiments, sensitivity analysis and uncertainty quantification will be necessary to evaluate the model uncertainty.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Analysis of Present Day and Future OH and Methane Lifetime in the ACCMIP Simulations

Results from simulations performed for the Atmospheric Chemistry and Climate Modeling Intercomparison Project (ACCMIP) are analysed to examine how OH and methane lifetime may change from present day to the future, under different climate and emissions scenarios. Present day (2000) mean tropospheric chemical lifetime derived from the ACCMIP multi-model mean is 9.8+/-1.6 yr (9.3+/-0.9 yr when only including selected models), lower than a recent observationally-based estimate, but with a similar range to previous multi-model estimates. Future model projections are based on the four Representative Concentration Pathways (RCPs), and the results also exhibit a large range. Decreases in global methane lifetime of 4.5 +/- 9.1% are simulated for the scenario with lowest radiative forcing by 2100 (RCP 2.6), while increases of 8.5+/-10.4% are simulated for the scenario with highest radiative forcing (RCP 8.5). In this scenario, the key driver of the evolution of OH and methane lifetime is methane itself, since its concentration more than doubles by 2100 and it consumes much of the OH that exists in the troposphere. Stratospheric ozone recovery, which drives tropospheric OH decreases through photolysis modifications, also plays a partial role. In the other scenarios, where methane changes are less drastic, the interplay between various competing drivers leads to smaller and more diverse OH and methane lifetime responses, which are difficult to attribute. For all scenarios, regional OH changes are even more variable, with the most robust feature being the large decreases over the remote oceans in RCP8.5. Through a regression analysis, we suggest that differences in emissions of non-methane volatile organic compounds and in the simulation of photolysis rates may be the main factors causing the differences in simulated present day OH and methane lifetime. Diversity in predicted changes between present day and future OH was found to be associated more strongly with differences in modelled temperature and stratospheric ozone changes. Finally, through perturbation experiments we calculated an OH feedback factor (F) of 1.24 from present day conditions (1.50 from 2100 RCP8.5 conditions) and a climate feedback on methane lifetime of 0.33+-0.13 yr/K, on average. Models that did not include interactive stratospheric ozone effects on photolysis showed a stronger sensitivity to climate, as they did not account for negative effects of climate-driven stratospheric ozone recovery on tropospheric OH, which would have partly offset the overall OH/methane lifetime response to climate change.

atmospheric composition↗

Carbon flux estimates are sensitive to data source: a comparison of field and lab temperature sensitivity data

Abstract A large literature exists on mechanisms driving soil production of the greenhouse gases CO 2 and CH 4 . Although it is common knowledge that measurements obtained through field studies vs. laboratory incubations can diverge because of the vastly different conditions of these environments, few studies have systematically examined these patterns. These data are used to parameterize and benchmark ecosystem- to global-scale models, which are then susceptible to the biases of the source data. Here, we examine how greenhouse gas measurements may be influenced by whether the measurement/incubation was conducted in the field vs. laboratory, focusing on CO 2 and CH 4 measurements. We use Q 10 of greenhouse gas flux (temperature sensitivity) for our analyses because this metric is commonly used in biological and Earth system sciences and is an important parameter in many modeling frameworks. We predicted that laboratory measurements would be less variable, but also less representative of true field conditions. However, there was greater variability in the Q 10 values calculated from lab-based measurements of CO 2 fluxes, because lab experiments explore extremes rarely seen in situ , and reflect the physical and chemical disturbances occurring during sampling, transport, and incubation. Overall, respiration Q 10 values were significantly greater in laboratory incubations (mean = 4.19) than field measurements (mean = 3.05), with strong influences of incubation temperature and climate region/biome. However, this was in part because field measurements typically represent total respiration (Rs), whereas lab incubations typically represent heterotrophic respiration (Rh), making direct comparisons difficult to interpret. Focusing only on Rh-derived Q 10 , these values showed almost identical distributions across laboratory ( n = 1110) and field ( n = 581) experiments, providing strong support for using the former as an experimental proxy for the latter, although we caution that geographic biases in the extant data make this conclusion tentative. Due to a smaller sample size of CH 4 Q 10 data, we were unable to perform a comparable robust analysis, but we expect similar interactions with soil temperature, moisture, and environmental/climatic variables. Our results here suggest the need for more concerted efforts to document and standardize these data, including sample and site metadata.

54 ENVIRONMENTAL SCIENCES↗

Aquarius Instrument Science Calibration During the Risk Reduction Phase

This final report presents the results of work performed under NASA Grant NAG512726 during the period 15 January 2003 through 30 June 2004. An analysis was performed of a possible vicarious calibration method for use by Aquarius to monitor and stabilize the absolute and relative calibration of its microwave radiometer. Stationary statistical properties of the brightness temperature (T(sub B)) measured by a low Earth orbiting radiometer operating at 1.4135 GHz are considered as a means of validating its absolute calibration. The global minimum, maximum, and average T(sub B) are considered, together with a vicarious cold reference method that detects the presence of a sharp lower bound on naturally occurring values for T(sub B). Of particular interest is the reliability with which these statistics can be extracted from a realistic distribution of T(sub B) measurements that would be observed by a typical sensor. Simulations of measurements are performed that include the effects of instrument noise and variable environmental factors such as the global water vapor and ocean surface temperature, salinity and wind distributions. Global minima can vary widely due to instrument noise and are not a reliable calibration reference. Global maxima are strongly influenced by several environmental factors as well as instrument noise and are even less stationary. Global averages are largely insensitive to instrument noise and, in most cases, to environmental conditions as well. The global average T(sub B) varies at only the 0.1 K RMS level except in cases of anomalously high winds, when it can increase considerably more. The vicarious cold reference is similarly insensitive to instrument effects and most environmental factors. It is not significantly affected by high wind conditions. The stability of the vicarious reference is, however, found to be somewhat sensitive (at the several tenths of Kelvins level) to variations in the background cold space brightness, T(sub c). The global average is much less sensitive to this parameter and so using two approaches together can be mutually beneficial.

Ruf, Christopher S.↗

Temporal variabilities of soil carbon dioxide fluxes from cornfield impacted by temperature and precipitation changes through high-frequent measurement and DAYCENT modelling

Soil carbon dioxide (CO 2 ) emissions from the field of corn (Zea mays L.) play an important role in global warming. This study investigated temporal variability of soil CO 2 fluxes (R s ) with soil temperature (T s ) and moisture (θ) and built DAYCENT models for predicting future impacts of climate changes on R s using the measured high-frequency data. R s trend was tested by Mann–Kendall and Sen Estimator. Predicted R s s under different climate scenarios were compared using Parallel-line Analysis. Herein, the findings indicated that daily R s exponentially increased with T s constrained by θ. During the θ of 27–31%, there was a strong exponential relationship between R s and T s , but the relationship was weaker for the θ of 38–41% and 22–26%. Soil environmental index (SEI, T s × θ) significantly impacted R s with linear regression R s 0.5 = 0.4599 + 0.002059 × SEI in 2008, 2009 and 2011. At the diurnal scale, there were different trends in R s s and relationships among R s and T s and θ in different years. Predicted yearly R s s, root R s s and corn yield in 2014–2049 increased with an increase in temperature scenarios, but the R s s significantly increased as temperature rose by 1°C or higher. Predicted yearly R s s, root R s s and yield reduced with precipitation scenario increase, and the root R s s and yield significantly diminished as precipitation increased by 15 and 30%. Predicted yearly R s from cornfields had a significantly increasing trend. Future research is needed to explore methods for mitigating cornfield R s and evaluating sensitivities of different cropland R s s to temperature changes.

54 ENVIRONMENTAL SCIENCES↗

Data-Driven Discovery of Linear Molecular Probes with Optimal Selective Affinity for PFAS in Water

Approaches to tackle the wide and growing variety of highly persistent per- and polyfluoroalkyl substances (PFAS) are of pressing global need because of their detrimental human health effects, such as cancer, birth defects, and hormone imbalance. Sensitive, selective, and easy-to-use real-time sensors to monitor and detect PFAS and sorbents to extract them are critical to meeting government-mandated environmental concentrations. In this work, we combine all-atom molecular dynamics simulations, enhanced sampling, deep representational learning, and Bayesian optimization to perform high-throughput virtual screening for highly sensitive and selective molecular probes. Our molecular design space consists of 3850 linear hydrocarbon chains with varying degrees of halogenation with and without amine- and phosphine-based headgroups. By employing a data-driven search process, we efficiently explore the molecular design space to optimize the sensitivity to perfluorooctanesulfonic acid (PFOS) as a prototypical PFAS analyte and selectivity relative to a sodium dodecyl sulfate (SDS) interferent. We calculate 504 Gibbs free energies of probe-analyte and probe-interferent interactions and identify probes with PFOS association free energies of up to (-ΔG PFOS ) = 9.8 ± 0.2 kJ/mol and selectivities relative to SDS of (-ΔΔG PFOS–SDS ) = 3.1 ± 1.5 kJ/mol. A C 11 Br 23 P(CH 3 ) 2 probe containing 11 backbone brominated carbons and a tertiary phosphine headgroup possesses the most sensitive binding constant to PFOS within the defined search space of K b PFOS = 177.4 ± 12.7, and a semibrominated probe C 5 H 11 C 7 Br 14 N(CH 3 ) 2 containing 12 backbone carbons and a tertiary amine headgroup possesses the highest selectivity relative to SDS of K b PFOS /K b SDS = 4.6 ± 1.7. A retrospective analysis of our data to extract interpretable design rules reveals that the sensitivity of linear hydrogenated probes increases by approximately 1 kJ/mol per C–C bond. The addition or removal of halogen atoms and amine or phosphine headgroups produces nonmonotonic changes in both sensitivity and selectivity with changes to the sensitivity of up to 2.5 kJ/mol. Finally, this work places empirical limitations on the performance of a wide range of linear probes for PFOS detection and offers a generic strategy for high-throughput computational screening to promote selective and sensitive binding.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimental validation of an organic rankine-vapor compression cooling cycle using low GWP refrigerant R1234ze(E)

There is a significant global opportunity to capture and utilize low grade waste heat to reduce fossil fuel consumption, greenhouse gas emissions, and improve energy efficiency across a wide range of industries. In this work, an advanced type of thermally activated cooling system, an organic Rankine-vapor compression cycle (ORVC) with novel heat integration strategies, was designed and tested at a relevant scale for industrial waste heat recovery (300 kW th cooling capacity). The ORVC linked an organic Rankine power cycle and a vapor compression cooling cycle using a turbine and compressor that shared a single shaft. The ORVC test facility absorbed waste heat from a liquid stream at 91 °C to simulate engine coolant in diesel generator sets, rejected heat to a glycol stream at 30 °C, and generated chilled water at 7 °C. A cooling capacity of 264 kW ± 3.5 kW was experimentally validated with a COP of 0.56 ± 0.01 during steady-state operation at the design temperatures. The thermal efficiency, accounting for pump work, of the Rankine cycle was 7.7% ± 0.22% and the COP of the vapor compression cycle was 5.25 ± 0.09. The centrifugal turbo-compressor operated at 31.5 kRPM ± 0.3 kRPM, with a turbine and compressor isentropic efficiencies of 76.7% ± 0.90% and 84.8% ± 0.54%, respectively, with near-perfect power transmission between these components. The pressure drop in the piping and heat exchangers were significantly larger than expected which had a detrimental impact on the performance of the ORVC. In addition, the condenser on the cooling cycle could not deliver the subcooling as specified from the design point modeling. Furthermore, the results from the sensitivity analysis showed that the higher condenser glycol outlet temperature had the largest impact on performance, which is consistent with other analytical models in the literature. When the ORVC simulations were updated with experimental values for isentropic efficiencies of the turbomachinery, the thermal COP was 0.66 which represents an estimate of the predicted performance if test facility limitations are overcome.

30 DIRECT ENERGY CONVERSION↗

Representing methane emissions from wet tropical forest soils using microbial functional groups constrained by soil diffusivity

Tropical ecosystems contribute significantly to global emissions of methane (CH 4 ), and landscape topography influences the rate of CH 4 emissions from wet tropical forest soils. However, extreme events such as drought can alter normal topographic patterns of emissions. Here we explain the dynamics of CH 4 emissions during normal and drought conditions across a catena in the Luquillo Experimental Forest, Puerto Rico. Valley soils served as the major source of CH 4 emissions in a normal precipitation year (2016), but drought recovery in 2015 resulted in dramatic pulses in CH 4 emissions from all topographic positions. Geochemical parameters including (i) dissolved organic carbon (C), acetate, and soil pH and (ii) hydrological parameters like soil moisture and oxygen (O 2 ) concentrations varied across the catena. During the drought, soil moisture decreased in the slope and ridge, and O 2 concentrations increased in the valley. We simulated the dynamics of CH 4 emissions with the Microbial Model for Methane Dynamics-Dual Arrhenius and Michaelis–Menten (M3D-DAMM), which couples a microbial functional group CH 4 model with a diffusivity module for solute and gas transport within soil microsites. Contrasting patterns of soil moisture, O 2 , acetate, and associated changes in soil pH with topography regulated simulated CH 4 emissions, but emissions were also altered by rate-limited diffusion in soil microsites. Changes in simulated available substrate for CH 4 production (acetate, CO 2 , and H 2 ) and oxidation (O 2 and CH 4 ) increased the predicted biomass of methanotrophs during the drought event and methanogens during drought recovery, which in turn affected net emissions of CH 4 . A variance-based sensitivity analysis suggested that parameters related to aceticlastic methanogenesis and methanotrophy were most critical to simulate net CH 4 emissions. This study enhanced the predictive capability for CH 4 emissions associated with complex topography and drought in wet tropical forest soils.

54 ENVIRONMENTAL SCIENCES↗

The Impact of Assimilating Large Volumes of GNSS Radio Occultation Observations from Spire’s Commercial Constellation into NASA’s Global Earth Observing System

The upcoming retrospective analysis for the 21st century (R21C) reanalysis product from NASA will include the full dataset of GNSS Radio Occultation (RO) observations collected by Spire with their constellation of smallsats. The Spire RO dataset, which was purchased by NASA for its Commercial Smallsat Data Acquisition (CSDA) archive, has global coverage and includes approximately 6 thousand RO profiles per day during 2019 increasing to nearly 20 thousand in 2022-2023. This more than doubles the volume of RO profiles from all other routinely assimilated RO missions combined during 2023. The increase in the number of RO profiles available is particularly important in the extratropics as the next largest RO constellation, the state-of-the-science Formosa Satellite Mission 7 (FORMOSAT-7)/Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC-2) mission (hereafter COSMIC-2), is focused on the Tropics and only provides approximately 4 thousand RO profiles per day. The large number of observations in the Spire RO dataset have great potential to improve analyses and forecasts of the Earth’s atmosphere produced by numerical weather prediction systems. Their impact is assessed through numerical experiments using NASA’s Global Earth Observing System (GEOS) Atmospheric Data Assimilation System (ADAS) with and without the observations from the COSMIC-2 mission. The ability of the Spire RO observations to make up for the omission of the RO observations from COSMIC-2 is directly assessed over the tropics and the impact of Spire over the extratropics is compared to that of the available RO missions that sample over this region. Finally, the quality of the Spire RO observations is compared to that of the routinely assimilated RO missions using Forecast Sensitivity-based Observation Impact (FSOI).

Michael J. Murphy, Jr.↗

The Impact of Assimilating Large Volumes of GNSS Radio Occultation Observations from Spire’s Commercial Constellation into NASA’s GEOS

The upcoming retrospective analysis for the 21st century (R21C) reanalysis product from NASA will include the full dataset of GNSS Radio Occultation (RO) observations collected by Spire with their constellation of smallsats. The Spire RO dataset, which was purchased by NASA for its Commercial Smallsat Data Acquisition (CSDA) archive, has global coverage and includes approximately 6 thousand RO profiles per day during 2019 increasing to nearly 20 thousand in 2022-2023. This more than doubles the volume of RO profiles from all other routinely assimilated RO missions combined during 2023. The increase in the number of RO profiles available is particularly important in the extratropics as the next largest RO constellation, the state-of-the-science Formosa Satellite Mission 7 (FORMOSAT-7)/Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC-2) mission (hereafter COSMIC-2), is focused on the Tropics and only provides approximately 4 thousand RO profiles per day. The large number of observations in the Spire RO dataset have great potential to improve analyses and forecasts of the Earth’s atmosphere produced by numerical weather prediction systems. Their impact is assessed through numerical experiments using NASA’s Global Earth Observing System (GEOS) Atmospheric Data Assimilation System (ADAS) with and without the observations from the COSMIC-2 mission. The ability of the Spire RO observations to make up for the omission of the RO observations from COSMIC-2 is directly assessed over the tropics and the impact of Spire over the extratropics is compared to that of the available RO missions that sample over this region. Finally, the quality of the Spire RO observations is compared to that of the routinely assimilated RO missions using Forecast Sensitivity-based Observation Impact (FSOI).

Michael J. Murphy↗

Technical and Economic Evaluation of the First Ever Polymer Flood Field Pilot to Enhance the Recovery of Heavy Oils on Alaska's North Slope via Machine Assisted History Matching

Polymer flooding has become globally established as a potential enhanced oil recovery method for heavy oils. To determine whether this technology may be useful in developing the substantial heavy oil resources on the Alaska North Slope, a polymer flood field pilot commenced at the Milne Point Unit in August 2018. This study seeks to evaluate the results of the field pilot on a technical and economic basis. A reservoir simulation model is constructed and calibrated to predict the oil recovery performance of the pilot through machine-assisted reservoir simulation techniques. To replicate the early water breakthrough observed during waterflooding, transmissibility contrasts are introduced into the simulation model, forcing viscous fingering effects. In the ensuing polymer flood, these transmissibility contrasts are reduced to replicate the restoration of injection conformance during polymer flooding. Transmissibility contrasts are later reinstated to replicate fracture overextension interpreted in one of the producing wells. The calibrated simulation models produced at each stage of the history matching process are used to forecast oil recovery. These forecasts are used as input for economic analysis, incremental to waterflooding expectations. The simulation forecasts indicate that polymer flooding significantly increases the heavy oil production for this field pilot compared to waterflooding alone, yielding attractive project economics. However, meaningful variations between simulation scenarios demonstrate that a simulation model is only valid for prediction if flow behavior in the reservoir remains consistent with that observed during the history matched period. Critically, this means that a simulation model calibrated for waterflooding may not fully capture the technical and economic benefits of an enhanced oil recovery process such as polymer flooding. Subsequently, the simulation model and economic model are used in conjunction to conduct a sensitivity analysis for polymer flood design parameters, from which recommendations are provided for both the continued operation of the current field pilot and future polymer flood designs. The results demonstrate that a higher polymer concentration can be injected due to the development of fractures in the reservoir. The throughput rate should remain high without exceeding operating constraints. A calculated point-forward polymer utilization parameter demonstrates the decreasing efficiency of the polymer flood at later times in the pattern life. Future projects will benefit from starting polymer injection earlier in the pattern life. A pattern with tighter horizontal well spacing will observe a greater incremental benefit from polymer flooding.

Keith, Cody↗

The Impact of Assimilating Large Volumes of GNSS Radio Occultation Observations from Spire’s Commercial Constellation into NASA’s GEOS

The upcoming retrospective analysis for the 21st century (R21C) reanalysis product from NASA will include the full dataset of GNSS Radio Occultation (RO) observations collected by Spire with their constellation of smallsats. The Spire RO dataset, which was purchased by NASA for its Commercial Smallsat Data Acquisition (CSDA) archive, has global coverage and increases in number from approximately 6 thousand RO profiles per day during 2019 to approximately 20 thousand in 2022-2023. This more than doubles the volume of RO profiles from all other routinely assimilated RO missions combined during 2022. The increase in the number of RO profiles available is particularly important in the extratropics as the next largest RO constellation, the state-of-the-science Formosa Satellite Mission 7 (FORMOSAT-7)/Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC-2) mission (hereafter COSMIC-2), is focused on the Tropics and only provides approximately 4 thousand RO profiles per day. The large number of observations in the Spire RO dataset have great potential to improve analyses and forecasts of the Earth’s atmosphere produced by numerical weather prediction systems. Their impact is assessed through numerical experiments using NASA’s Global Earth Observing System (GEOS) Atmospheric Data Assimilation System (ADAS) with and without the observations from the COSMIC-2 mission. The ability of the Spire RO observations to make up for the omission of the RO observations from COSMIC-2 is directly assessed over the tropics and the impact of Spire over the extratropics is compared to that of the available RO missions that sample over this region. Finally, the quality of the Spire RO observations is compared to that of the routinely assimilated RO missions using Forecast Sensitivity-based Observation Impact (FSOI).

Michael J. Murphy↗

Sensitivity of Climate to Changes in NDVI

The sensitivity of global and regional climate to changes in vegetation density is investigated using a coupled biosphere-atmosphere model. The magnitude of the vegetation changes and their spatial distribution are based on natural decadal variability of the normalized difference vegetation index (ndvi). Different scenarios using maximum and minimum vegetation cover were derived from satellite records spanning the period 1982-1990. Albedo decreased in the northern latitudes and increased in the tropics with increased ndvi. The increase in vegetation density revealed that the vegetation's physiological response was constrained by the limits of the available water resources. The difference between the maximum and minimum vegetation scenarios resulted in a 46% increase in absorbed visible solar radiation and a similar increase in gross photosynthetic C02 uptake on a global annual basis. This caused the canopy transpiration and interception fluxes to increase, and reduced those from the soil. The redistribution of the surface energy fluxes substantially reduced the Bowen ratio during the growing season, resulting in cooler and moister near-surface climate, except when soil moisture was limiting. Important effects of increased vegetation on climate are : (1) A cooling of about 1.8 K in the northern latitudes during the growing season and a slight warming during the winter, which is primarily due to the masking of high albedo of snow by a denser canopy. and (2) A year round cooling of 0.8 K in the tropics. These results suggest that increases in vegetation density could partially compensate for parallel increases in greenhouse warming . Increasing vegetation density globally caused both evapotranspiration and precipitation to increase. Evapotranspiration, however increased more than precipitation resulting in a global soil-water deficit of about 15 %. A spectral analysis on the simulated results showed that changes in the state of vegetation could affect the low-frequency modes of the precipitation distribution and might reduce its low frequency variability in the tropics while increasing it in northern latitudes.

Bounoua, L.↗