Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Accounting”

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 919 records · Page 51

How Equivalent Are Equivalent Porous Media?

Abstract Geoenergy and geoengineering applications usually involve fluid injection into and production from fractured media. Accounting for fractures is important because of the strong poromechanical coupling that ties pore pressure changes and deformation. A possible approach to the problem uses equivalent porous media to reduce the computational cost and model complexity instead of explicitly including fractures in the models. We investigate the validity of this simplification by comparing these two approaches. Simulation results show that pore pressure distribution significantly differs between the two approaches even when both are calibrated to predict identical values at the injection and production wells. Additionally, changes in fracture stability are not well captured with the equivalent porous medium. We conclude that explicitly accounting for fractures in numerical models may be necessary under some circumstances to perform reliable coupled thermohydromechanical simulations, which could be used in conjunction with other tools for induced seismicity forecasting.

58 GEOSCIENCES↗

Disentangling the Impact of the COVID‐19 Lockdowns on Urban NO 2 From Natural Variability

Abstract TROPOMI satellite data show substantial drops in nitrogen dioxide (NO 2 ) during COVID‐19 physical distancing. To attribute NO 2 changes to NO x emissions changes over short timescales, one must account for meteorology. We find that meteorological patterns were especially favorable for low NO 2 in much of the United States in spring 2020, complicating comparisons with spring 2019. Meteorological variations between years can cause column NO 2 differences of ~15% over monthly timescales. After accounting for solar angle and meteorological considerations, we calculate that NO 2 drops ranged between 9.2% and 43.4% among 20 cities in North America, with a median of 21.6%. Of the studied cities, largest NO 2 drops (>30%) were in San Jose, Los Angeles, and Toronto, and smallest drops (<12%) were in Miami, Minneapolis, and Dallas. These normalized NO 2 changes can be used to highlight locations with greater activity changes and better understand the sources contributing to adverse air quality in each city.

54 ENVIRONMENTAL SCIENCES↗

Effects of organized convection parameterization on the MJO and precipitation in E3SMv1. Part I: Mesoscale heating

Mesoscale organization of convection is typically not represented in global circulation models, and hence its influence on the global circulation is not accounted for. The heating component of a parameterization that represents the dynamical and physical effects of circulations associated with organized convection, referred to as the multiscale coherent structure parameterization (MCSP), is implemented in the Energy Exascale Earth System Model version 1 (E3SMv1). Numerical simulations are conducted to assess its impact on the simulated climate. Besides E3SMv1 simulations, we performed high-resolution (1 km) simulations using the Weather Research and Forecasting (WRF) Model to determine the temperature tendencies induced by mesoscale convective systems embedded in deep convection. We tuned the free parameters of the MCSP based on the WRF simulations. MCSP heating enhances Kevin wave spectra in E3SMv1, improves the representation of the Madden-Julian Oscillation, and impacts on mean precipitation that vary with locations and seasons, including a reduction of biases over the tropical Pacific.

Climate Modeling↗

Influences on Discharge Partitioning on a Large River Delta: Case Study of the Mississippi-Atchafalaya Diversion, 1926-1950

The modern Mississippi River (M.R.) Delta is plumbed by the Mississippi and Atchafalaya rivers, setting water and sediment dispersal pathways for Earth’s fifth-largest river system. The Atchafalaya River’s (A.R.) partial annexation of discharge from the M.R., particularly between 1926 and 1950, prompted warnings of a rapid river avulsion and the construction of the Old River Control Structure to regulate flow. Natural and anthropogenic causes of this discharge annexation are difficult to disentangle. Here, we develop and validate a hydrodynamic model of flow partitioning through the historic channel network. We then isolate how several key changes to the system affected discharge partitioning and stage at the A.R.-M.R. diversion. Simulations show that erosion of the upper A.R. can account for 73% of the water discharge increase. Dredging in the lower A.R. between 1932 and 1950 can account for 35% of the water discharge increase, and was also an important control on shear stress distribution. The lower M.R. was slightly erosional during this period, and therefore hindered the A.R. discharge increase. Significant lacustrine delta deposition in A.R. had little effect on partitioning. These findings highlight the importance of A.R. enlargement processes on avulsion dynamics at this site. Given the essential nature of this river junction to the society, transportation, and commerce of the United States, improved attribution of discharge increases may lead to future management strategies that are broadly impactful.

54 ENVIRONMENTAL SCIENCES↗

The Lithophile Element Budget of Earth's Core

Abstract The relative composition of Earth's core and mantle were set during core formation. By determining how elements partition between metal and silicate at high pressures and temperatures, measurements of the mantle composition and geophysical observations of the core can be used to understand the mechanisms by which Earth formed. Here we present the results of metal‐silicate partitioning experiments for a range of nominally lithophile elements (Al, Ca, K, Mg, O, Si, Th, and U) and S to 85 GPa and up to 5400 K. With our results and a compilation of literature data, we developed a parameterization for partitioning that accounts for compositional dependencies in both the metal and silicate phases. Using this parameterization in a range of planetary growth models, we find that, in general, lithophile element partitioning into the metallic phase is enhanced at high temperatures. The relative abundances of FeO, SiO 2 , and MgO in the mantle vary significantly between planetary growth models, and the mantle abundances of these elements can be used to provide important constraints on Earth's accretion. To match Earth's core mass and mantle composition, Earth's building blocks must have been enriched in Fe and depleted in Si compared with CI chondrites. Finally, too little Mg, Si, and O are partitioned into the core for precipitation of oxides to be a major source of energy for the geodynamo. In contrast, several ppb of U can be partitioned into the core at high temperatures, and this energy source must be accounted for in thermal evolution models.

Chidester, B. A.↗

Estimating the Acute Health Impacts of Fire-Originated PM 2.5 Exposure During the 2017 California Wildfires: Sensitivity to Choices of Inputs

Exposure to wildfire smoke increases the risk of respiratory and cardiovascular hospital admissions. Health impact assessments, used to inform decision-making processes, characterize the health impacts of environmental exposures by combining preexisting epidemiological concentration–response functions (CRFs) with estimates of exposure. These two key inputs influence the magnitude and uncertainty of the health impacts estimated, but for wildfire-related impact assessments the extent of their impact is largely unknown. We first estimated the number of respiratory, cardiovascular, and asthma hospital admissions attributable to fire-originated PM 2.5 exposure in central California during the October 2017 wildfires, using Monte Carlo simulations to quantify uncertainty with respect to the exposure and epidemiological inputs. We next conducted sensitivity analyses, comparing four estimates of fire-originated PM 2.5 and two CRFs, wildfire and nonwildfire specific, to understand their impact on the estimation of excess admissions and sources of uncertainty. We estimate the fires accounted for an excess 240 (95% CI: 114, 404) respiratory, 68 (95% CI: -10, 159) cardiovascular, and 45 (95% CI: 18, 81) asthma hospital admissions, with 56% of admissions occurring in the Bay Area. Although differences between impact assessment methods are not statistically significant, the admissions estimates' magnitude is particularly sensitive to the CRF specified while the uncertainty is most sensitive to estimates of fire-originated PM 2.5 . Not accounting for the exposure surface's uncertainty leads to an underestimation of the uncertainty of the health impacts estimated. Employing context-specific CRFs and using accurate exposure estimates that combine multiple data sets generates more certain estimates of the acute health impacts of wildfires.

58 GEOSCIENCES↗

Evidence for the 3D Radiative Effects of Boundary-Layer Clouds From Observations of Direct and Diffuse Surface Solar Fluxes

Numerical experiments have revealed the importance of horizontal transport of light in the presence of clouds (“3D effects”), with consequences for climate, weather, and solar resource availability predictions. Yet, analysis of 3D effects from observations remain sparse because of the difficulty to isolate the effect of horizontal transport in radiation measurements. In this study, we provide observational evidence for 3D effects based on the direct-diffuse partition of surface solar fluxes. It is compared to outputs from the ecRad radiative transfer scheme run on retrieved cloud profiles. Furthermore, the direct-beam calculation takes careful account of the field-of-view of the pyrheliometer to ensure consistency between observed and modeled direct fluxes. Only the solver that accounts for 3D effects is able to reproduce the observed mean direct-diffuse partition as a function of solar zenith angle and cloud cover, in particular at large solar zenith angles where cloud sides intercept most of the direct beam.

54 ENVIRONMENTAL SCIENCES↗

Dissipation and Bathymetric Sensitivities in an Unstructured Mesh Global Tidal Model

Abstract The mechanisms and geographic distribution of global tidal dissipation in barotropic tidal models are examined using a high resolution unstructured mesh finite element model. Mesh resolution varies between 2 and 25 km and is especially focused on inner shelves and steep bathymetric gradients. Tidal response sensitivities to bathymetric changes are examined to put into context response sensitivities to frictional processes. We confirm that the Ronne Ice Shelf dramatically affects Atlantic tides but also find that bathymetry in the Hudson Bay system is a critical control. We follow a sequential frictional parameter optimization process and use TPXO9 data‐assimilated tidal elevations as a reference solution. From simulated velocities and depths, dissipation within the global model is estimated and allows us to pinpoint dissipation at high resolution. Boundary layer dissipation is extremely focused with 1.4% of the ocean accounting for 90% of the total. Internal tide friction is much more distributed with 16.7% of the ocean accounting for 90% of the total. Often highly regional dissipation can impact basin‐scale and even ocean wide tides. Optimized boundary layer friction parameters correlate very well with the physical characteristics of the locality with high friction factors associated with energetic tidal regions, deep ocean island chains, and ice covered areas. Global complex M 2 tide errors are 1.94 cm in deep waters. Total global boundary layer and internal tide dissipation are estimated, respectively, at 1.83 and 1.49 TW. This continues the trend in the literature toward attributing more dissipation to internal tides.

54 ENVIRONMENTAL SCIENCES↗

The Impact of Sudden Commencements on Ground Magnetic Field Variability: Immediate and Delayed Consequences

Abstract We examine how Sudden Commencements (SCs) and Storm Sudden Commencements (SSCs) influence the occurrence of high rates of change of the magnetic field ( R ) as a function of geomagnetic latitude. These rapid, high amplitude variations in the ground‐level geomagnetic field pose a significant risk to ground infrastructure, such as power networks, as the drivers of geomagnetically induced currents. We find that rates of change of ∼30 nT min −1 at near‐equatorial stations are up to 700 times more likely in an SC than in any random interval. This factor decreases with geomagnetic latitude such that rates of change around 30 nT min −1 are only up to 10 times more likely by 65°. At equatorial latitudes we find that 25% of all R in excess of 50 nT min −1 occurs during SCs. This percentage also decreases with geomagnetic latitude, reaching ≤1% by 55°. However, the time period from the SC to 3 days afterward accounts for ≥90% of geomagnetic field fluctuations over 50 nT min −1 , up to ∼60° latitude. Above 60°, other phenomena such as isolated substorms account for the majority of large R . Furthermore, the elevated rates of change observed during and after SCs are solely due to those classified as SSCs. These results show that SSCs are the predominant risk events for large R at mid and low latitudes, but that the risk from the SC itself decreases with latitude.

Smith, Andrew W.↗

Process Interactions Can Change Process Ranking in a Coupled Complex System Under Process Model and Parametric Uncertainty

For a complex hydrologic system with multiple processes and process interactions, global sensitivity analysis is often used to identify important or influential parameters for model development and improvement. The identification is complicated by process model uncertainty, when a system process can be represented by multiple process models. This study develops a new total-effect process sensitivity index to identify influential processes under model uncertainty. This is done by extending Sobol's total-effect parameter sensitivity index for one system model to total-effect process sensitivity index for multiple system models to account for uncertainty in process models and model parameters. The total-effect process sensitivity index includes not only the first-order process sensitivity index for measuring the importance of individual processes but also higher-order indices that account for process interactions. The total-effect process sensitivity index can identify an influential process that itself and its interactions with other processes influence a model output. Here, the total-effect process sensitivity index is applied to two numerical examples: (a) Sobol's G*-functions with analytical solutions of first-order and total-effect process sensitivity indices, and (b) groundwater flow models with interactions between recharge, geology, and snowmelt processes. The second evaluation shows that, due to second-order and higher-order process interactions, the first-order and total-effect process sensitivity indices give different process ranking. It is thus necessary to estimate both first-order and total-effect process sensitivity indices to appreciate the difference between the first-order impact of a process alone and the overall total-effect impact of the process itself and its interactions with other processes on a model output.

54 ENVIRONMENTAL SCIENCES↗

Abundant Nitrate and Nitric Acid Aerosol in the Upper Troposphere and Lower Stratosphere

The tropospheric and stratospheric nitrate aerosol is simulated by a sectional aerosol model coupled to the Community Earth System Model. The simulated nitrate mass fractional contribution to aerosols is significantly higher in the upper troposphere and lower stratosphere (UTLS) than that at the surface. Both in situ measurements and simulations show that nitrate aerosol accounts for about 30%–40% of the aerosol mass at the tropopause of the Asian summer monsoon (ASM) region. Furthermore, simulated condensed nitric acid particles account for ~20% of the annual mean aerosol mass at the tropical tropopause, and over 95% in the UTLS at the South Pole in June-July-August. Here, our study suggests that the extremely cold ambient conditions in the UTLS of the tropics, ASM and polar regions thermodynamically favor the condensation of ammonia and nitric acid. The widely distributed nitrate aerosol in the global UTLS may be overlooked by climate models.

54 ENVIRONMENTAL SCIENCES↗

Wave and Particle Analysis of Z-Mode and O-Mode Emission in the Jovian Inner Magnetosphere

We report some of the most intense Z-mode and O-mode observations obtained by the Juno spacecraft while in orbit about Jupiter in a low to mid-latitude region near the inner edge of the Io torus. We have been able to estimate the density of the plasma in this region based on the lower frequency cutoff of the observed Z-mode emission. The results are compatible with the electron density measurements of the Jovian Auroral Distributions Experiment (JADE), on board the Juno spacecraft, if we account for unmeasured cold plasma. Direction-finding measurements indicate that the Z- and O-mode emission have distinct source regions. We have also used the measured phase space density of the JADE and the Jupiter energetic particle detector instruments to calculate estimated local growth rates of the observed O-mode and Z-mode emission assuming a loss cone instability and quasilinear analysis. The results suggest the emissions were observed near, but not within, a source region, and the free energy source is consistent with a loss cone. We have thus carried out the quasilinear wave analysis of the assumed remote Z- and O-mode wave growths. It is shown that the remotely generated waves, propagated through an inhomogeneous medium to the satellite location, may account for the observed wave characteristics. The importance of Z-mode in accelerating electrons in the inner Jovian magnetosphere makes these new wave mode confirmations at Jupiter of particular interest.

79 ASTRONOMY AND ASTROPHYSICS↗

Global Barotropic Tide Modeling Using Inline Self‐Attraction and Loading in MPAS‐Ocean

Abstract We examine ocean tides in the barotropic version of the Model for Prediction Across Scales (MPAS‐Ocean), the ocean component of the Department of Energy Earth system model. We focus on four factors that affect tidal accuracy: self‐attraction and loading (SAL), model resolution, details of the underlying bathymetry, and parameterized topographic wave drag. The SAL term accounts for the tidal loading of Earth's crust and the self‐gravitation of the ocean and the load‐deformed Earth. A common method for calculating SAL is to decompose mass anomalies into their spherical harmonic constituents. Here, we compare a scalar SAL approximation versus an inline SAL using a fast spherical harmonic transform package. Wave drag accounts for energy lost by breaking internal tides that are produced by barotropic tidal flow over topographic features. We compare a series of successively finer quasi‐uniform resolution meshes (62.9, 31.5, 15.7, and 7.87 km) to a variable resolution (45 to 5 km) configuration. We ran MPAS‐Ocean in a single‐layer barotropic mode forced by five tidal constituents. The 45 to 5 km variable resolution mesh obtained the best total root‐mean‐square error (5.4 cm) for the deep ocean ( 1,000 m) tide compared to TPXO8 and ran twice as fast as the quasi‐uniform 8 km mesh, which had an error of 5.8 cm. This error is comparable to those found in other forward (non‐assimilative) ocean tide models. In future work, we plan to use MPAS‐Ocean to study tidal interactions with other Earth system components, and the tidal response to climate change.

54 ENVIRONMENTAL SCIENCES↗

A Consistent Representation of Cloud Overlap and Cloud Subgrid Vertical Heterogeneity

Many global climate models underestimate the cloud cover and overestimate the cloud albedo, especially for low-level clouds. We determine how a correct representation of the vertical structure of clouds can fix part of this bias. We use the 1D McICA framework and focus on low-level clouds. Using Large Eddy Simulations results as reference, we propose a method based on exponential-random overlap that represents the cloud overlap between layers and the subgrid cloud properties over several vertical scales, with a single value of the overlap parameter. Starting from a coarse vertical grid, representative of atmospheric models, this algorithm is used to generate the vertical profile of the cloud fraction with a finer vertical resolution, or to generate it on the coarse grid but with subgrid heterogeneity and cloud overlap that ensures a correct cloud cover. Doing so we find decorrelation lengths are dependent on the vertical resolution, except if the vertical subgrid heterogeneity and interlayer overlap are taken into account coherently. We confirm that the frequently used maximum-random overlap leads to a significant error by underestimating the low-level cloud cover with a relative error of about 50%, that can lead to an error of SW cloud albedo as big as 70%. Not taking into account the subgrid vertical heterogeneity of clouds can cause a relative error of 20% in brightness, assuming the cloud cover is correct.

54 ENVIRONMENTAL SCIENCES↗

Modeling the Contribution of Precipitation Loss to a Radiation Belt Electron Dropout Observed by Van Allen Probes

A drift-diffusion model is used to simulate the low-altitude electron distribution, accounting for azimuthal drift, pitch angle diffusion, and atmospheric backscattering effects during a rapid electron dropout event on 21st August 2013, at L = 4.5. Additional external loss effects are introduced during times when the low-altitude electron distribution cannot be reproduced by diffusion alone. Here, the model utilizes low-altitude electron count rate data from five POES/MetOp satellites to quantify pitch angle diffusion rates. Low-altitude data provides critical constraint on the model because it includes the drift loss cone region where the electron distribution in longitude is highly dependent on the balance between azimuthal drift and pitch angle diffusion. Furthermore, a newly derived angular response function for the detectors onboard POES/MetOp is employed to accurately incorporate the bounce loss cone measurements, which have been previously contaminated by electrons from outside the nominal field-of-view. While constrained by low-altitude data, the model also shows reasonable agreement with high-altitude data. Pitch angle diffusion rates during the event are quantified and are faster at lower energies. Precipitation is determined to account for all of the total loss observed for 450 keV electrons, 88% for 600 keV and 38% for 900 keV. Predictions made in the MeV range are deemed unreliable as the integral energy channels E3 and P6 fail to provide the necessary constraint at relativistic energies.

79 ASTRONOMY AND ASTROPHYSICS↗

Subseasonal Clustering of Atmospheric Rivers Over the Western United States

Abstract The serial occurrence of atmospheric rivers (ARs) along the US West Coast can lead to prolonged and exacerbated hydrologic impacts, threatening flood‐control and water‐supply infrastructure due to soil saturation and diminished recovery time between storms. Here a statistical approach for quantifying subseasonal temporal clustering among extreme events is applied to a 41‐year (1979–2019) wintertime AR catalog across the western United States (US). Observed AR occurrence, compared against a randomly distributed AR timeseries with the same average event density, reveals temporal clustering at a greater‐than‐random rate across the western US with a distinct geographical pattern. Compared to the Pacific Northwest, significant AR clusters over the northern Coastal Range of California and Sierra Nevada are more frequent and occur over longer time periods. Clusters along the California Coastal Range typically persist for 2 weeks, are composed of 4–5 ARs per cluster, and account for over 85% of total AR occurrence. Across the northwest Coast‐Cascade Ranges, clusters account for ∼50% of total AR occurrence, typically last 8–10 days, and contain 3–4 individual AR events. Based on precipitation data from a high‐resolution dynamical downscaling of reanalysis, the fractions of total and extreme hourly precipitation attributable to AR clusters are largest along the northern California coast and in the Sierra Nevada. Interannual variability among clusters highlights their importance for determining whether a particular water year is anomalously wet or dry. The mechanisms behind this unusual clustering are unclear and require further research.

Meteorology & Atmospheric Sciences↗

Significance of Low‐Velocity Zones on Solute Retention in Rough Fractures

Natural fractures are characterized by high internal heterogeneity. This internal variability is the cause of flow channeling, which in turn leads to contaminant transport taking place primarily along the high-velocity channels. Mass exchange between the high-velocity channels and the low-velocity zones has the potential to enhance contaminant retention, due to solute diffusion into the low-velocity zones and subsequent exposure to additional surface area for diffusion into the bordering rock matrix. Here, we derive a random walk particle tracking method for heterogeneous fractures, which includes an additional term to account for the aperture gradient. The method takes into account advection, diffusion in the fracture and matrix diffusion. The developed numerical framework is applied to assess the effect of low-velocity zones in rough self-affine fractures. The results show that diffusion into low-velocity zones has a visible but modest impact on contaminant retention. The magnitude of this impact does not change considerably, regardless of whether diffusion into the rock matrix is considered in the model, and increases for a decreasing average Péclet number of the fracture.

58 GEOSCIENCES↗

The Mystery of Baselevel Controls in the Incision History of the Central Colorado Plateau

Erosion can remain active and changing in landscapes long after tectonic drivers have ceased, potentially due to local-geologic controls, climate changes, or geodynamics. We present new fluvial incision-rate histories and terrain analyses of the Colorado River system through the central Colorado Plateau to understand what has caused the variable erosion across this post-orogenic landscape. Results from new cosmogenic and luminescence dating of fluvial terrace and upland gravel deposits in Glen and Meander Canyons establish incision-rate histories that are marked by an Early-Middle Pleistocene erosion hiatus, followed by ~200 m of rapid incision over the last ~350 kyr. Projection of fluvial topography from above knickzones of the Colorado River drainage system roughly agree with the observed magnitude of recent incision and reflect a common baselevel fall from Pliocene river integration through Grand Canyon, which is still propagating through the drainage. A response-time model indicates that baselevel fall from integration likely took 2–4 Myr to reach the central Colorado Plateau and 100s kyr to travel across the study area, potentially accounting for incision rate changes in the fluvial terrace records of Meander and Glen Canyons. The upstream-migrating incision has likely been partitioned into multiple waves across the landscape due to the local geologic controls of lava damming, salt tectonics, and heterogenous bedrock. As baselevel fall from Pliocene Colorado River integration diffuses upstream, it can only account for perhaps a quarter of the total ~2 km of exhumation in the central Colorado Plateau, demanding an unknown driver for significant erosion in the Pliocene.

58 GEOSCIENCES↗