Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SAMPLERS”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

THE ORBIT AND MASS OF THE THIRD PLANET IN THE KEPLER-56 SYSTEM

While the vast majority of multiple-planet systems have orbital angular momentum axes that align with the spin axis of their host star, Kepler-56 is an exception: its two transiting planets are coplanar yet misaligned by at least 40° with respect to the rotation axis of their host star. Additional follow-up observations of Kepler-56 suggest the presence of a massive, non-transiting companion that may help explain this misalignment. We model the transit data along with Keck/HIRES and HARPS-N radial velocity data to update the masses of the two transiting planets and infer the physical properties of the third, non-transiting planet. We employ a Markov Chain Monte Carlo sampler to calculate the best-fitting orbital parameters and their uncertainties for each planet. We find the outer planet has a period of 1002 ± 5 days and minimum mass of 5.61 ± 0.38 M{sub Jup}. We also place a 95% upper limit of 0.80 m s{sup −1} yr{sup −1} on long-term trends caused by additional, more distant companions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

DEVELOPING ATMOSPHERIC RETRIEVAL METHODS FOR DIRECT IMAGING SPECTROSCOPY OF GAS GIANTS IN REFLECTED LIGHT. I. METHANE ABUNDANCES AND BASIC CLOUD PROPERTIES

Upcoming space-based coronagraphic instruments in the next decade will perform reflected light spectroscopy and photometry of cool directly imaged extrasolar giant planets. We are developing a new atmospheric retrieval methodology to help assess the science return and inform the instrument design for such future missions, and ultimately interpret the resulting observations. Our retrieval technique employs a geometric albedo model coupled with both a Markov chain Monte Carlo Ensemble Sampler ( emcee ) and a multimodal nested sampling algorithm ( MultiNest ) to map the posterior distribution. This combination makes the global evidence calculation more robust for any given model and highlights possible discrepancies in the likelihood maps. As a proof of concept, our current atmospheric model contains one or two cloud layers, methane as a major absorber, and a H{sub 2}–He background gas. This 6-to-9 parameter model is appropriate for Jupiter-like planets and can be easily expanded in the future. In addition to deriving the marginal likelihood distribution and confidence intervals for the model parameters, we perform model selection to determine the significance of methane and cloud detection as a function of expected signal-to-noise ratio in the presence of spectral noise correlations. After internal validation, the method is applied to realistic spectra of Jupiter, Saturn, and HD 99492c, a model observing target. We find that the presence or absence of clouds and methane can be determined with high confidence, while parameter uncertainties are model dependent and correlated. Such general methods will also be applicable to the interpretation of direct imaging spectra of cloudy terrestrial planets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Isotropic at the break? 3D kinematics of Milky Way halo stars in the foreground of M31

We present the line-of-sight (LOS) velocities for 13 distant main sequence Milky Way halo stars with published proper motions (PMs). The PMs were measured using long baseline (5–7 years) multi-epoch Hubble Space Telescope/Advanced Camera for Surveys photometry, and the LOS velocities were extracted from deep (5–6 hr integrations) Keck II/DEIMOS spectra. We estimate the parameters of the velocity ellipsoid of the stellar halo using a Markov chain Monte Carlo ensembler sampler method. The velocity second moments in the directions of the Galactic (l, b, LOS) coordinate system are 〈v{sub l}{sup 2}〉{sup 1/2}=138{sub −26}{sup +43} km s{sup −1}, 〈v{sub b}{sup 2}〉{sup 1/2}=88{sub −17}{sup +28} km s{sup −1}, and 〈v{sub LOS}{sup 2}〉{sup 1/2}=91{sub −14}{sup +27} km s{sup −1}. We use these ellipsoid parameters to constrain the velocity anisotropy of the stellar halo. Ours is the first measurement of the anisotropy parameter β using 3D kinematics outside of the solar neighborhood. We find β=−0.3{sub −0.9}{sup +0.4}, consistent with isotropy and lower than solar neighborhood β measurements by 2σ (β{sub SN} ∼ 0.5–0.7). We identify two stars in our sample that are likely members of the known TriAnd substructure, and excluding these objects from our sample increases our estimate of the anisotropy to β=0.1{sub −1.0}{sup +0.4}, which is still lower than solar neighborhood measurements by 1σ. The potential decrease in β with Galactocentric radius is inconsistent with theoretical predictions, though consistent with recent observational studies, and may indicate the presence of large, shell-type structure (or structures) at r ∼ 25 kpc. The methods described in this paper will be applied to a much larger sample of stars with 3D kinematics observed through the ongoing HALO7D program.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Cosmological parameters from CMB maps without likelihood approximation

We propose an efficient Bayesian Markov chain Monte Carlo (MCMC) algorithm for estimating cosmological parameters from cosmic microwave background (CMB) data without the use of likelihood approximations. It builds on a previously developed Gibbs sampling framework that allows for exploration of the joint CMB sky signal and power spectrum posterior, P(s,C{sub ℓ}|d), and addresses a long-standing problem of efficient parameter estimation simultaneously in regimes of high and low signal-to-noise ratio. To achieve this, our new algorithm introduces a joint Markov chain move in which both the signal map and power spectrum are synchronously modified, by rescaling the map according to the proposed power spectrum before evaluating the Metropolis–Hastings accept probability. Such a move was already introduced by Jewell et al., who used it to explore low signal-to-noise posteriors. However, they also found that the same algorithm is inefficient in the high signal-to-noise regime, since a brute-force rescaling operation does not account for phase information. This problem is mitigated in the new algorithm by subtracting the Wiener filter mean field from the proposed map prior to rescaling, leaving high signal-to-noise information invariant in the joint step, and effectively only rescaling the low signal-to-noise component. To explore the full posterior, the new joint move is then interleaved with a standard conditional Gibbs move for the sky map. We apply our new algorithm to simplified simulations for which we can evaluate the exact posterior to study both its accuracy and its performance, and find good agreement with the exact posterior; marginal means agree to ≲0.006σ and standard deviations to better than ∼3%. The Markov chain correlation length is of the same order of magnitude as those obtained by other standard samplers in the field.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Characterizing Earth Analogs in Reflected Light: Atmospheric Retrieval Studies for Future Space Telescopes

Space-based high-contrast imaging mission concepts for studying rocky exoplanets in reflected light are currently under community study. We develop an inverse modeling framework to estimate the science return of such missions given different instrument design considerations. By combining an exoplanet albedo model, instrument noise model, and ensemble Markov chain Monte Carlo sampler, we explore retrievals of atmospheric and planetary properties for Earth twins as a function of signal-to-noise ratio (S/N) and resolution (R). Our forward model includes Rayleigh-scattering, single-layer water clouds with patchy coverage, and pressure-dependent absorption due to water vapor, oxygen, and ozone. We simulate data at R = 70 and 140 from 0.4 to 1.0 μm with S/N = 5, 10, 15, and 20 at 550 nm (i.e., for HabEx/LUVOIR-type instruments). At these same S/Ns, we simulate data for WFIRST paired with a starshade, which includes two photometric points between 0.48 and 0.6 μm and R = 50 spectroscopy from 0.6 to 0.97 μm. Given our noise model for WFIRST-type detectors, we find that weak detections of water vapor, ozone, and oxygen can be achieved with observations with at least R = 70/S/N = 15 or R = 140/S/N = 10 for improved detections. Meaningful constraints are only achieved with R = 140/S/N = 20 data. The WFIRST data offer limited diagnostic information, needing at least S/N = 20 to weakly detect gases. Most scenarios place limits on planetary radius but cannot constrain surface gravity and, thus, planetary mass.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Detecting and Characterizing Mg ii Absorption in DESI Survey Validation Quasar Spectra

Abstract We present findings of the detection of Magnesium II (Mg ii , λ = 2796, 2803 Å) absorbers from the early data release of the Dark Energy Spectroscopic Instrument (DESI). DESI is projected to obtain spectroscopy of approximately 3 million quasars (QSOs), of which over 99% are anticipated to be at redshifts greater than z > 0.3, such that DESI would be able to observe an associated or intervening Mg ii absorber illuminated by the background QSO. We have developed an autonomous supplementary spectral pipeline that detects these systems through an initial line-fitting process and then confirms the line properties using a Markov Chain Monte Carlo sampler. Based upon a visual inspection of the resulting systems, we estimate that this sample has a purity greater than 99%. We have also investigated the completeness of our sample in regard to both the signal-to-noise properties of the input spectra and the rest-frame equivalent width ( W 0 ) of the absorber systems. From a parent catalog containing 83,207 quasars, we detect a total of 23,921 Mg ii absorption systems following a series of quality cuts. Extrapolating from this occurrence rate of 28.8% implies a catalog at the completion of the five-year DESI survey that will contain over eight hundred thousand Mg ii absorbers. The cataloging of these systems will enable significant further research because they carry information regarding circumgalactic medium environments, the distribution of intervening galaxies, and the growth of metallicity across the redshift range 0.3 ≤ z < 2.5.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Optimizing Optical Searches for Supermassive Black Hole Binaries in Active Galactic Nuclei Light Curves: Fourier versus Bayesian Periodicity Detection

Simulations predict that supermassive black hole binaries (SMBHBs) will exhibit periodic brightness variations that may exceed the stochastic variability intrinsic to active galactic nuclei (AGN). In this paper, we simulate SMBHBs with damped random walk (DRW) AGN variability and an added sinusoidal signal from the orbital motion, and test three methods—a generalized Lomb–Scargle periodogram (GLSP), a nested Bayesian sampler (NBS), and a weighted wavelet z-transform (or WWZ)—to determine which is best at recovering the periodicity. Our simulated light curves follow the properties of the Catalina Real-Time Transient Survey (or CRTS), Legacy Survey of Space and Time (LSST), and Zwicky Transient Facility (ZTF) to best inform current and future SMBHB searches. We map a broad range of parameter space and identify which DRW-only light curves best mimic periodicity and pass each method’s model selection. The NBS performs best at detecting periodicity and filtering out DRW-only light curves. Combined candidate selection with both the NBS and GLSP significantly reduces false-positive rates (FPRs) with marginal impact on true-positive rates (TPRs). With this joint model selection pipeline, we find the lowest FPRs in ZTF-like simulations and the highest detection rates in LSST-like simulations. Using a modified computation of the false-alarm probability with GLSP, we efficiently triage LSST AGN light curves (∼10 7 light curves in ∼10–30 hr) and achieve TPRs and FPRs of ∼40% and ∼0.5%, respectively.

Banaszak, Sebastian M. [Vanderbilt Univ., Nashvill↗
DOE OSTI2019

Estimation and uncertainty quantification for the output from quantum simulators

The problem of estimating certain distributions over {0, 1} d is considered here. The distribution represents a quantum system of d qubits, where there are non-trivial dependencies between the qubits. A maximum entropy approach is adopted to reconstruct the distribution from exact moments or observed empirical moments. The Robbins Monro algorithm is used to solve the intractable maximum entropy problem, by constructing an unbiased estimator of the un-normalized target with a sequential Monte Carlo sampler at each iteration. In the case of empirical moments, this coincides with a maximum likelihood estimator. A Bayesian formulation is also considered in order to quantify uncertainty a posteriori. Several approaches are proposed in order to tackle this challenging problem, according to recently developed methodologies. In particular, unbiased estimators of the gradient of the log posterior are constructed and used within a provably convergent Langevin-based Markov chain Monte Carlo method. The methods are illustrated on classically simulated output from quantum simulators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Impact of collection container material and holding times on sample integrity for mercury and methylmercury in water

Mercury (Hg) and methylmercury (MeHg) concentrations in streamwater can vary on short timescales (hourly or less) during storm flow and on a diel cycle; the frequency and timing of sampling required to accurately characterize these dynamics may be difficult to accomplish manually. Automated sampling can assist in sample collection; however use has been limited for Hg and MeHg analysis due to stability concerns of trace concentrations during extended storage times. We examined the viability of using automated samplers with disposable low-density polyethylene (LDPE) sample bags to collect industrially contaminated streamwater for unfiltered and filtered Hg and MeHg analysis. Specifically we investigated the effect of holding times ranging from hours to days on streamwater collected during baseflow and storm flow. Unfiltered and filtered Hg and MeHg concentrations decreased with increases in time prior to sample processing; holding times of 24 hours or less resulted in concentration changes (mean 11 7% different) similar to variability in duplicates collected manually during analogous field conditions (mean 7 10% different). Comparisons of samples collected with manual and automated techniques throughout a year for a wide range of stream conditions were also found to be similar to differences observed between duplicate grab samples. These results demonstrate automated sampling into LDPE bags with holding times of 24 hours or less can be effectively used to collect streamwater for Hg and MeHg analysis, and encourage the testing of these materials and methods for implementation in other aqueous systems where high-frequency sampling is warranted.

holding time↗
DOE OSTI2014

Ambient aromatic hydrocarbon measurements at Welgegund, South Africa

Aromatic hydrocarbons are associated with direct adverse human health effects and can have negative impacts on ecosystems due to their toxicity, as well as indirect negative effects through the formation of tropospheric ozone and secondary organic aerosol that affect human health, crop production and regional climate. Measurements were conducted at the Welgegund measurement station (South Africa) that is considered to be a regionally representative background site. However, the site is occasionally impacted by plumes from major anthropogenic source regions in the interior of South Africa, which include the western Bushveld Igneous Complex (e.g. platinum, base metal and ferrochrome smelters), the eastern Bushveld Igneous Complex (platinum and ferrochrome smelters), the Johannesburg-Pretoria metropolitan conurbation (>10 million people), the Vaal Triangle (e.g. petrochemical and industries), the Mpumalanga Highveld (e.g. coal-fired power plants and petrochemical industry) and also a region of anti-cyclonic recirculation of air mass over the interior of South Africa. The aromatic hydrocarbon measurements were conducted with an automated sampler on Tenax-TA and Carbopack-B adsorbent tubes with heated inlet for one year. Samples were collected twice a week for two hours during daytime and two hours 1 during night-time. A thermal desorption unit, connected to a gas chromatograph and a mass 2 selective detector was used for sample preparation and analysis. Results indicated that the 3 monthly median total aromatic hydrocarbon concentrations ranged between 0.01 to 3.1 ppb. 4 Benzene levels did not exceed local air quality standards. Toluene was the most abundant 5 species, with an annual median concentration of 0.63 ppb. No statistically significant 6 differences in the concentrations measured during daytime and night-time were found and no distinct seasonal patterns were observed. Air mass back trajectory analysis proved that the lack of seasonal cycles could be attributed to patterns determining the origin of the air masses sampled. Aromatic hydrocarbon concentrations were in general significantly higher in air masses that passed over anthropocentrically impacted regions. Interspecies correlations and ratios gave some indications of the possible sources for the different aromatic hydrocarbons in the source regions defined in the paper. The highest contribution of aromatic hydrocarbon concentrations to ozone formation potential was also observed in plumes passing over anthropocentrically impacted regions.

Jaars, K.↗
DOE OSTI2025

Seasonal investigation of ultrafine-particle organic composition in an eastern Amazonian rainforest

Abstract. Reports on the composition of ultrafine particles (<100 nm in diameter) in the Amazon are scarce, due in part to the fact that new-particle formation has rarely been observed near ground level. Ultrafine particles near the surface have nevertheless been observed, leaving open questions regarding the sources and chemistry of their formation and growth, particularly as these vary across seasons. Here, we present measurements of the composition of ultrafine particles collected in the Tapajós National Forest (2.857° S, 54.959° W) during three different seasonal periods: 10–30 September 2016 (SEP), 18 November–23 December 2016 (DEC), and 22 May–21 June 2017 (JUN). Size-selected (5–70 nm) particles were collected daily (for 22 h each day) using an offline sampler. Samples collected during the three time periods were compiled and analyzed using liquid chromatography coupled with Orbitrap high-resolution mass spectrometry. Our findings suggest a sustained influence of isoprene organosulfate chemistry on ultrafine particles from the different periods. We present chemical evidence that indicates that biological-spore fragmentation impacted ultrafine-particle composition during the late wet season (JUN), while chemical markers for biomass burning and secondary chemistry peaked during the dry season (SEP and DEC). Higher oxidation states and degrees of unsaturation were observed for organics in the dry season (SEP and DEC), suggesting greater extents of aerosol aging. Finally, applying a volatility parameterization to the observed compounds suggests that organic sulfur species are likely key drivers of new-particle growth in the region due to their low volatility compared to other species.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Quantitative chemical assay of nanogram-level particulate matter using aerosol mass spectrometry: characterization of particles collected from uncrewed atmospheric measurement platforms

Abstract. Aerosol generation techniques have expanded the utility of aerosol mass spectrometry (AMS) for offline chemical analysis of airborne particles and droplets. However, standard aerosolization techniques require relatively large liquid volumes (e.g., several milliliters) and high sample masses that limit their utility. Here we report the development and characterization of a micronebulization AMS (MN-AMS) technique that requires as low as 10 µL of sample and can provide the quantification of the nanogram level of organic and inorganic substances via the usage of an isotopically labeled internal standard (34SO42-). Using standard solutions, the detection limits for this technique were determined at 0.19, 0.75, and 2.2 ng for sulfate, nitrate, and organics, respectively. The analytical recoveries for these species are 104 %, 87 %, and 94 %, respectively. This MN-AMS technique was applied successfully to analyze filter and impactor samples collected using miniature particulate matter (PM) samplers deployable on uncrewed atmospheric measurement platforms, such as uncrewed aerial systems (UASs) and tethered balloon systems (TBSs). Chemical composition of PM samples collected from a UAS field campaign conducted at the Department of Energy's (DOE) Southern Great Plains (SGP) observatory was characterized. The offline MN-AMS data compared well with the in situ PM composition measured by a co-located aerosol chemical speciation monitor (ACSM). In addition, the MN-AMS and ion chromatography (IC) agreed well for measurements of sulfate and nitrate concentrations in the PM extracts. This study demonstrates the utility of combining MN-AMS with uncrewed measurement platforms to provide quantitative measurements of ambient PM composition.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Simulations of the collection of mesospheric dust particles with a rocket instrument

We investigate the collection of dust particles in the mesosphere with the MESS (MEteoric Smoke Sampler) instrument that is designed to fly on a sounding rocket. We assume that the ice particles that form in the polar mesosphere between 80 and 85 km altitude in summer contain meteoric smoke particles; and these should be collected with MESS. The instrument consists of a collection device with an opening and closure mechanism, as well as an attached conic funnel which increases the sampling area in comparison to the collection area. Dust particles are collected either directly after passing through the instrument or indirectly after colliding with and fragmenting on the funnel wall. We calculate the dust and fragment trajectories in the detector to determine the collection efficiency for different particle sizes, rocket velocities, and heights, and we find the final velocities and the temperatures of the particles. The considered design has a sampling area of 62.78 mm diameter and a collection area of 20 mm diameter. For the conditions at the rocket launch site in Andøya, Norway, we estimate the collection of meteoric smoke particles contained in the ice particles to be ~10 12 –10 14 amu mm –2 . The estimated temperatures suggest that the composition of these smoke particles is not affected by the collection. Our calculations also show that keeping the instrument open above 85 km altitude increases the amount of small smoke particles that are directly collected. The directly collected smoke particles are heated as they decelerate, which can affect their composition.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

A novel methodology for assessing the hygroscopicity of aerosol filter samples

Abstract. Due to US regulations, concentrations of hygroscopic inorganic sulfate and nitrate have declined in recent years, leading to an increased importance of the hygroscopic nature of organic matter (OM). The hygroscopicity of OM is poorly characterized because only a fraction of the multitude of organic compounds in the atmosphere is readily measured, and there is limited information on their hygroscopic behaviors. Hygroscopicity of aerosol is traditionally measured using a humidified tandem differential mobility analyzer (HTDMA) or electrodynamic balance (EDB). EDB measures water uptake by a single particle. For ambient and chamber studies, HTDMA measurements provide water uptake and particle size information but not chemical composition. To fill this information gap, we developed a novel methodology to assess the water uptake by particles collected on Teflon filters. This method uses the same filter sample for both hygroscopicity measurements and chemical characterization, thereby providing an opportunity to link the measured hygroscopicity with ambient particle composition. To test the method, hygroscopic measurements were conducted in the laboratory for ammonium sulfate, sodium chloride, glucose, and malonic acid, which were collected on 25 mm Teflon filters using an aerosol generator and sampler. Constant-humidity solutions (CHSs), including potassium chloride, barium chloride dihydrate, and potassium sulfate, were employed in a saturated form to maintain the relative humidity (RH) at approximately 84 %, 90 %, and 97 % in small chambers. Our preliminary experiments revealed that, without the pouch, water uptake measurements were not feasible due to rapid water loss during weighing. Additionally, we observed some absorption by the aluminum pouch itself. To account for this, concurrent measurements were conducted for both the loaded and the blank filters at each RH level. Thus, the dry loaded and blank Teflon filters were placed in aluminum pouches with one side open and in RH-controlled chambers for more than 24 h. The wet loaded samples and wet blanks were then weighed using an ultramicrobalance to determine the water uptake by the respective compound and the blank Teflon filter. The net amount of water absorbed by each compound was calculated by subtracting the water uptake of the blank filter from that of the wet loaded filter. Hygroscopic parameters, including the water-to-solute (W / S) ratio, molality, mass fraction solute (mfs), and growth factors (GFs), were calculated from the measurements. The results obtained are consistent with those reported by the Extended Aerosol Inorganics Model (E-AIM) and previous studies utilizing HTDMA and EDB for these compounds, highlighting the accuracy of this new methodology. This new approach enables the hygroscopicity and chemical composition of individual filter samples to be assessed so that in complex mixtures, such as chamber and ambient samples, the total water uptake can be parsed between the inorganic and organic components of the aerosol.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

A novel Whole Air Sample Profiler (WASP) for the quantification of volatile organic compounds in the boundary layer

The emission and fate of reactive VOCs is of inherent interest to those studying chemical biosphere-atmosphere interactions. In-canopy VOC observations are obtainable using tower-based samplers, but the lack of suitable sampling systems for the full boundary 5 layer has limited the data characterizing the vertical structure of such gases above the canopy height and still in the boundary layer. This is the important region where many reactive VOCs are oxidized or otherwise removed. Here we describe an airborne sampling system designed to collect a vertical profile of air into a 3/800 OD tube 150m in length. The inlet ram air pressure is used to flow sampled air through the 10 tube, which results in a varying flow rate based on aircraft speed and altitude. Since aircraft velocity decreases during ascent, it is necessary to account for the variable flow rate into the tube. This is accomplished using a reference gas that is pulsed into the air stream so that the precise altitude of the collected air can be reconstructed post-collection. The pulsed injections are also used to determine any significant effect 15 from diffusion/mixing within the sampling tube, either during collection or subsequent extraction for gas analysis. This system has been successfully deployed, and we show some measured vertical profiles of isoprene and its oxidation products methacrolein and methyl vinyl ketone from a mixed canopy near Columbia, Missouri.

Mak, J. E.↗
DOE OSTI2016

A flexible importance sampling method for integrating subgrid processes

Numerical models of weather and climate need to compute grid-box-averaged rates of physical processes such as microphysics. These averages are computed by integrating subgrid variability over a grid box. For this reason, an important aspect of atmospheric modeling is spatial integration over subgrid scales. The needed integrals can be estimated by Monte Carlo integration. Monte Carlo integration is simple and general but requires many evaluations of the physical process rate. To reduce the number of function evaluations, this paper describes a new, flexible method of importance sampling. It divides the domain of integration into eight categories, such as the portion that contains both precipitation and cloud, or the portion that contains precipitation but no cloud. It then allows the modeler to prescribe the density of sample points within each of the eight categories. The new method is incorporated into the Subgrid Importance Latin Hypercube Sampler (SILHS). The resulting method is tested on drizzling cumulus and stratocumulus cases. In the cumulus case, the sampling error can be considerably reduced by drawing more sample points from the region of rain evaporation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

High-frequency monitoring of water fluxes and nutrient loads to assess the effects of controlled drainage on water storage and nutrient transport

High nitrogen (N) and phosphorus (P) fluxes from upstream agriculture threaten aquatic ecosystems in surface waters and estuaries, especially in areas characterized by high agricultural N and P inputs and densely drained catchments like the Netherlands. Controlled drainage has been recognized as an effective option to optimize soil moisture conditions for agriculture and to reduce unnecessary losses of fresh water and nutrients. This is achieved by introducing control structures with adjustable overflow levels into subsurface tube drain systems. A small-scale (1 ha) field experiment was designed to investigate the hydrological and chemical changes after introducing controlled drainage. Precipitation rates and the response of water tables and drain fluxes were measured in the periods before the introduction of controlled drainage (2007–2008) and after (2009–2011). For the N and P concentration measurements, auto-analyzers for continuous records were combined with passive samplers for time-averaged concentrations at individual drain outlets. The experimental setup enabled the quantification of changes in the water and solute balance after introducing controlled drainage. The results showed that introducing controlled drainage reduced the drain discharge and increased the groundwater storage in the field. To achieve this, the overflow levels have to be elevated in early spring, before the drain discharge stops due to dryer conditions and falling groundwater levels. The groundwater storage in the field would have been larger if the water levels in the adjacent ditch were controlled as well by an adjustable weir. The N concentrations and loads increased, which was largely related to elevated concentrations in one of the three monitored tube drains. The P loads via the tube drains reduced due to the reduction in discharge after introducing controlled drainage. Furthermore, this may be counteracted by the higher groundwater levels and the larger contribution of N- and P-rich shallow groundwater and overland flow to the surface water.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Understanding the relative importance of vertical and horizontal flow in ice-wedge polygons

Abstract. Ice-wedge polygons are common Arctic landforms. The future of these landforms in a warming climate depends on the bidirectional feedback between the rate of ice-wedge degradation and changes in hydrological characteristics. This work aims to better understand the relative roles of vertical and horizontal water fluxes in the subsurface of polygonal landscapes, providing new insights and data to test and calibrate hydrological models. Field-scale investigations were conducted at an intensively instrumented location on the Barrow Environmental Observatory (BEO) near Utqiaġvik, AK, USA. Using a conservative tracer, we examined controls of microtopography and the frost table on subsurface flow and transport within a low-centered and a high-centered polygon. Bromide tracer was applied at both polygons in July 2015 and transport was monitored through two thaw seasons. Sampler arrays placed in polygon centers, rims, and troughs were used to monitor tracer concentrations. In both polygons, the tracer first infiltrated vertically until encountering the frost table and was then transported horizontally. Horizontal flow occurred in more locations and at higher velocities in the low-centered polygon than in the high-centered polygon. Preferential flow, influenced by frost table topography, was significant between polygon centers and troughs. Estimates of horizontal hydraulic conductivity were within the range of previous estimates of vertical conductivity, highlighting the importance of horizontal flow in these systems. This work forms a basis for understanding complexity of flow in polygonal landscapes.

54 ENVIRONMENTAL SCIENCES↗