Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “net present value”

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 343 records · Page 19

Seasonal-to-Interannual Variability in Antarctic Sea-Ice Dynamics, and Its Impact on Surface Fluxes and Water Mass Production

Strong seasonal and interannual signals in Antarctic bottom-water outflow remain unexplained yet are highly correlated with anomalies in net sea-ice growth in coastal polynyas. The mechanisms responsible for driving salination and replenishment and rejuvenation of the dense shelf "source" waters likely also generate pulses of bottom water outflow. The objective of this research is to investigate time-scales of variability in the dynamics of sea-ice in the Southern Ocean in order to determine the primary sites for production of dense shelf waters. We are using a merged satellite/buoy sea-ice motion data set for the period 1978-present day to compute the dynamics of opening and closing of coastal polynyas over the continental shelf. The Ocean Circulation and Climate Advanced Model (OCCAM) ocean general circulation model with coupled sea-ice dynamics is presently forced using National Center for Environmental Prediction (NCEP) data to simulate fluxes and the salination impact of the ocean shelf regions. This work is relevant in the context of measuring the influence of polar sea-ice dynamics upon polar ocean characteristics, and thereby upon global thermohaline ocean circulation. Interannual variability in simulated net freezing rate in the Southern Weddell Sea is shown for the period 1986-1993. There is a pronounced maximum of ice production in 1988 and minimum in 1991 in response to anomalies in equatorward meridional wind velocity. This follows a similar approximate 8-year interannual cycle in Sea Surface Temperature (SST) and satellite-derived ice-edge anomalies reported elsewhere as the "Antarctic Circumpolar Wave." The amplitude of interannual fluctuations in annual net ice production are about 40% of the mean value, implying significant interannual variance in brine rejection and upper ocean heat loss. Southward anomalies in wind stress induce negative anomalies in open water production, which are observed in passive microwave satellite images. Thus, cycles of enhanced poleward wind stress reduce ice growth by compacting the ice along the coastline and closing open water in leads and polynyas. Model simulations confirm that years of low ice production, such as 1991, coincide with years of lower than normal bottom water outflow. Future plans include the assimilation of satellite ice concentrations and ice drift dynamics to more accurately constrain boundary conditions in the model.

Drinkwater, Mark R.↗

Optimizing Direct Air Capture Solvents to Minimize Energy Consumption of CO 2 Release in a Carbonate Electrolyzer

Addressing climate change by carbon management is critical to achieving the goal of net zero carbon emissions by 2050. In this work, we examined the electrochemically-driven recovery of CO 2 during alkaline solvent regeneration for solvent-based direct air capture. A mathematical model was developed by incorporating carbonate chemistry with water electrolysis to predict the energy consumption per unit of CO 2 released. The predicted results were consistent with the experimental data, in which the experimental work was achieved by characterizing alkalinity and carbon loading values of solvent collected from a flow carbonate electrolyzer. Through this study, we learned that minimizing the energy expended on CO 2 release can be achieved by using an anolyte with a lower alkalinity, increasing the electric charge input to the electrolyzer, and reducing the ohmic resistance of the electrolyzer. Furthermore, using a supporting electrolyte, e.g., Na 2 SO 4 in the present work, effectively compensates for the higher ohmic resistance from using an anolyte with a lower alkalinity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

New Approaches to Low-Cost Scalable Doping of Interdigitated back Contact Silicon Solar Cells (Final Report)

The goal of this project was to develop novel approaches to patterning of dopants in the rear fingers of interdigitated back contact (IBC) Si solar cells to reduce the cost of manufacturing this high-efficiency-potential cell architecture. The work throughout this project can be divided into four categories: (a) development of dopant patterning technique using laser scribed Si contacts masks that are mechanically aligned with ~10 μm resolution to the underlying Si substrate; (b) measuring dopant spreading profiles in the isolation region between n- and p-type dopant fingers during plasma-enhanced chemical vapor deposition (PECVD) of doped hydrogenated amorphous silicon (a-Si:H) via shadow masks, and dopant desorption and re-adsorption during high-temperature annealing; (c) understanding the role of dopant compensation on the shunt resistance in contaminated isolation regions through analysis of defect-enhanced compensation; and (d) simulation and fabrication of passivated two-sided grid and back-contact solar cells to demonstrate the use of direct dopant patterning in cell fabrication. For the passivated two-sided grid solar cells, masked deposition was used to demonstrate an improvement in the blue response of the cell by creating a shallow front emitter. During development of the masked PECVD patterning process, we measured 3-D dopant profiles using secondary ion mass spectrometry. After deposition, in the masked region, the phosphorus dopant tail was >100 µm at concentrations >10 19 cm -3 while the boron dopant tail was shorter. During high-temperature crystallization of doped a-Si:H films to polycrystalline Si (poly-Si), phosphorus atoms spread by desorbing from the poly-Si surface and re-adsorbing onto intrinsic poly-Si on adjacent wafers that were separated by several millimeters. These contamination mechanisms resulted in a decrease in resistivity from ~10 5 Ω·cm for intrinsic poly-Si to ~10 -1 Ω·cm for contaminated poly-Si. Mitigation strategies for each contamination mechanism were developed to maintain a resistivity of ~10 5 Ω·cm between doped fingers. During fabrication of the 209 cells created during this project, it was found that despite contamination of the IBC gap through the abovementioned mechanisms, high shunt resistances and FF ~75% were still reached. Investigation into this led to the discovery of defect-enhanced compensation which exists within highly defective poly-Si when net doping concentrations reach the value of defect density (~10 18 cm -3 for many poly-Si films). Simulations guided us in the fabrication of IBCs and the cells fabricated were able to meet the year-end goals for BP 2 and 3 of 15% and 17% IBC cell efficiency, as well as the BP 4 goal of a 1% absolute increase in efficiency for PERC-like devices. However, the most efficient cell created during this project of 18.6% fell short of the 21% final project target. While the champion device fell short in V oc and J sc , many devices fabricated were able to reach the necessary goals of ~40 mA/cm 2 , ~700 mV, and ~75% FF required for a 21% device. The results generated from this project were disseminated through 11 conference presentations and proceedings. Presentations included oral talks at 2019 IEEE PVSC, 2019 MRS Fall Meeting, 2020 PVSEC-30 and 2021 IEEE PVSC, as well as poster presentations at 2020 IEEE PVSC and 2021 SiPV. The project resulted in 2 peer reviewed publications – one published in IEEE Journal of Photovoltaics and one in ACS Applied Energy Materials. The information gained in this project will aid in development of improved processes for fabrication of high-efficiency solar cells and other areas of the semiconductor device industry as well. The IBC cells will also pave the way for higher efficiency tandem devices. Through further manufacturing of high-efficiency solar cells, more of the world’s energy demands can be met through renewable sources, helping to stave off the worst effects that may come about from global climate change.

14 SOLAR ENERGY↗

An OpenStreetMaps based tool to study the energy demand and emissions impact of electrification of medium and heavy-duty freight trucks

In this paper, we present the mathematical formulation of an OpenStreetMaps (OSM) based tool that compares the costs and emissions of long-haul medium and heavy-duty (M&HD) electric and diesel freight trucks, and determines the spatial distribution of added energy demand due to M&HD EVs. The optimization utilizes a combination of information on routes from OSM, utility rate design data across the United States, and freight volume data, to determine these values. In order to deal with the computational complexity of this problem, we formulate the problem as a convex optimization problem that is scalable to a large geographic area. In our analysis, we further evaluate various scenarios of utility rate design (energy charges) and EV penetration rate across different geographic regions and their impact on the operating cost and emissions of the freight trucks. Our approach determines the net emissions reduction benefits of freight electrification by considering the primary energy source in different regions. Such analysis will provide insights to policy makers in designing utility rates for electric vehicle supply equipment (EVSE) operators depending upon the specific geographic region and to electric utilities in deciding infrastructure upgrades based on the spatial distribution of the added energy demand of M&HD EVs. To showcase the results, a case study for the U.S. state of Texas is conducted.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Cumulants and correlation functions of net-proton, proton, and antiproton multiplicity distributions in Au+Au collisions at energies available at the BNL Relativistic Heavy Ion Collider

Here, we report a systematic measurement of cumulants, C n , for net-proton, proton, and antiproton multiplicity distributions, and correlation functions, κ n , for proton and antiproton multiplicity distributions up to the fourth order in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7, 11.5, 14.5, 19.6, 27, 39, 54.4, 62.4, and 200 GeV. The C n and κ n are presented as a function of collision energy, centrality and kinematic acceptance in rapidity, y, and transverse momentum, p T . The data were taken during the first phase of the Beam Energy Scan (BES) program (2010–2017) at the BNL Relativistic Heavy Ion Collider (RHIC) facility. The measurements are carried out at midrapidity (|y| < 0.5) and transverse momentum 0.4 < p T < 2.0 GeV/c, using the STAR detector at RHIC. We observe a nonmonotonic energy dependence ($\sqrt{s_{NN}}$ = 7.7–62.4 GeV) of the net-proton C 4 /C 2 with the significance of 3.1σ for the 0–5% central Au+Au collisions. This is consistent with the expectations of critical fluctuations in a QCD-inspired model. Thermal and transport model calculations show a monotonic variation with $\sqrt{s_{NN}}$. For the multiparticle correlation functions, we observe significant negative values for a two-particle correlation function, κ 2 , of protons and antiprotons, which are mainly due to the effects of baryon number conservation. Furthermore, it is found that the four-particle correlation function, κ 4 , of protons plays a role in determining the energy dependence of proton C 4 /C 1 below 19.6 GeV, which cannot be understood by the effect of baryon number conservation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Impact of Magnetohydrodynamic (MHD) Induced Drag on High Ballistic Coefficient Aerocapture at Neptune

This paper presents analyses illustrating the benefits of magnetohydrodynamic (MHD) drag force control for a capsule performing aerocapture at Neptune relative to the same capsule utilizing conventional aerodynamic-only control across a range of design and flight conditions. MHD technology utilizes an onboard mag-net to create a magnetic field that interacts with the ionized plasma surrounding the capsule at hypersonic velocities, thus creating an electron current and inducing a Lorentz force used to increase capsule drag. The study assumes a 4.57-meterMars Science Laboratory (MSL) 70-degree sphere-cone capsule aerodynamic shape arrives hyperbolically at Neptune and aerocaptures into a synchronous Tri-ton orbit (3986 km x 430,000 km orbit inclined 157 degrees) that offers repeated flyby opportunities. Hyperbolic arrival velocity and the capsule entry mass (i.e., ballistic coefficient) are varied over a wide range of values to study the impact of these variables on key flight performance and sizing parameters. Results show reduced heat rates, heat loads, and thermal protection system mass for MHD cases at all conditions, with the benefit of MHD relative to conventional aerodynamic-only options increasing as ballistic coefficient increases, suggesting the payoff for MHD is significantly better for high-mass entry system design concepts with high ballistic coefficients.

Entry Descent and Landing↗

Plasma parameters and the reduction potential at a plasma-liquid interface

We report nonthermal plasmas in contact with liquids have been shown to generate a variety of reactive species capable of initiating reduction–oxidation (redox) reactions at the electrochemically active plasma–liquid interface. In conventional electrochemical cells, selective redox chemistry is achieved by controlling the reduction potential at the solid electrode–electrolyte interface by applying a bias via an external circuit. In the case of plasma–liquid systems, an analogous means of tuning the reduction potential near the interface has not clearly been identified. When treated as a floating surface, the liquid is expected to adopt a net negative charge to balance the flux of hot electrons and relatively cold positive ions. The reduction potential near the plasma–liquid interface is hypothesized to be proportional to the floating potential, which can be approximated using an analytical model provided the plasma parameters are known. Herein, we present a framework for correlating the electron density and electron temperature of a noble gas plasma jet to the reduction potential near the plasma–liquid interface. The plasma parameters were acquired for an argon atmospheric plasma jet in contact with an aqueous solution by means of laser Thomson scattering. The reduction potential was determined using identical reference electrodes to measure the potential difference between the plasma–liquid interface and bulk solution. Interestingly, the measured reduction potentials near the plasma–liquid interface were found to be in good agreement with the model-predicted values determined using the plasma parameters obtained from the Thomson scattering experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Algorithm for Compressing Time-Series Data

An algorithm based on Chebyshev polynomials effects lossy compression of time-series data or other one-dimensional data streams (e.g., spectral data) that are arranged in blocks for sequential transmission. The algorithm was developed for use in transmitting data from spacecraft scientific instruments to Earth stations. In spite of its lossy nature, the algorithm preserves the information needed for scientific analysis. The algorithm is computationally simple, yet compresses data streams by factors much greater than two. The algorithm is not restricted to spacecraft or scientific uses: it is applicable to time-series data in general. The algorithm can also be applied to general multidimensional data that have been converted to time-series data, a typical example being image data acquired by raster scanning. However, unlike most prior image-data-compression algorithms, this algorithm neither depends on nor exploits the two-dimensional spatial correlations that are generally present in images. In order to understand the essence of this compression algorithm, it is necessary to understand that the net effect of this algorithm and the associated decompression algorithm is to approximate the original stream of data as a sequence of finite series of Chebyshev polynomials. For the purpose of this algorithm, a block of data or interval of time for which a Chebyshev polynomial series is fitted to the original data is denoted a fitting interval. Chebyshev approximation has two properties that make it particularly effective for compressing serial data streams with minimal loss of scientific information: The errors associated with a Chebyshev approximation are nearly uniformly distributed over the fitting interval (this is known in the art as the "equal error property"); and the maximum deviations of the fitted Chebyshev polynomial from the original data have the smallest possible values (this is known in the art as the "min-max property").

Hawkins, S. Edward, III↗

An Integrated Paradigm for the Management of Delivery Risk in Electricity Markets: From Batteries to Insurance and Beyond [Slides]

In wholesale electricity markets today, flexibility from a limited number of distributed energy resources (DERs) is offered daily, and the value of flexibility is not yet recognized for economic hedging of delivery risk. Under a three-year project funded by the ARPA-E PERFORM program, a collaborative team is working towards developing an integrated risk management framework that will leverage flexibility from distributed and bulk resources to cost-effectively and reliably manage delivery risk of intermittent resources. Two concepts are at the core of the proposed integrated risk management framework: (A) flexibility options, which are a novel type of options and enable wholesale electricity market participants to hedge uncertainty by buying flexibility. (B) DER flexibility scores, which provide a way for utilities or aggregators to classify assets in groups with different likelihood of delivering contracted flexibility. This report presentation will focus on the proposed ISO-product "flexibility options," which is complementary to ramp and other products being introduced by ISOs/RTOs to manage net load uncertainties. Participating resources with imbalance risk can buy flexibility options to hedge their production, whereas grid-connected resources that can provide physical flexibility can offer flexibility options. We will present basics of the formulation for a day-ahead ISO market that matches buyers and sellers of this hedge in coordination with existing capabilities to schedule energy and ancillary services, and outline how their settlements mitigate the impact of imbalance risk.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Mapping the dark matter in the NGC 5044 group with ROSAT: Evidence for a nearly homogeneous cooling flow with a cooling wake

The NGC 5044 group of galaxies was observed by the ROSAT Position Sensitive Proportional Counter (PSPC) for 30 ks during its reduced pointed phase (1991 July). Due to the relatively cool gas temperature in the group (kT = 0.98 +/- 0.02 keV) and the excellent photon statistics (65,000 net counts), we are able to determine precisely a number of fundamental properties of the group within 250 kpc of the central galaxy. In particular, we present model-independent measurements of the total gravitating mass, the temperature and abundance profiles of the gas, and the mass accretion rate. Between 60 and 250 kpc, the gas is nearly isothermal with T varies as r(exp (-0.13 +/- 0.03)). The total gravitating mass of the group can be unambiguously determined from the observed density and temperature profiles of the gas using the equation of hydrostatic equilibrium. Within 250 kpc, the gravitating mass is 1.6 x 10(exp 13) solar mass, yielding a mass-to-light ratio of 130 solar mass/solar luminosity. The baryons (gas and stars) comprise 12% of the total mass within this radius. At small radii, the temperature clearly increases outward and attains a maximum value at 60 kpc. The positive temperature gradient in the center of the group confirms the existence of a cooling flow. The cooling flow region extends well beyond the temperature maximum with a cooling radius between 100 and 150 kpc. There are two distinct regions in the cooling flow separated by the temperature maximum. In the outer region, the gas is nearly isothermal with a unifor m Fe abundance of approximately 80% solar, the flow is nearly homogeneous with dot-M= 20 to 25 solar mass/year, the X-ray contours are spherically symmetric, and rho(sub gas) varies as r(exp -1.6). In the inner region, the temperature profile has a positive gradient, the mass accretion rate decreases rapidly inward, the gas density profile is steeper, and the X-ray image shows some substrucutre. NGC 5044 is offset from the centroid of the outer X-ray contours indicating that the central galaxy may have a residual velocity with respect to the center of the group potential. There is also a linear X-ray feature with an extent of approximately 30 kpc with one end coincident with NGC 5044. The X-ray emission from this feature is softer than the ambient gas. We interpret this feature as a 'cooling wake' formed by the accreting gas as it is gravitationally focused into the wake of NGC 5044. One of the most surprising results of our PSPC observation is the discovery of a nearly homogeneous cooling flow. Prior results concerning the mass accretion profile in cooling flows indicate that dot-M varies as r. This relation implies that significant mass deposition occurs at large radii which generates an inhomogeneous flow. The mass accretion rate in the NGC 5044 group is essentially a constant beyond 40 kpc (well within the cooling radius). Significant mass deposition (a declining dot-M) does not commence until the gas accretes to within 40 kpc of the group center where the radiative cooling time is approximately equals 10(exp 9) year. Th is radius also corresponds to the temperature maximum, the break in gas density profile, and the onset of structure in the X-ray image. A Hubble constant of H(sub 0) = 50 km/sec/Mpc is used throughout the paper.

David, Laurence P.↗

Building Energy Modeling Enhancements to Identify Least-Cost Pathways to Net-Zero Carbon Homes: Preprint

Residential grid-interactive efficient buildings (GEBs) can utilize high levels of energy efficiency and demand flexibility to deliver value to occupants, the grid, and society. However, without the ability to analyze, design, and optimize residential GEBs there is a risk that homes will not realize their full potential value as efficient and dynamic resources, which could ultimately lead to higher than necessary energy costs for U.S. households and lower realization of potential energy efficiency and demand flexibility benefits such as increased convenience/automation, thermal comfort, durability, indoor air quality, and resilience. In 2018, the National Renewable Energy Laboratory (NREL) Residential Modeling Team developed a vision for an open source Residential GEB Analytics Platform built within DOE's EnergyPlus and OpenStudio modeling environment that would enable the design and optimization of residential GEBs, including the identification of least-cost pathways to highly energy-efficient and energy-flexible homes (e.g., net zero carbon homes). We created a detailed workplan for the development of new and enhanced residential GEB component-level models for EnergyPlus and OpenStudio needed to progress toward the vision for the analytics platform while delivering near-term benefits to industry. This paper 1) presents the vision for the platform, summarizing the workplan for residential GEB modeling enhancements; 2) highlights significant advancements that have been achieved between 2018 and 2022, including stochastic residential occupancy modeling, flexible water heater modeling, detailed lithium-ion stationary battery modeling, and realistic residential HVAC modeling; and 3) outlines ongoing efforts and next steps toward the full vision for the platform.

building energy modeling↗

Permafrost Active Layer Microbes From Ny Ålesund, Svalbard (79°N) Show Autotrophic and Heterotrophic Metabolisms With Diverse Carbon-Degrading Enzymes

The active layer of permafrost in Ny Ålesund, Svalbard (79°N) around the Bayelva River in the Leirhaugen glacier moraine is measured as a small net carbon sink at the brink of becoming a carbon source. In many permafrost-dominating ecosystems, microbes in the active layers have been shown to drive organic matter degradation and greenhouse gas production, creating positive feedback on climate change. However, the microbial metabolisms linking the environmental geochemical processes and the populations that perform them have not been fully characterized. In this paper, we present geochemical, enzymatic, and isotopic data paired with 10 Pseudomonas sp. cultures and metagenomic libraries of two active layer soil cores (BPF1 and BPF2) from Ny Ålesund, Svalbard, (79°N). Relative to BPF1, BPF2 had statistically higher C/N ratios (15 ± 1 for BPF1 vs. 29 ± 10 for BPF2; n = 30, p < 10 –5 ), statistically lower organic carbon (2% ± 0.6% for BPF1 vs. 1.6% ± 0.4% for BPF2, p < 0.02), statistically lower nitrogen (0.1% ± 0.03% for BPF1 vs. 0.07% ± 0.02% for BPF2, p < 10 –6 ). The d 13 C values for inorganic carbon did not correlate with those of organic carbon in BPF2, suggesting lower heterotrophic respiration. An increase in the δ 13 C of inorganic carbon with depth either reflects an autotrophic signal or mixing between a heterotrophic source at the surface and a lithotrophic source at depth. Potential enzyme activity of xylosidase and N-acetyl-β-D-glucosaminidase increases twofold at 15°C, relative to 25°C, indicating cold adaptation in the cultures and bulk soil. Potential enzyme activity of leucine aminopeptidase across soils and cultures was two orders of magnitude higher than other tested enzymes, implying that organisms use leucine as a nitrogen and carbon source in this nutrient-limited environment. Besides demonstrating large variability in carbon compositions of permafrost active layer soils only ∼84 m apart, results suggest that the Svalbard active layer microbes are often limited by organic carbon or nitrogen availability and have adaptations to the current environment, and metabolic flexibility to adapt to the warming climate.

54 ENVIRONMENTAL SCIENCES↗

TPSAS-NF1676L-26925-DND

Observationally-based atmospheric energy balance is analyzed using Clouds and the Earth's Radiant Energy System (CERES)-derived TOA and surface irradiance, Global Precipitation Climatology Project (GPCP)-derived precipitation, dry static and kinetic energy tendency and divergence estimated from ERA-Interim, and surface sensible heat flux from SeaFlux. The residual tends to be negative over tropics and positive over mid-latitudes. A negative residual implies that precipitation rate is too small, divergence is too large, or radiative cooling is too large. The residual of atmospheric energy is spatially and temporally correlated with cloud objects to identify cloud types associated with the residual. Spatially, shallow cumulus, cirrostratus, and deep convective cloud object occurrence are positively correlated with the absolute value of the residual. The temporal correlation coefficient between the number of deep convective cloud objects and individual energy components, net atmospheric irradiance, precipitation rate, and the sum of dry static and kinetic energy divergence and their tendency over western Pacific are, respectively, 0.84, 0.95, and 0.93. However, when all energy components are added, the atmospheric energy residual over tropical Pacific temporally correlates well with the number of shallow cumulus cloud objects over tropical Pacific. Because shallow cumulus does not alter enough atmospheric energy compared to the residual, these suggest 1) if retrieval errors associated with deep convective clouds are causing the column integrated atmospheric energy residual, the errors vary among individual deep convective clouds, 2) it is possible that the residual is associated with processes in which shallow cumulus clouds affect deep convective clouds and hence atmospheric energy budget over tropical western Pacific, and 3) it is also possible that cumulus congestus present cirrostratus that missed in massive sensor retrieval contribute the residual. Cloud objects used in the study that leads to the above conclusions are derived from Moderate Resolution Imaging Spectroradiometer (MODIS). In this investigation, we use Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) and CloudSat data to evaluate and verify above findings. Specifically, we focus on 3) above and analyze frequency of occurrence of cumulus congestus present under cirrostratus and correlation to the energy balance residual.

Seiji Kato↗

Thermodynamics and Cloud Radiative Effect from the First Year of GoAmazon

Deforestation is an ongoing concern for the Amazon Rainforest of Brazil and associated changes to the land surface have been hypothesized to alter the climate in the region. A comprehensive set of meteorological observations at the surface and within the lower troposphere above Manacapuru, Brazil and data from the Modern Era Retrospective Analysis for Research and Applications Version 2 (MERRA-2) are used to evaluate the seasonal cycle of cloudiness, thermodynamics, and the radiation budget. While ample moisture is present in the Amazon Rainforest year round, the northward progression of the Hadley circulation during the dry season contributes to a drying of the middle troposphere and inhibits the formation of deep convection. This results in a reduction in cloudiness and precipitation as well as an increase in the height of the lifting condensation level, which is shown to have a negative correlation to the fraction of low clouds. Frequent cloudiness prevents solar radiation from reaching the surface and clouds are often reflective with high values of shortwave cloud radiative effect at the surface and top of the atmosphere. Cloud radiative effect is reduced during the dry season however the dry season surface shortwave cloud radiative effect is still double what is observed during the wet season in other tropical locations. Within the column, the impact of clouds on the radiation budget is more prevalent in the longwave part of the spectrum, with a net warming in the wet season.

Amazon Rainforest↗

Modelled land use and land cover change emissions – a spatio-temporal comparison of different approaches

Quantifying the net carbon flux from land use and land cover changes (fLULCC) is critical for understanding the global carbon cycle and, hence, to support climate change mitigation. However, large-scale fLULCC is not directly measurable and has to be inferred from models instead, such as semi-empirical bookkeeping models and process-based dynamic global vegetation models (DGVMs). By definition, fLULCC estimates are not directly comparable between these two different model types. As an important example, DGVM-based fLULCC in the annual global carbon budgets is estimated under transient environmental forcing and includes the so-called loss of additional sink capacity (LASC). The LASC results from the impact of environmental changes on land carbon storage potential of managed land compared to potential vegetation and accumulates over time, which is not captured in bookkeeping models. The fLULCC from transient DGVM simulations, thus, strongly depends on the timing of land use and land cover changes mainly because LASC accumulation is cut off at the end of the simulated period. To estimate the LASC, the fLULCC from pre-industrial DGVM simulations, which is independent of changing environmental conditions, can be used. Additionally, DGVMs using constant present-day environmental forcing enable an approximation of bookkeeping estimates. Here, we analyse these three DGVM-derived fLULCC estimations (under transient, pre-industrial, and present-day forcing) for 12 models within 18 regions and quantify their differences as well as climate- and CO2-induced components and compare them to bookkeeping estimates. Averaged across the models, we find a global fLULCC (under transient conditions) of 2.0±0.6 PgC/yr for 2009–2018, of which ∼40 % are attributable to the LASC (0.8±0.3 PgC/yr). From 1850 onward, the fLULCC accumulated to 189±56 PgC with 40±15 PgC from the LASC. Around 1960, the accumulating nature of the LASC causes global transient fLULCC estimates to exceed estimates under present-day conditions, despite generally increased carbon stocks in the latter. Regional hotspots of high cumulative and annual LASC values are found in the USA, China, Brazil, equatorial Africa, and Southeast Asia, mainly due to deforestation for cropland. Distinct negative LASC estimates in Europe (early reforestation) and from 2000 onward in the Ukraine (recultivation of post-Soviet abandoned agricultural land), indicate that fLULCC estimates in these regions are lower in transient DGVM compared to bookkeeping approaches. Our study unravels the strong dependence of fLULCC estimates on the time a certain land use and land cover change event happened to occur and on the chosen time period for the forcing of environmental conditions in the underlying simulations. We argue for an approach that provides an accounting of the fLULCC that is more robust against these choices, for example by estimating a mean DGVM ensemble fLULCC and LASC for a defined reference period and homogeneous environmental changes (CO2 only).

land use and land cover change↗

The Value of Information From Horizontal Distributed Acoustic Sensing Compared to Multicomponent Geophones Via Machine Learning

Abstract Faults play an important role in recharging many geothermal reservoirs, and seismic information can image the locations of these faults. The value of information (VOI) metric is used to objectively quantify and compare the value of two types of seismic receiver data via a machine learning approach. The demonstrated VOI methodology is novel by including spatial models from seismic data and obtaining the information statistics from machine learning. Our two-dimensional numerical experiments compare images created from sparsely spaced (80 m), two-component geophone sampling to high spatial resolution (1 m), single-component DAS. We used a three-fold cross validation of a U-Net convolutional neural networks to achieve average classification statistics. The results suggest that when horizontal sources are utilized, geophones and DAS identify reflectors and non-reflectors at roughly the same rate. The average F1 score for horizontal DAS is 0.939 and 0.931 for geophones. For images created from a vertical source, DAS performed marginally better (F1 = 0.919) than geophones (F1 = 0.877). Our transferrable methodology can provide guidance on which acquisition scenarios can improve images of important structures in the subsurface and present an efficient method for obtaining reliability statistics from high-dimensional, spatial data.

Energy & Fuels↗

Coupled primary production and respiration in a large river contrasts with smaller rivers and streams

Abstract Although time series in ecosystem metabolism are well characterized in small and medium rivers, patterns in the world's largest rivers are almost unknown. Large rivers present technical difficulties, including depth measurements, gas exchange (, ) estimates, and the presence of large dams, which can supersaturate gases. We estimated reach‐scale metabolism for the Hanford Reach of the Columbia River (Washington state, USA), a free‐flowing stretch with an average discharge of 3173 . We calculated from semi‐empirical models and directly estimated it from tracer measurements. We fixed at the median value from these calculations (0.5 ), and used maximum likelihood to estimate reach‐scale, open‐channel metabolism. Both gross primary production (GPP) and ecosystem respiration (ER) were high (GPP range: 0.3–30.8 g , ER range: 0.8–30.6 g ), with peak GPP and ER occurring in the late summer or early fall. GPP increased exponentially with temperature, consistent with metabolic theory, while light was seasonally saturating. Annual average GPP, estimated at 1500 g carbon , was in the top 2% of estimates for other rivers. GPP and ER were tightly coupled and 90% of GPP was immediately respired, resulting in net ecosystem production near 0. Patterns in the Hanford Reach contrast with those in small‐medium rivers, suggesting that metabolism magnitudes and patterns in large rivers may not be simply scaled from knowledge of smaller rivers.

54 ENVIRONMENTAL SCIENCES↗

Dust Trails of SP/Tuttle and the Unusual Outbursts of the Ursid Shower

Halley-type comets tend to have a series of dust trails that remain spatially correlated for extended periods of time, each dating from a specific return of the comet. Encounters with 1 - 9 revolution old individual dust trails of 55P/Tempel-Tuttle have led to well recognized Leonid shower maxim, the peak time of which was well predicted by recent models. Now. we used the same model to calculate the position of dust trails of comet Shuttle, a Halley-type comet in an (approximately) 13.6 year orbit passing just outside of Earth's orbit. We discovered that the meteoroids tend to be trapped in the 14:12 mean motion resonance with Jupiter, while the comet librates in a slightly shorter period orbit around the 13:15 resonance. It takes six centuries to change the orbit enough to intersect Earth's orbit. During that time, the meteoroids and comet separate in mean anomaly by six years. thus explaining the unusual aphelion occurrences of Ursid outbursts. The resonances also prevent dispersion, so that the dust trail encounters (specifically, from dust trails of AD 1378 - 1405) occur only in one year in each orbit. We predicted enhanced activity on December 22, 2000, at around 7:29 and 8:35 UT (universal time) from dust trails dating to the 1405 and 1392 return, respectively. This event was observed from California using video and photographic techniques. At the same time, five Global-MS-Net stations in Finland, Japan and Belgium counted meteors using forward meteor scatter. The outburst peaked at 8:06:07 UT, December 22, at Zenith Hourly Rate (approx.) 90 per hour. The Ursid rates were above half peak intensity during 4.2 hours. This is only the second Halley type comet for which a meteor outburst can be dated to a specific return of the parent comet, and traces their presence back form 9 to at least 45 revolutions of the comet. New orbital elements of Ursid meteoroids are presented. We find that most orbits do scatter around the anticipated positions, confirming the link with comet Shuttle and the epoch of ejection. The 1405 and.1392 dust trails appear to have contributed similar amounts to the activity profile. Some orbits provide a hint of much older debris being present as well. Some of the dispersion in the radiant position may reflect a true variation in inclinations, with two groupings at low and high values, which is not understood at present.

Jenniskens, Peter↗