Search NASA⌕ Search

SEARCH · Search NASA

Results for “surface exchange”

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 19 records

Influence of Heterointerfaces on the Kinetics of Oxygen Surface Exchange on Epitaxial La 1.85 Sr 0.15 CuO 4 Thin Films

Considerable attention has been directed to understanding the influence of heterointerfaces between Ruddlesden–Popper (RP) phases and ABO 3 perovskites on the kinetics of oxygen electrocatalysis at elevated temperatures. Here, we report the effect of heterointerfaces on the oxygen surface exchange kinetics by employing heteroepitaxial oxide thin films formed by decorating LaNiO 3 (LNO) on La 1.85 Sr 0.15 CuO 4 (LSCO) thin films. Regardless of LNO decoration, tensile in-plane strain on LSCO films does not change. The oxygen surface exchange coefficients (kchem) of LSCO films extracted from electrical conductivity relaxation curves significantly increase with partial decorations of LNO, whereas full LNO coverage leads to the reduction in the kchem of LSCO films. The activation energy for oxygen exchange in LSCO films significantly decreases with partial LNO decorations in contrast with the full coverage of LNO. Optical spectroscopy reveals the increased oxygen vacancies in the partially covered LSCO films relative to the undecorated LSCO film. We attribute the enhanced oxygen surface exchange kinetics of LSCO to the increased oxygen vacancies by creating the heterointerface between LSCO and LNO.

25 ENERGY STORAGE↗

Task Sharing of Proton Incorporation in Vertically Aligned Nanocomposite Triple Conductors: Growth, Structure, and Surface Exchange Kinetics

As protonic ceramic electrolysis cells emerge for efficient H 2 production, there is a need to develop air electrode materials enabling fast, durable steam splitting and proton incorporation. Single-phase triple conductors may fail to satisfy the myriad performance/stability requirements, and their critical charge-carriers (holes, oxygen vacancies, and protons) are in competition, limiting their concentrations. Instead, we propose task-sharing, vertically aligned nanocomposites (VANs), comprising a proton conductor (BaZr 0.9 Y 0.1 O 3-δ ) and a redox-active mixed ionic electronic conductor (Ce 0.9 Pr 0.1 O 2-δ ), that may enable rapid proton surface exchange at the solid–gas interface and transport along the solid–solid heterointerfaces. We grew VANs by pulsed laser deposition and investigated the interplay between their processing conditions, structure, and proton and oxygen surface exchange kinetics. We varied the substrate temperature, laser repetition rate, laser fluence, and processing oxygen pressure. The crystallinity and phases were characterized by grazing-incidence X-ray diffraction, and the strain and structural order as a function of depth were evaluated by angle-dependent synchrotron X-ray pair distribution function analysis. To evaluate the potential for interdiffusion, the formation energies of substitutional defects were simulated with density functional theory. Corresponding structural analysis and elemental mapping were performed by scanning/transmission electron microscopy, energy-dispersive X-ray spectroscopy, and electron energy-loss spectroscopy, indicating distinct nanoscale compositional regions with a hierarchical structure embedded in individual VANs columns and minimal interdiffusion across a bilayer film. Proton and oxygen surface exchange coefficients (k H , k O ) and polarization resistances were evaluated by electrical and optical relaxations and impedance spectroscopy of VAN-incorporated protonic ceramic electrochemical cells, respectively, at 400–500 °C, demonstrating values comparable to some of the best-known triple and mixed conductors.

36 MATERIALS SCIENCE↗

Proton surface exchange kinetics of perovskite triple conducting thin films for protonic ceramic electrolysis cells: BaPr 0.9 Y 0.1 O 3–δ (BPY) vs. Ba 1–x Co 0.4 Fe 0.4 Zr 0.1 Y 0.1 O 3–δ (BCFZY)

Protonic ceramic electrolysis cells (PCECs) are an attractive green H 2 production technology, given their intermediate-temperature operating range and ability to produce dry H 2 . However, PCECs will benefit from development of more efficient and durable “triple conducting” anodes where steam is split, H incorporated, and oxygen evolved. In this work, we evaluated the kinetics of the steam-splitting/H incorporation reaction on BaPr 0.9 Y 0.1 O 3–δ (BPY) in comparison to the benchmark Ba 1–x Co 0.4 Fe 0.4 Z r0.1 Y 0.1 O 3–δ (BCFZY) composition, replacing most of the transition metal elements (Co, Fe, Zr) with the lanthanide Pr. We prepared geometrically well-defined perovskite BPY and BCFZY thin films by pulsed laser deposition and performed simultaneous optical transmission relaxation and electrical conductivity relaxation measurements at 400–500 °C in 0.21 atm O 2 during switching of the steam partial pressure to isolate and compare their proton surface exchange coefficients (k). The k values of BPY were comparable to those of BCFZY and more stable over time. According to angle-resolved XPS and STEM-EDS mapping of FIB cross-sections, the surface of BPY exhibited Ba enrichment, Pr deficiency, and Si contamination. In contrast, BCFZY exhibited Ba deficiency throughout, no obvious surface segregation, and less Si contamination. The Ba segregation on the BPY film appears to have promoted steam splitting/H incorporation kinetics even though the more basic surface reacted with the acidic environmental SiO x H y . Faster kinetics observed on stoichiometric BCFZY vs. Ba-deficient BCFZY confirmed the benefit of a high A-site Ba concentration. This result contrasts with most work on perovskites applied in solid oxide electrolysis cell anodes, in which A-site segregation is considered deleterious for surface reaction kinetics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Power Plant Testing of Ultrasonic Measurements of Temperature Distributions and Heat Fluxes to Heat Exchange Surfaces

Extreme environmental conditions are common in energy conversion applications, in which the most hardened insertion sensors often do not perform reliably for long. On the other hand, ultrasonic measurements can be acquired noninvasively, with sensitive components kept away from the damaging environment. We have previously developed the ultrasonic method for measuring the spatial distribution of temperatures in solid materials applicable when large thermal gradients are present. In the developed approach, we use the echogenically segmented ultrasound propagation path, structured to contain engineered or naturally occurring echogenic features, to produce a train of echoes in response to an external pulse of ultrasonic excitation. This paper outlines the testing results of this approach applied to the temperature distribution measurements to a waterwall of a 500 MW power-plant utility boiler. The validation results obtained over the prolonged power plant test show that the estimated temperature profile is correctly captured, and the measurement accuracy can be comparable with traditional insertion sensors, such as thermocouples. Overall, the testing has confirmed that the developed approach has matured to become an attractive alternative to conventional sensing in solving challenging problems of long-term temperature measurements in extreme environments. Heat fluxes and thermal stresses in the structure can then be characterized noninvasively using the measured temperature distribution as the basis.

Walton, Kenneth↗

Riverbed Temperature and 4D ERT Monitoring Reveals Heterogenous Horizontal and Vertical Groundwater-Surface Water Exchange Flows Under Dynamic Stage Conditions

Groundwater surface water exchange plays a critical role in physical, biological, and geochemical function of coastal and riverine systems. Observing exchange flow behavior in heterogeneous systems is a primary challenge, particularly when flows are governed by dynamic river stage or tidal variations. In this paper we demonstrate a novel application of time-lapse 3D electrical resistivity tomography and temperature monitoring where an array of thermistors installed beneath a riverbed double as resistivity electrodes. We use the array to monitor stage driven exchange flows over a 6-day period in a dynamic, stage-driven high order stream. We present a method for addressing the otherwise confounding effects of the moving river-surface boundary on the raw resistivity data, thereby enabling successful tomographic imaging. Temperature time-series at each thermistor location and time-lapse 3D images of changes in bulk electrical conductivity together provide a detailed description of exchange dynamics over a 10-meter by 45-meter section of the riverbed, to a depth of approximately 5 m. Results reveal highly variable flux behavior throughout the monitoring domain including both horizontal and vertical exchange flows.

54 ENVIRONMENTAL SCIENCES↗

Surface oxygen exchange kinetics of mixed conducting oxides: Dilatometric vs electrical conductivity relaxation study

The surface oxygen exchange kinetics of mixed ionic-electronic conducting (MIEC) oxides play a crucial role in a variety of applications, including solid oxide fuel/electrolysis cells (SOFCs/SOECs), permeation membranes and sensors. To date, a variety of methods, including isotope exchange, electrical conductivity, optical absorptivity and thermogravimetry relaxation, and impedance spectroscopy have been used to measure the oxygen exchange coefficient of MIEC oxides. Each of these methods have their advantages and limitations, depending on the physical and chemical properties of the investigated materials and their sample dimensions and densities. Here, we demonstrate the ability to precisely measure the surface oxygen exchange coefficient (k chem ) of MIEC model material Pr 0.1 Ce 0.9 O 2-δ (PCO) by use of dilatometric relaxation measurements, particularly of interest for materials that exhibit larger chemical expansion coefficients. This is achieved by use of porous bulk specimens to ensure that the contribution of oxygen exchange to the overall kinetics is dominant. We demonstrate that k chem values extracted from chemical expansion relaxation measurements on PCO are nearly identical to those derived from electrical conductivity relaxation measurements. As a result, this provides the opportunity to precisely investigate the oxygen exchange and chemical expansion kinetics of a wide range of materials used in high-temperature applications, particularly where more conventional methods are difficult or inappropriate to apply.

(Pr,Ce)O2-x↗

Effect of H + Exchange and Surface Impurities on Bulk and Interfacial Electrochemistry of Garnet Solid Electrolytes

Contact loss and current constriction pose significant challenges at the Li metal interface of solid-state batteries. For garnet-structured Li 7 La 3 Zr 2 O 12 (LLZO), these effects are amplified by Li + /H + exchange and surface contamination reactions, which lead to conductivity losses and poor Li wetting. In this study, we utilize a variety of surface treatment processes across 37 cells to selectively induce proton exchange and contamination reactions in LLZO. The resulting bulk and surface chemistry is systematically characterized and correlated to changes in electrochemical properties. Additionally, we combine impedance analysis and finite element method modeling to deconvolute sources of impedance contributions at the Li metal interface. Specifically, we show that constriction impedance at the Li metal interface arises not solely from voids, but also from ionically-resistive surface contaminants. Further, these findings emphasize the connection between ionic conductivity and constriction, demonstrating that micron-scale ionically-resistive components increase constriction even with identical contact geometries. Finally, we leverage our comprehensive dataset to highlight unstable overpotential growth as a failure mechanism, additionally showing that the phase of a cell’s impedance is a sensitive indicator for the onset of interfacial instability. Overall, this study clarifies the impacts of proton exchange and surface contamination on electrochemical properties at the Li|solid electrolyte interface and elucidates insights that are generalizable to other solid-state battery systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-Temperature 2D Optical Relaxation Visualizes Enhanced Oxygen Exchange Kinetics at Metal-Mixed Conducting Oxide Interfaces

Solid-state heterointerfaces are of interest for emergent local behavior that is distinct from either bulk parent compound. One technologically relevant example is the case of mixed ionic/electronic conductor (MIEC)–metal interfaces, which play an important role in electrochemistry. Metal–MIEC composite electrodes can demonstrate improved catalytic activity vs single-phase MIECs, improving fuel cell efficiency. Similarly, MIEC surface reaction kinetics are often evaluated using techniques that place metal current collectors in contact with the surface under evaluation, potentially altering the response vs the native surface. Here, techniques enabling direct and local in situ observation of the behavior at and around such heterointerfaces are needed. Here, we develop a spatially resolved optical transmission relaxation (2D-OTR) method providing continuous evaluation of local, high-temperature, controlled atmosphere defect kinetics across a ~1 cm2 sample area simultaneously in a contact-free manner. We apply it to observe the spatial variance of oxygen incorporation and evolution rates at ~525–620 °C, in response to step changes in oxygen partial pressure, on MIEC SrTi 0.65 Fe 0.35 O 3–x films as a function of distance from porous Pt and Au layers. Using this model geometry, we find significant enhancements in kinetics adjacent to the metals that decay over a few millimeter distance. To extract kinetic parameters, we fit the short-term optical data (initial portion of relaxations) with an exponential decay function appropriate for surface-exchange-limited kinetics, yielding apparent surface exchange coefficients (k chem ) with spatial resolution, decreasing with distance from the metal. To understand the kinetic processes governing the complete (long-term) optical relaxations, we performed COMSOL simulations, which demonstrated that a combination of laterally varying k chem and in-plane diffusion controls the observed kinetics over the full time range. Further support for spatially varying kchem comes from demonstrations of changing surface and bulk chemistry vs distance from the metal–MIEC interface, by X-ray photoelectron and optical absorption spectroscopies, respectively. Although microporous Pt and Au are not excellent electrodes in isolation, both metals exert a synergistic effect on the oxygen surface exchange rate in the presence of the mixed conducting film.

36 MATERIALS SCIENCE↗

Interactions Between Clouds and Wind-Driven Surface Heat Exchanges over Land

Earth system model experiments show that increasing horizontal resolution fundamentally alters the simulated soil-moisture-precipitation feedback. Kilometer-scale simulations often produce weaker or even negative feedback compared to coarse-resolution models. A key difference of kilometer-scale models is that they resolve mesoscale secondary circulations, including boundary layer horizontal rolls and cellular structures, in addition to cold pools and downdrafts associated with convective precipitation. However, because the relevant processes occur on yet-smaller scales, these circulations are often poorly resolved. This project demonstrated that boundary layer secondary circulations significantly affect surface heat exchanges and wind gusts, and that current model parameterizations can misrepresent these processes at kilometer-scale resolution. Using DOE Atmospheric Radiation Measurement (ARM) observations and targeted experiments with the DOE Energy Exascale Earth System Model (E3SM), we identified physically unrealistic wind gust and surface flux responses to secondary circulations, diagnosed a systematic overestimation of wind shear in convective cold pools, and uncovered a multivariate relationship between land surface fluxes and the scales of updrafts that form shallow cumulus clouds. These findings provide observation-based recommendations for improving parameterizations of surface fluxes and wind gusts in high-resolution Earth system models, thereby reducing uncertainty in convective storm prediction and land-atmosphere feedbacks.

54 ENVIRONMENTAL SCIENCES↗

Surface Complexation/Ion Exchange Hybrid Model for Radionuclide Sorption to Clay Minerals (M4SF-23LL010301062)

This progress report (Level 4 Milestone Number M4SF-23LL010301062) summarizes research conducted at Lawrence Livermore National Laboratory (LLNL) within the Argillite International Collaborations Activity Number SF-23LL01030106. The activity is focused on our long-term commitment to engaging our partners in international nuclear waste repository research. The focus of this milestone is the establishment of international collaborations for surface complexation modeling and the associated impacts of unlocking larger, community-based datasets. More specifically, we are developing a database framework for Spent Fuel and Waste and Science Technology (SFWST) that is aligned with the Helmholtz Zentrum Dresden Rossendorf (HZDR) sorption database development group in support of the database needs of the SFWST program. In our FY22 effort, we described a detailed analysis of U(VI) sorption to quartz through both traditional surface complexation modeling and through a hybrid ML framework. In FY23, effort was placed on publication of these results and expansion of the LLNL surface complexation and ion exchange database (L-SCIE) in order to assess mineral-based radionuclide retardation under a wider variety of geochemical conditions (e.g., ionic strength, varying electrolyte compositions). Efforts were initiated to expand L-SCIE to include radionuclide surface complexation and ion exchange to clays that are relevant to subsurface geochemical processes occurring at nuclear waste repositories. In particular, a large source of sorption data for clays resides at the Paul Scherrer Institute (PSI) (work primarily by Bradbury and Baeyens) and we initiated discussions on how to retrieve those data and apply FAIR principles to those datasets. In addition to L-SCIE development, two hybrid models that incorporate AI/ML were investigated and compared to discern the most promising approaches for accurate and precise estimations of radionuclide retardation. Key considerations for future model development include (1) the ability to reduce computational burden on determining retardation coefficients for PA and (2) the ability to quantify and predict radionuclide-mineral partitioning at a more efficient, rapid pace due to automated workflows. Upon the careful consideration of the most effective modeling approaches, we are identifying ways to implement these approaches into PA. Ultimately, the data science-based workflows will provide a major incentive for other institutions to adopt a FAIR-formatted, interoperable database. LLNL will play a key role in disseminating sorption data and acting as good data stewards by updating the database in a consistent format and assessing the quality of the newly assimilated data in an organized fashion. To this end, all data and workflows are open access and made available on the LLNL Seaborg research website (https://seaborg.llnl.gov/resources/geochemical-databases-modeling-codes).

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Unlocking bulk and surface oxygen transport properties of mixed oxide-ion and electron conducting membranes with combined oxygen permeation cell and oxygen probe method

Surface exchange kinetics and bulk diffusion of oxygen are of paramount importance to the activity of oxygen electrocatalysis and performance of electrochemical devices such as fuel cell, metal-air batteries, and oxygen separation membranes. Conventional approaches to obtaining these transport properties are often limited to single property under a specific non-operation related condition. Here we use a combined oxygen permeation cell and oxygen probe methodology to simultaneously attain rates of oxygen surface exchange and bulk conductivity/chemical diffusivity of three representative mixed oxide-ion and electron conductors, namely SrCo 0.9 Ta 0.1 O 3-δ (SCT), La 0.6 Sr 0.4 CoO 3-δ (LSC) and La 0.6 Sr 0.4 FeO 3-δ (LSF), operated under a steady-state oxygen flux. The results explicitly show that SCT exhibit the highest oxide-ion conductivity/chemical diffusivity, fastest rates of surface oxygen exchange kinetics, thus promising to be the best oxygen electrocatalyst. Here, we have also mapped out the distribution of oxygen chemical potential gradient across the membranes and applied B-transport number concept to illustrate the rate-limiting steps in the overall oxygen permeation process.

42 ENGINEERING↗

Optimization of Triply Periodic Minimal Surface Heat Exchanger to Achieve Compactness, High Efficiency, and Low-Pressure Drop

With advancements in additive manufacturing (AM) techniques, high-quality triply periodic minimal surface (TPMS) structures can now be produced. TPMS walled heat exchangers (HX) hold significant potential for industrial applications and are receiving increasing attention. This paper explores the impact of various TPMS design variables on flow and thermal performance to optimize TPMS heat exchangers for compactness, high efficiency, and low pressure drop. The design variables examined include the type of TPMS lattice, unit cell size, wall thickness, aspect ratio, TPMS orientation, and equivalent thickness. The study reveals that the flow and heat transfer performance of TPMS structures are significantly affected by these design variables. For the Gyroid, Diamond, and SplitP lattices, performance is nearly identical when the surface-to-volume ratio is kept constant. The average velocity of the fluid in the TPMS HX should be 0.3 m/s. The corresponding Re is between 300~800. Thin wall thickness, small equivalent thickness, and flat lattice configurations can significantly reduce pressure drop while maintaining the overall heat transfer coefficient. Additionally, the angle between the flow direction and TPMS orientation can increase pressure drop. Three aluminum heat exchangers were successfully printed using an AM machine, and testing results are comparable with theoretical prediction.

42 ENGINEERING↗

Influence of Concentration Gradients on Electroconvection at a Cation-Exchange Membrane Surface

Membrane-based systems, such as electrodialysis, play an important role in desalination and industrial separation processes. Electrodialysis uses alternating anion- and cation-exchange membranes with a perpendicular electric field to generate concentrated and diluate streams from a feed solution. It is known that under overlimiting current conditions, reduced charge and mass transfer at the membrane interface leads to regions of high ion depletion generating instability and vortices termed electroconvection. While electroconvective mixing is known to directly impact the separation efficiency of electrodialysis, the influence of ion concentration gradients across the membrane experienced in a functional electrodialysis system is not known. Here, we report the influence of ion concentration gradients across a cation exchange membrane (Nafion) that is both aligned with and opposed to the applied electric field. Experiments were conducted by coflowing NaCl solutions of different concentrations (0.1–100 mM) on each side of the membrane, and electroconvection was visualized with a fluorescence dye (Rhodamine 6G). Here, we obtained concentration profiles from fluorescence image data and systematically measured the thickness of the depletion boundary layer d BL under different conditions. We found smaller d BL values at a higher flow rate both with and without concentration gradients. Our results show that electroconvection is enhanced when the electric field is opposite to the direction of the concentration gradient.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Adsorption of REEs to Kaolinite via Ion Exchange and Surface Complexation as a Function of Water Chemistry

Rare earth elements (REEs) are critical components of modern technology behind renewable energy, transportation, and electronics but have a limited current supply. A substantial portion of global REE production relies on ion adsorption deposits. A high abundance of kaolinite in REE enrichment zones within these deposits suggests that kaolinite controls the subsurface migration of REEs. This study aimed to improve the current understanding of REE binding to kaolinite under varying water chemistry conditions. We conducted batch experiments with kaolinite (KGa-2) and three REEs (Nd, Dy, and Yb) at varying pH, electrolyte concentration, dissolved inorganic carbon (DIC), low molecular weight organic acids (citric and oxalic acids), and total REE concentration conditions. Increasing electrolyte concentration inhibits REE adsorption at pH < 7, suggesting that ion exchange contributes to adsorption at these pH values. DIC affects adsorption above pH 7–8 by forming strong aqueous complexes with heavy REEs. Citric acid decreases REE adsorption via aqueous complexation of REEs at pH > 5 but does not affect adsorption at pH < 5. The surface complexation model captures the main adsorption trends with two mechanisms: ion exchange on basal planes at pH < ∼6 and inner-sphere surface complexation to edge sites at pH > ∼6. Equilibrium constants for surface complexation increase in the order of Yb > Dy > Nd, indicating a higher strength of adsorption for heavy REEs. This study demonstrates how water chemistry conditions control the adsorption mechanisms that may determine the mobility of REEs in subsurface environments rich in kaolinite.

58 GEOSCIENCES↗

Airborne flux measurements of ammonia over the southern Great Plains using chemical ionization mass spectrometry

Abstract. Ammonia (NH3) is an abundant trace gas in the atmosphere and an important player in atmospheric chemistry, aerosol formation and the atmosphere–surface exchange of nitrogen. The accurate determination of NH3 emission rates remains a challenge, partly due to the propensity of NH3 to interact with instrument surfaces, leading to high detection limits and slow response times. In this paper, we present a new method for quantifying ambient NH3, using chemical ionization mass spectrometry (CIMS) with deuterated benzene cations as reagents. The setup aimed at limiting sample–surface interactions and achieved a 1σ precision of 10–20 pptv and an immediate 1/e response rate of < 0.4 s, which compares favorably to the existing state of the art. The sensitivity exhibited an inverse humidity dependence, in particular in relatively dry conditions. Background of up to 10 % of the total signal required consideration as well, as it responded on the order of a few minutes. To showcase the method's capabilities, we quantified NH3 mixing ratios from measurements obtained during deployment on a Gulfstream I aircraft during the HI-SCALE (Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems) field campaign in rural Oklahoma during May 2016. Typical mixing ratios were 1–10 parts per billion by volume (ppbv) for the boundary layer and 0.1–1 ppbv in the lower free troposphere. Sharp plumes of up to tens of ppbv of NH3 were encountered as well. We identified two of their sources as a large fertilizer plant and a cattle farm, and our mixing ratio measurements yielded upper bounds of 350 ± 50 and 0.6 kg NH3 h−1 for their respective momentary source rates. The fast response of the CIMS also allowed us to derive vertical NH3 fluxes within the turbulent boundary layer via eddy covariance, for which we chiefly used the continuous wavelet transform technique. As expected for a region dominated by agriculture, we observed predominantly upward fluxes, implying net NH3 emissions from the surface. The corresponding analysis focused on the most suitable flight, which contained two straight-and-level legs at ∼ 300 m above ground. We derived NH3 fluxes between 1 and 11 mol km−2 h−1 for these legs, at an effective spatial resolution of 1–2 km. The analysis demonstrated how flux measurements benefit from suitably arranged flight tracks with sufficiently long straight-and-level legs, and it explores the detrimental effect of measurement discontinuities. Following flux footprint estimations, comparison to the NH3 area emissions inventory provided by the U.S. Environmental Protection Agency indicated overall agreement but also the absence of some sources, for instance the identified cattle farm. Our study concludes that high-precision CIMS measurements are a powerful tool for in situ measurements of ambient NH3 mixing ratios, and even allow for the airborne mapping of the air–surface exchange of NH3.

58 GEOSCIENCES↗