Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion Couple Experiments”

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 109 records · Page 6

Resolving ECRH deposition broadening due to edge turbulence in DIII-D

Microwave heat pulse propagation experiments have demonstrated a correlation between millimeter-scale turbulence and deposition profile broadening of electron cyclotron (EC) waves on the DIII-D tokamak. In a set of discharges in DIII-D, a variation in edge density fluctuations on the mm-scale is associated with 40%–150% broader deposition profiles, expressed in terms of normalized minor radius, as compared with equilibrium ray tracing. The 1D power profile is determined from transport analysis of the electron temperature response to EC power modulation using perturbative analysis with a square wave power modulation at 20–70 Hz, producing a series of Fourier harmonics that are fit collectively to resolve transport. Fitting an integrated heat flux expressed in the Fourier basis of the modulation to diffusive, convective, and coupled transport terms in a linear model can resolve the broadened EC deposition width from the power perturbation to resolve a broadening in each case. The best fit degree of beam broadening observed scales approximately linearly with the Doppler backscattering measured fluctuation level in the steep gradient region. Quantifying the effect of edge fluctuation broadening on EC current drive power needs of future devices will require 3D full-wave codes that can be validated on the current generation of machines. These DIII-D experiments provide a quantitative measure of fluctuation effects and a dataset to benchmark full-wave simulations that can model and eventually predict nonlinear effects neglected by 1D equilibrium beam and ray tracing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Polarizabilities of neutral atoms and atomic ions with a noble gas electron configuration

Atomic polarizabilities play an important role in the development of force fields for molecular simulations, as well as for the development of qualitative concepts of atomic and molecular behavior. Coupled cluster theory at the coupled cluster singles doubles triples level with very large correlation-consistent basis sets with extended diffuse functions has been used to predict the polarizabilities of the atomic neutrals, mono-cations and mono-anions with a noble gas configuration. Additional corrections for scalar relativistic and spin–orbit effects were also included for the electron configurations of Kr, Xe, and Rn. The results are in excellent agreement with experiment or with other high level calculations where available. The current results for most of these species represent the best available values for the polarizabilities. The results show that the polarizability of H - is very difficult to calculate without extremely diffuse functions. The polarizability of H - is the largest value, 34.05 Å 3 , calculated for all species in the current study. The polarizabilities of the remaining halogen anions are also the best available values. The polarizabilities of the halogen anions (excluding F - ) and H - have a linear correlation with the electron affinity of the neutral atom. Spin– orbit effects, even for closed shell species, cannot be ignored for quantitative accuracy, and the inclusion of spin–orbit effects for Fr + , Rn, and At - increases the polarizability by 4%, 6%, and 15%, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Combined Experimental and Modeling Study of the Interactions of Acid Gas with Common Spacecraft Surfaces for Fire Safety Applications

A fire in a spacecraft poses detrimental consequences and risks mission success in addition to crew safety. This is compounded during long-duration missions when the crew has limited options to recover from a fire. A common spacecraft fire concern is the smoldering of wire insulation, typically made from Polyvinyl chloride (PVC) or Polytetrafluoroethylene (PTFE). This creates acid gases such as Hydrogen Chloride (HCl), Hydrogen Fluoride (HF) and Hydrogen Cyanide (HCN). These poisonous gases are hazardous to the crew. They also interact with common surfaces within the spacecraft more than dominant combustion products such as CO2 and H2O. This makes them more difficult to track for potential fire detection techniques, or for postfire clean-up. It is imperative to be able to understand and predict the fate of these poisonous species in a microgravity environment in order to design a safe vehicle. HCl interacts with a number of materials inside a spacecraft. Primary among these materials is aluminum, which is abundantly used due to its strong and light weight nature. Aluminum has a natural oxide layer that protects it from corrosion but is typically treated to enhance this oxide layer. Among these treatments is a chromate conversion coating (CCC), which provides a thin enough protective oxide layer to still conduct electricity, and a traditional anodized material that has a thicker oxide layer that does not conduct electricity. Nomex is another common material found inside a spacecraft. It is a flame-resistant woven polymer that is related to nylon. This commercially available material is used for cargo storage bags and as a fire barrier. Physics-based models were developed to predict the uptake of HCl by these materials. The ultimate objective of these models is to predict the fate of HCl within the spacecraft so that sensors can be placed in meaningful locations in future missions based on the model predictions. To support these modeling efforts, experiments were performed in a cast acrylic test cell that measured the difference between the inlet and outlet concentration of HCl after inserting a sample rod of the test material. Different uptake capacities were realized for each type of sample tested. A computational fluid dynamics model (CFD) model of the reactor was then constructed that used a one-step global reaction rate with calibratable reaction (or kinetic) constants. These constants were calibrated to match the HCl uptake on the CCC aluminum samples, and the same kinetic constants were then tested for the stock and anodized aluminum samples. Model predictions matched the experimental data for the stock aluminum, and to a much lesser extent, the anodized aluminum. The model was additionally validated at different flow rates, sample surface areas, and inlet concentrations, and showed good agreement for all stock and CCC samples. The model did not accurately predict the HCl uptake in the anodized samples compared to the other two types of aluminum. Adjusting the kinetic constants and transport properties did little to improve the prediction. X-Ray Photoelectron Spectroscopy (XPS) was used to determine that the oxide layer thickness of anodized aluminum is approximately 5,000 nm, compared to 250 nm for CCC and 50 nm for stock. XPS also revealed presence of chlorine further down in the aluminum oxide layer in anodized samples than CCC and stock samples after the samples were saturated with HCl, indicating that accounting for diffusion of HCl into the oxide layer is important for accurate prediction of HCl uptake onto anodized aluminum. Consequently, a multi-scale model was developed and tested. First, a single pore inside the anodized aluminum oxide layer was modeled and is referred to as the pore-scale model. In this model, HCl diffused through the pore and reacted with the aluminum oxide pore wall to create aluminum chloride. The sample was then saturated when the mass transfer resistance through the growing aluminum chloride layer became too large for the HCl to reach the aluminum oxide wall and continue the reaction. This pore-scale model was coupled to the reactor-scale model using a concentration-dependent diffusion coefficient, resulting in much more accurate predictions (approximately half the sum square error of the aforementioned reactor-scale model that produced good agreement for stock and CCC) for a variety of operating conditions. The amount of water vapor or relative humidity (RH) in the flow during a reactor experiment was determined to influence HCl uptake. Experiments were performed to understand the interaction of gaseous HCl with aluminum surfaces in the presence of water vapor. The results show that increasing levels of RH increased the capacity of aluminum to adsorb HCl but decreased the capacity of Nomex to uptake HCl. A series of tests were performed on individual aluminum samples after they had been saturated with a fixed concentration of HCl in dry air conditions with the goal of determining how their HCl uptake capacity changes after various treatments with water relative to the original saturation tests. HCl-saturated aluminum samples subjected to a second dry air flow at the same HCl concentration as the original test had an uptake of 23.5% of the original sample with no treatment in between. Saturated aluminum samples subjected to an air flow with a RH of 90% in between tests had an uptake of 35.6% of the original. Saturated aluminum samples submerged in distilled water for 12 hours in between tests had an uptake of 82.2% of the original sample. Previously saturated aluminum tested with HCl and a 50% RH air flow resulted in similar uptake characteristics in multiple repeated tests. The experiments show the profound effect water vapor has on HCl uptake onto aluminum surfaces. In the samples subjected to water vapor or liquid water, capillary condensation and capillary diffusion alters the transport of HCl significantly. A model was proposed that developed a relationship between RH and the coefficient of HCl diffusion in aluminum chloride. This produced an “S-shaped” curve with diffusion coefficient as a function of RH, with 45% RH represented as the point where the diffusion coefficient is halfway between no water saturation and 100% water saturation in the aluminum chloride product layer. No difference in uptake characteristics for the experiment or model were realized between 50% and 62% RH. The results from the large-scale microgravity experiment, Saffire, are discussed as they pertain to the fate of HCl throughout a spacecraft. HCl was released, both as a standalone event, and in concurrence with the burning of a structured cloth. These events only produced a small response in the far field HCl sensor, while a PMMA burn that did not produce HCl had a significantly greater response. A ground-based large-scale facility was constructed to flow acid gas at the scale and configuration realized in the Saffire experiments. A CFD model of this duct was constructed to test kinetic parameters developed in this work at a larger scale and different geometric configuration and to predict the results of the large-scale facility. The models developed in this work were used to interpret the results of the microgravity tests and lead the discussion on what further experiments and models are needed in order to predict the fate of acid gas in a spacecraft environment. To summarize, the major contributions of this work are as follows: the capacity to uptake HCl, with and without the presence of water vapor, was measured for a variety of real spacecraft surfaces. Several different models (single reactor-scale, multiscale, spacecraft-scale) were developed and with the aid of modeling, the rate of uptake for those surfaces was also predicted and validated. The kinetic parameters determined from the small-scale reactor experiments and models were used to predict large-scale and microgravity tests. Conclusions from this research will be used in the design of spacecraft vehicles and large-scale microgravity fire safety experiments. The models built by this work will aid designers in sensor placement and could be used to predict acid gas transport from fires in partial gravity, as would be seen in Lunar and Martian habitats.

fire safety↗

Effect of penetrant–polymer interactions and shape on the motion of molecular penetrants in dense polymer networks

The diffusion of dilute molecular penetrants within polymers plays a crucial role in the advancement of material engineering for applications such as coatings and membrane separations. The potential of highly cross-linked polymer networks in these applications stems from their capacity to adjust the size and shape selectivity through subtle changes in network structures. In this paper, we use molecular dynamics simulation to understand the role of penetrant shape (aspect ratios) and its interaction with polymer networks on its diffusivity. We characterize both local penetrant hopping and the long-time diffusive motion for penetrants and consider different aspect ratios and penetrant–network interaction strengths at a variety of cross-link densities and temperatures. The shape affects the coupling of penetrant motion to the cross-link density- and temperature-dependent structural relaxation of networks and also affects the way a penetrant experiences the confinement from the network meshes. The attractive interaction between the penetrant and network primarily affects the former since only the system of dilute limit is of present interest. These results offer fundamental insights into the intricate interplay between penetrant characteristics and polymer network properties and also suggest future directions for manipulating polymer design to enhance the separation efficiency.

Chemistry↗

Rain Reevaporation, Boundary Layer Convection Interactions, and Pacific Rainfall Patterns in an AGCM

Sensitivity experiments with an atmospheric general circulation model (AGCM) show that parameterized rain re-evaporation has a large impact on simulated precipitation patterns in the tropical Pacific, especially on the configuration of the model s intertropical convergence zone (ITCZ). Weak re-evaporation leads t o the formation of a "double ITCZ" during the northern warm season. The double ITCZ is accompanied by strong coupling between precipitation and high-frequency vertical motion in the planetary boundary layer (PBL). Strong reevaporation leads to a better overall agreement of simulated precipitation with observations. The model s double ITCZ bias is reduced. At the same time, correlation between high-frequency vertical motion in the PBL and precipitation is reduced. Experiments with modified physics suggest that evaporative cooling by rain near the PBL top weakens the coupling between precipitation and vertical motion. This may reduce the model s tendency to form double ITCZs. The strength of high-frequency vertical motions in the PBL was also reduced directly through the introduction of a diffusive cumulus momentum transport (DCMT) parameterization. The DCMT had a visible impact on simulated precipitation in the tropics, but did not reduce the model s double bias in all cases.

Bacmeister, Julio T.↗

Understanding grain boundary segregation under irradiation: A mesoscale modeling study

Grain boundaries (GBs) are an important component in the design of irradiation resistant materials as they play a crucial role as point defect sinks for non-equilibrium vacancies and self-interstitials. However, in the presence of sustained point defect fluxes, the coupled transport of atomic components with inherently different diffusivities often leads to radiation-induced segregation (RIS), with undesired effects on properties such as intergranular corrosion and stress-assisted cracking. In this talk, we investigate the fundamental mechanisms leading to GB segregation in BCC FeCr-based alloys using multiscale modeling and experiments. A phase-field model is employed to account for both thermodynamic and kinetic segregation mechanisms. The full Onsager transport matrix is derived using kinetic Monte Carlo, molecular dynamics, and kinetic cluster expansion (KineCluE) calculations. The thermodynamics of GB segregation is described via an atomic density-based modification to the CALPHAD free energy. The modeling predictions are then compared against segregation data from atom probe characterization for various temperatures, compositions, and neutron irradiation conditions. We will conclude by discussing the utility of the multiscale modeling approach to shed light on the fundamental mechanisms leading to solute-GB interaction.

36 MATERIALS SCIENCE↗

Dark Photons from Perturbative Decay of a Misaligned Higgs Field

We reconsider the production of dark photons $A'$ as dark matter, from the perturbative decay of a dark Higgs field $h$, that is stochastically misaligned from the minimum of its potential during inflation. This is a simple and predictive framework for generating the $A'$ relic abundance. It is constrained by structure formation, since the $A'$ are initially boosted, and inflationary isocurvature fluctuations, which require small quartic couplings $λh^4$. We identify $A'$ masses between 100 eV and 1 GeV and gauge couplings $g\sim 10^{-15}-10^{-10}$ that are consistent in this scenario, and which become more tightly constrained if a generic level of kinetic mixing is present. The favored parameter region could be tested through future CMB or Lyman-$α$ observations, and, in the presence of kinetic mixing, by direct detection experiments or diffuse soft gamma-ray searches.

Cline, James M. [McGill U., Montreal (main)]↗

Multiphonon diffuse scattering in solids from first principles: Application to layered crystals and two-dimensional materials

Time-resolved diffuse scattering experiments have gained increasing attention due to their potential to reveal nonequilibrium dynamics of crystal lattice vibrations with full momentum resolution. Although progress has been made in interpreting experimental data on the basis of one-phonon scattering, understanding the role of individual phonons can be sometimes hindered by multiphonon excitations. In Ref. [M. Zacharias, H. Seiler, F. Caruso, D. Zahn, F. Giustino, P. C. Kelires, and R. Ernstorfer, Phys. Rev. Lett. 127, 207401 (2021)], we have introduced a rigorous approach for the calculation of the all-phonon inelastic scattering intensity of solids from first-principles. In the present work, we describe our implementation in detail and show that multiphonon interactions are captured efficiently by exploiting translational and time-reversal symmetries of the crystal. We demonstrate its predictive power by calculating the scattering patterns of monolayer molybdenum disulfide (MoS 2 ), bulk MoS 2 , and black phosphorus (bP), and we obtain excellent agreement with our measurements of thermal electron diffuse scattering. Remarkably, our results show that multiphonon excitations dominate in bP across multiple Brillouin zones, while in MoS 2 they play a less pronounced role. We expand our analysis for each system and examine the effect of individual atomic and interatomic vibrational motion on the diffuse scattering signals. We further demonstrate the high-throughput capability of our approach by reporting all-phonon scattering maps of two-dimensional MoSe 2 , WSe 2 , WS 2 , graphene, and CdI 2 , rationalizing in each case the effect of multiphonon processes. As a side point, we show that the special displacement method reproduces the thermally distorted configuration that generates precisely the all-phonon diffuse pattern. In conclusion, the present methodology opens the way for systematic calculations of the scattering intensity in crystals and the accurate interpretation of static and time-resolved inelastic scattering experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A simple model for the formation of 'reflected', 'intermediate', and 'diffuse' ion distributions upstream of earth's bow shock

Detailed energetic ion and low-frequency wave experiments on ISEE 1, 2, and 3, have established correlations between the types of ion populations and waves observed to exist in the earth's foreshock region. The present calculation has the objective to quantitatively test the picture obtained for the earth's foreshock region. An analysis is conducted of the coupled spatial and temporal evolution of upstream protons and hydromagnetic waves by means of a simple model in which the protons are represented by counterstreaming beams which resonate with hydromagnetic waves, exciting the outward propagating modes and damping the inward propagating modes at a single wave number.

Lee, M. A.↗

Loss and Isotopic Fractionation of Alkali Elements during Diffusion-Limited Evaporation from Molten Silicate: Theory and Experiments

Moderately volatile elements (MVEs) are variably depleted in planetary bodies, reflecting the imprints of nebular and planetary processes. Among MVEs, Na, K, and Rb are excellent tracers for unraveling the history of MVE depletion in planetary bodies because they have similar geochemical behaviors but can be chemically fractionated by evaporation and condensation processes. Furthermore, K and Rb are amenable to high-precision isotopic analyses, which can help constrain the conditions of evaporation and condensation. To quantitatively understand why Na, K, and Rb are depleted in planetary bodies, we have carried out vacuum evaporation experiments from basaltic melt at 1200 and 1400 °C to study their evaporation kinetics and isotopic fractionations. We chose this composition because it is relevant to evaporation from small differentiated planetesimals. The Rb isotopic compositions of the evaporation residues were measured by multicollector inductively coupled plasma mass spectrometry (MC-ICPMS), and the K isotopic compositions were measured along profiles across the residues by secondary ion mass spectrometry (SIMS). In the 1400 °C run products, we found that the concentrations of both K and Rb in the run products decreased from core to rim, which was accompanied by a heavy K isotope enrichment near the surface. This indicates that, in this run, evaporation was limited by diffusion. To use those data quantitatively, we derive analytical equations that describe the evaporation rate and isotopic fractionation associated with diffusion-limited evaporation from a sphere, slab, and cylinder in transient and quasi-steady state regimes. This model is used to tease out the roles that diffusive transport in the melt and evaporation at the melt/gas interface play in setting the elemental depletion and isotopic composition of the residue. Under our experimental conditions, volatility decreases in the order of Na, Rb, and K. Using our experimental results in a thermodynamic model, we have estimated the product γΓ of activity coefficients × evaporation coefficients of Na, Rb, and K. The measured isotopic compositions of the residues are well explained using Rayleigh distillations, whereby the relative volatilities of K and Rb isotopes are given by the square root of their masses. We use our results and previously published data to predict how K and Rb could have been lost as a function of temperature, melt composition, oxygen fugacity, and saturation degree relevant to Vesta’s building blocks. We find that the K and Rb depletions, K/Rb elemental fractionation, and δ 41 K and δ 87 Rb isotopic fractionations of Vesta (as sampled by howardite-eucrite-diogenite (HED) meteorites) are best explained by evaporation of submillimeter size objects for 0.1-10 years at moderate temperatures (~1050 °C) in a medium ~98.8% saturated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interface Pattern Selection in Directional Solidification

The central focus of this research is to establish key scientific concepts that govern the selection of cellular and dendritic patterns during the directional solidification of alloys. Ground-based studies have established that the conditions under which cellular and dendritic microstructures form are precisely where convection effects are dominant in bulk samples. Thus, experimental data can not be obtained terrestrially under pure diffusive regime. Furthermore, reliable theoretical models are not yet possible which can quantitatively incorporate fluid flow in the pattern selection criterion. Consequently, microgravity experiments on cellular and dendritic growth are designed to obtain benchmark data under diffusive growth conditions that can be quantitatively analyzed and compared with the rigorous theoretical model to establish the fundamental principles that govern the selection of specific microstructure and its length scales. In the cellular structure, different cells in an array are strongly coupled so that the cellular pattern evolution is controlled by complex interactions between thermal diffusion, solute diffusion and interface effects. These interactions give infinity of solutions, and the system selects only a narrow band of solutions. The aim of this investigation is to obtain benchmark data and develop a rigorous theoretical model that will allow us to quantitatively establish the physics of this selection process.

Trivedi, Rohit↗

Modeling of the amplified optical stochastic cooling experiment at IOTA

Optical Stochastic Cooling (OSC), a beam cooling technique based on Stochastic Cooling, is in the early stages of experimental development. It uses radiation produced by the beam in an undulator magnet (the pickup) to correct the momentum deviation of particles downstream in another undulator (the kicker). The OSC mechanism was recently demonstrated at Fermilab’s IOTA ring using a passive configuration. However, the cooling rate of OSC can be dramatically increased by first amplifying the undulator radiation before applying the corrective kick. In collaboration with the IOTA experiment, we developed a computational model of the OSC mechanism. This paper presents beam-dynamics simulations of the amplified-OSC configuration. We implement a model of intrabeam scattering and study the effects on beam equilibrium and diffusion rate as a function of bunch charge. Finally, we investigate the phase-space dynamics with various coupling configurations between the transverse and longitudinal planes.

43 PARTICLE ACCELERATORS↗

Data and Modeling of Dendrites Subject to A Step Change in Pressure (TDSE)

There is considerable interest in dendritic solidification because of the influence dendrites have in the determination of microstructure, and thereby in the physical properties of cast metals and alloys. Current theories and models of dendritic growth generally couple diffusion effects in the melt with the physics of the interface. Data and subsequent analysis prior of the tip growth speed and radii of thermal succinonitrile dendrites in the near-convection free, on-orbit, free-fall environment demonstrate that these theories yield predictions that are reasonably in agreement with the results of experiment. However, data and analysis for assessing the interfacial physics component of theory are not sufficiently detailed or definitive. To study fundamental aspects of dendritic interface stability, we are measuring and modeling the kinetics and morphology of dendrites as they evolve from one well-defined steady state at a pre-set supercooling, through a transient stage, to a different well-defined steady state. More specifically, we subject succinonitrile dendrites, growing under steady-state conditions, to a rapid change in pressure. This leads to a rapid change in thermal driving force from the corresponding change in both the equilibrium melting temperature due to the Clapeyron effect, and a change in the far-field temperature due to adiabatic temperature changes in the bulk liquid and solid. Subsequently, we observe transformations from a well-characterized initial state into a new steady-state. Initial data reveal that the dendrite tip velocity changes almost as fast as the pressure charges, while the tip radius changes occur more slowly, taking from 10 60 seconds depending on the size of the step change and the final supercooling. Computer modeling of this process shows both agreements and disagreements with the experimental data. In making these observations and measurements, we are gaining new understandings of interfacial dynamics and state-selection physics.

Koss, Matthew B.↗

Transport and Growth Kinetics in Microgravity Protein Crystal Growth

The dynamic coupling between mass transport and incorporation of growth units into the surface of a crystal growing from solution in microgravity is used to derive quantitative information on the crystal growth kinetics. To this end, new procedures for experiment preparation, interferometric data processing and model fitting have been developed. The use of experimental data from the bulk diffusive maw transport together with a model for steady state stagnant crystal growth allows the detailed quantitative understanding of the kinetics of both the concentration depletion zone around the crystal and the growth of the crystal interface. The protein crystal used in the experiment is shown to be growing in the mixed kinetic regime (0.2 x 10(exp -6) centimeters per second less than beta R/D less than 0.9 x 10(exp -6) centimeters per second).

Otalora, F.↗

Analysis of the High-Pressure High-Temperature (HPHT) growth of single crystal diamond

A multi-scale, computational model is developed to describe the growth characteristics of single-crystal diamond in the High-Pressure, High-Temperature (HPHT) process. This model is the first to connect phase-change kinetics governing crystal growth to the continuum transport of carbon through the growth cell. Results show the importance of convective transport driven by buoyant flow in the metallic solvent, which increases the growth rate by nearly an order of magnitude over that obtained under diffusion alone. Parametric studies show how crystal growth may be kinetically-limited or transport-limited, depending on the value of the macroscopic kinetic coefficient. Estimating this kinetic coefficient from growth experiments yields a phase-change Damköhler number of unity, indicating a mixed regime where phase-change kinetics and transport are comparable and strongly coupled in this system. Mechanisms responsible for slowing growth as the crystal size increases are explained. Finally, supersaturation inhomogeneities along the facets of larger crystals are predicted, which may be relevant to solvent inclusion formation during growth.

36 MATERIALS SCIENCE↗

The response of the Amazon ecosystem to the photosynthetically active radiation fields: integrating impacts of biomass burning aerosol and clouds in the NASA GEOS Earth system model

The Amazon experiences fires every year, and the resulting biomass burning aerosols, together with cloud particles, influence the penetration of sunlight through the atmosphere, increasing the ratio of diffuse to direct photosynthetically active radiation (PAR) reaching the vegetation canopy and thereby potentially increasing ecosystem productivity. In this study, we use the NASA Goddard Earth Observing System (GEOS) model with coupled aerosol, cloud, radiation, and ecosystem modules to investigate the impact of Amazon biomass burning aerosols on ecosystem productivity, as well as the role of the Amazon's clouds in tempering this impact. The study focuses on a 7-year period (2010–2016) during which the Amazon experienced a variety of dynamic environments (e.g., La Niña, normal years, and El Niño). The direct radiative impact of biomass burning aerosols on ecosystem productivity – called here the aerosol diffuse radiation fertilization effect – is found to increase Amazonian gross primary production (GPP) by 2.6 % via a 3.8 % increase in diffuse PAR (DFPAR) despite a 5.4 % decrease in direct PAR (DRPAR) on multiyear average during burning seasons. On a monthly basis, this increase in GPP can be as large as 9.9 % (occurring in August 2010). Consequently, the net primary production (NPP) in the Amazon is increased by 1.5 %, or ~92 Tg C yr -1 – equivalent to ~37 % of the average carbon lost due to Amazon fires over the 7 years considered. Clouds, however, strongly regulate the effectiveness of the aerosol diffuse radiation fertilization effect. The efficiency of this fertilization effect is the highest in cloud-free conditions and linearly decreases with increasing cloud amount until the cloud fraction reaches ~0.8, at which point the aerosol-influenced light changes from being a stimulator to an inhibitor of plant growth. Nevertheless, interannual changes in the overall strength of the aerosol diffuse radiation fertilization effect are primarily controlled by the large interannual changes in biomass burning aerosols rather than by changes in cloudiness during the studied period.

54 ENVIRONMENTAL SCIENCES↗

The Response of the Amazon Ecosystem to the Photosynthetically Active Radiation Fields: Integrating Impacts of Biomass Burning Aerosol and Clouds in the NASA GEOS Earth System Model

The Amazon experiences fires every year, and the resulting biomass burning aerosols, together with cloud particles, influence the penetration of sunlight through the atmosphere, increasing the ratio of diffuse to direct photosynthetically active radiation (PAR) reaching the vegetation canopy and thereby potentially increasing ecosystem productivity. In this study, we use the NASA Goddard Earth Observing System (GEOS) model with coupled aerosol, cloud, radiation, and ecosystem modules to investigate the impact of Amazon biomass burning aerosols on ecosystem productivity, as well as the role of the Amazon’s clouds in tempering this impact. The study focuses on a seven-year period (2010-2016) during which the Amazon experienced a variety of dynamic environments (e.g., La Niña, normal years, and El Niño). The direct radiative impact of biomass burning aerosols on ecosystem productivity—called here the aerosol diffuse radiation fertilization effect —is found to increase Amazonian Gross Primary Production (GPP) by 2.6% via a 3.8% increase in diffuse PAR (DFPAR) despite a 5.4% decrease in direct PAR (DRPAR) on multiyear average during burning seasons. On a monthly basis, this increase in GPP can be as large as 9.9% (occurring in August 2010). Consequently, the net primary production (NPP) in Amazon is increased by 1.5%, or ~92 Tg C a-1– equivalent to ~37% of the average carbon lost due to Amazon fires over the seven years considered. Clouds, however, strongly regulate the effectiveness of the aerosol diffuse radiation fertilization effect. The efficiency of this fertilization effect is the highest in cloud-free conditions and linearly decreases with increasing cloud amount until the cloud fraction reaches ~0.8, at which point the aerosol-influenced light changes from being a stimulator to an inhibitor of plant growth. Nevertheless, interannual changes in the overall strength of the aerosol diffuse radiation fertilization effect are primarily controlled by the large interannual changes in biomass burning aerosols rather than by changes in cloudiness during the studied period.

biomass burning aerosols on ecosystem productivity↗

Coupling Surface Flow with High-performance Subsurface Reactive Flow and Transport Code PFLOTRAN

Water exchange between the surface and subsurface is important for both water resource management and environmental protection. In this paper, we develop coupled surface and subsurface flow simulation capability in a parallel subsurface flow and reactive transport code PFLOTRAN. We sequentially couple the diffusion wave-based surface flow with the subsurface flow governedby the Richards equation in PFLOTRAN. These two flow domains are linked with a boundary condition switching method that ensures continuity of pressure and flux at the surface-subsurface interface. We verify the coupled code against other existing hydrologic models and observation data using a number of numerical experiments. The coupled hydrological model exhibits good performance in strong parallel scaling tests. The new coupled surface and subsurface simulator significantly advance community simulation capability towards improving integrated hydrologic and biogeochemical understanding of complex systems such as watersheds and river corridors. Keywords: Surface flow, Integrated hydrological modeling, Boundary condition switching, Parallel computing

Wu, Runjian↗