The occurrence and geochemistry of arsenic in groundwaters of the Newark basin of Pennsylvania

Elevated As concentrations in groundwater in the eastern United States have been recognized predominantly in the accretionary geologic terranes of northern New England. A retrospective examination of more than 18,000 existing groundwater samples from the Pennsylvania Department of Environmental Protection (PA DEP) Drinking Water and Sampling Information System database indicates that elevated groundwater As concentrations occur throughout the northern half of the Piedmont Province of Pennsylvania. Chemical analyses of 53 samples collected in 2005 from drinking water wells in this area all had detectable As, and 23% of these samples contained elevated (>133 nmol/L or >10 μg/L) concentrations of As. Elevated concentrations of As in the groundwater samples were most common in the Mesozoic sedimentary strata composed of sandstone and red mudstone with interbedded gray shale, and gray to black siltstone and shale. Arsenic was typically not elevated in groundwater of diabase intrusions of the Newark Basin or in crystalline and calcareous aquifers to the north of the Newark Basin. Geochemical parameters such as pH and oxidation–reduction potential can indicate mobility mechanisms of As in some regions. In this area, measured groundwater conditions were predominantly oxidizing (Eh > +50 mV), and more than 85% of samples contained arsenate as the dominant As species. Variations in pH were strongly correlated to the As concentration, with highest As concentrations observed at pH values greater than 6.4. The original source of As is most likely the black and gray shales that contain some arsenian pyrite with groundwater concentrations likely to be controlled by adsorption/desorption reactions with Fe oxides in the red mudstone aquifer materials.     

Slowness-backazimuth weighted migration: A new array approach to a high-resolution image

We have developed a new array method combining conventional migration with a slowness-backazimuth deviation weighting scheme. All seismic traces are shifted based on the theoretical traveltime of the scattered wave from specific gridpoints in a 3-D volume. Observed slowness and backazimuth are calculated for each raypath and compared with theoretical values in order to estimate slowness and backazimuth deviations. Subsequently, stacked energy calculated by a conventional migration method is weighted by the slowness and backazimuth deviations to suppress any arrival energy whose slowness and backazimuth are inconsistent with the expected theoretical values. This new method was applied to two P- wave data sets which comprise (1) underside reflections at the 410 and 660 km mantle discontinuities and (2) D" reflections as well as their corresponding synthetic data sets. The results show that the weighting scheme dramatically increases the resolution of the migrated images and enables us to obtain well-constrained, focused images, making upper-mantle discontinuities and D" reflections more distinct by reducing their surrounding energy.     

New equipment and processing for magnetotelluric remote reference observations

Robust estimates of magnetotelluric and geomagnetic response functions are determined using the coherency and expected uniformity of the magnetic source field as quality criteria. The method is applied on data sets of three simultaneously recording sites. For the data acquisition we used a new generation of geophysical equipment (S.P.A.M. MkIII), which comprises novel concepts of parallel computing and networked, digital data transmission. The data-processing results show that the amount of noise on the horizontal components of the magnetic field varies considerably in time, between sites and over the frequency range. The removal of such contaminated data beforehand is essential for most data-processing schemes, as the magnetic channels are usually assumed to be free of noise. The standard remote reference method is aimed at reducing bias in response function estimates. However, this does not necessarily improve their precision as our results clearly show. With our method, on the other hand, we can filter out source field irregularities, thereby providing suitable working conditions for the robust algorithm, and eventually obtain considerably improved results. Contrary to previous concepts, we suggest rejecting as much data as feasible in order to concentrate on the remaining parts of high-quality observations.     

The angular momentum content of dwarf galaxies: new challenges for the theory of galaxy formation

We compute the specific angular momentum distributions for a sample of low-mass disc galaxies observed by Swaters. We compare these distributions to those of dark matter haloes obtained by Bullock et al. from high-resolution N -body simulations of structure formation in a ΛCDM universe. We find that although the disc mass fractions are significantly smaller than the universal baryon fraction, the total specific angular momenta of the discs are in good agreement with those of dark matter haloes. This suggests that discs form out of only a small fraction of the available baryons, but yet manage to draw most of the available angular momentum. In addition we find that the angular momentum distributions of discs are clearly distinct from those of the dark matter; discs lack predominantly both low and high specific angular momenta. Understanding these findings in terms of a coherent picture for disc formation is challenging. Cooling, feedback and stripping, which are the main mechanisms to explain the small disc mass fractions found, seem unable to simultaneously explain the angular momentum distributions of the discs. In fact, it seems that the baryons that make up the discs must have been born out of angular momentum distributions that are clearly distinct from those of ΛCDM haloes. However, the dark and baryonic mass components experience the same tidal forces, and it is therefore expected that they should have similar angular momentum distributions. Therefore, understanding the angular momentum content of disc galaxies remains an important challenge for our picture of galaxy formation.     

Detection and extraction of signals from the epoch of reionization using higher-order one-point statistics

Detecting redshifted 21-cm emission from neutral hydrogen in the early Universe promises to give direct constraints on the epoch of reionization (EoR). It will, though, be very challenging to extract the cosmological signal (CS) from foregrounds and noise which are orders of magnitude larger. Fortunately, the signal has some characteristics which differentiate it from the foregrounds and noise, and we suggest that using the correct statistics may tease out signatures of reionization. We generate mock data cubes simulating the output of the Low Frequency Array (LOFAR) EoR experiment. These cubes combine realistic models for Galactic and extragalactic foregrounds and the noise with three different simulations of the CS. We fit out the foregrounds, which are smooth in the frequency direction, to produce residual images in each frequency band. We denoise these images and study the skewness of the one-point distribution in the images as a function of frequency. We find that, under sufficiently optimistic assumptions, we can recover the main features of the redshift evolution of the skewness in the 21-cm signal. We argue that some of these features – such as a dip at the onset of reionization, followed by a rise towards its later stages – may be generic, and give us a promising route to a statistical detection of reionization.     

iVINE – Ionization in the parallel tree/sph code VINE: first results on the observed age-spread around O-stars

We present a three-dimensional, fully parallelized, efficient implementation of ionizing ultraviolet (UV) radiation for smoothed particle hydrodynamics ( sph ) including self-gravity. Our method is based on the sph / tree code vine . We therefore call it iVINE (for Ionization + VINE). This approach allows detailed high-resolution studies of the effects of ionizing radiation from, for example, young massive stars on their turbulent parental molecular clouds. In this paper, we describe the concept and the numerical implementation of the radiative transfer for a plane-parallel geometry and we discuss several test cases demonstrating the efficiency and accuracy of the new method. As a first application, we study the radiatively driven implosion of marginally stable molecular clouds at various distances of a strong UV source and show that they are driven into gravitational collapse. The resulting cores are very compact and dense exactly as it is observed in clustered environments. Our simulations indicate that the time of triggered collapse depends on the distance of the core from the UV source. Clouds closer to the source collapse several 10⁵ yr earlier than more distant clouds. This effect can explain the observed age spread in OB associations where stars closer to the source are found to be younger. We discuss possible uncertainties in the observational derivation of shock front velocities due to early stripping of protostellar envelopes by ionizing radiation.     

Ionizing radiation in smoothed particle hydrodynamics

A new method for the inclusion of ionizing radiation from uniform radiation fields into 3D smoothed particle hydrodynamics (SPHI) simulations is presented. We calculate the optical depth for the Lyman continuum radiation from the source towards the SPHI particles by ray-tracing integration. The time-dependent ionization rate equation is then solved locally for the particles within the ionizing radiation field. Using test calculations, we explore the numerical behaviour of the code with respect to the implementation of the time-dependent ionization rate equation. We also test the coupling of the heating caused by the ionization to the hydrodynamical part of the SPHI code.