Grain‐scale pressure variations and chemical equilibrium in high‐grade metamorphic rocks

In the classical view of metamorphic microstructures, fast viscous relaxation (and so constant pressure) is assumed, with diffusion being the limiting factor in equilibration. This contribution is focused on the only other possible scenario – fast diffusion and slow viscous relaxation – and brings an alternative interpretation of microstructures typical of high‐grade metamorphic rocks. In contrast to the pressure vessel mechanical model applied to pressure variation associated with coesite inclusions in various host minerals, a multi‐anvil mechanical model is proposed in which strong single crystals and weak grain boundaries can maintain pressure variation at geological time‐scales in a polycrystalline material. In such a mechanical context, exsolution lamellae in feldspar are used to show that feldspar can sustain large differential stresses (>10 kbar) at geological time‐scales. Furthermore, it is argued that the existence of grain‐scale pressure gradients combined with diffusional equilibrium may explain chemical zoning preserved in reaction rims. Assuming zero net flux across the microstructure, an equilibrium thermodynamic method is introduced for inferring pressure variation corresponding to the chemical zoning. This new barometric method is applied to plagioclase rims around kyanite in felsic granulite (Bohemian Massif, Czech Republic), yielding a grain‐scale pressure variation of 8 kbar. In this approach, kinetic factors are not invoked to account for mineral composition zoning preserved in rocks metamorphosed at high grade.     

Corundum- and kyanite-bearing anatexites from the Precambrian of tanzania

In the Precambrian near Morogoro (East Tanzania) there are anatectic gneisses which are composed of two different domains: (1) Medium grained gneiss with the assemblage (+ rutile + baddeleyite): albite (ab₉₅an₀₃or₀₂) + muscovite + phlogopite_ss + corundum, and albite (ab₈₈an₀₉or₀₃ + kyanite or sillimanite + phlogopite_ss. (2) Coarse grained nests of corundum and antiperthite with minor albite, muscovite, phlogopite, rutile, baddeleyite, and tourmaline. The bulk composition of the antiperthite id ab₇₅an₀₂or₂₃, the composition of the Na-phase ab₉₅an₀₃or₀₂, that of the K-phase ab₁₃an₀₀or₈₇. Phase relationship of the rock can be modelled in the Na-rich, Si-undersaturated part of the KNASH system. P-T conditions are close to the intersection of the melting curve muscovite_ss + H₂O = corundum + liquid with the kyanite-sillimanite transition, i.e. 7.7 kb/695°C assuming a_H2O = 1.     

Hydrologic and climatic implications of a multidisciplinary study of late Holocene sediment from Kenosee Lake, southeastern Saskatchewan, Canada

Sediment lithology and mineralogy, as well as ostracode, plant macrofossil and stable isotope stratigraphies of lake sediment cores, are used to reconstruct late Holocene hydrologic changes at Kenosee Lake, a relatively large, hyposaline lake in southeastern Saskatchewan. Chronological control is provided by AMS radiocarbon ages of upland and shoreline plant macrofossils. All indicators outline an early, low-water, saline phase of lake history (4100–3000 BP), when the basin was occupied by a series of small, interconnected, sulfate-rich brine pools, as opposed to the single, topographically-closed lake that exists today. A rapid rise in lake-level (3000–2300 BP) led to the establishment of carbonate-rich, hyposaline lake conditions like those today. Lithostratigraphic data and ostracode assemblages indicate peak salinities were attained early in this period of lake infilling, suggesting that the lake-level rise was initially driven by an influx of saline groundwater. Lake-level and water chemistry have remained relatively stable over the last 2000 years, compared to earlier events. Because of a lack of datable organic material in sediments deposited during the last 2000 years, the chronology of recent events is not well resolved. Plant macrofossil, lithostratigraphic and ostracode evidence suggests that lake draw-down, accompanied by slightly higher than present salinites, occurred sometime prior to 600 BP, followed by peak lake-level and freshwater conditions. This most recent high lake stand, indicative of a high water table on the surrounding upland, may also have led to the establishment of an extensive cover of Betula in the watershed, possibly in response to paludification. Ostracode assemblages indicate that peak freshwater conditions occurred within the last 100 years. Since historically documented lake-level fluctuations correlate with decadal scale climatic fluctuations in the meteorological record, and late-Holocene hydrologic dynamics correspond to well documented climatic excursions of the Neoglacial and Little Ice Age, Kenosee Lake dynamics offer insight into the susceptibility of the region's water resources to climate change.     

Distribution and interplay of geologic processes on Titan from Cassini radar data

The Cassini Titan Radar Mapper is providing an unprecedented view of Titan’s surface geology. Here we use Synthetic Aperture Radar (SAR) image swaths (Ta-T30) obtained from October 2004 to December 2007 to infer the geologic processes that have shaped Titan’s surface. These SAR swaths cover about 20% of the surface, at a spatial resolution ranging from ∼350 m to ∼2 km. The SAR data are distributed over a wide latitudinal and longitudinal range, enabling some conclusions to be drawn about the global distribution of processes. They reveal a geologically complex surface that has been modified by all the major geologic processes seen on Earth - volcanism, tectonism, impact cratering, and erosion and deposition by fluvial and aeolian activity. In this paper, we map geomorphological units from SAR data and analyze their areal distribution and relative ages of modification in order to infer the geologic evolution of Titan’s surface. We find that dunes and hummocky and mountainous terrains are more widespread than lakes, putative cryovolcanic features, mottled plains, and craters and crateriform structures that may be due to impact. Undifferentiated plains are the largest areal unit; their origin is uncertain. In terms of latitudinal distribution, dunes and hummocky and mountainous terrains are located mostly at low latitudes (less than 30°), with no dunes being present above 60°. Channels formed by fluvial activity are present at all latitudes, but lakes are at high latitudes only. Crateriform structures that may have been formed by impact appear to be uniformly distributed with latitude, but the well-preserved impact craters are all located at low latitudes, possibly indicating that more resurfacing has occurred at higher latitudes. Cryovolcanic features are not ubiquitous, and are mostly located between 30° and 60° north. We examine temporal relationships between units wherever possible, and conclude that aeolian and fluvial/pluvial/lacustrine processes are the most recent, while tectonic processes that led to the formation of mountains and Xanadu are likely the most ancient.     

Recent progress in estimating uncertainty in geomagnetic field modeling

Abstract

Recent work pertaining to estimating error and accuracies in geomagnetic field modeling is reviewed from a unified viewpoint and illustrated with examples. The formulation of a finite dimensional approximation to the underlying infinite dimensional problem is developed. Central to the formulation is an inner product and norm in the solution space through which a priori information can be brought to bear on the problem. Such information is crucial to estimation of the effects of higher degree fields at the Core-Mantle boundary (CMB) because the behavior of higher degree fields is masked in our measurements by the presence of the field from the Earth's crust. Contributions to the errors in predicting geophysical quantities based on the approximate model are separated into three categories: (1) the usual error from the measurement noise; (2) the error from unmodeled fields, i.e. from sources in the crust, ionosphere, etc.; and (3) the error from truncating to a finite dimensioned solution and prediction space. The combination of the first two is termed low degree error while the third is referred to as truncation error.

The error analysis problem consists of “characterizing” the difference δz = z—z, where z is some quantity depending on the magnetic field and z is the estimate of z resulting from our model. Two approaches are discussed. The method of Confidence Set Inference (CSI) seeks to find an upper bound for |z—?|. Statistical methods, i.e. Bayesian or Stochastic Estimation, seek to estimate E(δz² ), where E is the expectation value. Estimation of both the truncation error and low degree error is discussed for both approaches. Expressions are found for an upper bound for |δz| and for E(δz² ). Of particular interest is the computation of the radial field, B., at the CMB for which error estimates are made as examples of the methods. Estimated accuracies of the Gauss coefficients are given for the various methods. In general, the lowest error estimates result when the greatest amount of a priori information is available and, indeed, the estimates for truncation error are completely dependent upon the nature of the a priori information assumed. For the most conservative approach, the error in computing point values of B_r at the CMB is unbounded and one must be content with, e.g., averages over some large area. The various assumptions about a priori information are reviewed. Work is needed to extend and develop this information. In particular, information regarding the truncated fields is needed to determine if the pessimistic bounds presently available are realistic or if there is a real physical basis for lower error estimates. Characterization of crustal fields for degree greater than 50 is needed as is more rigorous characterization of the external fields.     

A summary of sulfur dioxide emission rate measuremnts from Guatemalan volcanoes

Measurements of the sulfur dioxide (SO₂) emission rate from three Guatemalan volcanoes provide data which are consistent with theoretical and laboratory studies of eruptive and shallow magma chamber processes. In particular, unerupted magma makes a major contribution to the measured SO₂ emission rates at Santiaguito, a continuously erupting dacitic volcanic dome. Varying shallow magma convection rates can explain the variations in SO₂ emission rates at Santiaguito. At Fuego, a basaltic volcano currently in repose, SO₂ emission rate measurements are consistent with a high level magma body that is crystallizing and releasing volatiles. At Pacaya, a continuously erupting basaltic volcano, recent SO₂ emission rate measurements support laboratory simulation studies of strombolian eruptions; these studies indicate that the majority of gas escapes during eruptions and little gas escapes between eruptions.Average SO₂ emission rates over the last 20 years for Santiaguito, Fuego and Pacaya are 80, 160 and 260 Mg/d, respectively. On a global scale, these three volcanoes account for 1% of the annual global volcanic output of SO₂. Santiaguito and Pacaya, together, emit 6% of the total annual SO₂ emitted by continuously erupting volcanoes.Even though SO₂ measurements at these volcanoes have been made infrequently and by different investigators, the collective data help to establish a useful baseline by which to judge future changes. A more complete record of SO₂ emission rates from these volcanoes could lead to a better understanding of their eruption mechanisms and reduce the impact of their future eruptions on Guatemalan society.     

Resolving the 180-degree ambiguity in vector magnetic field measurements: The ‘minimum’ energy solution

I present a robust algorithm that resolves the 180-deg ambiguity in measurements of the solar vector magnetic field. The technique simultaneously minimizes both the divergence of the magnetic field and the electric current density using a simulated annealing algorithm. This results in the field orientation with approximately minimum free energy. The technique is well-founded physically and is simple to implement.     

Stochastic finite fault modelling of M w 4.8 earthquake in Kachchh, Gujarat, India

The modified stochastic finite fault modelling technique based on dynamic corner frequency has been used to simulate the strong ground motions of M _w 4.8 earthquake in the Kachchh region of Gujarat, India. The accelerograms have been simulated for 14 strong motion accelerographs sites (11 sites in Kachchh and three sites in Saurashtra) where the earthquake has been recorded. The region-specific source, attenuation and generic site parameters, which are derived from recordings of small to moderate earthquakes, have been used for the simulations. The main characteristics of the simulated accelerograms, comprised of peak ground acceleration (pga), duration, Fourier and response spectra, predominant period, are in general in good agreement with those of observed ones at most of the sites. The rate of decay of simulated pga values with distance is found to be similar with that of observed values. The successful modelling of the empirical accelerograms indicates that the method can be used to prepare wide range of scenarios based on simulation which provide the information useful for evaluating and mitigating the seismic hazard in the region.