First-order asymptotics for the eigenfrequencies of the Earth and application to the retrieval of large-scale lateral variations of structure

Summary . We present a variety of examples, showing systematic fluctuations as a function of angular order of measured eigenfrequencies for given normal modes of the Earth. The data are single station measurements from the GEOSCOPE network. Such fluctuations are attributed to departures from the lowest order asymptotic expression of the geometrical optics approximation. We derive first-order asymptotic expressions for the location parameter for all three components of the Earth's motion, by a method based on the stationary phase approximation and geometric relations on the unit sphere.
We illustrate the sensitivity of the fluctuations to the different parameters involved (source parameters, epicentral distance, laterally heterogeneous earth model) with synthetic examples corresponding to GEOSCOPE observations. Finally, we show the results of first attempts at inversion, which indicate that, when the fluctuations are taken into account, more accurate estimates of the great circle average eigenfrequencies can be obtained, and additional constraints put on the structure in the neighbourhood of the great circle.     

Rapid magmatic processes accompany arc–continent collision: the Western Bismarck arc, Papua New Guinea

New U–Th–Ra, major and trace element, and Sr–Nd–Pb isotope data are presented for young lavas from the New Britain and Western Bismarck arcs in Papua New Guinea. New Britain is an oceanic arc, whereas the latter is the site of an arc–continent collision. Building on a recent study of the Manus Basin, contrasts between the two arcs are used to evaluate the processes and timescales of magma generation accompanying arc–continent collision and possible slab detachment. All three suites share many attributes characteristic of arc lavas that can be ascribed to the addition of a regionally uniform subduction component derived from the subducting altered oceanic crust and sediment followed by dynamic melting of the modified mantle. However, the Western Bismarck arc lavas diverge from the Pb isotope mixing array formed by the New Britain and the Manus Basin lavas toward elevated ²⁰⁸Pb/²⁰⁴Pb. We interpret this to reflect a second and subsequent addition of sediment melt at crustal depth during collision. ²³⁸U and ²²⁶Ra excesses are preserved in all of the lavas and are greatest in the Western Bismarck arc. High-Mg andesites with high Sr/Y ratios in the westernmost arc are attributed to recent shallow mantle flux melting at the slab edge. Data for two historical rhyolites are also presented. Although these rhyolites formed in quite different tectonic settings and display different geochemical and isotopic compositions, both formed from mafic parents within millennia.     

Impact of rock-water interactions and recharge on water resources quality of the Agadir-Essaouira basin,southwestern Morocco

The Agadir-Essaouira area in the occidental High Atlas Mountains of Morocco is characterized by a semi-arid climate. The scarcity and quality of water resources, exacerbated by long drought periods, constitute a major problem for a sustainable development of this region. Groundwater resources of carbonate units within Jurassic and Cretaceous aquifers are requested for drinking and irrigation purposes. In this study, we collected 84 samples from wells, boreholes, springs, and rivers. Hydrochemical and isotopic data were used to examine the mineralization and origin of water, which control groundwater quality. The chemical composition of water seems to be controlled by water-rock interactions, such as dissolution of carbonates (calcite and dolomite), weathering of gypsum, as well as ion exchange processes, which explain the observed variability. Stable isotopes results show that groundwater from the mainly marly Cretaceous aquifer are submitted to an evaporation effect, while samples from the chiefly calcareous Jurassic aquifer indicate a meteoric origin, due to a rapid infiltration of recharge runoff through the karstic outcrops. The low values of δ¹⁸O and δ²H suggest a local recharge from areas with elevations ranging from 400 to 1200 m for the Cretaceous aquifer and from 800 to 1500 m for the Jurassic units.     

Sampling Trace-Level Organic Solutes with Polymeric Tubing Part 2. Dynamic Studies

This is the second part of a study conducted to determine whether polymeric sampling tubing can affect organic analyte concentrations during a sampling event. In this part of the study, we looked for sorption and desorption of tricholoroethylene (TCE) and leaching of organic constituents in water pumped through five types of polymeric tubing. The materials tested were a rigid fluoropolymer, a flexible fluoropolymer, low-density polyethylene (LDPE), and two plasticized polypropylene tubings. The effects of tubing length and flow rate were examined.
The least sorptive tubings, both initially and at equilibrium, were the fluoropolymers. However, in some instances the LDPE tubing had little effect on TCE concentrations. This was when a slow flow rate was used to sample relatively shallow wells (50 feet [15 m] or less) or when a faster flow rate (1 L/min) was used to sample wells that are less than 500 feet (152 m). Further testing is recommended using more sorptive analytes.
Using high performance liquid chromatography (HPLC), we were unable to detect any constituents leaching from any of the tubings used in these studies, even when a slow flow rate was used. However, desorption of sorbed analytes is a concern for all the tubings tested, including the rigid fluoropolymer.     

Monazite U-Pb dating of staurolite grade metamorphism in pelitic schists

A study of the occurrence of and relations between rare-earth element (REE) minerals in pelitic schists indicates that monazite forms at or near the P and T of the staurolite isograd. Samples at staurolite grade from the Silurian Perry Mountain Formation in the Rumford quadrangle of Maine yield monazite in sufficient quantities to permit accurate dating of the metamorphic events forming the monazites. The bulk chemistry of the metapelites, as seen in the major element abundances and REE patterns, does not vary significantly across the study area. Thus the appearance and disappearance of REE phases is assumed to reflect changes in metamorphic grade. In a sample from the biotite zone, scanning electron microscope and microprobe studies show allanite and monazite intimately associated on a 10 m scale. The texture suggest that metastable detrital monazite breaks down, distributing its REE components to allanite. From samples below staurolite grade in which monazite is not present, our observations suggest that REEs are partitioned into allanite. At or near the staurolite isograd monazite forms as a metamorphic mineral, initiating its role as a geochronometer. Garnet-biotite geothermometry on samples at this grade from this and other studies places constraints on the minimum temperature necessary to form monazite: 525° C±25°C at 3.1±0.25 kbar. A total of 15 separates from nine schist samples ranging up to sillimanite grade have been dated. Each date is remarkably concordant, even though petrologic and textural studies by previous workers have shown that the rocks in the area have been affected by at least three metamorphic episodes. Calculations indicate insignificant Th disequilibrium in these monazites. The conditions associated with the metamorphic events suggest that monazite remains closed to lead loss provided that subsequent metamorphisms are at or below sillimanite grade. Two distinct metamorphic events are resolved, one at around 400 Ma and one at about 370 Ma. The latter was due to thermal effects of a nearby pluton that yields concordant monazite ages of 363 Ma. This work suggests that in addition to dating plutonism and high-grade metamorphism, monazite should be viewed as a reliable geochronometer for moderate metamorphism of pelitic schists.     

Calcium isotope (δCa) fractionation along hydrothermal pathways, Logatchev field (Mid-Atlantic Ridge, 14°45′N)

We investigate the Logatchev Hydrothermal Field at the Mid-Atlantic Ridge, 14°45′N to constrain the calcium isotope hydrothermal flux into the ocean. During the transformation of seawater to a hydrothermal solution, the Ca concentration of pristine seawater ([Ca]_SW) increases from about 10 mM to about 32 mM in the hydrothermal fluid endmember ([Ca]_HydEnd) and thereby adopts a δ^44/40Ca_HydEnd of −0.95 ± 0.07‰ relative to seawater (SW) and a ⁸⁷Sr/⁸⁶Sr isotope ratio of 0.7034(4). We demonstrate that δ^44/40Ca_HydEnd is higher than that of the bedrock at the Logatchev field. From mass balance calculations, we deduce a δ^44/40Ca of −1.17 ± 0.04‰ (SW) for the host-rocks in the reaction zone and −1.45 ± 0.05‰ (SW) for the isotopic composition of the entire hydrothermal cell of the Logatchev field. The values are isotopically lighter than the currently assumed δ^44/40Ca for Bulk Earth of −0.92 ± 0.18‰ (SW) [Skulan J., DePaolo D. J. and Owens T. L. (1997) Biological control of calcium isotopic abundances in the global calcium cycle. Geochim. Cosmochim. Acta61,(12) 2505-2510] and challenge previous assumptions of no Ca isotope fractionation between hydrothermal fluid and the oceanic crust [Zhu P. and Macdougall J. D. (1998) Calcium isotopes in the marine environment and the oceanic calcium cycle. Geochim. Cosmochim. Acta62,(10) 1691-1698; Schmitt A. -D., Chabeaux F. and Stille P. (2003) The calcium riverine and hydrothermal isotopic fluxes and the oceanic calcium mass balance. Earth Planet. Sci. Lett. 6731, 1-16]. Here we propose that Ca isotope fractionation along the fluid flow pathway of the Logatchev field occurs during the precipitation of anhydrite. Two anhydrite samples from the Logatchev Hydrothermal Field show an average fractionation of about Δ^44/40Ca = −0.5‰ relative to their assumed parental solutions. Ca isotope ratios in aragonites from carbonate veins from ODP drill cores indicate aragonite precipitation directly from seawater at low temperatures with an average δ^44/40Ca of −1.54 ± 0.08‰ (SW). The relatively large fractionation between the aragonite precipitates and seawater in combination with their frequent abundance in weathered mafic and ultramafic rocks suggest a reconsideration of the marine Ca isotope budget, in particular with regard to ocean crust alteration.     

Man and climate in the Maya lowlands

A 15-m sedimentary core from Lake Salpeten provides the first complete Holocene sequence for the lowlying Peten District, Guatemala. Today, Lake Salpeten is a brackish, calcium sulfate lake near saturation surrounded by tropical semievergreen forest. The basal pollen record depicts sparse juniper scrub surrounding a lake basin that held ephermal pools and halophytic marshes. The lake rapidly deepened to > 27 m in the early Holocene and may have been meromictic, because nearly 2 m of gypsum “mush” was deposited. Mesic forests were quickly established and persisted until the Maya entered the district 3000 yr ago and caused extensive deforestation. Any climatic information contained in the pollen record of the Maya period is thus masked, but a regional pollen sequence linked to the archaeological record is substantiated because environmental disturbance was pervasive. Local intensification of occupation and population growth are seen as an increased deposition of pollen of agricultural weeds and colluviation into the lake, while the Classic Maya collapse is marked by a temporary decline in Compositae pollen. Effects of perturbations induced by the Maya persist in the pollen and limnetic record 400 yr after the Spanish conquest.     

Landscape variations in stream water SO4 and δSSO4 in a boreal stream network

Despite reduced anthropogenic deposition during the last decades, deposition sulphate may still play an important role in the biogeochemical cycles of S and many catchments may act as net sources of S that may remain for several decades. The aim of this study is to elucidate the temporal and spatial dynamics of both SO₄²⁻ and δ³⁴S_SO4 in stream water from catchments with varying percentage of wetland and forest coverage and to determine their relative importance for catchment losses of S. Stream water samples were collected from 15 subcatchments ranging in size from 3 to 6780 ha, in a boreal stream network, northern Sweden. In forested catchments (<2% wetland cover) S-SO₄²⁻ concentrations in stream water averaged 1.7 mg L⁻¹ whereas in wetland dominated catchments (>30% wetland cover) the concentrations averaged 0.3 mg L⁻¹. A significant negative relationship was observed between S-SO₄²⁻ and percentage wetland coverage (r² = 0.77, p < 0.001) and the annual export of stream water SO₄²⁻ and wetland coverage (r² = 0.76, p < 0.001). The percentage forest coverage was on the other hand positively related to stream water SO₄²⁻ concentrations and the annual export of stream water SO₄²⁻ (r² = 0.77 and r² = 0.79, respectively). The annual average δ³⁴S_SO4 value in wetland dominated streams was +7.6‰ and in streams of forested catchments +6.7‰. At spring flood the δ³⁴S_SO4 values decreased in all streams by 1‰ to 5‰. The δ³⁴S_SO4 values in all streams were higher than the δ³⁴S_SO4 value of +4.7‰ in precipitation (snow). The export of S ranged from 0.5 kg S ha⁻¹ yr⁻¹ (wetland headwater stream) to 3.8 kg S ha⁻¹ yr⁻¹ (forested headwater stream). With an average S deposition in open field of 1.3 kg S ha⁻¹ yr⁻¹ (2002-2006) the mass balance results in a net export of S from all catchments, except in catchments with >30% wetland. The high temporal and spatial resolution of this study demonstrates that the reducing environments of wetlands play a key role for the biogeochemistry of S in boreal landscapes and are net sinks of S. Forested areas, on the other hand were net sources of S.     

Localization and characterization of an active fault in an urbanized area in central Guatemala by means of geoelectrical imaging

The Polochic and Motagua faults define the active plate boundary between the North American and Caribbean plates in central Guatemala. A splay of the Polochic Fault traverses the rapidly growing city of San Miguel Uspantán that is periodically affected by destructive earthquakes. This fault splay was located using a 2D electrical resistivity tomography (ERT) survey that also characterized the fault damage zone and evaluated the thickness and nature of recent deposits upon which most of the city is built. ERT images show the fault as a ~50 m wide, near-vertical low-resistivity anomaly, bounded within a few meters by high resistivity anomalies. Forward modeling reproduces the key aspects of the observed electrical resistivity data with remarkable fidelity thus defining the overall location, geometry, and internal structure of the fault zone as well as the affected lithologies. Our results indicate that the city is constructed on a ~20 m thick surficial layer consisting of poorly consolidated, highly porous, water-logged pumice. This soft layer is likely to amplify seismic waves and to liquefy upon moderate to strong ground shaking. The electrical conductivity as well as the major element chemistry of the groundwater provides evidence to suggest that the local aquifer might, at least in part, be fed by water rising along the fault. Therefore, the potential threat posed by this fault splay may not be limited to its seismic activity per se, but could be compounded its potential propensity to enhance seismic site effects by injecting water into the soft surficial sediments. The results of this study provide the basis for a rigorous analysis of seismic hazard and sustainable development of San Miguel Uspantán and illustrate the potential of ERT surveying for paleoseismic studies.