Model for the formation of arsenic contamination in groundwater. 1. Datong Basin,China

Mineral equilibria were analyzed in the system As-bearing rock-meteoric water. It was shown that carbonate rocks are the most probable source of As and Sr in the waters of the Datong Basin (People's Republic of China). The reason for groundwater enrichment in As is the shift of the equilibrium FeCO₃ (siderite) + H₂O = FeOOH(goethite) + CO₂(g) + H₂(g) to the left (toward siderite formation) owing to organic matter oxidation by atmospheric oxygen and an increase in the equilibrium partial pressure of CO₂, while the Eh of the system remains below ?0.30 ± 0.06 V.     

Friedrichbeckeite, K (□0.5Na0.5)2 (Mg0.8Mn0.1Fe0.1)2 (Be0.6 Mg0.4)3 [Si12O30], a new milarite-type mineral from the Bellerberg volcano, Eifel area, Germany

Friedrichbeckeite is a new milarite-type mineral. It was found in a single silicate-rich xenolith from a quarry at the Bellerberg volcano near Ettringen, eastern Eifel volcanic area, Germany. It forms thin tabular crystals flattened on {0001}, with a maximum diameter of 0.6 mm and a maximum thickness of 0.1 mm. It is associated with quartz, tridymite, augite, sanidine, magnesiohornblende, enstatite, pyrope, fluorapatite, hematite, braunite and roedderite. Friedrichbeckeite is light yellow, with white to light cream streak and vitreous lustre. It is brittle with irregular fracture and no cleavage, Mohs hardness of 6, calculated density is 2.686 gcm^?3. Optically, it is uniaxial positive with n_ω = 1.552(2) and n_ε = 1.561(2) at 589.3 nm and a distinct pleochroism from yellow (//ω) to light blue (//ε). Electron microprobe analyses yielded (wt.%): Na₂O 2.73, K₂O 4.16, BeO 4.67, MgO 11.24, MnO 2.05, FeO 1.76, Al₂O₃ 0.15, SiO₂ 73.51, (Σ CaO, TiO₂ = 0.06) sum 100.33 (BeO determined by LA-ICP-MS). The empirical formula based on Si = 12 is K_0.87 Na_0.86 (Mg_1.57Mn_0.28Fe_0.24)_Σ2.09 (Be_1.83?Mg_1.17)_Σ3.00 [Si₁₂O₃₀], and the simplified formula can be given as K (□_0.5Na_0.5)₂ (Mg_0.8Mn_0.1Fe_0.1)₂ (Be_0.6?Mg_0.4)₃ [Si₁₂O₃₀]. Friedrichbeckeite is hexagonal, space-group P6/mcc, with a = 9.970(1), c = 14.130(3) Å, V = 1216.4(3) ų, and Z = 2. The strongest lines in the X-ray powder diffraction pattern are (d in Å / I _obs / hkl): 3.180 / 100 / 121, 2.885 / 70 / 114, 4.993 / 30 / 110, 4.081 / 30 / 112, 3.690 / 30 / 022. A single-crystal structure refinement (R1 = 3.62 %) confirmed that the structure is isotypic with milarite and related ^[12] C ^[9] B ₂ ^[6] A ₂ ^[4] T2₃ [^[4] T1₁₂O₃₀] compounds. The C-site is dominated by potassium, the B-site is almost half occupied by sodium, and the A-site is dominated by Mg. The site-scattering at the T2-site can be refined to a Be/(Be?+?Mg) value close to 0.61; the T1-site is occupied by Si. Micro-Raman spectroscopy reveals an increasing splitting of scattering bands around 550 cm^?1 for friedrichbeckeite. The mineral can be classified as an unbranched ring silicate or as a beryllo-magnesiosilicate. With respect to the end-member formula K (□0.5Na0.5)₂ Mg₂ Be₃ [Si₁₂O₃₀] friedrichbeckeite represents the Mg-dominant analogue of almarudite, milarite or oftedalite. The mineral and its paragenesis were formed during pyrometamorphic modifications of the silicate-rich xenoliths enclosed in Quaternary leucite-tephritic lava of the Bellerberg volcano. Holotype material of friedrichbeckeite has been deposited at the mineral collection of the Naturhistorisches Museum Wien, Austria. The mineral is named friedrichbeckeite in honour of the Austrian mineralogist and petrographer Friedrich Johann Karl Becke (1855–1931).     

PREDICTION ERROR FILTERS,WHITE NOISE AND ORTHOGONAL COORDINATES*

Optimum filters can be computed using orthogonal coordinates obtained from the eigenvalues and eigenvectors of the autocorrelation matrix. The method is used to obtain unit distance prediction error filters. The output of a unit distance prediction error filter when applied to the input wavelet is an impulse at zero time. The effect on the output of added white noise is easily obtained using the approach through the orthogonal coordinates. The added white noise results in output wavelets which are no longer impulses at zero time. The decrease in time resolution gives a filter that does not increase undesirable high frequency noise as much as filters computed without white noise. Orthogonal coordinates with little signal energy can be omitted from the filter computation resulting in output wavelets resembling those computed using added white noise.     

TUTORIAL MINIMUM PHASE FOR CONTINUOUS TIME AND DISCRETE TIME FUNCTIONS*

The concept of minimum phase is clarified for geophysicists by collecting in one place the properties of minimum phase functions. The “earliest energy arrival” property in the time domain, the “minimum phase-slope property” in the frequency domain, and of the role of causal all-pass filters are demonstrated. The emphasis is placed on keeping the mathematics within the realm familiar to geophysicists and on making clear the somewhat imperfect match between physical continuous time functions and their associated discrete time representations.     

Application of kriging with external drift to estimate hydraulic conductivity from electrical-resistivity data in unconsolidated deposits near Montalto Uffugo, Italy

In the framework of various studies to characterize the aquifer at the groundwater experimental field near Montalto Uffugo, Italy, the present work estimates the spatial distribution of hydraulic conductivity of the unconsolidated deposits that underlie the area, by applying the geostatistical technique of kriging with external drift to electrical-resistivity and hydraulic-conductivity data. The reliability of the estimation method was tested by implementing a model, based on the method of cells, that simulates groundwater flow, with the estimated values of hydraulic conductivity. The results obtained indicate that the estimation method used has an acceptable degree of reliability.     

Calculating the two-dimensional magnetotelluric Jacobian in finite elements using reciprocity

To speed up the calculation of the field Jacobian for 2-D magnetoteliuric inversion using finite elements, the principle of electromagnetic reciprocity is applied. The governing relationship for the Jacobian of the field along strike is obtained by differentiating the Helmholtz equation with respect to the resistivity of each region in the finite-element mesh. The result is a similar Helmholtz equation for the Jacobian, with new sources distributed over all nodes within the parameter medium. However, according to the principle of electromagnetic reciprocity, the roles of sources and receivers are interchangeable. Utilizing reciprocity, the field values obtained from the original forward problem and for new unit sources imposed at the receivers are then utilized in the calculation of the Jacobian by simple multiplication and summation with finite-element terms at each rectangle in the mesh. For the auxiliary (across-strike) fields, the Jacobian terms are obtained by solving source vectors loaded with parabola coefficients used in the approximation to Maxwell's equations. Jacobian terms for the apparent resistivity ( p _a), the impedance phase (φ) and the vertical magnetic field ( K _zy) are then calculated utilizing the parallel- and auxiliary-field Jacobians. Comparison of Jacobian values obtained from reciprocity calculations and by differencing two forward solutions show that the reciprocity method is accurate and can be used to decrease the number of calculations required to obtain sensitivities by one to two orders of magnitude.     

On the climate response of the low-latitude Pacific Ocean to changes in the global freshwater cycle

Under global warming, the predicted intensification of the global freshwater cycle will modify the net freshwater flux at the ocean surface. Since the freshwater flux maintains ocean salinity structures, changes to the density-driven ocean circulation are likely. A modified ocean circulation could further alter the climate, potentially allowing rapid changes, as seen in the past. The relevant feedback mechanisms and timescales are poorly understood in detail, however, especially at low latitudes where the effects of salinity are relatively subtle. In an attempt to resolve some of these outstanding issues, we present an investigation of the climate response of the low-latitude Pacific region to changes in freshwater forcing. Initiated from the present-day thermohaline structure, a control run of a coupled ocean–atmosphere general circulation model is compared with a perturbation run in which the net freshwater flux is prescribed to be zero over the ocean. Such an extreme experiment helps to elucidate the general adjustment mechanisms and their timescales. The atmospheric greenhouse gas concentrations are held constant, and we restrict our attention to the adjustment of the upper 1,000 m of the Pacific Ocean between 40°N and 40°S, over 100 years. In the perturbation run, changes to the surface buoyancy, near-surface vertical mixing and mixed-layer depth are established within 1 year. Subsequently, relative to the control run, the surface of the low-latitude Pacific Ocean in the perturbation run warms by an average of 0.6°C, and the interior cools by up to 1.1°C, after a few decades. This vertical re-arrangement of the ocean heat content is shown to be achieved by a gradual shutdown of the heat flux due to isopycnal (i.e. along surfaces of constant density) mixing, the vertical component of which is downwards at low latitudes. This heat transfer depends crucially upon the existence of density-compensating temperature and salinity gradients on isopycnal surfaces. The timescale of the thermal changes in the perturbation run is therefore set by the timescale for the decay of isopycnal salinity gradients in response to the eliminated freshwater forcing, which we demonstrate to be around 10–20 years. Such isopycnal heat flux changes may play a role in the response of the low-latitude climate to a future accelerated freshwater cycle. Specifically, the mechanism appears to represent a weak negative sea surface temperature feedback, which we speculate might partially shield from view the anthropogenically-forced global warming signal at low latitudes. Furthermore, since the surface freshwater flux is shown to play a role in determining the ocean’s thermal structure, it follows that evaporation and/or precipitation biases in general circulation models are likely to cause sea surface temperature biases.     

A study of three massive star-forming regions in CS lines

We have mapped three star-forming regions (G265.14+1.45, G269.16?1.14, G291.27?0.71) in the CS(3–2) and C³⁴S(2–1) lines using the 15 m SEST telescope (Chile), and analyzed the relative positions of methanol and H₂O masers, IRAS sources, and emission maxima in the CS lines. In most cases, the maser positions are close to those of the IRAS sources. We compared the radial velocities of the maser sources and high-density CS cores, and estimated the CS column densities assuming LTE. The sizes, densities, and masses of the dense core are estimated; the masses obtained in the LTE approximation agree with the virial masses.