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 379 records · Page 21

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↗

Influence of halide flux on the crystallinity, microstructure and thermoluminescence properties of CdSiO{sub 3}:Co{sup 2+} nanophosphor

Graphical abstract: TL glow curves of CdSiO{sub 3}:Co{sup 2+} different alkali flux (inset without adding flux). Display Omitted Highlights: ► CdSiO{sub 3}:Co{sup 2+} (1–7 mol%) nanocrystalline phosphors synthesized by combustion route. ► Flux effect on thermoluminescence behavior of CdSiO{sub 3}:Co{sup 2+} reported for first time. ► Addition of 2 wt% of flux would drastically enhance the TL properties. ► Well resolved single glow peak at ∼170 °C was recorded for all the samples. ► Among all the alkali flux, NaCl shows highest TL peak intensity. -- Abstract: CdSiO{sub 3}:Co{sup 2+} (1–7 mol %) nanophosphors have been prepared via solution combustion method with post calcination at 800 °C for 2 h for the first time. The formation of expected monoclinic phase was investigated by Powder X-ray diffraction (PXRD) measurements. The effect of different fluxes like NaF, NaCl, NH{sub 4}F and NH{sub 4}Cl on the crystallinity, phase and morphology of CdSiO{sub 3} was investigated in detail. The crystallinity of the samples can be greatly enhanced by using fluxes rather than increasing the calcination temperature. Scanning electronic micrograph (SEM) image shows that the powder morphologies are highly influenced by flux addition. The addition of 2 wt% of fluxes would drastically enhance the crystallinity when NaCl, NH{sub 4}F and NH{sub 4}Cl fluxes are used. A well resolved single thermoluminescent glow peak at ∼170 °C was recorded for all the samples. Among all the halide fluxes, NaCl flux was found to be the potential one in enhancing the TL peak intensity along with crystallinity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

The Effects of Disturbance and Climate on Carbon Storage and the Exchanges of CO 2 Water Vapor and Energy Exchange of Evergreen Coniferous Forests in the Pacific Northwest: Integration of Eddy Flux, Plant and Soil Measurements at a Cluster of Supersites. Final report

This is the final technical report containing a summary of all findings with regard to the following objectives of the project: (1) To quantify and understand the effects of wildfire on carbon storage and the exchanges of energy, CO2, and water vapor in a chronosequence of ponderosa pine (disturbance gradient); (2) To investigate the effects of seasonal and interannual variation in climate on carbon storage and the exchanges of energy, CO2, and water vapor in mature conifer forests in two climate zones: mesic 40-yr old Douglas-fir and semi-arid 60-yr old ponderosa pine (climate gradient); (3) To reduce uncertainty in estimates of CO2 feedbacks to the atmosphere by providing an improved model formulation for existing biosphere-atmosphere models; and (4) To provide high quality data for AmeriFlux and the NACP on micrometeorology, meteorology, and biology of these systems. Objective (1): A study integrating satellite remote sensing, AmeriFlux data, and field surveys in a simulation modeling framework estimated that the pyrogenic carbon emissions, tree mortality, and net carbon exchange associated with four large wildfires that burned ~50,000 hectares in 2002-2003 were equivalent to 2.4% of Oregon statewide anthropogenic carbon emissions over the same two-year period. Most emissions were from the combustion of the forest floor and understory vegetation, and only about 1% of live tree mass was combusted on average. Objective (2): A study of multi-year flux records across a chronosequence of ponderosa pine forests yielded that the net carbon uptake is over three times greater at a mature pine forest compared with young pine. The larger leaf area and wetter and cooler soils of the mature forest mainly caused this effect. A study analyzing seven years of carbon and water dynamics showed that interannual and seasonal variability of net carbon exchange was primarily related to variability in growing season length, which was a linear function of plant-available soil moisture in spring and early summer. A multi-year drought (2001-2003) led to a significant reduction of net ecosystem exchange due to carry-over effects in soil moisture and carbohydrate reserves in plant-tissue. In the same forest, the interannual variability in the rate carbon is lost from the soil and forest floor is considerable and related to the variability in tree growth as much as it is to variability in soil climatic conditions. Objective (3): Flux data from the mature ponderosa pine site support a physical basis for filtering nighttime data with friction velocity above the canopy. An analysis of wind fields and heat transport in the subcanopy at the mesic 40-year old Douglas site yielded that the non-linear structure and behavior of spatial temperature gradients and the flow field require enhanced sensor networks to estimate advective fluxes in the subcanopy of forest to close the surface energy balance in forests. Reliable estimates for flux uncertainties are needed to improve model validation and data assimilation in process-based carbon models, inverse modeling studies and model-data synthesis, where the uncertainties may be as important as the fluxes themselves. An analysis of the time scale dependence of the random and flux sampling error yielded that the additional flux obtained by increasing the perturbation timescale beyond about 10 minutes is dominated by random sampling error, and therefore little confidence can be placed in its value. Artificial correlation between gross ecosystem productivity (GEP) and ecosystem respiration (Re) is a consequence of flux partitioning of eddy covariance flux data when GEP is computed as the difference between NEE and computed daytime Re (e.g. using nighttime Re extrapolated into daytime using soil or air temperatures). Tower-data must be adequately spatially averaged before comparison to gridded model output as the time variability of both is inherently different. The eddy-covariance data collected at the mature ponderosa pine site and the mesic Douglas fir site were used to develop and evaluate a new method to extract the signal of ecosystem respiration directly from daytime net ecosystem exchange. This approach may help reducing uncertainty in carbon budgets by providing direct measurements of ecosystem respiration during daylight conditions by replacing modeled estimates. Objective (4): We submitted our flux and biological and ancillary data to the AmeriFlux web site and to Fluxnet. This includes atmospheric carbon, water, and heat fluxes, soil fluxes, NPP, carbon stocks, LAI, and disturbance history. Fluxnet is updating the original La Thuile files and will include the more recent years of data. They will be using a new approach to compute GPP, following discussions within the network about the need to improve GPP methodology

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.↗

Constraints on hydrothermal heat flux through the oceanic lithosphere from global heat flow

A significant discrepancy exists between the heat flow measured at the seafloor and the higher values predicted by thermal models of the cooling lithosphere. This discrepancy is generally interpreted as indicating that the upper oceanic crust is cooled significantly by hydrothermal circulation. The magnitude of this heat flow discrepancy is the primary datum used to estimate the volume of hydrothermal flow, and the variation in the discrepancy with lithospheric age is the primary constraint on how the hydrothermal flux is divided between near-ridge and off-ridge environments. The resulting estimates are important for investigation of both the thermal structure of the lithosphere and the chemistry of the oceans. We reevaluate the magnitude and age variation of the discrepancy using a global heat flow data set substantially larger than in earlier studies, and the GDHI (Global Depth and Heat Flow) model that better predicts the heat flow. We estimate that of the predicted global oceanic heat flux of 32 x 10(exp 12) W, 34% (11 x 10(exp 12) W) occurs by hydrothermal flow. Approximately 30% of the hydrothermal heat flux occurs in crust younger than 1 Ma, so the majority of this flux is off-ridge. These hydrothermal heat flux estimates are upper bounds, because heat flow measurements require sediment at the site and so are made preferentially at topographic lows, where heat flow may be depressed. Because the water temperature for the near-ridge flow exceeds that for the off-ridge flow, the near-ridge water flow will be even a smaller fraction of the total water flow. As a result, in estimating fluxes from geochemical data, use of the high water temperatures appropriate for the ridge axis may significantly overestimate the heat flux for an assumed water flux or underestimate the water flux for an assumed heat flux. Our data also permit improved estimates of the 'sealing' age, defined as the age where the observed heat flow approximately equals that predicted, suggesting that hydrothermal heat transfer has largely ceased. Although earlier studies suggested major differences in sealing ages for different ocean basins, we find that the sealing ages for the Atlantic, Pacific, and Indian oceans are similar and consistent with the sealing age for the entire data set, 65 +/- 10 Ma. The previous inference of a young (approximately 20 Ma) sealing age for the Pacific appears to have biased downward several previous estimates of the global hydrothermal flux. The heat flow data also provide indirect evidence for the mechanism by which the hydrothermal heat flux becomes small, which has often been ascribed to isolation of the igneous crust from seawater due to the hydraulic conductivity of the intervening sediment. We find, however, that even the least sedimented sites show the systematic increase of the ratio of observed to predicted heat flow with age, although the more sedimented sites have a younger sealing age. Moreover, the heat flow discrepancy persists at heavily sedimented sites until approximately 50 Ma. It thus appears that approximately 100-200 m of sediment is neither necessary nor sufficient to stop hydrothermal heat transfer. We therefore conclude that the age of the crust is the primary control on the fraction of heat transported by hydrothermal flow and that sediment thickness has a lesser effect. This inference is consistent with models in which hydrothermal flow decreases with age due to reduced crustal porosity and hence permeability.

Stein, Carol A.↗

Asymptotic domination of cold relativistic MHD winds by kinetic energy flux

We study the conditions which lead to the conversion of most Poynting flux into kinetic energy flux in cold, relativistic hydromagnetic winds. It is shown that plasma acceleration along a precisely radial flow is extremely inefficient due to the near cancellation of the toroidal magnetic pressure and tension forces. However, if the flux tubes in a flow diverge even slightly faster than radially, the fast magnetosonic point moves inward from infinity to a few times the light cylinder radius. Once the flow becomes supermagnetosonic, further divergence of the flux tubes beyond the fast point can accelerate the flow via the 'magnetic nozzle' effect, thereby further converting Poynting flux to kinetic energy flux. We show that the Grad-Shafranov equation admits a generic family of kinetic energy-dominated asymptotic wind solutions with finite total magnetic flux. The Poynting flux in these solutions vanishes logarithmically with distance. The way in which the flux surfaces are nested within the flow depends only on the ratio of angular velocity to poliodal 4-velocity as a function of magnetic flux. Radial variations in flow structure can be expressed in terms of a pressure boundary condition on the outermost flux surface, provided that no external toriodal field surrounds the flow. For a special case, we show explicitly how the flux surfaces merge gradually to their asymptotes. For flows confined by an external medium of pressure decreasing to zero at infinity we show that, depending on how fast the ambient pressure declines, the final flow state could be either a collimated jet or a wind that fills the entire space. We discuss the astrophysical implications of our results for jets from active galactic nuclei and for free pulsar winds such as that believed to power the Crab Nebula.

Begelman, Mitchell C.↗

Calculation of surface and top of atmosphere radiative fluxes from physical quantities based on ISCCP data sets. 1: Method and sensitivity to input data uncertainties

The largest uncertainty in upwelling shortwave (SW) fluxes (approximately equal 10-15 W/m(exp 2), regional daily mean) is caused by uncertainties in land surface albedo, whereas the largest uncertainty in downwelling SW at the surface (approximately equal 5-10 W/m(exp 2), regional daily mean) is related to cloud detection errors. The uncertainty of upwelling longwave (LW) fluxes (approximately 10-20 W/m(exp 2), regional daily mean) depends on the accuracy of the surface temperature for the surface LW fluxes and the atmospheric temperature for the top of atmosphere LW fluxes. The dominant source of uncertainty is downwelling LW fluxes at the surface (approximately equal 10-15 W/m(exp 2)) is uncertainty in atmospheric temperature and, secondarily, atmospheric humidity; clouds play little role except in the polar regions. The uncertainties of the individual flux components and the total net fluxes are largest over land (15-20 W/m(exp 2)) because of uncertainties in surface albedo (especially its spectral dependence) and surface temperature and emissivity (including its spectral dependence). Clouds are the most important modulator of the SW fluxes, but over land areas, uncertainties in net SW at the surface depend almost as much on uncertainties in surface albedo. Although atmospheric and surface temperature variations cause larger LW flux variations, the most notable feature of the net LW fluxes is the changing relative importance of clouds and water vapor with latitude. Uncertainty in individual flux values is dominated by sampling effects because of large natrual variations, but uncertainty in monthly mean fluxes is dominated by bias errors in the input quantities.

Zhang, Y.-C.↗