Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed storage”

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 1,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

CyDER: A Cyber Physical Co-simulation Platform for Distributed Energy Resources in Smartgrids

The CyDER project aimed at developing an open-source, modular and scalable co-simulation platform for power grids with large shares of Distributed Energy Resources (DERs). The project partners are the Lawrence Berkeley National Lab (LBNL), Lawrence Livermore National Lab (LLNL), PG&E, SolarCity, and ChargePoint. The prime recipient is LBNL; SolarCity and ChargePoint were partners for the project’s first two years. Increased DER integration introduces a number of challenges in power grid operation including a more dynamic interaction between the transmission grid and distribution grids, and increased modeling complexity. Although specialized software exists to precisely model different components of the power system, it is far from trivial to integrate all various models and perform a holistic simulation. Instead of replicating all models in a common simulation program, a commonly accepted approach to tackle this model diversity is to couple third-party simulators and models through a co-simulation platform that coordinates information exchange among the various components. Following this line of research, this project’s objective was to develop a co-simulation platform based on a widely accepted industrial standard called Functional Mockup Interface (FMI). Within this process, the project developed models compliant with the FMI standard, called Functional Mockup Units (FMUs), and used them to perform various operational and planning power system analyses. Relying and building upon an industrial standard is the main differentiation of this project compared with previous or parallel efforts in the co-simulation area. Particular emphasis was put on delivering software utilities to facilitate setting up and running co-simulations by end-users. Furthermore, a strong aspect of this project is demonstrating that co-simulation techniques can be used to perform Hardware-in-the-Loop (HIL) simulations that couple software components (e.g., simulated models) with hardware components (e.g., real devices such PV systems and batteries). The long-term goal of CyDER project is to help establish FMI as a powerful standard for co-simulation and promote adoption by electric utilities and other interested stakeholders. The main accomplishments of the project include the development of several FMUs including distribution and transmission grid models, PV inverters with Volt/Var/Watt controllers, batteries, and predictive optimal controllers. Additionally, a unique software package was developed, called SimulatorToFMU, which is capable of exporting any Python-driven simulator or Python script as an FMU. This is an important contribution towards establishing FMI as one of the main co-simulation standards, because more and more third-party programs for sub-system modeling and simulation are delivered with Python APIs. The CyDER platform was used to perform PV hosting capacity analyses in real utility feeders with and without smart inverter controls, battery storage, and EV charging. Smart inverter controls include conventional Volt/Var/Watt controls for reactive power support and active power curtailment, but also predictive controls that optimize the charging and discharging profile of the battery connected on the DC side in order to minimize the customer’s economic benefit. Finally, an important result of this project is delivering an experimental setup that consists of residential-scale PV inverters with battery storage, a real-time grid simulator with an ideal voltage source as grid emulator, and micro Phasor Measurement Units (PMUs). All these components and additional software modules are coupled to one another using the FMI standard and can be co-simulated with the CyDER platform.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Tundra SGS – Critical Importance of Site-Specific Injectivity Testing for Carbon Dioxide Injection Well Simulation Calibration

Paper presented at 16th International Conference on Greenhouse Gas Control Technologies (GHGT-16), October 23–27, 2022, Lyon, France. Brine injectivity tests were conducted on both the Broom Creek (BC) and Deadwood (DW) storage formations in one of the cored stratigraphic test wells prior to finalized simulations. The injectivity tests entailed injection of brine into the formations, followed by shutting in the wells. The shut-in duration represents the fall-off period which was used to carry out pressure transient analyses of the tests. The analysis of injectivity tests indicated that core-indicated and modeled permeability distributions should be conservatively multiplied by 2.5× and 5× in the BC and DW geologic models, respectively, prior to numerical simulations of CO 2 injection using Computer Modelling Group’s (CMG’s) GEM equation-of-state reservoir simulator. The scaler permeability multipliers allowed the fluid flow simulations to fully represent well test-measured injection well rates in the respective formations.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2023

Effect of Formation Heterogeneity on CO 2 Dissolution in Subsurface Porous Media

Dissolution trapping is one of the most dominant mechanisms for the secure storage of CO 2 injected in porous subsurface formations saturated with brine. This trapping mechanism is enhanced by convective mixing, which occurs due to the gravitational instability between the different fluid layers in the aquifer. The reservoir permeability also plays a crucial role in the dissolution rate and overall fluid flow dynamics during the density-driven convection in porous media. Herein this study investigates the role of complex heterogeneity, i.e., irregular permeability distribution in CO 2 dissolution, using a novel experimental approach to create medium permeability heterogeneity in Hele–Shaw cells. Complex subsurface transport phenomena such as a preferential dissolution path, CO 2 sweep efficiency, changes in finger morphology, and CO 2 concentration distribution are visualized by creating heterogeneous media. Experimental results showed that reservoir permeability heterogeneity causes significant channeling effects and poor sweep efficiency. A scaling relationship between average finger growth rate (Gr) and permeability (k) was obtained as Gr [m s –1 ] = 266.8k [m 2 ] + 1.20 × 10 –6 . Furthermore, the mass of CO 2 dissolved is calculated using the spectrophotometric method to characterize the convective instability. The convective flux was analyzed by comparing the experimental dissolution flux with the theoretical diffusion flux, calculating a maximum Sherwood number of 6.8. The study’s findings improve the current understanding of the CO 2 convection morphology in heterogeneous media, allowing better assessment of long-term CO 2 storage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Stable all-solid-state sodium-sulfur batteries for low-temperature operation enabled by sodium alloy anode and confined sulfur cathode

All-solid-state sodium-sulfur (Na-S) batteries are promising for stationary energy storage devices because of their low operating temperatures (less than 100 °C), improved safety, and low-cost fabrication. Using Na alloy instead of Na metal as an anode in Na-S batteries can prevent dendrite growth and improve interfacial stability between the anode and solid electrolytes to achieve long-cycling stability. A high-sulfur content cathode possessing high sulfur utilization is also important to enable an energy-dense Na-S battery. In this work, we studied Na-Sb and Na- Sn alloy anodes and demonstrated the superiority of Na 3 Sb alloy undergoing a stable Na alloying/dealloying process at 0.04 mA cm -2 for over 500 hours. Combining the optimized Na3Sb alloy anode with sulfur-carbon composites prepared by the vapor deposition approach, the full cell shows a high sulfur specific capacity and improved rate performance. Moreover, the all-solid-state Na alloy-S battery can deliver a high initial discharge specific capacity of 1377 mAh g -1 and maintain good capacity retention of 70 % after 180 cycles at 60 °C. Further, post-cycle characterizations show that both the anode and cathode perform a reversible discharge/charge process after the 1st cycle, and the cathode undergoes significantly rearranged distributions of carbon and solid-state electrolytes after 180 cycles due to severe volume change induced by repeated sodiation/desodiation process.

25 ENERGY STORAGE↗
DOE OSTI2022

Dimensionality effect of conductive carbon fillers in LiNi 1/3 Mn 1/3 Co 1/3 O 2 cathode

Developing advanced electrode architectures through modifying active materials, conductive fillers, binders, and electrolytes as well as processing methods has drawn significant research interest. Due to the insufficient electrical conductivity of many active materials, adding conductive carbon fillers to composite electrodes provides the necessary electrical conductivity. The dimensionality effect among different conductive fillers has a significant impact on electrochemistry, which can be associated with morphological and chemical heterogeneities of electrodes. In this study, synchrotron X-ray mosaic nanotomography and X-ray spectroscopy nanoimaging provided direct three-dimensional (3D) visualization and quantification capabilities to investigate the dimensionality effects of Super P (SP) and single-walled carbon nanotube (SWCNT) fillers on the capacity retention of LiNi 1/3 Mn 1/3 Co 1/3 O 2 (NMC111). The results indicate that NMC/SWCNT electrodes, with a wrapping effect from the SWCNTs, exhibited more homogeneous particle size distributions, morphological changes, and chemical states than NMC/SP electrodes, without the wrapping effect. This work developed a framework of 3D quantification methods to study the capacity fading behavior associated with morphological and chemical heterogeneities and paved the way toward designing electrodes for high rate energy storage applications.

25 ENERGY STORAGE↗
DOE OSTI2021

Nanostructured Block Polymer Electrolytes: Tailoring Self-Assembly to Unlock the Potential in Lithium-Ion Batteries

Ion-containing solid block polymer (BP) electrolytes can self-assemble into microphase-separated domains to facilitate the independent optimization of ion conduction and mechanical stability; this assembly behavior has the potential to improve the functionality and safety of lithium-ion batteries over liquid electrolytes to meet future demands (e.g., large capacities and long lifetimes) in various applications. However, significant enhancements in the ionic conductivity and processability of BPs must be realized for BP-based electrolytes to become robust alternatives in commercial devices. Toward this end, the controlled modification of BP electrolytes’ intra-domain (nanometer-scale) and multi-grain (micrometer-scale) structure is one viable approach; intra-domain ion transport and segmental compatibility (related to the effective Flory–Huggins parameter, χ eff ) can be increased by tuning the ion and monomer-segment distributions, and the morphology can be selected such that the multi-grain transport is less sensitive to grain size and orientation.

25 ENERGY STORAGE↗
DOE OSTI2020

X-ray Nanocomputed Tomography in Zernike Phase Contrast for Studying 3D Morphology of Li–O 2 Battery Electrode

In this work, we investigated the lithium peroxide (Li 2 O 2 ) and pore size distribution in lithium-O 2 battery electrodes at different states of charge using transmission X-ray microscopy coupled with Zernike phase contrast to carry out nano-Computed Tomography. We report that such technique enables, at the nanoscale, to distinguish light elements such as carbon and Li 2 O 2 in Li-O 2 battery cathode electrodes. We verified by wave-propagation simulation that this approach efficiently improves the contrast of images in comparison with the pure absorption. The Li 2 O 2 distribution and thickness, interphases, and pore network are visualized and quantified, which give a valuable insight into our cathode architecture. From this 3D analysis, we highlight modifications of the air-cathode morphology and the Li 2 O 2 spatial organization as well as their potential implication in terms of carbon surface passivation and pore-clogging. After the full recharge process, this technique can also reveal the spatial distribution of the residual Li 2 O 2 and other byproducts.

25 ENERGY STORAGE↗
DOE OSTI2023

Resonance Scattering Treatment with the Windowed Multipole Formalism

A new method for directly sampling the resonance upscattering effect is presented. Alternatives have relied on inefficient rejection sampling techniques or large tabular storage of relative velocities. None of these approaches, which require pointwise energy data, are particularly well suited to the windowed multipole cross-section representation. The new method, called multipole analytic resonance scattering, overcomes these limitations by inverse transform sampling from the target relative velocity distribution where the cross section is expressed in the multipole formalism. The closed-form relative speed distribution contains a novel special function we deem the incomplete Faddeeva function, and we present the first results on its efficient numerical evaluation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Improving SOFC Anode Electrocatalytic Activity Using Nanoparticle Infiltration into MIEC Compositions

Mitigating activation polarization in the anode is one of the major challenges in intermediate-temperature operation of solid oxide fuel cells (SOFCs). Liquid phase infiltration of nanoscale electrocatalysts has been shown to result in significant reductions in activation polarization in SOFC anodes. In this study, we explore liquid-phase infiltration of nickel, gadolinium doped ceria (GDC), and Ni/GDC electrocatalysts into two different types of cermet anodes: one with a conventional Ni-YSZ composition, and the other with a Ni-MIEC cermet anode where the YSZ has been doped with 3 mol% TiO 2 to impart electronic conductivity. The principal goal of this study is to explore the role of electronic transport in the MIEC phase in effective utilization of the infiltrated nanoscale electrocatalysts. The role of temperature, infiltration cycles and the type of electrocatalysts have been experimentally studied in symmetric cells using electrochemical impedance spectroscopy (EIS). Distribution of relaxation times (DRT) modeling has been used to elucidate the contributions of various charge transfer processes.

25 ENERGY STORAGE↗
DOE OSTI2026

Machine Protection Studies of Fast Instability-Driven Beam Losses in the EIC Electron Storage Ring

Fast transverse beam instabilities represent one of the most challenging protection scenarios for the Electron-Ion Collider (EIC) electron storage ring (ESR). Coherent betatron oscillations driven by collective effects can grow rapidly and lead to substantial beam losses within only a few tens of turns, posing a threat to collimators, superconducting (SC) magnets, cryogenic systems, and detector components. A dedicated simulation framework based on distributed transverse excitations was developed to evaluate instability-driven beam losses in the ESR. A conservative instability rise time on the order of 10 turns was used to define a bounding protection scenario. Simulations show that the original two-collimator configuration is insufficient to localize losses under all operating conditions, whereas an optimized four-collimator configuration reduces loss leakage into protected regions by more than two orders of magnitude. Nevertheless, the most severe instability scenarios produce catastrophic beam impacts, with nearly half of the stored beam being intercepted by the collimators within a single turn, exceeding the preliminary estimated survivability limits of the baseline collimator design. The results demonstrate that collimation alone is insufficient to guarantee safe ESR operation and must be complemented by a fast machine-protection system. Dedicated turn-by-turn beam-position monitors and fast beam-loss monitors provide sufficient advance warning to detect instability growth and initiate a beam abort before equipment-damaging loss levels are reached. These findings establish quantitative requirements for the ESR machine-protection architecture and support the implementation of a multilayer protection strategy for both accelerator and detector systems.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Simulation-Guided Decision-Making for Enhancing Energy Resilience in a Remote Alaskan Community

The increasing frequency and severity of extreme weather events underscore the need to bolster the resilience of energy infrastructure in remote coastal communities exposed to climate hazards. In recent years, considerable effort has been made to harden the grid infrastructure of remote communities through investment in energy storage, advanced metering, renewable generation and energy-efficient loads. However, simulation-based studies are still needed to identify appropriate locations for investment and determine the adequacy of existing infrastructure in supporting new resources. Preparing the required models may often be challenging due to insufficient metering, disparate data sources and workforce limitations. This paper describes modeling efforts undertaken to represent, within a unified co-simulation platform: (a) the electric distribution network and (b) the thermal behavior of a community medical center, in the remote community of Cordova, Alaska. With the help of case-studies, it is shown how the developed simulation platform can help the local utility in making decisions regarding capital investment aimed at enhancing community resilience.

Microgrid, energy storage, resilience↗
DOE OSTI2021

Soil carbon stocks in temperate grasslands differ strongly across sites but are insensitive to decade-long fertilization

We report enhancing soil carbon (C) storage has the potential to offset human-caused increases in atmospheric CO 2 . Rising CO 2 has occurred concurrently with increasing supply rates of biologically limiting nutrients such as nitrogen (N) and phosphorus (P). However, it is unclear how increased supplies of N and P will alter soil C sequestration, particularly in grasslands, which make up nearly a third of non-agricultural land world-wide. Here, we leverage a globally distributed nutrient addition experiment (the Nutrient Network) to examine how a decade of N and P fertilization (alone and in combination) influenced soil C and N stocks at nine grassland sites spanning the continental United States. We measured changes in bulk soil C and N stocks and in three soil C fractions (light and heavy particulate organic matter, and mineral-associated organic matter fractions). Nutrient amendment had variable effects on soil C and N pools that ranged from strongly positive to strongly negative, while soil C and N pool sizes varied by more than an order of magnitude across sites. Piecewise SEM clarified that small increases in plant C inputs with fertilization did not translate to greater soil C storage. Nevertheless, peak season aboveground plant biomass (but not root biomass or production) was strongly positively related to soil C storage at seven of the nine sites, and across all nine sites, soil C covaried with moisture index and soil mineralogy, regardless of fertilization. Overall, we show that site factors such as moisture index, plant productivity, soil texture, and mineralogy were key predictors of cross-site soil C, while nutrient amendment had weaker and site-specific effects on C sequestration. This suggests that prioritizing the protection of highly productive temperate grasslands is critical for reducing future greenhouse gas losses arising from land use change.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Optimized V1G and V2G Electric Vehicle Fleet Management and Grid Transaction at Marine Corps Air Station Miramar in San Diego, CA

The overall technical goal of the project was to demonstrate an all-electric bi-directional non-tactical fleet at Marine Corps Air Station (MCAS) Miramar that was integrated and controlled with other distributed energy resources (DERs) (i.e., PV, stationary battery, and building loads) to provide resilience to critical electric loads in the event of grid outages, to minimize charging costs, and to provide economic energy resources to electricity markets. In this project, the specific, technical objectives were: 1. Demonstrate that bi-directional electric vehicles can provide critical complementary services to fixed storage batteries in microgrid applications while performing function as non-tactical vehicles. 2. Demonstrate participation of bi-directional (V2G) and unidirectional (V1G) PEVs for demand management and minimization of charging costs. 3. Demonstrate integration of multiple DERs for grid service participation. US Marine Corps Air Station (MCAS) Miramar in San Diego was the site of this electric vehicle-to-microgrid-utility grid test and demonstration project. Existing microgrid assets in this study included (1) a public works building; (2) a 30-kW rooftop photovoltaic (PV) system and (3) a separate 250 kW carport PV system. In this project, six bi-directional V2G vans were located at the MCAS Miramar’s showcase building-scale microgrid to develop and test technical capabilities that V2G can provide in microgrid applications (e.g., cost reduction and resiliency). These resources provided aggregated demand management and simulated participation in current retail DR programs. The vehicles used in this demonstration were selected because they provided functionality that MCAS Miramar needed, 15 passenger transport and facilities work cargo carrying capacity, and bi-directional charging capability that the research project required. All vehicles in this study were manufactured and distributed by VIA Motors, Inc. There were six vehicles total and each was VIA’s VTRUX eREV V2G model, a modified General Motors Chevrolet 2500 2WD van. Three of the vans were configured as passenger vans and the other three were configured as cargo vans. Each van had an on-board bi-direcrtional inverter/charger, Bel Power Solutions model 350INVCHGT150-120-240-8G nominally rated at +/-15 kW. The VIA van’s charging connector follows the J1772 charging protocol. The bi-directional EVSEs demonstrated in this study were manufactured by Coritech, Inc. Each VGI-80-AC charging station enabled enhanced V2G charging capability to a Clipper Creek CS-100 charging module. The enhanced capabilities included ethernet communication following the SEP2.0 protocol with a distributed energy resource function set and an operator screen displaying real-time SOC, voltage, and current. The VGI-80-AC charging stations are classified as level 2 with a maximum current output of 80 A or effectively 19 kW. The VIA van’s onboard charger limited the charging and discharging power to 15 kW in each direction. A control computer was installed in the EWOC and connected to an existing monitor. The V2G control communication network was a completely stand-alone closed system that did not have any connection to any other networks on the base. A cybersecure remote communication connection was created with a cellular modem, firewall hardware, and a virtual private network configuration.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Factors controlling injection-induced rupture of intersecting faults during geological sequestration of CO 2

This study addresses coupled multiphase fluid flow and geomechanics effects on potential fault activation associated with subsurface CO 2 injection around intersecting faults. An enhanced fault-representation model is used to capture geomechanical responses of two intersecting faults with finite length during CO 2 injection. The faults are embedded in a strike-slip stress regime of a caprock-reservoir-basement system with the faults represented by zero-thickness interfaces with adjacent finite-thickness damage zones. A sensitivity analysis is conducted to study the effect of fault permeability, slip-weakening behavior, well location relative to the orientation of faults, and well placement (the number and location of injection wells). Five metrics (pressure, CO 2 plume, shear state on the fault, as well as shear displacement and stress path at selected fault monitoring points) are selected to assess CO 2 migration and reactivation of intersecting faults. The results show that induced ruptures are favored by low permeability faults due to high pressure buildup and by slip-weakening behavior resulting from fault strength reduction. The location of one injection well relative to fault orientation determines the magnitude of changes in effective normal stress and shear stress, affecting the location of induced ruptures. Well placement (two injection wells used in the paper) dominates pressure diffusion around the intersection and tips of faults. This redistributes changes in effective normal stress caused by each injection well, influencing the spatial distribution of ruptures along faults. A larger injection volume induces far-field ruptures that are controlled by stress transfer within the injection layer. The findings presented here can provide valuable insights into engineering operations for a long-term, safe, and reliable geologic CO 2 storage.

Fault permeability↗
DOE OSTI2024

Nonstationarity in the global terrestrial water cycle and its interlinkages in the Anthropocene

Climate change and human activities alter the global freshwater cycle, causing nonstationary processes as its distribution shifting over time, yet a comprehensive understanding of these changes remains elusive. Here, we develop a remote sensing–informed terrestrial reanalysis and assess the nonstationarity of and interconnections among global water cycle components from 2003 to 2020. We highlight 20 hotspot regions where terrestrial water storage exhibits strong nonstationarity, impacting 35% of the global population and 45% of the area covered by irrigated agriculture. Emerging long-term trends dominate the most often (48.2%), followed by seasonal shifts (32.8%) and changes in extremes (19%). Notably, in mid-latitudes, this encompasses 34% of Asia and 27% of North America. The patterns of nonstationarity and their dominant types differ across other water cycle components, including precipitation, evapotranspiration, runoff, and gross primary production. These differences also manifest uniquely across hotspot regions, illustrating the intricate ways in which each component responds to climate change and human water management. Our findings emphasize the importance of considering nonstationarity when assessing water cycle information toward the development of strategies for sustainable water resource usage, enhancing resilience to extreme events, and effectively addressing other challenges associated with climate change.

Science & Technology - Other Topics↗
DOE OSTI2025

Resolving the Solvation Structure and Transport Properties of Aqueous Zinc Electrolytes from Salt-in-Water to Water-in-Salt Using Neural Network Potential

Zn Cl 2 solutions are promising electrolytes for aqueous zinc-ion batteries. Here, we report a joint computational and experimental study of the structural and dynamic properties of aqueous Zn Cl 2 electrolytes with concentrations ranging from salt-in-water to water-in-salt (WIS). By developing a neural network potential (NNP) model, we perform molecular dynamics (MD) simulations with accuracy but at much larger lengths and longer timescales. The NNP predicted structures are validated by the structure factors measured by X-ray total scattering experiments. The MD trajectories provide a comprehensive and quantitative picture of the Zn 2 + solvation shell structures. Additionally, we find that the O − H covalent bonds in water are strengthened with increasing salt concentration, thus expanding the electrochemical stability window of aqueous electrolytes. In terms of dynamic properties, the calculated and experimentally measured conductivities are in good agreement. Through the analysis of the calculated cation transference number, we propose a three-stage charge carrier transport mechanism with increasing concentration: independent ion transport, strongly correlated ion transport, and small positive charge carrier diffusion through negatively charged polymeric clusters. Our study provides fundamental atomic scale insights into the structure and transport properties of the Zn Cl 2 electrolyte that can aid the optimization and development of WIS electrolytes. Published by the American Physical Society 2025

25 ENERGY STORAGE↗
DOE OSTI2020

Sandia National Laboratories Annual Site Environmental Report, 2019: New Mexico

Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the United States Department of Energy (DOE) National Nuclear Security Administration. The National Nuclear Security Administration’s Sandia Field Office administers the contract and oversees contractor operations at Sandia National Laboratories, New Mexico. Activities at the site support research and development programs with a wide variety of national security missions, resulting in technologies for nonproliferation, homeland security, energy and infrastructure, and defense systems and assessments. DOE and its management and operating contractor are committed to safeguarding the environment, assessing sustainability practices, and ensuring the validity and accuracy of the monitoring data presented in this Annual Site Environmental Report. This report summarizes the environmental protection and monitoring programs in place at Sandia National Laboratories, New Mexico, during calendar year 2019. Environmental topics include air quality, ecology, environmental restoration, oil storage, site sustainability, terrestrial surveillance, waste management, water quality, and implementation of the National Environmental Policy Act. This report is prepared in accordance with and as required by DOE O 231.1B, Admin Change 1, Environment, Safety, and Health Reporting, and has been approved for public distribution.

54 ENVIRONMENTAL SCIENCES↗