Internally consistent gahnitic spinel-cordierite-garnet equilibria in the FMASHZn system: geothermobarometry and applications

The equilibrium (Mg, Fe, Zn)₃Al₂Si₃O₁₂+2Al₂SiO₅=3(Mg, Fe, Zn)Al₂O₄+5SiO₂ garnet + sillimanite/kyanitc = spinel + quartz was calibrated in the piston-cylinder apparatus between 11 and 30 kbar, and over the temperature range of 950 to 1200°C. Three experimental mixes of Mg no. [100^*MgO/(MgO+FeO)] 40, 47 and 60, in the FeO –MgO–Al₂O₃–SiO₂–ZnO (FMASZn) system were used under low oxygen fugacities and anhydrous conditions. We derive a ternary Fe–Mg–Zn symmetric mixing model for aluminous spinels in equilibrium with garnet, to quantify the increase in gahnitic end-member of spinel with increasing pressure and descreasing temperature. Further experiments in the spinel-cordieritequartz-sillimanite field were combined with garnet-cordierite data from the literature to produce a consistent set of equations describing the exchange reactions in FMASHZn relevant to quartz-sillimanite bearing rocks at granulite facies conditions. As spinel is an important mineral participating in many rocks of aluminous composition at granulite-facies conditions, and as zinc contributes to an enlargement of spinel's stability field towards higher pressures and lower temperatures, the thermobarometric calibrations presented here will be most significant in delineating the prograde and retrograde trajectory of P-T paths.     

William Smith's map brought alive by digital remastering

The geological mapping carried out by William Smith, which resulted in the publication of his famous map in 1815, was remarkable in many respects, not least because it relied on him being able to make consistent and accurate observations on the rock types he encountered during his fieldwork. This ability, gained from his many years studying rocks, allowed him to observe features with his own eyes (or at the very least, with the aid of a simple magnifying device) that others could not. We take a new look at William Smith's original stratigraphical sequences, and with samples collected from his classic field areas (many of which are around the city of Bath, Somerset, UK), demonstrate how spatial mineralogy mapping can be incorporated into the modern age of digital mapping.     

Avoiding numerical stability problems of long duration DGPS/INS Kalman filters

A current pursuit of the geodetic community is the optimal integration of differential GPS (DGPS) and inertial navigation system (INS) data streams for precise and efficient position and gravity vector surveying. Therein a complete INS and multiple-antenna GPS receiver payload, mounted on a moving platform, is used in conjunction with a network of ground-fixed single antenna GPS receivers. This paper presents a complete, GPS-based, external updating measurement model for the applicable Kalman filter. The model utilizes four external observation types for every GPS satellite in-view: DGPS range differences, single phase differences, and single phase-rate differences; as well as the mobile, multipleantenna GPS receiver's measurement of theerrors in the INS's estimate of the phase difference between any two vehicle-borne GPS antennae. Although not widely conveyed in the geodetic world, the inertial navigation community has long known that traditional Kalman filter covariance propagation recurrences are inherently unstable when such highly accurate external updates are repeatedly applied (every 1 second) over long time durations. A hybrid square root covariance/U — D covariance factorization approach is a numerically stable alternative and is reviewed herein. The hybrid makeup of the algorithm is necessitated by the correlated nature of the fourth type of GPS external measurement listed above (each vehicle-borne GPS antenna formstwo baselines). Such measurement correlations require a functional transformation of the overall external updating model to permit the multiple updates (simultaneously available at each updating epoch) to be sequentially (and efficiently) processed. An appropriate transformation is given. Stable covariance propagation relationships are presented and the transformed Kalman gain is also furnished and its use in the determination of the externally updated error states is discussed. Specific DGPS/INS instabilities produced by the traditional recurrences are displayed. The stable alternative method requires about 25% more CPU time than the traditional Kalman recurrences. With the ever-increasing computational speeds of microprocessors, this added CPU time is of no real concern.     

http://www.sciencedirect.com/science/article/pii/S167498711400108X

We present three 3D numerical models of deep subduction where buoyant material from an oceanic plateau and a plume interact with the overriding plate to assess the influence on subduction dynamics,trench geometry,and mechanisms for plateau accretion and continental growth.Transient instabilities of the convergent margin are produced,resulting in:contorted trench geometry;trench migration parallel with the plate margin;folding of the subducting slab and orocline development at the convergent margin;and transfer of the plateau to the overriding plate.The presence of plume material beneath the oceanic plateau causes flat subduction above the plume,resulting in a "bowed" shaped subducting slab.In plateau-only models,plateau accretion at the edge of the overriding plate results in trench migration around the edge of the plateau before subduction is re-established directly behind the trailing edge of the plateau.The plateau shortens and some plateau material subducts.The presence of buoyant plume material beneath the oceanic plateau has a profound influence on the behaviour of the convergent margin.In the plateau + plume model,plateau accretion causes rapid trench advance.Plate convergence is accommodated by shearing at the base of the plateau and shortening in the overriding plate.The trench migrates around the edge of the plateau and subduction is re-established well behind the trailing edge of the plateau,effectively embedding the plateau into the overriding plate.A slab window forms beneath the accreted plateau and plume material is transferred from the subducting plate to the overriding plate through the window.In all of the models,the subduction zone maintains a relatively stable configuration away from the buoyancy anomalies within the downgoing plate.The models provide a dynamic context for plateau and plume accretion in Phanerozoic accretionary orogenic systems such as the East China Orogen and the Central Asian Orogen(Altiads),which are characterised by accreted ophiolite complexes with diverse geochemical affinities,and a protracted evolution of accretion of exotic terranes including oceanic plateau and terranes with plume origins.     

An experiment on reservoir representation schemes to improve hydrologic prediction: coupling the national water model with the HEC-ResSim

ABSTRACT

This study experiments with reservoir representation schemes to improve the ability to model active water management in the National Water Model (NWM). For this purpose, we developed an integrated water management model, NWM-ResSim, by coupling the NWM with HEC-ResSim, and two reservoir representation schemes are tested: simulation of reservoir operations and retrieval of scheduled operations. The experiments focus on a pilot reservoir domain in the Russian River basin – Lake Mendocino, California – and its contributing watershed. The evaluation results suggest that the NWM-ResSim improves the simulation performance of reservoir outflow from this managed reservoir over the NWM default level pool routing scheme. The degree of this improvement depends on the suitability of the operation guidance; the reservoir operations simulation scheme could have acceptable errors for the purposes of water resources management, but not for flood operations. Results of the retrieval scheme of scheduled operations demonstrated better performance for sub-daily flood operations.     

Helium isotope geochemistry of mid-ocean ridge basalts from the South Atlantic

We report new helium isotope results for 49 basalt glass samples from the Mid-Atlantic Ridge between 1°N and 47°S.³He/⁴He in South Atlantic mid-ocean ridge basalts (MORB) varies between 6.5 and 9.0 R_A (R_A is the atmospheric ratio of1.39 × 10⁻⁶), encompassing the range of previously reported values for MORB erupted away from high³He/⁴He hotspots such as Iceland. He, Sr and Pb isotopes show systematic relationships along the ridge axis. The ridge axis is segmented with respect to geochemical variations, and local spike-like anomalies in³He/⁴He, Pb and Sr isotopes, and trace element ratios such as(La/Sm)_N are prevalent at the latitudes of the islands of St. Helena, Tristan da Cunha and Gough to the east of the ridge. The isotope systematics are consistent with injection beneath the ridge of mantle “blobs” enriched in radiogenic He, Pb and Sr, derived from off-axis hotspot sources. The variability in³He/⁴He along the ridge can be used to refine the hotspot source-migrating-ridge sink model.

MORB from the 2–7°S segment are systematically the least radiogenic samples found along the mid-ocean ridge system to date. Here the depleted mantle source is characterized by⁸⁷Sr/⁸⁶Sr of 0.7022, Pb isotopes close to the geochron and with²⁰⁶Pb/²⁰⁴Pb of 17.7, and³He/⁴He of 8.6–8.9 R_A. The “background contamination” of the subridge mantle, by radiogenic helium derived from off-ridge hotspots, displays a maximum between 20 and 24°S. The HePb and HeSr isotope relations along the ridge indicate that the³He/⁴He ratios are lower for the hotspot sources of St. Helena, Tristan da Cunha and Gough than for the MORB source, consistent with direct measurements of³He/⁴He ratios in the island lavas. Details of the HeSrPb isotope systematics between 12 and 22°S are consistent with early, widespread dispersion of the St. Helena plume into the asthenosphere, probably during flattening of the plume head beneath the thick lithosphere prior to continental breakup. The geographical variation in theHe/Pbratio deduced from the isotope systematics suggests only minor degassing of the plume during this stage. Subsequently, it appears that the plume component reaching the mid-Atlantic ridge was partially outgassed of He during off-ridge hotspot volcanism and related melting activity.

Overall, the similar behavior of He and Pb isotopes along the ridge indicates that the respective mantle sources have evolved under conditions which produced related He and Pb isotope variations.     

Turbulent kinetic energy budget in a gravel-bed channel flow

The present experimental investigation focuses on the characteristics of near bed turbulence in a fully rough, uniform open-channel flow over a gravel-type bed. Due to bed topography small scale heterogeneity, the flow is not uniform locally in the near bed region and a double averaging methodology is applied over a length scale much larger than the gravel size. The double-averaged Turbulent Kinetic Energy (TKE) budget derived in the context of the present flow over a gravel bed differs from the TKE budget written for flow over a vegetation canopy. The non-constant shape of the roughness function measured in our gravel bed leads to an additional bed-induced production term which is null for vertical roughness elements, such as simplified vegetation elements. The experimental estimation of the terms of the TKE budget reveals that the maximum turbulent activity takes place away from the reference plane, near the roughness crests. However, within the interface sublayer the work of the bed induced velocity fluctuations against the Reynolds stress is of the same magnitude as the main turbulence production term. Consequently, the characteristics of the TKE budget have similarities with uniform flows over canopies and strongly differ from uniform flows over smooth and transitionally rough flows over sedimentlike beds.     

Contraction of Radiating Shell in Classical and Post-NewtonianApproximation

A semi-numeric method was used to study the evolution of a general-relativistic radiating shell in its Post-Newtonian and Newtonian approximations. The solutions, where the main parameters are given reasonable values, show that a relativistic effect can be found in X-ray-burster-like objects. This revised version was published online in July 2006 with corrections to the Cover Date.     

Microwave brightness of Saturn's rings

It is shown that a lower limit exists on the microwave brightness of the rings of Saturn, if they are assumed to be composed of Mie scatterers of geological composition. The lower limit (about 15°K) is due to scattering of planetary microwave emission. Significant variation of brightness with azimuth along the rings is expected if the particles are typically of 2–3cm radius. Implications for the multiple-scattering hypothesis of the radar cross section of the rings are noted.