Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spatial Partitioning”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

High-resolution Y-axis readout for delay-line microchannel anodes

We have devised a method to obtain a high-resolution Y-axis event position determination from microchannel plate delay-line detectors. The method is based on the double-delay-line wedge-wedge charge partition principle of Lampton et al. (1990), where the X axis is read out by a pair of side-by-side delay lines and the Y-axis coordinate is determined by comparing the charges on the two parallel delay lines. However, our new method abandons the common field return path of the two delay lines, and splits them into two adjacent but independent transformer coupled timing circuits having no common ground connection except via the charge-measurement amplifiers. Thus we eliminate the large electrostatic capacitance of the delay lines from the input of the charge-measuring system. Over its limited working field, the spatial resolution in the Y coordinate can easily equal the resolution provided in the X coordinate by the delay line. Applications to extreme and FUV photon-counting spectroscopy are envisioned.

Raffanti, R.↗
NASA NTRS1997

An Adaptive Flow Solver for Air-Borne Vehicles Undergoing Time-Dependent Motions/Deformations

This report describes a concurrent Euler flow solver for flows around complex 3-D bodies. The solver is based on a cell-centered finite volume methodology on 3-D unstructured tetrahedral grids. In this algorithm, spatial discretization for the inviscid convective term is accomplished using an upwind scheme. A localized reconstruction is done for flow variables which is second order accurate. Evolution in time is accomplished using an explicit three-stage Runge-Kutta method which has second order temporal accuracy. This is adapted for concurrent execution using another proven methodology based on concurrent graph abstraction. This solver operates on heterogeneous network architectures. These architectures may include a broad variety of UNIX workstations and PCs running Windows NT, symmetric multiprocessors and distributed-memory multi-computers. The unstructured grid is generated using commercial grid generation tools. The grid is automatically partitioned using a concurrent algorithm based on heat diffusion. This results in memory requirements that are inversely proportional to the number of processors. The solver uses automatic granularity control and resource management techniques both to balance load and communication requirements, and deal with differing memory constraints. These ideas are again based on heat diffusion. Results are subsequently combined for visualization and analysis using commercial CFD tools. Flow simulation results are demonstrated for a constant section wing at subsonic, transonic, and a supersonic case. These results are compared with experimental data and numerical results of other researchers. Performance results are under way for a variety of network topologies.

Singh, Jatinder↗
NASA NTRS1987

Energetic particle acceleration in spherically symmetric accretion flows and shocks

Steady state, spherically symmetric solutions of the cosmic-ray transport equation describing the acceleration of energetic particles in galactic accretion flows onto neutron stars, black holes, white dwarfs, and protostars are studied. The results indicate that astrophysical accretion flows can be partitioned into distinct classes depending upon whether the accretion rate lies above or below a critical value of a few times 10 to the -7th stellar masses/yr. When the particle transport is convection-dominated, both classes of accretion flows exhibit a spectral index appropriate for first-order Fermi acceleration at a plane shock in the absence of losses. As the particle transport becomes diffusion-dominated, both classes show a break and precipitous falloff in the particle spectrum due to the escape of these particles from the accretion flow. The precise nature of the spectrum depends on the relationship between the particle momentum and the spatial diffusion coefficient.

Webb, G. M.↗
NASA NTRS2022

Historically inconsistent productivity and respiration fluxes in the global terrestrial carbon cycle

The terrestrial carbon cycle is a major source of uncertainty in climate projections. Its dominant fluxes, gross primary productivity (GPP), and respiration (in particular soil respiration, RS), are typically estimated from independent satellite-driven models and upscaled in situ measurements, respectively. We combine carbon-cycle flux estimates and partitioning coefficients to show that historical estimates of global GPP and RS are irreconcilable. When we estimate GPP based on RS measurements and some assumptions about RS:GPP ratios, we found the resulted global GPP values (bootstrap mean 149^(+29)_(−23) Pg C/yr) are significantly higher than most GPP estimates reported in the literature (113^(+18)_(−18) Pg C/yr). Similarly, historical GPP estimates imply a soil respiration flux (Rs(_GPP), bootstrap mean of 68^(+10)_(−8) Pg C/yr) statistically inconsistent with most published RS values (87^(+9)_(−8) Pg C/yr), although recent, higher, GPP estimates are narrowing this gap. Furthermore, global R_(S):GPP ratios are inconsistent with spatial averages of this ratio calculated from individual sites as well as CMIP6 model results. This discrepancy has implications for our understanding of carbon turnover times and the terrestrial sensitivity to climate change. Future efforts should reconcile the discrepancies associated with calculations for GPP and Rs to improve estimates of the global carbon budget.

Jinshi Jian↗
NASA NTRS1989

Estimation of regional surface resistance to evapotranspiration from NDVI and thermal-IR AVHRR data

Infrared surface temperatures from satellite sensors have been used to infer evaporation and soil moisture distribution over large areas. However, surface energy partitioning to latent versus sensible heat changes with surface vegetation cover and water availability. The hypothesis that the relationship between surface temperature and canopy density is sensitivite to seasonal changes in canopy resistance of conifer forests is presently tested. Surface temperature and canopy density were computed for a 20 x 25 km forested region in Montana, from the NOAA/AVHRR for 8 days during the summer of 1985. A forest ecosystem model, FOREST-BGC, simulated canopy resistance for the same period. For all eight days, surface temperatures had high association with canopy density, measured as Normalized Difference Vegetation Index, implying that latent heat exchange is the major cause of spatial variations in surface radiant tmeperatures.

Nemani, Ramakrishna R.↗
NASA NTRS2019

Applicability of Spectral Decomposition by Varimax-Rotated, Principal Component Analysis to the Surface Biology and Geology (SBG) VNIR Mission Concept

Cyanobacterial and Harmful Algal Blooms (CyanoHABs) are a growing concern in coastal and inland waters. But, spectral interference from multiple constituents in optically complex waters can hamper application of remote sensing using traditional image processing methods. The Kent State University (KSU) spectral decomposition method can be applied to multispectral and hyperspectral remote sensing images (e.g. HICO and the NASA Glenn HSI2) to partition and identify signals related to cyanobacteria, algae, pigment degradation products and suspended sediment in each pixel. Fundamental to the use of remote sensing data is the ability to extract independent signals from correlated hyperspectral VNIR data cubes. The Kent State University varimax-rotated, principal component analysis method (VPCA) is important to integrate into the SBG VNIR mission concept because it provides greater specificity, a software-based SNR boost relative to hardware performance, and can assist with Cal/Val, Modeling and Applications. We present examples of the hyperspectral application of the KSU VPCA method with relevance to SBG. The information extracted by VPCA can be validated spectrally or spatially with laboratory and/or in situ sensors, which capture spatial or time series of information at discrete points within remote sensing images. Comparisons show hyperspectral sensors extract more components than multispectral ones, but more independent information can be extracted from multispectral sensors by VPCA than traditional band ratio approaches. The spectral decomposition method is capable of enhancing the signal to noise ratio (SNR) of the NASA Glenn, second-generation hyperspectral imager by a factor of 7x to 20x, with a spectral reproducibility of ±3%. The spectral decomposition method, when compared against existing remote sensing monitoring methods exhibits both greater specificity and a lower detection limit. The method has been validated with multispectral images in Lake Erie to quantify the Microcystis CyanoHAB and from the Indian River Lagoon, Florida to quantify the Brown Tide resulting from A. lagunesnsis. Field operations in the Western Basin of Lake Erie were conducted using a bbe Fluoroprobe to collect vertical profiles and horizontal tows along a transect from the Toledo to the Detroit Lighthouse during coincident satellite overpasses. Extraction of pixel values from the MODIS Aqua sensor yields agreement between in situ field and lab-based measures of cyanobacterial, cryptophyte, diatoms and green algae, suspended sediment and pigment degradation products with R2>0.8.

hyperspectral↗
DOE OSTI2024

Partition between supercooled liquid droplets and ice crystals in mixed-phase clouds based on airborne in situ observations

Abstract. The onset of ice nucleation in mixed-phase clouds determines the lifetime and microphysical properties of ice clouds. In this work, we develop a novel method that differentiates between various phases of mixed-phase clouds, such as clouds dominated by pure liquid or pure ice segments, compared with those having ice crystals surrounded by supercooled liquid water droplets or vice versa. Using this method, we examine the relationship between the macrophysical and microphysical properties of Southern Ocean mixed-phase clouds at −40 to 0 °C (e.g. stratiform and cumuliform clouds) based on the in situ aircraft-based observations during the US National Science Foundation Southern Ocean Clouds, Radiation, Aerosol Transport Experimental Study (SOCRATES) flight campaign. The results show that the exchange between supercooled liquid water and ice crystals from a macrophysical perspective, represented by the increasing spatial ratio of regions containing ice crystals relative to the total in-cloud region (defined as ice spatial ratio), is positively correlated with the phase exchange from a microphysical perspective, represented by the increasing ice water content (IWC), decreasing liquid water content (LWC), increasing ice mass fraction, and increasing ice particle number fraction (IPNF). The mass exchange between liquid and ice becomes more significant during phase 3 when pure ice cloud regions (ICRs) start to appear. Occurrence frequencies of cloud thermodynamic phases show a significant phase change from liquid to ice at a similar temperature (i.e. −17.5 °C) among three types of definitions of mixed-phase clouds based on ice spatial ratio, ice mass fraction, or IPNF. Aerosol indirect effects are quantified for different phases using number concentrations of aerosols greater than 100 or 500 nm (N>100 and N>500, respectively). N>500 shows stronger positive correlations with ice spatial ratios compared with N>100. This result indicates that larger aerosols potentially contain ice-nucleating particles (INPs), which facilitate the formation of ice crystals in mixed-phase clouds. The impact of N>500 is also more significant in phase 2 when ice crystals just start to appear in the mixed phase compared with phase 3 when pure ICRs have formed, possibly due to the competing aerosol indirect effects on primary and secondary ice production in phase 3. The thermodynamic and dynamic conditions are quantified for each phase. The results show stronger in-cloud turbulence and higher updraughts in phases 2 and 3 when liquid and ice coexist compared with pure liquid or ice (phases 1 and 4, respectively). The highest updraughts and turbulence are seen in phase 3 when supercooled liquid droplets are surrounded by ice crystals. These results indicate both updraughts and turbulence support the maintenance of supercooled liquid water amongst ice crystals. Overall, these results illustrate the varying effects of aerosols, thermodynamics, and dynamics through various stages of mixed-phase cloud evolution based on this new method that categorizes cloud phases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The Colorado East River Community Observatory Data Collection

Abstract The U.S. Department of Energy's (DOE) Colorado East River Community Observatory (ER) in the Upper Colorado River Basin was established in 2015 as a representative mountainous, snow‐dominated watershed to study hydrobiogeochemical responses to hydrological perturbations in headwater systems. The ER is characterized by steep elevation, geologic, hydrologic and vegetation gradients along floodplain, montane, subalpine, and alpine life zones, which makes it an ideal location for researchers to understand how different mountain subsystems contribute to overall watershed behaviour. The ER has both long‐term and spatially‐extensive observations and experimental campaigns carried out by the Watershed Function Scientific Focus Area (SFA), led by Lawrence Berkeley National Laboratory, and researchers from over 30 organizations who conduct cross‐disciplinary process‐based investigations and modelling of watershed behaviour. The heterogeneous data generated at the ER include hydrological, genomic, biogeochemical, climate, vegetation, geological, and remote sensing data, which combined with model inputs and outputs comprise a collection of datasets and value‐added products within a mountainous watershed that span multiple spatiotemporal scales, compartments, and life zones. Within 5 years of collection, these datasets have revealed insights into numerous aspects of watershed function such as factors influencing snow accumulation and melt timing, water balance partitioning, and impacts of floodplain biogeochemistry and hillslope ecohydrology on riverine geochemical exports. Data generated by the SFA are managed and curated through its Data Management Framework. The SFA has an open data policy, and over 70 ER datasets are publicly available through relevant data repositories. A public interactive map of data collection sites run by the SFA is available to inform the broader community about SFA field activities. Here, we describe the ER and the SFA measurement network, present the public data collection generated by the SFA and partner institutions, and highlight the value of collecting multidisciplinary multiscale measurements in representative catchment observatories.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRSDate not supplied

Evaluating Subgrid-Scale Properties of Low Clouds over the Beaufort Sea in Arctic System Reanalysis using ARISE Airborne In Situ Observations

Arctic low clouds influence the evolution of the Arctic system through their effects on radiative fluxes, boundary layer mixing, stability, turbulence, humidity, and precipitation. Unfortunately, atmospheric models and retrospective analysis (reanalysis) products struggle to accurately simulate the occurrence and properties of low clouds in the Arctic. To better understand the magnitude and possible causes of this problem, we analyze the disagreements between the Arctic System Reanalysis version 2 (ASR) and data taken during the September 2014 Arctic Radiation-IceBridge Sea and Ice Experiment (ARISE) airborne campaign conducted over the Beaufort Sea. Previous analysis of grid scale-averaged meteorological variables revealed a lack of cloud water produced in ASR, related to warm/dry biases, a failure to produce sufficient cloud water in conditions that are observed to be favorable for large cloud water values, and sampling issues with the ARISE flight paths. However, ARISE collected meteorological observations at one second intervals, offering spatial data resolution comparable to the subgrid scale of ASR. Subgrid scale properties and processes must be parameterized in conventional atmospheric models and reanalyses such as ASR, and often rely on assumed distributions of meteorological variables about the grid box mean state. To enable an examination of the realism of ASR’s subgrid scale properties, we use high resolution data from ARISE determine the “subgrid scale” distributions of cloud water and thermodynamic variables in the real atmosphere. Preliminary results show that the subgrid scale distribution of thermodynamic variables about the grid scale mean resembles the triangular distribution assumed in many bulk microphysical parameterization schemes, but also with tails resembling Gaussian distributions. The main goal is to understand the partitioning of the total available atmospheric water into water vapor and cloud water. Models typically covert some percentage of total water in a grid box to cloud water based on the estimated amount of supersaturation from the assumed vapor distribution. The ARISE measurements can reveal the relationship between subgrid scale supersaturation and the partition between cloud water and water vapor.

Jason Brant Dodson↗
DOE OSTI2022

BioRT-Flux-PIHM v1.0: a biogeochemical reactive transport model at the watershed scale

Watersheds are the fundamental Earth surface functioning units that connect the land to aquatic systems. Many watershed-scale models represent hydrological processes but not biogeochemical reactive transport processes. This has limited our capability to understand and predict solute export, water chemistry and quality, and Earth system response to changing climate and anthropogenic conditions. Here we present a recently developed BioRT-Flux-PIHM (BioRT hereafter) v1.0, a watershed-scale biogeochemical reactive transport model. The model augments the previously developed RT-Flux-PIHM that integrates land-surface interactions, surface hydrology, and abiotic geochemical reactions. It enables the simulation of (1) shallow and deep-water partitioning to represent surface runoff, shallow soil water, and deeper groundwater and of (2) biotic processes including plant uptake, soil respiration, and nutrient transformation. The reactive transport part of the code has been verified against the widely used reactive transport code CrunchTope. BioRT-Flux-PIHM v1.0 has recently been applied in multiple watersheds under diverse climate, vegetation, and geological conditions. This paper briefly introduces the governing equations and model structure with a focus on new aspects of the model. It also showcases one hydrology example that simulates shallow and deep-water interactions and two biogeochemical examples relevant to nitrate and dissolved organic carbon (DOC). These examples are illustrated in two simulation modes of complexity. One is the spatially lumped mode (i.e., two land cells connected by one river segment) that focuses on processes and average behavior of a watershed. Another is the spatially distributed mode (i.e., hundreds of cells) that includes details of topography, land cover, and soil properties. Whereas the spatially lumped mode represents averaged properties and processes and temporal variations, the spatially distributed mode can be used to understand the impacts of spatial structure and identify hot spots of biogeochemical reactions. The model can be used to mechanistically understand coupled hydrological and biogeochemical processes under gradients of climate, vegetation, geology, and land use conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Historically inconsistent productivity and respiration fluxes in the global terrestrial carbon cycle

The terrestrial carbon cycle is a major source of uncertainty in climate projections. Its dominant fluxes, gross primary productivity (GPP), and respiration (in particular soil respiration, R S ), are typically estimated from independent satellite-driven models and upscaled in situ measurements, respectively. We combine carbon-cycle flux estimates and partitioning coefficients to show that historical estimates of global GPP and R S are irreconcilable. When we estimate GPP based on R S measurements and some assumptions about R S :GPP ratios, we found the resulted global GPP values (bootstrap mean ${149}^{+29}_{-23}$ Pg C yr ⁻1 ) are significantly higher than most GPP estimates reported in the literature (${113}^{+18}_{-18}$ Pg C yr ⁻1 ). Similarly, historical GPP estimates imply a soil respiration flux (Rs GPP , bootstrap mean of ${68}^{+10}_{-8}$ Pg C yr ⁻1 ) statistically inconsistent with most published R S values (${87}^{+9}_{-8}$ Pg C yr ⁻1 ), although recent, higher, GPP estimates are narrowing this gap. Furthermore, global R S :GPP ratios are inconsistent with spatial averages of this ratio calculated from individual sites as well as CMIP6 model results. This discrepancy has implications for our understanding of carbon turnover times and the terrestrial sensitivity to climate change. Future efforts should reconcile the discrepancies associated with calculations for GPP and Rs to improve estimates of the global carbon budget.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2010

Changes in the Mass Balance of the Greenland Ice Sheet in a Warming Climate During 2003-2009

Mass changes of the Greenland ice sheet (GIS) derived from ICESat and GRACE data both show that the net mass loss from GIS during 2003-2009 is about 175 Gt/year, which contributes 0.5mm/yr global sea-level rise. The rate of mass loss has increased significantly since the 1990's when the GIS was close to mass balance. Even though the GIS was close to mass balance during the 1990's, it was already showing characteristics of responding to8 warmer climate, specifically thinning at the margins and thickening inland at higher elevations. During 2003-2009, increased ice thinning due to increases in melting and acceleration of outlet glaciers began to strongly exceed the inland thickening from increases in accumulation. Over the entire GIS, the mass loss between the two periods, from increased melting and ice dynamics, increased by about 190 Gt/year while the mass gain, from increased precipitation and accumulation, increased by only about 15Gt/year. These ice changes occurred during a time when the temperature on GIS changed at rate of about 2K/decade. The distribution of elevation and mass changes derived from ICESat have high spatial resolution showing details over outlet glaciers, by drainage systems, and by elevation. However, information on the seasonal cycle of changes from ICESat data is limited, because the ICESat lasers were only operated during two to three campaigns per year of about 35 days duration each. In contrast, the temporal resolution of GRACE data, provided by the continuous data collection, is much better showing details of the seasonal cycle and the inter-annual variability. The differing sensitivity of the ICESat altimetry and the GRACE gravity methods to motion of the underlying bedrock from glacial isostatic adjustment (GIA) is used to evaluate the GIA corrections provided by models. The two data types are also combined to make estimates of the partitioning of the mass gains and losses among accumulation, melting, and ice discharge from outlet glaciers.

Zwally, H. Jay↗
NASA NTRS2019

Effects of Volcanic Emissions on Clouds During Kilauea Degassing Events

Aerosols influence Earth’s radiative balance directly by scattering and absorbing solar radiation, and indirectly by modifying cloud properties. Current scientific consensus indicates that these effects may offset as much as 50% of the warming due to greenhouse gas emissions. Over the last two decades dramatic volcanic events in Hawaii have produced localized aerosol emissions in otherwise clean environments. These are “natural experiments" where the aerosol effects on clouds and climate can be partitioned from other effects like meteorology and industrial emissions. Therefore, these events provide a unique opportunity to learn about possible effects of aerosol pollution on climate through cloud modification. In this work we use the version 5 of the NASA Goddard Earth Observing System (GEOS-5) and satellite retrievals to analyze and evaluate the strength of the aerosol indirect effect on liquid and ice clouds during the 2008 and 2018 Kilauea degassing events using different emissions scenarios (0, 1, and 5 actual emissions). Our results suggested that the 2018 event was stronger and more regionally significant with respect to cloud formation process for both liquid and ice clouds, while the 2008 affected local liquid clouds only. GEOS-5 predictions reproduced spatial patterns for all parameters, however better precision could be gained by using more accurate plume parameters for height and ash concentration.

Breen, K. H.↗
DOE OSTI2020

LLNL Kimberlina 1.2 NUFT Simulations June 2018 (v2)

This dataset contains the output 6,000, 3-dimensional reactive multi-phase flow and transport aquifer simulations of brine and CO2 leakage into a protective aquiver in California’s San Joaquin Valley and input data files detailing the geologic mesh, aquifer physical properties and CO2 and brine injection rates. This data set was generated as an ongoing effort with the US DOE National Risk Assessment Partnership (NRAP) to evaluate the effectiveness of monitoring techniques to detect brine and CO2 leakage from legacy wells into underground sources of drinking water overlaying a CO2 storage reservoir. Each simulation contains a unique set of input parameters, generated stochastically. The outputs consist of these upper three geologic layers (from top): the Etchegoin, Macoma-Chanac, Santa Margarita-McLure formations. These simulations span the several distances (1, 3 and 6 km or wells W31-0.2, W31-0.5 and W31-1.0, respectively) from the CO2 injector, initiated from bottom hole pressure and saturation to calculate wellbore leakage from the storage reservoir, with low and high regional groundwater gradients and wellbore leakage into 5 leaky nodes. The dataset includes 1,000 unique simulations for each distance, which each contain a unique aquifer heterogeneity, aquifer and caprock permeability, and two model generations are included with a high permeability (prod07) and hybrid permeability (prod09). The range of permeability distributions is listed in Table 1. Each model generation consists of 3,000 simulations. Included in the dataset are the leakage rates determined from 2D wellbore models which utilize the pressure and CO2 saturation from LBL's reservoir simulations, NUFT mesh files with distributed lithology, NUFT rocktab files which describe the material properties for the geologic layers and the NUFT input files and post-processed output 'ntab' files. Each ntab file contains spatial (rows) and temporal (columns) model output tables for each model cell, the locations (x,y,z) and dimensions for each cells (dx, dy, dz). Table 1. Permeability distribution ranges for prod07 and prod09 model generations Geologic Layer: Permeability Range (log10 m^2) prod07 prod09 Etchegoin -12.92 to -10.92 -13.70 to -11.44 Macoma-Chanac -12.72 to -10.72 -13.50 to -11.24 Santa Margarita-McLure -12.70 to -10.70 -13.48 to -11.22 The input files used to generate the model include which are included in the dataset are: Time series of CO2 leakage input into the model (ex: Q_brn.W31-0.2.sim1000.layers123.tab) Time series of CO2 leakage input into the model (ex: Q_CO2.W31-0.2.sim1000.layers123.tab) Physical properties of the aquifer materials detailing the aquifer porosity, solid density, partitioning coefficients, permeabilities and van-Genuchten parameters detailed in a NUFT rocktab file: (ex: sim1000.usnt.rocktab) Numerical mesh and geologic data assigned to each model cell detailed in a NUFT genmsh format (ex: sim1000.mesh_k16.prod07.trans.genmsh) The primary output parameters are: pH (use absolute value) Change in TDS (mg/kg) Change in Pressure (Pa) Change CO2 gas saturation (fraction range 0.0-1.0) for example, the directory /p/lscratchh/mansoor1/nrap/kimberlina/prod09/mainfiles/sim1000/W31- 0.2 contains: sim1000.W31-0.2.trans.pH.red.ntab sim1000.W31-0.2.no_bg.trans.TDS.red.ntab sim1000.W31-0.2.usnt.P.deltabg.red.ntab sim1000.W31-0.2.usnt.CO2_sat.deltabg.red.ntab Each row in the NTAB files consist of model output per numerical grid cell. Each output file contains 33 columns (variables), including the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the ith grid of y-axis 4. k: the ith grid of z-axis 5. element_ref: element reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The remainder (14, 15, 16...) variables are the simulated parameter values over time, take Pressure as an example, are: 14. 0.0y: initial pressure per cell. 15. 10.0y: simulated pressure at the end of the 10th year. 16. 20.0y: simulated pressure at the end of the 20th year. ... (repeated for every 10 years until 200 years)... The model extends 10,000 m, 5,000 m and 1,411 m in the x,y and z dimensions, respectively. The mesh consists of 164,832 cells with mesh dimensions of 101 x 51 x 32 (nx, ny, nz), with cell dimensions ranging from 100 m laterally (along x and y-axis) and model layers are as designated in the z-axis: Layer 1: atmosphere (1e-30 m thick) Layer 2: upper caprock (10 m thick) Layers 3-13: Etchegoin (536.23 m thck) Layers 14-27: Macoma-Chanac (679.04 m thick) Layers 28-32: Santa Margarita-McLure (185.94 m thick) The wellbore is placed along node i=51, j=26, and extends vertically along 5 nodes from the top to the bottom of the model. Special instructions when extracting files: Each Gzip archive (ex: prod07.sim1000-sim00099.tar.gz) contains 100 simulations. Gzip archives should be transferred into base directories (ie. In Linux: mkdir prod07; mv prod07.*.tar.gz prod07/.) before extracting, or files will be overwritten. Each sub-simulation tree should have the following file structure pattern (using the linux 'tree' command): |-- prod07 | |-- sim0001 | |-- W31-0.2 | | |-- Q_brn.W31-0.2.sim0001.layers123.tab | | |-- Q_co2.W31-0.2.sim0001.layers123.tab | | |-- sim0001.W31-0.2.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.2.trans.pH.red.ntab | | |-- sim0001.W31-0.2.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.2.usnt.P.deltabg.red.ntab | |-- W31-0.5 | | |-- Q_brn.W31-0.5.sim0001.layers123.tab | | |-- Q_co2.W31-0.5.sim0001.layers123.tab | | |-- sim0001.W31-0.5.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.5.trans.pH.red.ntab | | |-- sim0001.W31-0.5.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.5.usnt.P.deltabg.red.ntab | |-- W31-1.0 | | |-- Q_brn.W31-1.0.sim0001.layers123.tab | | |-- Q_co2.W31-1.0.sim0001.layers123.tab | | |-- sim0001.W31-1.0.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-1.0.trans.pH.red.ntab | | |-- sim0001.W31-1.0.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-1.0.usnt.P.deltabg.red.ntab | |-- sim0001.mesh_k16.prod07.trans.genmsh Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. This report was reviewed and released as LLNL-MI-753464.

aquifer↗
NASA NTRS2007

Estimation of Forest Fuel Load from Radar Remote Sensing

Understanding fire behavior characteristics and planning for fire management require maps showing the distribution of wildfire fuel loads at medium to fine spatial resolution across large landscapes. Radar sensors from airborne or spaceborne platforms have the potential of providing quantitative information about the forest structure and biomass components that can be readily translated to meaningful fuel load estimates for fire management. In this paper, we used multifrequency polarimetric synthetic aperture radar imagery acquired over a large area of the Yellowstone National Park (YNP) by the AIRSAR sensor, to estimate the distribution of forest biomass and canopy fuel loads. Semi-empirical algorithms were developed to estimate crown and stem biomass and three major fuel load parameters, canopy fuel weight, canopy bulk density, and foliage moisture content. These estimates when compared directly to measurements made at plot and stand levels, provided more than 70% accuracy, and when partitioned into fuel load classes, provided more than 85% accuracy. Specifically, the radar generated fuel parameters were in good agreement with the field-based fuel measurements, resulting in coefficients of determination of R(sup 2) = 85 for the canopy fuel weight, R(sup 2)=.84 for canopy bulk density and R(sup 2) = 0.78 for the foliage biomass.

Saatchi, Sassan↗
DOE OSTI2020

Time of Emergence and Large Ensemble Intercomparison for Ocean Biogeochemical Trends

Anthropogenically forced changes in ocean biogeochemistry are underway and critical for the ocean carbon sink and marine habitat. Detecting such changes in ocean biogeochemistry will require quantification of the magnitude of the change (anthropogenic signal) and the natural variability inherent to the climate system (noise). Here we use Large Ensemble (LE) experiments from four Earth system models (ESMs) with multiple emissions scenarios to estimate Time of Emergence (ToE) and partition projection uncertainty for anthropogenic signals in five biogeochemically important upper-ocean variables. We find ToEs are robust across ESMs for sea surface temperature and the invasion of anthropogenic carbon; emergence time scales are 20–30 yr. For the biological carbon pump, and sea surface chlorophyll and salinity, emergence time scales are longer (50+ yr), less robust across the ESMs, and more sensitive to the forcing scenario considered. We find internal variability uncertainty, and model differences in the internal variability uncertainty, can be consequential sources of uncertainty for projecting regional changes in ocean biogeochemistry over the coming decades. In combining structural, scenario, and internal variability uncertainty, this study represents the most comprehensive characterization of biogeochemical emergence time scales and uncertainty to date. Our findings delineate critical spatial and duration requirements for marine observing systems to robustly detect anthropogenic change.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2015

Impact of Asian Aerosols on Precipitation Over California: An Observational and Model Based Approach

Dust and pollution emissions from Asia are often transported across the Pacific Ocean to over the western United States. Therefore, it is essential to fully understand the impact of these aerosols on clouds and precipitation forming over the eastern Pacific and western United States, especially during atmospheric river events that account for up to half of California's annual precipitation and can lead to widespread flooding. In order for numerical modeling simulations to accurately represent the present and future regional climate of the western United States, we must account for the aerosol-cloud-precipitation interactions associated with Asian dust and pollution aerosols. Therefore, we have constructed a detailed study utilizing multi-sensor satellite observations, NOAA-led field campaign measurements, and targeted numerical modeling studies where Asian aerosols interacted with cloud and precipitation processes over the western United States. In particular, we utilize aerosol optical depth retrievals from the NASA Moderate Resolution Imaging Spectroradiometer (MODIS), NOAA Geostationary Operational Environmental Satellite (GOES-11), and Japan Meteorological Agency (JMA) Multi-functional Transport Satellite (MTSAT) to effectively detect and monitor the trans-Pacific transport of Asian dust and pollution. The aerosol optical depth (AOD) retrievals are used in assimilating the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem) in order to provide the model with an accurate representation of the aerosol spatial distribution across the Pacific. We conduct WRF-Chem model simulations of several cold-season atmospheric river events that interacted with Asian aerosols and brought significant precipitation over California during February-March 2011 when the NOAA CalWater field campaign was ongoing. The CalWater field campaign consisted of aircraft and surface measurements of aerosol and precipitation processes that help extensively validate our WRF-Chem model simulations. After validating the capability of the WRF-Chem in realistically simulating the aerosol-cloud precipitation interactions, we conduct sensitivity studies where the AOD is doubled to diagnose whether an increasing concentration of Asian aerosols over the western United States will lead to further impacts on the cloud and precipitation processes over California. We also perform sensitivity studies where the aerosols will be partitioned into dust-only and pollution-only in order to separate the impacts of the differing Asian aerosol species. The results of our WRF-Chem model simulations aim to show that the trans-Pacific transport of Asian aerosols influence the precipitation associated with atmospheric river events that can ultimately impact the regional climate of the western United States. 1 University

Naeger, Aaron R.↗
DOE OSTI2020

Data-driven surrogates for high dimensional models using Gaussian process regression on the Grassmann manifold

This paper introduces a surrogate modeling scheme based on Grassmannian manifold learning to be used for cost-efficient predictions of high-dimensional stochastic systems. The method exploits subspace-structured features of each solution by projecting it onto a Grassmann manifold. This point-wise linear dimensionality reduction harnesses the structural information to assess the similarity between solutions at different points in the input parameter space. The method utilizes a solution clustering approach in order to identify regions of the parameter space over which solutions are sufficiently similarly such that they can be interpolated on the Grassmannian. In this clustering, the reduced-order solutions are partitioned into disjoint clusters on the Grassmann manifold using the eigen-structure of properly defined Grassmannian kernels and, the Karcher mean of each cluster is estimated. Then, the points in each cluster are projected onto the tangent space with origin at the corresponding Karcher mean using the exponential mapping. For each cluster, a Gaussian process regression model is trained that maps the input parameters of the system to the reduced solution points of the corresponding cluster projected onto the tangent space. Using this Gaussian process model, the full-field solution can be efficiently predicted at any new point in the parameter space. In certain cases, the solution clusters will span disjoint regions of the parameter space. In such cases, for each of the solution clusters we utilize a second, density-based spatial clustering to group their corresponding input parameter points in the Euclidean space. The proposed method is applied to two numerical examples. Here, the first is a nonlinear stochastic ordinary differential equation with uncertain initial conditions where the surrogate is used to predict the time history solution. The second involves modeling of plastic deformation in a model amorphous solid using the Shear Transformation Zone theory of plasticity, where the proposed surrogate is used to predict the full strain field of a material specimen under large shear strains.

42 ENGINEERING↗