Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FLUX”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19

GEM-CEDAR Challenge: Poynting Flux at DMSP and Modeled Joule Heat

Poynting flux into the ionosphere measures the electromagnetic energy coming from the magnetosphere. This energy flux can vary greatly between quiet times and geomagnetic active times. As part of the Geospace Environment Modeling-coupling energetics and dynamics of atmospheric regions modeling challenge, physics-based models of the 3-D ionosphere and ionospheric electrodynamics solvers of magnetosphere models that specify Joule heat and empirical models specifying Poynting flux were run for six geomagnetic storm events of varying intensity. We compared model results with Poynting flux values along the DMSP-15 satellite track computed from ion drift meter and magnetic field observations. Although being a different quantity, Joule heat can in practice be correlated to incoming Poynting flux because the energy is dissipated primarily in high latitudes where Poynting flux is being deposited. Within the physics-based model group, we find mixed results with some models overestimating Joule heat and some models agreeing better with observed Poynting flux rates as integrated over auroral passes. In contrast, empirical models tend to underestimate integrated Poynting flux values. Modeled Joule heat or Poynting flux patterns often resemble the observed Poynting flux patterns on a large scale, but amplitudes can differ by a factor of 2 or larger due to the highly localized nature of observed Poynting flux deposition that is not captured by the models. In addition, the positioning of modeled patterns appear to be randomly shifted against the observed Poynting flux energy input. This study is the first to compare Poynting flux and Joule heat in a large variety of models of the ionosphere.

Rastaetter, Lutz↗

Prior biosphere model impact on global terrestrial CO2 fluxes estimated from OCO-2 retrievals

This study assesses the impact of different state of the art global biospheric CO2 flux models, when applied as prior information, on inverse model “top-down” estimates of terrestrial CO2 fluxes obtained when assimilating Orbiting Carbon Observatory 2 (OCO-2) observations. This is done with a series of observing system simulation experiments (OSSEs) using synthetic CO2 column-average dry air mole fraction (XCO2) retrievals sampled at the OCO-2 satellite spatiotemporal frequency. The OSSEs utilized a 4-D variational (4D-Var) assimilation system with the GEOS-Chem global chemical transport model (CTM) to estimate CO2 net ecosystem exchange (NEE) fluxes using synthetic OCO-2 observations. The impact of biosphere models in inverse model estimates of NEE is quantified by conducting OSSEs using the NASA-CASA, CASA-GFED, SiB-4, and LPJ models as prior estimates and using NEE from the multi-model ensemble mean of the Multiscale Synthesis and Terrestrial Model Intercomparison Project as the “truth”. Results show that the assimilation of simulated XCO2 retrievals at OCO-2 observing modes over land results in posterior NEE estimates which generally reproduce “true” NEE globally and over terrestrial TransCom-3 regions that are well-sampled. However, we find larger spread among posterior NEE estimates, when using different prior NEE fluxes, in regions and seasons that have limited OCO-2 observational coverage and a large range in “bottom-up” NEE fluxes. Seasonally averaged posterior NEE estimates had standard deviations (SD) of ∼10 % to ∼50 % of the multi-model-mean NEE for different TransCom-3 land regions with significant NEE fluxes (regions/seasons with a NEE flux ≥0.5 PgC/yr). On a global average, the seasonally averaged residual impact of the prior model NEE assumption on the posterior NEE spread is ∼10 %–20 % of the posterior NEE mean. Additional OCO-2 OSSE simulations demonstrate that posterior NEE estimates are also sensitive to the assumed prior NEE flux uncertainty statistics, with spread in posterior NEE estimates similar to those when using variable prior model NEE fluxes. In fact, the sensitivity of posterior NEE estimates to prior error statistics was larger than prior flux values in some regions/times in the tropics and Southern Hemisphere where sufficient OCO-2 data were available and large differences between the prior and truth were evident. Overall, even with the availability of spatiotemporally dense OCO-2 data, noticeable residual differences (up to ∼20 %–30 % globally and 50 % regionally) in posterior NEE flux estimates remain that were caused by the choice of prior model flux values and the specification of prior flux uncertainties.

Sajeev Philip↗

Towards a More Realistic Representation of NASA CERES-derived Surface Radiative Fluxes during Polar Night: A Comparison with the MOSAiC Field Campaign

The Arctic is one of the most sensitive regions of Earth to climate change, and yet remains one of the more difficult regions to both observe and simulate. The remoteness of the Arctic from most of civilization makes large volumes of in situ observations difficult to collect, and so satellite observations are a critical tool for observing the region, such as those from Clouds and the Earth’s Radiant Energy System (CERES). But validating the satellite radiative flux estimates is difficult because of the lack of in situ measurements. A useful remedy is utilizing measurements from field campaigns, such as the year-long Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition during 2019-2020. The extensive high-quality surface radiative flux and meteorological measurements collected from MOSAiC provide a useful check on flux retrievals from CERES instruments. We compare MOSAiC and CERES surface radiative fluxes during October through March of 2019-2020. As this time period is mostly during polar night and twilight, the focus of the study is on longwave surface fluxes. To better understand the reasons for errors in CERES fluxes, we use the large set of meteorological measurements also collected by MOSAiC, including surface temperature, thermodynamic vertical profiles, surface turbulent fluxes, and cloud properties. Previous work identified a significant source of error in the CERES estimate of surface downwelling longwave flux as the estimate of low level cloud amount when compared with MOSAiC W-band radar measurements. An overestimate of low cloud amount yielded an underestimate of the downwelling flux on the order of tens of W m2. Continuing this work, we test other meteorological properties such as the vertical thermodynamic profile and surface turbulent fluxes. In addition, we examine different estimates of surface fluxes derived from CERES measurements. Previous work used the SYN1deg product, which had the limitation of having a 1º x 1 º horizontal resolution. A newer product estimates the surface fluxes at the CERES footprint resolution, which allows a more precise colocation, and thus a representative flux estimate, with the MOSAiC site.

J Brant Dodson↗

Radon Fluxes from an Earthen Barrier Over Uranium Mill Tailings After Two Decades of Service - 17234

Radon (Rn) fluxes were measured at a disposal facility for uranium mill tailings after two decades of service. The facility was closed with an earthen cover vegetated with grasses. Measurements were made on the surface of the Rn barrier and directly on the surface of the tailings. Fluxes were measured using extra small (area = 0.018 m{sup 2}), small (0.071 m{sup 2}), medium (0.59 m{sup 2}), and large (2.32 m{sup 2}) flux chambers to evaluate the impact of measurement scale on Rn flux. Activated carbon (AC) passive collectors and electric radon detectors (RAD7) were used to measure Rn concentrations. Tests were conducted at various locations on each cover representing conditions that can lead to different levels of soil structural development and different water content. Rn fluxes at the surface of the Rn barrier were much lower than Rn fluxes measured at the surface of the tailings, indicating the barrier remained effective for Rn containment. Geometric mean Rn fluxes measured at the surface of the Rn barrier in each test pit were below the regulatory requirement (0.74 Bq/m{sup 2}- s). Rn concentrations measured using AC samplers in the flux chambers were 60% of concentrations measured using the RAD7, on average, and Rn fluxes computed with the AC data were 9% of those computed with the RAD7 data, on average. Size of the flux chamber had no systematic effect on Rn flux, indicating that 20 y of soil forming processes had not created a pore network causing scale-dependent Rn flux. Geometric mean Rn fluxes in the test pits were similar regardless of differences in surface conditions (vegetation, thickness of soil cover) known to influence soil forming processes and soil structure, but were lower in areas where surface conditions promoted higher water content in the Rn barrier. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Ocean Surface Flux Algorithm Effects on Tropical Indo-Pacific Intraseasonal Precipitation

Surface latent heat fluxes help maintain tropical intraseasonal precipitation. We develop a latent heat flux diagnostic that depicts how latent heat fluxes vary with the near-surface specific humidity vertical gradient (Δq) and surface wind speed (|V|). Compared to fluxes estimated from |V| and Δq measured at tropical moorings and the Coupled Ocean Atmosphere Response Experiment 3.0 (COARE3.0) algorithm, tropical latent heat fluxes in the National Center for Atmospheric Research CEMS2 and Department of Energy E3SMv1 models are significantly overestimated at |V| and Δq extrema. Madden–Julian oscillation (MJO) sensitivity to surface flux algorithm is tested with offline and inline flux corrections. The offline correction adjusts model output fluxes toward mooring-estimated fluxes; the inline correction replaces the original bulk flux algorithm with the COARE3.0 algorithm in atmosphere-only simulations of each model. Both corrections indicate reduced latent heat flux feedback to intraseasonal precipitation, in better agreement with observations, suggesting that model-simulated fluxes are overly supportive for maintaining MJO convection.

54 ENVIRONMENTAL SCIENCES↗

Frozen flux violation, electron demagnetization and magnetic reconnection

We argue that the analogue in collisionless plasma of the collisional diffusion region of magnetic reconnection is properly defined in terms of the demagnetization of the plasma electrons that enable “frozen flux” slippage to occur. This condition differs from the violation of the “frozen-in” condition, which only implies that two fluid effects are involved, rather than the necessary slippage of magnetic flux as viewed in the electron frame. Using 2D Particle In Cell (PIC) simulations, this approach properly finds the saddle point region of the flux function. Our demagnetization conditions are the dimensionless guiding center approximation expansion parameters for electrons which we show are observable and determined locally by the ratio of non-ideal electric to magnetic field strengths. Proxies for frozen flux slippage are developed that (a) are measurable on a single spacecraft, (b) are dimensionless with theoretically justified threshold values of significance, and (c) are shown in 2D simulations to recover distinctions theoretically possible with the (unmeasurable) flux function. A new potentially observable dimensionless frozen flux rate, Λ{sub Φ}, differentiates significant from anecdotal frozen flux slippage. A single spacecraft observable, ϒ, is shown with PIC simulations to be essentially proportional to the unobservable local Maxwell frozen flux rate. This relationship theoretically establishes electron demagnetization in 3D as the general cause of frozen flux slippage. In simple 2D cases with an isolated central diffusion region surrounded by separatrices, these diagnostics uniquely identify the traditional diffusion region (without confusing it with the two fluid “ion-diffusion” region) and clarify the role of the separatrices where frozen flux violations do occur but are not substantial. In the more complicated guide and asymmetric 2D cases, substantial flux slippage regions extend out along, but inside of, the preferred separatrices, demonstrating that Λ{sub Φ} ≠ 0 violations are present over significant distances (in ion inertial units) from the separator identified by the 2D flux function; these violations are, however, generally weaker than seen at known separators in 2D simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

FLUX EMERGENCE IN A MAGNETIZED CONVECTION ZONE

We study the influence of a dynamo magnetic field on the buoyant rise and emergence of twisted magnetic flux ropes and their influence on the global external magnetic field. We ran three-dimensional MHD numerical simulations using the ASH code (anelastic spherical harmonics) and analyzed the dynamical evolution of such buoyant flux ropes from the bottom of the convection zone until the post-emergence phases. The global nature of this model can only very crudely and inaccurately represent the local dynamics of the buoyant rise of the implanted magnetic structure, but nonetheless allows us to study the influence of global effects, such as self-consistently generated differential rotation and meridional circulation, and of Coriolis forces. Although motivated by the solar context, this model cannot be thought of as a realistic model of the rise of magnetic structures and their emergence in the Sun, where the local dynamics are completely different. The properties of initial phases of the buoyant rise are determined essentially by the flux-rope's properties and the convective flows and consequently are in good agreement with previous studies. However, the effects of the interaction of the background dynamo field become increasingly strong as the flux ropes evolve. During the buoyant rise across the convection zone, the flux-rope's magnetic field strength scales as B{proportional_to}{rho}{sup {alpha}}, with {alpha} {approx}< 1. An increase of radial velocity, density, and current density is observed to precede flux emergence at all longitudes. The geometry, latitude, and relative orientation of the flux ropes with respect to the background magnetic field influences the resulting rise speeds, zonal flow amplitudes (which develop within the flux ropes), and the corresponding surface signatures. This influences the morphology, duration and amplitude of the surface shearing, and the Poynting flux associated with magnetic flux-rope emergence. The emerged magnetic flux influences the system's global polarity, leading in some cases to a polarity reversal while inhibiting the background dynamo from doing so in others. The emerged magnetic flux is slowly advected poleward while being diffused and assimilated by the background dynamo field.

79 ASTRONOMY AND ASTROPHYSICS↗

A Hybrid Bulk Algorithm to Predict Turbulent Fluxes over Dry and Wet Bare Soils

Measurements made in the Columbia River basin (Oregon) in an area of irregular terrain during the second Wind Forecast Improvement Project (WFIP2) field campaign are used to develop an optimized hybrid bulk algorithm to predict the surface turbulent fluxes from readily measured or modeled quantities over dry and wet bare or lightly vegetated soil surfaces. The hybrid (synthetic) algorithm combines (i) an aerodynamic method for turbulent flow, which is based on the transfer coefficients (drag coefficient and Stanton number), roughness lengths, and Monin–Obukhov similarity; and (ii) a modified Priestley–Taylor (P-T) algorithm with physically based ecophysiological constraints, which is essentially based on the surface energy budget (SEB) equation. Soil heat flux in the latter case was estimated from measurements of soil temperature and soil moisture. In the framework of the hybrid algorithm, bulk estimates of the momentum flux and the sensible heat flux are derived from a traditional aerodynamic approach, whereas the latent heat flux (or moisture flux) is evaluated from a modified P-T model. Direct measurements of the surface fluxes (turbulent and radiative) and other ancillary atmospheric/soil parameters made during WFIP2 for different soil conditions (dry and wet) are used to optimize and tune the hybrid bulk algorithm. The bulk flux estimates are validated against the measured eddy-covariance fluxes. We also discuss the SEB closure over dry and wet surfaces at various time scales based on the modeled and measured fluxes. Although this bulk flux algorithm is optimized for the data collected during the WFIP2, a hybrid approach can be used for similar flux-tower sites and field campaigns.

54 ENVIRONMENTAL SCIENCES↗

Hall Effect–Mediated Magnetic Flux Transport in Protoplanetary Disks

The global evolution of protoplanetary disks (PPDs) has recently been shown to be largely controlled by the amount of poloidal magnetic flux threading the disk. The amount of magnetic flux must also coevolve with the disk, as a result of magnetic flux transport, a process that is poorly understood. In weakly ionized gas as in PPDs, magnetic flux is largely frozen in the electron fluid, except when resistivity is large. When the disk is largely laminar, we show that the relative drift between the electrons and ions (the Hall drift), and the ions and neutral fluids (ambipolar drift) can play a dominant role on the transport of magnetic flux. Using two-dimensional simulations that incorporate the Hall effect and ambipolar diffusion (AD) with prescribed diffusivities, we show that when large-scale poloidal field is aligned with disk rotation, the Hall effect rapidly drags magnetic flux inward at the midplane region, while it slowly pushes flux outward above/below the midplane. This leads to a highly radially elongated field configuration as a global manifestation of the Hall-shear instability. This field configuration further promotes rapid outward flux transport by AD at the midplane, leading to instability saturation. In quasi-steady state, magnetic flux is transported outward at approximately the same rate at all heights, and the rate is comparable to the Hall-free case. For anti-aligned field polarity, the Hall effect consistently transports magnetic flux outward, leading to a largely vertical field configuration in the midplane region. The field lines in the upper layer first bend radially inward and then outward to launch a disk wind. Overall, the net rate of outward flux transport is about twice as fast as that of the aligned case. In addition, the rate of flux transport increases with increasing disk magnetization. The absolute rate of transport is sensitive to disk microphysics, which remains to be explored in future studies.

79 ASTRONOMY AND ASTROPHYSICS↗

Constraining Genome-Scale Models to Represent the Bow Tie Structure of Metabolism for 13 C Metabolic Flux Analysis

Determination of internal metabolic fluxes is crucial for fundamental and applied biology because they map how carbon and electrons flow through metabolism to enable cell function. 13 C Metabolic Flux Analysis ( 13 C MFA) and Two-Scale 13 C Metabolic Flux Analysis (2S- 13 C MFA) are two techniques used to determine such fluxes. Both operate on the simplifying approximation that metabolic flux from peripheral metabolism into central “core” carbon metabolism is minimal, and can be omitted when modeling isotopic labeling in core metabolism. The validity of this “two-scale” or “bow tie” approximation is supported both by the ability to accurately model experimental isotopic labeling data, and by experimentally verified metabolic engineering predictions using these methods. However, the boundaries of core metabolism that satisfy this approximation can vary across species, and across cell culture conditions. Here, we present a set of algorithms that (1) systematically calculate flux bounds for any specified “core” of a genome-scale model so as to satisfy the bow tie approximation and (2) automatically identify an updated set of core reactions that can satisfy this approximation more efficiently. First, we leverage linear programming to simultaneously identify the lowest fluxes from peripheral metabolism into core metabolism compatible with the observed growth rate and extracellular metabolite exchange fluxes. Second, we use Simulated Annealing to identify an updated set of core reactions that allow for a minimum of fluxes into core metabolism to satisfy these experimental constraints. Together, these methods accelerate and automate the identification of a biologically reasonable set of core reactions for use with 13 C MFA or 2S- 13 C MFA, as well as provide for a substantially lower set of flux bounds for fluxes into the core as compared with previous methods. We provide an open source Python implementation of these algorithms at https://github.com/JBEI/limitfluxtocore.

13C metabolic flux analysis↗

Coupling between radiative flux divergence and turbulence near the surface

Near‐surface turbulent flux and radiation divergence field observations are analysed over a grass‐covered surface located at the Wageningen observatory, the Netherlands. Net radiative flux divergence appears to be a large component of the energy budget near the surface, accounting for a cooling rate of several tens of degrees per day. Long‐wave radiation divergence dominates this radiation divergence flux. We show here that long‐wave flux divergence near the surface is strongly coupled to sensible heat flux, except in high relative humidity (>90%) and foggy conditions. The net long‐wave radiative flux divergence exhibits a sharp gradient within the first 20 m above the surface. This flux divergence is itself strongly coupled to sensible heat flux through adjustments in surface‐layer profiles. Nonetheless, no systematic effect of radiation is witnessed on the turbulent temperature spectrum so that the main effect of near‐surface radiation is on the mean heat budget. As typical flux‐gradient relationships are derived based on observations taken in the near‐surface boundary‐layer region where sharp long‐wave divergence is present, we suspect that those relationships implicitly represent some of the long‐wave divergence term. A modification of Monin–Obukhov Similarity Theory to include the effect of radiative divergence is proposed and discussed. This calls for independent measurements of turbulent fluxes and radiative flux divergence near the surface to re‐derive turbulent flux‐gradient relationships independently of radiative effects.

Gentine, Pierre↗

The importance of sub-meter-scale snow roughness on conductive heat flux of Arctic sea ice

Abstract The conductive heat flux through the snow and ice is a critical component of the mass and energy budgets in the Arctic sea ice system. We use high horizontal resolution (3–15 cm) measurements of snow topography to explore the impacts of sub-meter-scale snow surface roughness on heat flux as simulated by the Finite Element method. Simulating horizontal heat flux in a variable snow cover modestly increases the total simulated heat flux. With horizontal heat flux, as opposed to simple 1D-vertical heat flux modeling, the simulated heat flux is 10% greater than that for uniform snow with the same mean snow thickness for a 31.5 × 21 m region of sea ice (the largest region we studied). Vertical-only (1D) heat flux simulates just a 6% increase for the same region. However, this is highly dependent on observation resolution. Had we measured the snow cover at 1 m horizontal spacing or greater, simulating horizontal heat flux would not have changed the net heat flux from that simulated with vertical-only heat flux. These findings suggest that measuring and modeling snow roughness at sub-meter horizontal scales may be necessary to accurately represent horizontal heat flux on level Arctic sea ice.

Geology↗

Statistical Survey of Magnetic Flux Integral Quantities in Saturn's Magnetosphere

Magnetic flux integral quantities (e.g., flux tube entropy, flux tube content) are conserved quantities under the frozen-in assumption. The change of these quantities often indicates the violation of the frozen-in condition (e.g., interchange instability). In this study, we combine the Cassini CAPS and CHEMS moments data with a steady-state magnetic field model (i.e., the Caudal model) to estimate the flux tube mass and flux tube entropy in Saturn's magnetosphere. Our statistical survey found that the flux tube mass rapidly decreases with radial distance away from Saturn in the inner magnetosphere and roughly levels out in the middle magnetosphere. This indicates that the radial transport processes could occur via a double-reconnection process in the inner magnetosphere, while the radial transport processes appear to be more advective in the middle magnetosphere. Notice that Saturn's magnetosphere is stabilized by a radially increasing profile of flux tube entropy and destabilized by a radially decreasing profile of flux tube content. In this study, we also estimate the expected penetration location by using the flux tube interchange stability formalism developed by Southwood and Kivelson (1987, https://doi.org/10.1029/ja092ia01p00109). The results show that flux tube entropy can play a crucial role in braking the injections, while the flux tube content has a relatively smaller influence on the injected flux tube, being consistent with our previous case study by Wing et al. (2022, https://doi.org/10.3847/1538-4357/ac85b2).

Ma, Xuanye [Embry-Riddle Aeronautical University, ↗

Toward a more realistic representation of surface albedo in NASA CERES-derived surface radiative fluxes

Accurate multidecadal radiative flux records are vital to understand Arctic amplification and constrain climate model uncertainties. Uncertainty in the NASA Clouds and the Earth’s Radiant Energy System (CERES)-derived irradiances is larger over sea ice than any other surface type and comes from several sources. The year-long Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition in the central Arctic provides a rare opportunity to explore uncertainty in CERES-derived radiative fluxes. First, a systematic and statistically robust assessment of surface shortwave and longwave fluxes was conducted using in situ measurements from MOSAiC flux stations. The CERES Synoptic 1degree (SYN1deg) product overestimates the downwelling shortwave flux by +11.40 Wm –2 and underestimates the upwelling shortwave flux by –15.70 Wm –2 and downwelling longwave fluxes by –12.58 Wm –2 at the surface during summer. In addition, large differences are found in the upwelling longwave flux when the surface approaches the melting point (approximately 0°C). The biases in downwelling shortwave and longwave fluxes suggest that the atmosphere represented in CERES is too optically thin. The large negative bias in upwelling shortwave flux can be attributed in large part to lower surface albedo (–0.15) in satellite footprint relative to surface sensors. Additionally, the results show that the spectral surface albedo used in SYN1deg overestimates albedo in visible and mid-infrared bands. A series of radiative transfer model perturbation experiments are performed to quantify the factors contributing to the differences. The CERES-MOSAiC broadband albedo differences (approximately 20 Wm –2 ) explain a larger portion of the upwelling shortwave flux difference than the spectral albedo shape differences (approximately 3 Wm –2 ). In addition, the differences between perturbation experiments using hourly and monthly MOSAiC surface albedo suggest that approximately 25% of the sea ice surface albedo variability is explained by factors not correlated with daily sea ice concentration variability. Biases in net shortwave and longwave flux can be reduced to less than half by adjusting both albedo and cloud inputs toward observed values. The results indicate that improvements in the surface albedo and cloud data would substantially reduce the uncertainty in the Arctic surface radiation budget derived from CERES data products.

54 ENVIRONMENTAL SCIENCES↗

Conceptual Model Update of Mercury Sources and Flux at Y-12 and Upper East Fork Poplar Creek, Oak Ridge,Tennessee

The interactions between the key processes that control the fate and transport of mercury (Hg) in the sediments, groundwater, and stream water within the Y-12 National Security Complex (Y-12) in Oak Ridge, Tennessee are complex, and many aspects are not well understood. In contaminated environments such as Y-12, conceptual models are frequently developed to aid in visualizing and understanding the dynamic nature of the hydrologic, geochemical, and physical environment. These conceptual models integrate data in an internally consistent manner to understand processes that control the fate and transport of contaminants. Over the past few decades of environmental investigation at Y-12, a number of conceptual models have been developed to identify and define various technical processes at various scales. The main purpose of the activities described in this report is to update the 2011 conceptual model for Upper East Fork Poplar Creek using the most recent Hg concentration and flux data and informed scientific interpretation. Emphasis was placed on improving the previous model by implementing a more detailed spatial and temporal approach to visualize transport pathways in the watershed and trends in flux and concentration over time, and to compare baseflow and stormflow system dynamics. Detailed descriptions of historical and current Hg sources and transport pathways are also provided for the east and west ends of the facility. This conceptual model will allow the US Department of Energy to evaluate past and present remedial activities and provide a strong technical basis for prioritizing and optimizing remedial responses in a cost-effective and efficient manner. Twenty years ago, Hg flux from Outfall (OF) 200 at the headwaters of East Fork Poplar Creek represented approximately 20% of the overall flux leaving Y-12. By the time of the 2011 conceptual model report, that percentage had increased to 70%–80%. Flux estimates compiled for this report using data from 2009 to 2018 show that OF200 accounted for roughly 60% of the flux leaving Y-12. The relative role of the four storm drain conduits to OF200 appeared similar to 2011, with OF163 being the greatest contributor to downstream flux. An increase in Hg concentrations and flux occurred throughout and downstream of west end storm drains in 2011 in response to a storm drain cleanout. Leading up to the 2011 storm drain cleanout, annual baseflow flux was higher at OF200A6 than at Station 17. From 2011 to 2018, the trend reversed with concentrations and flux being higher at Station 17 than at OF200A6. A spike in concentrations occurred in July 2018, potentially due to a one-time influx of Hg that occurred during COLEX decontamination and decommissioning activities at Y-12. This report includes recommendations for improving our conceptual understanding of Hg sources, transport pathways, and flux at Y-12. Important recommendations include further evaluation of transport pathways, concurrent measurement of Hg concentration and flux under baseflow and stormflow conditions at multiple sites, and research to understand the connections among periphyton and methylmercury concentration, bioavailability, and bioaccumulation. Ongoing facility decontamination and decommissioning efforts, site characterization, remedial actions, and research are facilitating the collection of new Hg data at Y-12. The conceptual model update provided herein takes advantage of the extensive recent Hg sampling to provide an updated assessment of Hg mass balance and support a refined understanding of Hg behavior at or near Y-12. This assessment should assist in future environmental management decisions and in mitigating the impacts of Hg on the surrounding environment.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Upscaling Methane Flux From Plot Level to Eddy Covariance Tower Domains in Five Alaskan Tundra Ecosystems

Spatial heterogeneity in methane (CH 4 ) flux requires a reliable upscaling approach to reach accurate regional CH 4 budgets in the Arctic tundra. In this study, we combined the CLM-Microbe model with three footprint algorithms to scale up CH 4 flux from a plot level to eddy covariance (EC) tower domains (200 m × 200 m) in the Alaska North Slope, for three sites in Utqiaġvik (US-Beo, US-Bes, and US-Brw), one in Atqasuk (US-Atq) and one in Ivotuk (US-Ivo), for a period of 2013–2015. Three footprint algorithms were the homogenous footprint (HF) that assumes even contribution of all grid cells, the gradient footprint (GF) that assumes gradually declining contribution from center grid cells to edges, and the dynamic footprint (DF) that considers the impacts of wind and heterogeneity of land surface. Simulated annual CH 4 flux was highly consistent with the EC measurements at US-Beo and US-Bes. In contrast, flux was overestimated at US-Brw, US-Atq, and US-Ivo due to the higher simulated CH 4 flux in early growing seasons. The simulated monthly CH 4 flux was consistent with EC measurements but with different accuracies among footprint algorithms. At US-Bes in September 2013, RMSE and NNSE were 0.002 μmol m −2 s −1 and 0.782 using the DF algorithm, but 0.007 μmol m −2 s −1 and 0.758 using HF and 0.007 μmol m −2 s −1 and 0.765 using GF, respectively. DF algorithm performed better than the HF and GF algorithms in capturing the temporal variation in daily CH 4 flux each month, while the model accuracy was similar among the three algorithms due to flat landscapes. Temporal variations in CH 4 flux during 2013–2015 were predominately explained by air temperature (67–74%), followed by precipitation (22–36%). Spatial heterogeneities in vegetation fraction and elevation dominated the spatial variations in CH 4 flux for all five tower domains despite relatively weak differences in simulated CH 4 flux among three footprint algorithms. The CLM-Microbe model can simulate CH 4 flux at both plot and landscape scales at a high temporal resolution, which should be applied to other landscapes. Integrating land surface models with an appropriate algorithm provides a powerful tool for upscaling CH 4 flux in terrestrial ecosystems.

54 ENVIRONMENTAL SCIENCES↗

Long term microparticle impact fluxes on LDEF determined from optical survey of Interplanetary Dust Experiment (IDE) sensors

Many of the IDE metal-oxide-silicon (MOS) capacitor-discharge impact sensors remained active during the entire Long Duration Exposure Facility (LDEF) mission. An optical survey of impact sites on the active surfaces of these sensors has been extended to include all sensors from the low-flux sides of LDEF (i.e. the west or trailing side, the earth end, and the space end) and 5-7 active sensors from each LDEF's high-flux sides (i.e. the east or leading side, the south side, and the north side). This survey was facilitated by the presence of a relatively large (greater than 50 micron diameter) optical signature associated with each impact site on the active sensor surfaces. Of the approximately 4700 impacts in the optical survey data set, 84% were from particles in the 0.5 to 3 micron size range. An estimate of the total number of hypervelocity impacts on LDEF from particles greater than 0.5 micron diameter yields a value of approximately 7 x 10(exp 6). Impact feature dimensions for several dozen large craters on MOS sensors and germanium witness plates are also presented. Impact fluxes calculated from the IDE survey data closely matched surveys of similar size impacts (greater than or equal to 3 micron diameter craters in Al, or marginal penetrations of a 2.4 micron thick Al foil) by other LDEF investigators. Since the first year IDE data were electronically recorded, the flux data could be divided into three long term time periods: the first year, the entire 5.8 year mission, and the intervening 4.8 years (by difference). The IDE data show that there was an order of magnitude decrease in the long term microparticle impact flux on the trailing side of LDEF, from 1.01 to 0.098 x 10(exp -4) m(exp 2)/s, from the first year in orbit compared to years 2-6. The long term flux on the leading edge showed an increase from 8.6 to 11.2 x 10(exp -4) m(exp -2)/s over this same time period. (Short term flux increases up to 10,000 times the background rate were recorded on the leading side during LDEF's first year in orbit.) The overall east/west ratio was 44, but during LDEF's first year in orbit the ratio was 8.5, and during years 2-6 the ratio was 114. Long term microparticle impact fluxes on the space end decreased from 1.12 to 0.55 x 10(exp -4) m(exp -2)/s from the first year in orbit compared to years 2-6. The earth end showed the opposite trend with an increase from 0.16 to 0.38 x 10(exp -4) m(exp -2)/s. Fluxes on rows 6 and 12 decreased from 6.1 to 3.4 and 6.7 to 3.7 x 10(exp -4) m(exp -2)/s, respectively, over the same time periods. This resulted in space/earth microparticle impact flux ratios of 7.1 during the first year and 1.5 during years 2-6, while the south/north, space/north and space/south ratios remained constant at 1.1, 0.16 and 0.17, respectively, during the entire mission. This information indicates the possible identification of long term changes in discrete microparticle orbital debris component contributions to the total impact flux experienced by LDEF. A dramatic decrease in the debris population capable of striking the trailing side was detected that could possibly be attributed to the hiatus of western launch activity experienced from 1986-1989. A significant increase in the debris population that preferentially struck the leading side was also observed and could possibly be attributed to a single breakup event that occurred in September of 1986. A substantial increase in the microparticle debris population that struck the earth end of LDEF, but not the space end, was also detected and could possibly be the result of a single breakup event at low altitude. These results point to the importance of including discrete orbital debris component contribution changes in flux models in order to achieve accurate predictions of the microparticle environment that a particular spacecraft will experience in earth orbit. The only reliable, verified empirical measurements of these changes are reported in this paper. Further time-resolved in-situ measurements of these debris populations are needed to accurately assess model predictions and mitigation practices.

Simon, C. G.↗

Carbon dioxide fluxes over a raised open bog at the Kinosheo Lake tower site during the Northern Wetlands Study (NOWES)

Measurements of carbon dioxide concentration and flux were made above a raised open bog at Lake Kinosheo in the southern Hudson Bay lowlands during the Northern Wetlands Study (NOWES) experiment in 1990. The flux measurements were made using micrometeorological techniques. They provide the first nondisturbing, larger-scale CO2 flux measurements for this ecosystem and are the first to integrate the exchange over the whole 24 hours of the day. Continuous concentration measurements by infrared gas analyzers (IRGA) and spot flask samples were taken over the period July 1 to July 29. Afternoon CO2 values were only 5 to 7 parts per million by volume (ppmv) lower than measurements over the same period at Canadian background monitoring stations. This suggested that there was little draw-down by local photosynthetic sinks. CO2 fluxes were measured at 8 and 18 m by Bowen ratio and eddy correlation methods, respectively. The methods produced comparable results on averaged data but often diverged considerably on individual half-hour results. Fluxes were small. Daytime values averaged to -0.068 mg/sq m/s by eddy correlation and -0.077 mg/sq m/s by Bowen ratio over the period June 25 to July 28 (negative denotes downward flux), while at night, flux densities were +0.062 mg/sq m/s and +0.085 mg/sq m/s. Integration of the mean diurnal curve gave a net flux of -1.7 g/sq m/d. Comparable data for this type of ecosystem were not found. However, Coyne and Kelley (1975), measuring near Barrow, Alaska, over wet meadow tundra dominated by sedges and grasses, found net fluxes of -7.2 g/sq m/d. Typical net CO2 fluxes from other active temperature ecosystems have been found to be -10 to -20 g/sq m/d (Monteith, 1976). Mean half hourly fluxes were almost constant at +0.06 mg/sq m/s through the nighttime hours. About one half-hour after sunrise the flux reversed direction. Uptake peaked about 0900 eastern daylight time (EDT) and then gradually declined but remained downward until near sunset. The early peak was interpreted to signify that the many plants in the bog experienced water stress during the day as evaporative demand increased and nighttime dew was evaporated.

Neumann, H. H.↗