Search NASA⌕ Search

SEARCH · Search NASA

Results for “parallel planes”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

HUBBLE FRONTIER FIELDS FIRST COMPLETE CLUSTER DATA: FAINT GALAXIES AT z ∼ 5-10 FOR UV LUMINOSITY FUNCTIONS AND COSMIC REIONIZATION

We present comprehensive analyses of faint dropout galaxies up to z ∼ 10 with the first full-depth data set of the A2744 lensing cluster and parallel fields observed by the Hubble Frontier Fields (HFF) program. We identify 54 dropouts at z ∼ 5-10 in the HFF fields and enlarge the size of the z ∼ 9 galaxy sample obtained to date. Although the number of highly magnified (μ ∼ 10) galaxies is small because of the tiny survey volume of strong lensing, our study reaches the galaxies' intrinsic luminosities comparable to the deepest-field HUDF studies. We derive UV luminosity functions with these faint dropouts, carefully evaluating by intensive simulations the combination of observational incompleteness and lensing effects in the image plane, including magnification, distortion, and multiplication of images, with the evaluation of mass model dependencies. Our results confirm that the faint-end slope, α, is as steep as –2 at z ∼ 6-8 and strengthen the evidence for the rapid decrease of UV luminosity densities, ρ{sub UV}, at z > 8 from the large z ∼ 9 sample. We examine whether the rapid ρ{sub UV} decrease trend can be reconciled with the large Thomson scattering optical depth, τ{sub e}, measured by cosmic microwave background experiments, allowing a large space of free parameters, such as an average ionizing photon escape fraction and a stellar-population-dependent conversion factor. No parameter set can reproduce both the rapid ρ{sub UV} decrease and the large τ {sub e}. It is possible that the ρ{sub UV} decrease moderates at z ≳ 11, that the free parameters significantly evolve toward high z, or that there exist additional sources of reionization such as X-ray binaries and faint active galactic nuclei.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Magnetic Fields in the Massive Dense Cores of the DR21 Filament: Weakly Magnetized Cores in a Strongly Magnetized Filament

We present Submillimeter Array 880 μ m dust polarization observations of six massive dense cores in the DR21 filament. The dust polarization shows complex magnetic field structures in the massive dense cores with sizes of 0.1 pc, in contrast to the ordered magnetic fields of the parsec-scale filament. The major axes of the massive dense cores appear to be aligned either parallel or perpendicular to the magnetic fields of the filament, indicating that the parsec-scale magnetic fields play an important role in the formation of the massive dense cores. However, the correlation between the major axes of the cores and the magnetic fields of the cores is less significant, suggesting that during the core formation, the magnetic fields below 0.1 pc scales become less important than the magnetic fields above 0.1 pc scales in supporting a core against gravity. Our analysis of the angular dispersion functions of the observed polarization segments yields a plane-of-sky magnetic field strength of 0.4–1.7 mG for the massive dense cores. We estimate the kinematic, magnetic, and gravitational virial parameters of the filament and the cores. The virial parameters show that the gravitational energy in the filament dominates magnetic and kinematic energies, while the kinematic energy dominates in the cores. Our work suggests that although magnetic fields may play an important role in a collapsing filament, the kinematics arising from gravitational collapse must become more important than magnetic fields during the evolution from filaments to massive dense cores.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Poster - 38: On the physical and dosimetric properties of 3D printed electron bolus fabricated using polylactic acid

Purpose: 3D printing technology could simplify and improve electron bolus fabrication. The purpose of this study was to characterize the density, dimensional accuracy, uniformity, and attenuation of PLA boluses fabricated with a low-cost 3D printer. Methods: Several solid square slabs were printed with specific requested dimensions and 100% infill using different fill patterns. These pieces were imaged using an x-ray flat panel imager in order to check for uniformity of the prints. Percentage depth doses (PDDs) were measured downstream of the slabs in solid water using a parallel plate chamber and compared to measurements in water in order to characterize attenuation. The dimensions of the PLA slabs were measured using digital calipers. The slabs were also weighed to find their density. Results: The fill pattern used to create boluses can affect the attenuation of the bolus. Fill patterns should be chosen carefully and quality assurance should be done for each printed piece. PLA causes the electron PDD to shift towards shallower depths, compared to water, by 1.7 mm for each centimeter of PLA. Agreement with design dimensions was within 1 mm in the plane of the printer bed, and within 1/3 of a millimeter (roughly the thickness of a single layer), perpendicular to the printer bed. Average density was in the range 1.20 – 1.22. Conclusions: 3D printing shows great promise for use in fabricating electron bolus. This work indicates that printed PLA can be a suitable material provided the increased attenuation is properly accounted for.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Daylight performance of a microstructured prismatic window film in deep open plan offices

Daylight redirecting systems with vertical windows have the potential to offset lighting energy use in deep perimeter zones. Here, a microstructured prismatic film designed for such use was characterized using goniophotometric measurements and ray tracing simulations. The synthetically-generated bidirectional scattering distribution function (BSDF) data were shown to have good agreement with limited measured data for normal incident angles (0–60°). Measured data indicated that the prismatic film was most efficient when vertical angles of incidence were between 18 and 35° and within ±45° of normal incidence to the plane of the window so maximum energy savings across the full depth of the zone occurred over the equinox to winter solstice period. Annual lighting energy use and visual comfort in a deep open plan office zone were evaluated using the Radiance three-phase method in several climates and for south and east-facing window orientations. Lighting energy savings were 39–43% for a 12 m (40 ft) deep south-facing perimeter zone compared to the same zone with no lighting controls. The prismatic film with and without a diffuser controlled glare for views parallel to the window but produced glare for seated viewpoints looking toward the window. At mature market costs, the system was projected to have a simple payback of 2–6 years. Technical challenges encountered throughout the evaluation led to improvements in measurement and modeling tools and stressed the importance of having accurate input data for product development.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2016

Quantization and anomalous structures in the conductance of Si/SiGe quantum point contacts

Quantum point contacts (QPCs) are fabricated on modulation-doped Si/SiGe heterostructures and ballistic transport is studied at low temperatures. We observe quantized conductance with subband separations up to 4 meV and anomalies in the first conductance plateau at 4e{sup 2}/h. At a temperature of T = 22 mK in the linear transport regime, a weak anomalous kink structure arises close to 0.5(4e{sup 2}/h), which develops into a distinct plateau-like structure as temperature is raised up to T = 4 K. Under magnetic field parallel to the wire up to B = 14 T, the anomaly evolves into the Zeeman spin-split level at 0.5(4e{sup 2}/h), resembling the '0.7 anomaly' in GaAs/AlGaAs QPCs. Additionally, a zero-bias anomaly (ZBA) is observed in nonlinear transport spectroscopy. At T = 22 mK, a parallel magnetic field splits the ZBA peak up into two peaks. At B = 0, elevated temperatures lead to similar splitting, which differs from the behavior of ZBAs in GaAs/AlGaAs QPCs. Under finite dc bias, the differential resistance exhibits additional plateaus approximately at 0.8(4e{sup 2}/h) and 0.2(4e{sup 2}/h) known as '0.85 anomaly' and '0.25 anomaly' in GaAs/AlGaAs QPCs. Unlike the first regular plateau at 4e{sup 2}/h, the 0.2(4e{sup 2}/h) plateau is insensitive to dc bias voltage up to at least V{sub DS} = 80 mV, in-plane magnetic fields up to B = 15 T, and to elevated temperatures up to T = 25 K. We interpret this effect as due to pinching off one of the reservoirs close to the QPC. We do not see any indication of lifting of the valley degeneracy in our samples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Planar 2D wireframe DNA origami

Two-dimensional (2D) DNA origami is widely used for applications ranging from excitonics to single-molecule biophysics. Conventional, single-layer 2D DNA origami exhibits flexibility and curvature in solution; however, that may limit its suitability as a 2D structural template. In contrast, 2D wireframe DNA origami rendered with six-helix bundle edges offers local control over duplex orientations with enhanced in-plane rigidity. Here, we investigate the 3D structure of these assemblies using cryo–electron microscopy (cryo-EM). 3D reconstructions reveal a high degree of planarity and homogeneity in solution for polygonal objects with and without internal mesh, enabling 10-Å resolution for a triangle. Coarse-grained simulations were in agreement with cryo-EM data, offering molecular structural insight into this class of 2D DNA origami. Our results suggest that these assemblies may be valuable for 2D material applications and geometries that require high structural fidelity together with local control over duplex orientations, rather than parallel duplex assembly.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

Nonlinear Fourier analysis for discontinuous conductivities: Computational results

Two reconstruction methods of Electrical Impedance Tomography (EIT) are numerically compared for nonsmooth conductivities in the plane based on the use of complex geometrical optics (CGO) solutions to D-bar equations involving the global uniqueness proofs for Calderón problem exposed in Nachman (1996) [43] and Astala and Päivärinta (2006) [6]: the Astala–Päivärinta theory-based low-pass transport matrix method implemented in Astala et al. (2011) [3] and the shortcut method which considers ingredients of both theories. The latter method is formally similar to the Nachman theory-based regularized EIT reconstruction algorithm studied in Knudsen et al. (2009) [34] and several references from there. New numerical results are presented using parallel computation with size parameters larger than ever, leading mainly to two conclusions as follows. First, both methods can approximate piecewise constant conductivities better and better as the cutoff frequency increases, and there seems to be a Gibbs-like phenomenon producing ringing artifacts. Second, the transport matrix method loses accuracy away from a (freely chosen) pivot point located outside of the object to be studied, whereas the shortcut method produces reconstructions with more uniform quality.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Atomic-resolution lithography with an on-chip scanning tunneling microscope

In this work, atomic-resolution lithography with a Microelectromechanical-System (MEMS) based Scanning Tunneling Microscope (STM) is demonstrated for the first time. The microscope consists of a commercial UltraHigh-Vacuum (UHV) STM whose regular tip is replaced with a 1-Degree-of-Freedom (1-DOF) MEMS nanopositioner. This results in a hybrid STM system where XY-plane motions are provided by the piezotube of the original system and Z-axis motion by the MEMS with a higher bandwidth. Sharp tips made of Pt or W are added to the MEMS devices with postfabrication techniques. With this hybrid system, STM-based lithography is demonstrated on an H-passivated Si (100)-2×1 sample under UHV condition. Results prove the capability of the hybrid STM system for atomic-scale lithography. This capability, paired with the small footprint of the MEMS device, makes this approach a candidate for building a high-throughput parallel STM lithography platform by incorporating an array of 1-DOF MEMS devices that perform lithography in parallel.

Alipour, Afshin (ORCID:0000000263252067)↗
DOE OSTI2021

Equatorial Electronic Structure in the Uranyl Ion: Cs 2 UO 2 Cl 4 and Cs 2 UO 2 Br 4

Electric field gradient (EFG) tensors in the equatorial plane of the linear UO 2 2+ ion have been measured by nuclear magnetic resonance (NMR) and nuclear quadrupole resonance (NQR) experiments and computed by relativistic Kohn–Sham methods with and without environment embedding for Cs 2 UO 2 Cl 4 and Cs 2 UO 2 Br 4 . This approach expands the possibilities for probing the electronic structure in uranyl complexes beyond the strongly covalent U–O bonds. The combined analyses find that one of the two largest principal EFG tensor components at the halogen sites points along the U–X bond (X = Cl, Br), and the second is parallel to the UO 2 2+ ion; in Cs 2 UO 2 Cl 4 , the components are nearly equal in magnitude, whereas in Cs 2 UO 2 Br 4 , due to short-range bromide–cesium interactions, the equatorial component is dominant for one pair of Br sites and the axial component is larger for the second pair. Additionally, the directions and relative magnitudes of the field gradient principal axes are found to be sensitive to the σ and π electron donation by the ligands and the model of the environment. Chlorine-35 NQR spectra of 235 U-depleted and 235 U-enriched Cs 2 UO 2 Cl 4 exhibited no uranium-isotope-dependent shift, but the resonance of the depleted sample displayed a 58% broader line width.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Novel Layered Supercell Structure from Bi 2 AlMnO 6 for Multifunctionalities

Layered materials, e.g., graphene and transition metal (di)chalcogenides, holding great promises in nanoscale device applications have been extensively studied in fundamental chemistry, solid state physics and materials research areas. In parallel, layered oxides (e.g., Aurivillius and Ruddlesden–Popper phases) present an attractive class of materials both because of their rich physics behind and potential device applications. In this work, we report a novel layered oxide material with self-assembled layered supercell structure consisting of two mismatch-layered sublattices of [Bi 3 O 3+δ ] and [MO 2 ] 1.84 (M = Al/Mn, simply named BAMO), i.e., alternative layered stacking of two mutually incommensurate sublattices made of a three-layer-thick Bi–O slab and a one-layer-thick Al/Mn–O octahedra slab in the out-of-plane direction. Strong room-temperature ferromagnetic and piezoelectric responses as well as anisotropic optical property have been demonstrated with great potentials in various device applications. Furthermore, the realization of the novel BAMO layered supercell structure in this work has paved an avenue toward exploring and designing new materials with multifunctionalities.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Crystallographic orientation of orthorhombic aragonite using reflection generalized ellipsometry

Here, the 2-modulator generalized ellipsometry microscope (2-MGEM) has been used to study a natural crystal of aragonite. Like its polymorph calcite, aragonite has a large refractive index difference between light polarized parallel to the c-axis and light polarized perpendicular to the c-axis. Unlike calcite, aragonite is orthorhombic, so there is also a very small difference between the refractive indices polarized along the a- and b-directions. As a result, it is not possible to use the 2-MGEM to obtain a definitive map of the optic axis directions of a sample as was possible with calcite, but it is possible to determine approximately the orientation of the c-axis with respect to the surface normal. If the c-axis is in the sample surface plane, it is possible to measure very small deviations of the c-axis direction with an accuracy of ~0.2°. If the c-axis is oriented normal to the sample surface, 2-MGEM data can be used to identify different crystallites due to rotations about the c-axis. For comparison, the orientations of some of the crystallites have also been measured using X-ray Laue and electron beam backscatter diffraction. In addition, spectroscopic generalized ellipsometry measurements have been used to determine the refractive indices of aragonite.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Tuning the structure of the Josephson vortex lattice in Bi 2 Sr 2 CaCu 2 O 8+δ single crystals with pancake vortices

In extremely anisotropic cuprate superconductors a lattice of stacks of pancake vortices nucleates when a magnetic field is applied perpendicular to the copper oxide layers, while an orthogonal lattice of highly elliptical Josephson vortices forms when the applied field is parallel to the layers. Under tilted magnetic fields these sublattices can interact in complex ways to form systems of vortex chains and composite vortex lattices. Here we have used high-resolution scanning Hall microscopy to map the rich tilted-field vortex phase diagram in an underdoped Bi 2 Sr 2 CaCu 2 O 8+δ single crystal. We find that the Josephson vortex lattice spacing has an unexpected non-monotonic dependence on the pancake vortex density reflecting the delicate balance between attractive and repulsive vortex interactions, and actually undergoes a field-driven structural transformation at high out-of-plane fields. We also identify particularly stable composite structures composed of vortex chains separated by an integer number of rows of interstitial pancake vortex stacks and are able to establish the precise evolution of vortex chain phases as the out-of-plane field is increased at small in-plane fields. Lastly, our results are in good semi-quantitative agreement with theoretical models and could enable the development of vortex ratchets and lenses based on the interactions between Josephson and pancake vortices.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Umbrella sampling: a powerful method to sample tails of distributions

We present the umbrella sampling (US) technique and show that it can be used to sample extremely low-probability areas of the posterior distribution that may be required in statistical analyses of data. In this approach sampling of the target likelihood is split into sampling of multiple biased likelihoods confined within individual umbrella windows. We show that the US algorithm is efficient and highly parallel and that it can be easily used with other existing Markov Chain Monte Carlo (MCMC) samplers. The method allows the user to capitalize on their intuition and define umbrella windows and increase sampling accuracy along specific directions in the parameter space. Alternatively, one can define umbrella windows using an approach similar to parallel tempering. We provide a public code that implements US as a standalone PYTHON package. We present a number of tests illustrating the power of the US method in sampling low-probability areas of the posterior and show that this ability allows a considerably more robust sampling of multimodal distributions compared to standard sampling methods. We also present an application of the method in a real world example of deriving cosmological constraints using the supernova type Ia data. We show that US can sample the posterior accurately down to the approximate to 15 sigma credible region in the Omega(m) - Omega(Lambda) plane, while for the same computational effort the affine-invariant MCMC sampling implemented in the emcee code samples the posterior reliably only to approximate to 3 sigma.

cosmology: cosmological parameters↗
DOE OSTI2017

Extreme scale unstructured adaptive CFD for aerodynamic flow control

Understanding the flow of fluid, either liquid or gas, through and around solid bodies has challenged man since the dawn of scientific inquiry. Many of the great minds of science and mathematics have progressively built up a hierarchy of fluid models since fluid flow impacts our lives in so many fundamental ways - from the early days of flow of water in viaducts to cities, to today's flight of planes. This chapter is concerned with the computational modeling of turbulent flow around aerodynamic bodies such as planes and wind turbines. In this case, viscous effects near the solid bodies create very thin boundary layers that yield highly anisotropic (gradients normal to the surface may be 10 6 larger than gradients along the surface) solutions to the governing non-linear partial differential equations (PDE); the Navier-Stokes equations. Furthermore, turbulent flows develop extremely broad ranges of length and time scales motivating the use of discretization methods capable of employing adaptivity and implicit time integration. The combination of these features - non-linear, anisotropy, adaptivity, and implicit - dramatically raise the complexity of the discretization, posing large challenges to efficient scalable parallel implementation. However, through careful design, the more complex algorithms can provide great reductions in computational cost relative to simpler methods (e.g., Cartesian grids with explicit time integration) that are easier to mate efficiently to hardware. In this chapter, we not only describe our approach but we also address the fact that while complex algorithms may never be as efficient flop-for-flop as simple methods, in the important measure of science-per-core-hour, they can still win big by making complex features like adaptivity and implicit methods as efficient and scalable as possible.

42 ENGINEERING↗
DOE OSTI2015

Variational approach for static mirror structures

Anisotropic static plasma equilibria where the parallel and perpendicular pressures are only functions of the amplitude of the local magnetic field are shown to be amenable to a variational principle with a free energy density given by the parallel tension. This approach is used to demonstrate that two-dimensional small-amplitude static magnetic holes constructed from a Grad-Shafranov type equation slightly below the (subcritical) mirror instability threshold identify with lump solitons of KPII equation, but turn out to be unstable. Differently, large-amplitude magnetic structures, which are stable as they realize a minimum of the free energy, are computed using a gradient method within two-dimensional numerical simulations where the regularizing effect of finite Larmor radius corrections is retained. Interestingly, these structures transform from stripes to bubbles when the angle of the magnetic field with the coordinate plane is increased.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Extreme Scale Unstructured Adaptive CFD: From Multiphase Flow to Aerodynamic Flow Control

Understanding the flow of fluid, either liquid or gas, through and around solid bodies has challenged man since the dawn of scientific inquiry. Many of the great minds of science and math have progressively built up a hierarchy of fluid models. This report is concerned with the computational modeling of turbulent flow around aerodynamic bodies such as planes and wind turbines. In this case, viscous effects near the solid bodies create very thin boundary layers that yield highly anisotropic (gradients normal to the surface may be 10^6 times larger than gradients along the surface) solutions to the governing non-linear partial differential equations (PDE); the Navier-Stokes equations. Furthermore, turbulent flows develop extremely broad ranges of length and time scales.This disparity motivates the use of discretization methods capable of employing adaptivity and implicit time integration. The combination of these features (non-linear, anisotropy, adaptivity, andimplicit) dramatically raise the complexity of the discretization, posing large challenges to efficient scalable parallel implementation. However, through careful design, the more complex algorithms can provide great reductions in computational cost relative to simpler methods (e.g., Cartesian grids with explicit time integration) that are easier to mate efficiently to hardware. In this report, we not only describe our approach but we also address the fact that while complex algorithms may never be as efficient flop-for-flop as simple methods, in the important measure of science-per-core-hour, they can still win big by making complex features like adaptivity and implicit methods as efficient and scalable as possible.

42 ENGINEERING↗
DOE OSTI2014

Parallel-scanning tomosynthesis using a slot scanning technique: Fixed-focus reconstruction and the resulting image quality

Purpose: Parallel-scanning tomosynthesis (PS-TS) is a novel technique that fuses the slot scanning technique and the conventional tomosynthesis (TS) technique. This approach allows one to obtain long-view tomosynthesis images in addition to normally sized tomosynthesis images, even when using a system that has no linear tomographic scanning function. The reconstruction technique and an evaluation of the resulting image quality for PS-TS are described in this paper. Methods: The PS-TS image-reconstruction technique consists of several steps (1) the projection images are divided into strips, (2) the strips are stitched together to construct images corresponding to the reconstruction plane, (3) the stitched images are filtered, and (4) the filtered stitched images are back-projected. In the case of PS-TS using the fixed-focus reconstruction method (PS-TS-F), one set of stitched images is used for the reconstruction planes at all heights, thus avoiding the necessity of repeating steps (1)–(3). A physical evaluation of the image quality of PS-TS-F compared with that of the conventional linear TS was performed using a R/F table (Sonialvision safire, Shimadzu Corp., Kyoto, Japan). The tomographic plane with the best theoretical spatial resolution (the in-focus plane, IFP) was set at a height of 100 mm from the table top by adjusting the reconstruction program. First, the spatial frequency response was evaluated at heights of −100, −50, 0, 50, 100, and 150 mm from the IFP using the edge of a 0.3-mm-thick copper plate. Second, the spatial resolution at each height was visually evaluated using an x-ray test pattern (Model No. 38, PTW Freiburg, Germany). Third, the slice sensitivity at each height was evaluated via the wire method using a 0.1-mm-diameter tungsten wire. Phantom studies using a knee phantom and a whole-body phantom were also performed. Results: The spatial frequency response of PS-TS-F yielded the best results at the IFP and degraded slightly as the distance from the IFP increased. A visual evaluation of the spatial resolution using the x-ray test pattern indicated that the resolution was 1.8 lp/mm at the IFP and 1.2 lp/mm at heights of −100 and 100 mm from the IFP. The authors demonstrated that a spatial resolution of 1.2–1.8 lp/mm could be obtained within heights of 200 mm of the IFP. The slice sensitivity varied between 11.1 and 13.8 mm for heights between −50 and 100 mm, and there was no critical change in the slice sensitivity within a height range of 150 mm around the IFP. The phantom results demonstrated that tomosynthesis and long-view images could be reconstructed. Conclusions: PS-TS-F provides tomosynthesis images while using low-cost systems that have no tomographic scanning function, such as tableside-controlled universal R/F systems or universal radiographic systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Integrated hydrogeophysical modelling and data assimilation for geoelectrical leak detection

Time-lapse electrical resistivity tomography (ERT) measurements provide indirect observations of hydrological processes in the Earth's shallow subsurface at high spatial and temporal resolution. ERT has been used in the past decades to detect leaks and monitor the evolution of associated contaminant plumes. Specifically, inverted resistivity images allow visualization of the dynamic changes in the structure of the plume. However, existing methods do not allow the direct estimation of leak parameters (e.g. leak rate, location, etc.) and their uncertainties. We propose an ensemble-based data assimilation framework that evaluates proposed hydrological models against observed time-lapse ERT measurements without directly inverting for the resistivities. Each proposed hydrological model is run through the parallel coupled hydro-geophysical simulation code PFLOTRAN-E4D to obtain simulated ERT measurements. The ensemble of model proposals is then updated using an iterative ensemble smoother. In this paper, we demonstrate the proposed framework on synthetic and field ERT data from controlled tracer injection experiments. Our results show that the approach allows joint identification of contaminant source location, initial release time, and solute loading from the cross-borehole time-lapse ERT data, alongside with an assessment of uncertainties in these estimates. We demonstrate a reduction in site-wide uncertainty by comparing the prior and posterior plume mass discharges at a selected image plane. This framework is particularly attractive to sites that have previously undergone extensive geological investigation (e.g., nuclear sites). It is well suited to complement ERT imaging and we discuss practical issues in its application to field problems.

58 GEOSCIENCES↗