Petrogenetic grid for ultrahigh-pressure metamorphism in the model system CaO-MgO-SiO2-CO2-H2O

Abstract Petrogenetic grids for ultrahigh-pressure (UHP) metamorphism were calculated at different X_co2 conditions in the model system CaO-MgO-SiO₂-CO₂-H₂O involving coesite (Co), diopside (Di), dolomite (Do), enstatite (En), forsterite (Fo), magnesite (Ms), quartz (Qz), talc (Tc), tremolite (Tr) using a published internally consistent thermodynamic data set. Two P-T grids at X_co2= 0.01 and 0.5 are described. In the calculated P-T grid at X_co2= 0.01, four out of 10 stable invariant points, Co-En-Ms-Tc, Co-Di-En-Tc-Tr, Co-Di-Ms-Tc-Tr and Di-En-Ms-Tc-Tr lie within the stability field of coesite. If the fluid phase has X_co2= 0.5, no invariant point is stable under UHP conditions. Some magnesite-bearing assemblages are stabilized by the following three reactions: Di + Ms = Do + Fo + CO₂, Ms + Tr = Do + Fo + CO₂+ H₂O and Ms + Tc = Fo+ CO₂+ H₂O at X_co2= 0.01 and by reaction Ms + Tc = Fo + CO₂+ H₂O together with these three at X_co2= 0.5. Ten possible UHP assemblages for mafic and ultramafic compositions at very low X_co2 conditions include the following: Co-Do-Ms, Co-Di-Ms, Co-Di-Tc, Di-Ms-Tc, Di-En-Tc-, Di-En-Ms, Co-Di-En, Di-En-Fo, Di-Fo-Ms, Di-Do-Fo. Among them, talc-bearing assemblages are restricted to X_co2 < 0.02 and their high-P limit is 31.7 kb (749°C) at X_co2= 0.01. Dolomite-magnesite-silica assemblages have large P-T stability fields even if X_co2 is as low as 0.1, and could occur in cold subduction zones with very low geothermal gradients. Reported UHP coesite-dolomite assemblage is restricted only to a calc-silicate rock interlayered with marble where X_co2 is relatively higher; no such assemblage appears for mafic and ultramafic rocks with low X_co2 evidenced by the occurrence of diopside (or omphacite) at the expense of dolomite + coesite. The effect of X_co2 on the stability of coesite-dolomite-magnesite, diopside-enstatite-magnesite, diopside-talc assemblages is examined and the occurrence of coesite-dolomite, magnesite-bearing and talc-bearing assemblages in the Dabie UHP rocks are interpreted by employing the calculated P-T grids.     

230Th/234U and 14C dating of a lowstand coral reef beneath the insular shelf off Irabu Island, Ryukyus, southwestern Japan

Abstract   High-resolution seismic reflection profiles delineated the distribution of mound-shaped reflections, which were interpreted as reefs, beneath the insular shelf western off Irabu Island, Ryukyus, southwestern Japan. A sediment core through one of the mounded structures was recovered from the sea floor at a depth of −118.2 m by offshore drilling and was dated by radiometric methods. The lithology and coral fauna of the core indicate that the mounded structure was composed of coral–algal boundstone suggesting a small-scaled coral reef. High-precision α-spectrometric ²³⁰Th/²³⁴U dating coupled with calibrated accelerator mass spectrometric ¹⁴C ages of corals obtained reliable ages of this reef ranging from 22.18 ± 0.63 to 30.47 ± 0.98 ka. This proves that such a submerged reef was formed during the lowstand stage of marine oxygen isotope stages 3–2. The existence of low-Mg calcite in the aragonitic coral skeleton of 22.18 ± 0.63 ka provides evidence that the reef had once been exposed by lowering of the relative sealevel to at least −126 m during the last glacial maximum in the study area. There is no room for doubt that a coral reef grew during the last glacial period on the shelf off Irabu Island of Ryukyus in the subtropical region of western Pacific.     

GRAVITY GRADIENT TENSORS DUE TO A POLYHEDRON WITH POLYGONAL FACETS1

Using the conjugate complex variables formulation, closed-form formulae for the gravity gradient tensors of the gravitational potential due to a homogeneous polyhedral body composed of polygonal facets are derived. The treatise considers the cases of the observation point being inside the polyhedron, on the surface of a facet, or outside the polyhedron.     

A Geostatistically Based Ground Water Monitoring Study of Nonpoint Source No3--N Concentrations

Geostatistical interpretations of ground water monitoring data are presented to define the spatial distributions of NO₃^--N in the ground water at two demonstration test sites in the Idaho Snake River Plain. Sequential Gaussian simulation was used to delineate monthly ground water NO₃^--N changes during and after implementation of a prescribed crop rotation at test site 1. Trend surface analyses were used to illustrate monthly ground water NO₃^--N changes during and after a prescribed irrigation practice was implemented at test site 2. These evaluations suggest that geostatistically based ground water monitoring can be effective in the delineation of changes in ground water quality in shallow, unconfined aquifers in agricultural areas such as those in southern Idaho. Geostatistical methods showed spatial and temporal changes in ground water NO₃^--N inferred to be a result of the agricultural practices implemented.     

40Ar/39Ar ages of blueschist facies pelitic schists from Qingshuigou in the Northern Qilian Mountains, western China

Abstract   Pelitic schists from Qingshuigou in the Northern Qilian Mountains of China contain mainly glaucophane, garnet, white mica, clinozoisite, chlorite and piemontite. Isotopic age dating of these schists provides new constraints on the formation of the high-grade blueschists at Qingshuigou. White mica ⁴⁰Ar/³⁹Ar ages range from 442.1 to 447.5 Ma (total fusion age of single grain) and from 445.7 to 453.9 Ma (integrated age of white mica concentrates). These ages (442.1–453.9 Ma) represent the peak metamorphic ages or cooling ages of the blueschists during exhumation shortly after peak metamorphism. The ⁴⁰Ar/³⁹Ar dates in the present study are similar to ages previously reported for eclogites and blueschists in the area; this suggests that both the eclogites and pelitic sediments underwent high-grade metamorphism during the same subduction event. From this chronological evidence and the presence of well-developed Silurian remnant-sea flysch and Devonian molasse, it is concluded that the Northern Qilian Ocean had closed by the end of the Ordovician, and rapid orogenic uplift followed in the Devonian.     

Tectonothermal evolution of the Lesser Himalaya, Nepal: Constraints from 40Ar/39Ar ages from the Kathmandu Nappe

Abstract The Himalaya is a fold-and-thrust wedge formed along the northern margin of the Indian continent, and consists of three thrust-bounded lithotectonic units; the Sub-Himalaya, the Lesser Himalaya, and the Higher Himalaya with the overlying Tethys Himalaya from south to north, respectively. The orogen-scale, intracrustal thrusts which bound the above lithotectonic units are splays off an underlying subhorizontal dkcollement, and show a southward propagating piggy-back sequence with an out-of-sequence thrust. Among these thrusts, the Main Central Thrust zone (MCT zone) has played a major role in Himalayan tectonics. The MCT zone represents a shear zone which has accommodated southward thrusting of the Higher Himalayan crystalline thrust sheet over the Lesser Himalayan sequence for ～140 km. The Kathmandu Nappe in central Nepal has been transported over the Lesser Himalayan metasediments along the MCT zone, and is locally separated from the Higher Himalayan thrust sheet in the north by an out-of-sequence thrust. ⁴⁰Ar/³⁹Ar ages have been determined for one whole-rock phyllite and six muscovite concentrates from metasedimenta-ry rocks and variably deformed granites in the Kathmandu Nappe. These ages range from 44 Ma to 14 Ma, and suggest a record of both Eo-Himalayan (Eocene) and Neo-Himalayan (Miocene) tectonothermal events in the Tertiary Himalayan orogeny. The Miocene event was associated with translation along the MCT zone. No tectonothermal event of the Late Miocene to Early Pliocene ages have been reported near the MCT zone in southern Lesser Himalayan crystalline nappe or klippe, although such events have been documented within and around the MCT zone in the northern root zone of the Higher Himalaya. This suggests that out-of-sequence thrusting may have occurred between 14 Ma and 5 Ma, probably during the period 10-7.5 Ma. Since then the frontal MCT zone below the Kathmandu Nappe has been inactive, but the MCT zone in the northern root zone has remained active. The rapid increase in denudation rates of the Higher Himalaya since the Late Miocene may have been caused by ramping along the out-of-sequence thrust at depth.     

Adaptive subtraction of multiples using the L1-norm

A strategy for multiple removal consists of estimating a model of the multiples and then adaptively subtracting this model from the data by estimating shaping filters. A possible and efficient way of computing these filters is by minimizing the difference or misfit between the input data and the filtered multiples in a least‐squares sense. Therefore, the signal is assumed to have minimum energy and to be orthogonal to the noise. Some problems arise when these conditions are not met. For instance, for strong primaries with weak multiples, we might fit the multiple model to the signal (primaries) and not to the noise (multiples). Consequently, when the signal does not exhibit minimum energy, we propose using the L₁‐norm, as opposed to the L₂‐norm, for the filter estimation step. This choice comes from the well‐known fact that the L₁‐norm is robust to ‘large’ amplitude differences when measuring data misfit. The L₁‐norm is approximated by a hybrid L₁/L₂‐norm minimized with an iteratively reweighted least‐squares (IRLS) method. The hybrid norm is obtained by applying a simple weight to the data residual. This technique is an excellent approximation to the L₁‐norm. We illustrate our method with synthetic and field data where internal multiples are attenuated. We show that the L₁‐norm leads to much improved attenuation of the multiples when the minimum energy assumption is violated. In particular, the multiple model is fitted to the multiples in the data only, while preserving the primaries.     

Inversion of normal moveout for monoclinic media1

Multiple vertical fracture sets, possibly combined with horizontal fine layering, produce an equivalent medium of monoclinic symmetry with a horizontal symmetry plane. Although monoclinic models may be rather common for fractured formations, they have hardly been used in seismic methods of fracture detection due to the large number of independent elements in the stiffness tensor. Here, we show that multicomponent wide-azimuth reflection data (combined with known vertical velocity or reflector depth) or multi-azimuth walkaway VSP surveys provide enough information to invert for all but one anisotropic parameters of monoclinic media. In order to facilitate the inversion procedure, we introduce a Thomsen-style parametrization for monoclinic media that includes the vertical velocities of the P-wave and one of the split S-waves and a set of dimensionless anisotropic coefficients. Our notation, defined for the coordinate frame associated with the polarization directions of the vertically propagating shear waves, captures the combinations of the stiffnesses responsible for the normal-moveout (NMO) ellipses of all three pure modes. The first group of the anisotropic parameters contains seven coefficients (ε^(1,2), δ^(1,2,3) and γ^(1,2)) analogous to those defined by Tsvankin for the higher-symmetry orthorhombic model. The parameters ε^(1,2), δ^(1,2) and γ^(1,2) are primarily responsible for the pure-mode NMO velocities along the coordinate axes x₁ and x₂ (i.e. in the shear-wave polarization directions). The remaining coefficient δ⁽³⁾ is not constrained by conventional-spread reflection traveltimes in a horizontal monoclinic layer. The second parameter group consists of the newly introduced coefficients ζ^(1,2,3) which control the rotation of the P-, S₁- and S₂-wave NMO ellipses with respect to the horizontal coordinate axes. Misalignment of the P-wave NMO ellipse and shear-wave polarization directions was recently observed on field data by Pérez et al. Our parameter-estimation algorithm, based on NMO equations valid for any strength of the anisotropy, is designed to obtain anisotropic parameters of monoclinic media by inverting the vertical velocities and NMO ellipses of the P-, S₁- and S₂-waves. A Dix-type representation of the NMO velocity of mode-converted waves makes it possible to replace the pure shear modes in reflection surveys with the PS₁- and PS₂-waves. Numerical tests show that our method yields stable estimates of all relevant parameters for both a single layer and a horizontally stratified monoclinic medium.     

40Ar–39Ar analysis of phlogopite in the Horoman Peridotite Complex, Hokkaido, Japan and implications for its origin

Abstract ⁴⁰Ar–³⁹Ar analysis of phlogopite separated from a plagioclase lherzolite of the Horoman Peridotite Complex, Hokkaido, Japan, has yielded a plateau age of 20.6 ± 0.5 Ma in an environment where the metamorphic fluid was characterized by an almost atmospheric Ar isotopic ratio. The age spectrum is slightly saddle-shaped, implying some incorporation of excess ⁴⁰Ar during the formation of the phlogopite at a depth. As the phlogopite has been inferred to have formed in veins and/or interstitials during exhumation of the peridotite body, metasomatic fluids, to which ground- and sea water might have contributed, were probably involved in the formation of phlogopite in the crustal environment. A total ⁴⁰Ar–³⁹Ar age of 129 Ma of a whole rock sample of the plagioclase lherzolite, from which the phlogopite was separated and is representative of the main lithology of the Horoman Peridotite Complex, indicates the occurrence of excess ⁴⁰Ar. Hence, the age has no geological meaning.     

K–Ar and 40Ar/39Ar phengite ages of Sanbagawa schist clasts from the Kuma Group, central Shikoku, southwest Japan

Conglomerates of the Kuma Group, central Shikoku, southwest Japan contain Sanbagawa schist clasts with a variety of metamorphic grades and lithologies. K–Ar and ⁴⁰Ar/³⁹Ar dating of phengite show all the pelitic schist clasts from low- to high-grade zones have similar phengite ages (82–84 Ma) that are significantly older than those from the in situ Sanbagawa sequence of central Shikoku. This is because the Kuma–Sanbagawa sequence was exhumed earlier than the in situ Asemi sequence with an exhumation process intermediate between those for the Kanto Mountains and the in situ Asemi sequences. ⁴⁰A/³⁹Ar plateau ages (103 and 117 Ma) of phengite in amphibolites indicate the timing of the early stage of the exhumation of the metamorphic pile, probably close to the peak metamorphic age.     

Vitamin B12-Catalyzed Dechlorination of Perchloroethylene Present as Residual DNAPL

The goal of this study was the cleanup of residual solvents in the saturated zone using an in situ biochemical treatment. Perchloroethylene (PCE) was chosen as a model compound because it is the most commonly found organic ground water contaminant. A mixture of vitamin B₁₂ with titanium citrate was pumped as the remedial solution through a column containing 100 μL of PCE residual. The rate of reaction was found to be first order with respect 10 the concentration of PCE and to the concentration of vitamin B₂. At 10 ppm B₁₂, more than 85 percent PCM was degraded to trichloroelhylene (TCE) and dichloroelhylene (DCE) in two hours. The presence of low to moderate concentrations of organic carbon had no significant effect on the reaction. Vitamin B₁₂ reduced by titanium citrate was found lo be compatible with the survival of anaerobic bacteria. The four major advantages of the biochemical system over the use of anaerobic bacteria are that (1) the rate is faster: (2) there is no need for the careful balance of nutrients or the addition of an extraneous carbon source: (3) there is no restriction in the concentration range of the compound to be treated; and (4) the remedial solution is mobile, even in the presence of organic carbon.     

40Ar/ 39Ar and K–Ar geochronological age constraints for the inception and early evolution of the Izu–Bonin – Mariana arc system

K–Ar and ⁴⁰Ar/³⁹Ar dates are presented for locations in the Izu–Bonin – Mariana (IBM) forearc (Ocean Drilling Program (ODP) sites 786 & 782, Chichijima, Deep Sea Drilling Program (DSDP) sites 458 & 459, Saipan), and Palau on the remnant arc of the Kyushu–Palau Ridge. For a number of these locations, the ⁴⁰Ar/³⁹Ar plateau and ³⁶Ar/⁴⁰Ar versus ³⁹Ar/⁴⁰Ar isochrons give older ages than the K–Ar results. The most important results are: (i) at site 786, initial construction of the proto-IBM (now forearc) basement occurred at least by ca 47–45 Ma, consistent with the age of the immediately overlying sediments (middle Eocene nannofossil Zone CP13c); the younger pulse of construction dated at ca 35 Ma by K–Ar could not be confirmed by ⁴⁰Ar/³⁹Ar analysis; (ii) ⁴⁰Ar/³⁹Ar ages for the initial construction of the Mariana portion of the IBM system are as old as those of the Izu–Bonin portion, for example at site 458, initial construction commenced at least by ca 49 Ma and at ca 47 Ma at Saipan (Sankakayuma Formation); and (iii) a combination of K–Ar and ⁴⁰Ar/³⁹Ar ages indicate continued boninite magmatism in the Izu–Bonin forearc (and remnant arc at Palau) until ca 35 Ma. Subduction inception including boninite series rocks along most of the exposed length of the IBM system, clearly preceded by some 5 million years the Middle Eocene (ca 43.5 Ma) change in Pacific plate motion. Boninitic series magmatism persisted at locations now exposed in the forearc for ～ 15 million years after arc inception concurrently with low-K tholeiitic series eruptions from a subaerial arc system, established at ≥ 40 Ma, on the Kyushu–Palau Ridge. For the Mariana portion of the IBM system, reconstruction of the proto-arc places this activity adjacent to the concurrent but orthogonally spreading Central Basin Ridge of the West Philippine Basin. It is possible that a combination of subduction of a young North New Guinea Plate beneath newly created back-arc basin crust may account for some of the features of the Mariana system. It is clear, however, that the understanding of the processes of subduction initiation and early IBM arc development is incomplete.