Cokriging nonstationary data

Universal cokriging is used to obtain predictions when dealing with multivariate random functions. An important type of nonstationarity is defined in terms of multivariate random functions with increments which are stationary of orderk. The covariance between increments of different variables is modeled by means of the pseudo-cross-covariance function. Criteria are formulated to which the parameters of pseudo-cross-covariance functions must comply so as to ensure positive-definiteness. Cokriging equations and the induced cokriging equations are given. The study is illustrated by an example from soil science.     

The interannual leading modes of the extratropical variability in the Southern Hemisphere simulated by the ECHAM-4 atmospheric model

 An ensemble of twenty-three 14-year experiments conducted with the ECHAM-4 GCM has been examined to test the model's capability to simulate the principal modes of interannual variability. The integrations were performed under specified monthly SST between 1979–1993. The analysis was focused on the Southern Hemisphere (SH) extratropics. Empirical orthogonal functions analysis (EOF) using seasonal anomaly fields has been performed to isolate the principal modes that dominate the southern extratropical variability at the interannual time scale. Leading patterns of 500 hPa geopotential height (z500) have been compared with those estimated from the ECMWF re-analysis dataset. The model is able to adequately reproduce the spatial pattern of the annular mode, but it represents the temporal variations of the oscillation less satisfactorily. The model simulation of the Pacific South American (PSA) pattern is better, both in the shape of the pattern and in the temporal evolution. To verify if the capability of the model to adequately simulate the temporal oscillation of the propagating patterns is related to the increased influence of the tropical external forcing, covarying SST-atmospheric modes have been identified by singular value decomposition (SVD). In winter (July-August-September, JAS) the tropical SST variability is highly correlated with the ENSO mode. In summer (January-February-March, JFM) the strength of the teleconnections is related to strong westerly anomalies, disrupted by a meridional out of phase relation near to South America. The large size of the ensemble was exploited by comparing the time-varying model spread and degrees of freedom of the simulated extratropical circulation. Results show that when the extratropical circulation has a few degrees of freedom, the reproducibility is relatively low and the ensemble is governed by a fairly robust zonally symmetric structure of dispersion. Received: 9 May 2000 / Accepted: 30 January 2001     

Baryte and calcite cements in the 'breccias' of Ogmore-by-Sea, South Wales

Each year many hundreds, if not thousands, of Earth-science students visit the coastline at Ogmore-by-Sea, mid-Glamorgan. They study the excellent exposures of Lower Carboniferous Limestone, with their well-preserved coral and other shelly faunas, together with the presumed late Triassic 'breccias'* that rest unconformably on top. A study of the cements within the breccias may give an insight into syndepositional mineralization and the formation of sedimentary basins linked to the opening of the Atlantic in early Mesozoic times.     

Recharge and salinization in the Madrid Basin

The purpose of this work was to reinvestigate the existing hydrogeological conceptual model of the basin of Madrid, Spain. A cumulative chemical isotopic diagram which enabled the distinction between different groups of water as well as calculation of the mode of their blending was used for this investigation. It was found that the groups of discharge were lighter in their isotopic composition than that of recharge. The previous explanation of this fact, backed by carbon-14 dating, was the long residence time due to flow lines going down to depths of more than 1000 m. This flow model assumes homogenous conditions to these depths. This assumption can not be supported by evidence from deep wells. Thus a modified model is suggested which maintains homogenous conditions only to about 300 m and a deep confined aquifer below containing paleowater. The higher degree of depletion of this water has been explained by a colder climate on top of an altitude effect. Another interesting observation was the correlation between the isotopic composition of the rains, the month of the rain event and the composition of the recharge group groundwater. It could be seen that the winter rains resemble the groundwater composition, which shows that practically all the spring and summer rains were evapotranspirated.     

THE EUROPEAN MACROSEISMIC SCALE (MSK-92)

A new updated version of the MSK macroseismic intensity scale has been prepared by a Working Group of the European Seismological Commission and has been published in April 1993 (European Macroseismic Scale 1992: updated MSK scale, 1993, ed. by G. Grünthal, Cahiers du Centre Européen de Geodynamique et de Séismologie, no. 7).     

The stable isotope signature of kilometre-scale fracturedominated metamorphic fluid pathways, Mary Kathleen, Australia

Abstract Large calcite veins and pods in the Proterozoic Corella Formation of the Mount Isa Inlier provide evidence for kilometre-scale fluid transport during amphibolite facies metamorphism. These 10- to 100-m-scale podiform veins and their surrounding alteration zones have similar oxygen and carbon isotopic ratios throughout the 200 × 10-km Mary Kathleen Fold Belt, despite the isotopic heterogeneity of the surrounding wallrocks. The fluids that formed the pods and veins were not in isotopic equilibrium with the immediately adjacent rocks. The pods have δ¹³C_calcite values of –2 to –7% and δ¹⁸O_calcite values of 10.5 to 12.5%. Away from the pods, metadolerite wallrocks have δ¹⁸O_whole-rock values of 3.5 to 7%. and unaltered banded calc-silicate and marble wallrocks have δ¹³C_calcite of –1.6 to –0.6%, and δ¹⁸O_calcite of 18 to 21%. In the alteration zones adjacent to the pods, the δ¹⁸O values of both metadolerite and calc-silicate rocks approach those of the pods. Large calcite pods hosted entirely in calc-silicates show little difference in isotopic composition from pods hosted entirely in metadolerite. Thus, 100- to 500-m-scale isotopic exchange with the surrounding metadolerites and calc-silicates does not explain the observation that the δ¹⁸O values of the pods are intermediate between these two rock types. Pods hosted in felsic metavolcanics and metasiltstones are also isotopically indistinguishable from those hosted in the dominant metadolerites and calc-silicates. These data suggest the veins are the product of infiltration of isotopically homogeneous fluids that were not derived from within the Corella Formation at the presently exposed crustal level, although some of the spread in the data may be due to a relatively small contribution from devolatilization reactions in the calc-silicates, or thermal fluctuations attending deformation and metamorphism. The overall L-shaped trend of the data on plots of δ¹³C vs. δ¹⁸O is most consistent with mixing of large volumes of externally derived fluids with small volumes of locally derived fluid produced by devolatilization of calc-silicate rocks. Localization of the vein systems in dilatant sites around metadolerite/calc-silicate boundaries indicates a strong structural control on fluid flow, and the stable isotope data suggest fluid migration must have occurred at scales greater than at least 1 km. The ultimate source for the external fluid is uncertain, but is probably fluid released from crystallizing melts derived from the lower crust or upper mantle. Intrusion of magmas below the exposed crustal level would also explain the high geothermal gradient calculated for the regional metamorphism.     

Monsoon-induced sea surface temperature changes recorded in Indian corals

The South Asian Summer Monsoon induces vertical mixing in the surface Arabian Sea, leading to a reduction in the Sea Surface Temperature (SST) of the order of 3–4°C. This reduction in temperature is recorded by modern corals (Porites) that grow in the Lakshadweep Islands (coralline islands located at about 350 km off the south-west coast of India) in their stable oxygen isotope ratios (denoted by δ¹⁸O). As large coral colonies of this genus are available, our results show that palaeomonsoon records for a few centuries back in time, a crucial input for climatic models aimed at predicting the Asian Monsoon, can be obtained from these corals. We also show that two corals separated by ? 60 km show similar (518O variations as does a giant clam (Tridacna maximus) that grew near one end of the coral colonies. As this clam is known to precipitate CaCO3 in isotopic equilibrium with the ambient sea-water, it is possible to estimate the isotopic offset of coralline δ¹⁸O from that of the CaCO3 precipitated in isotopic equilibrium. This ‘disequilibrium effect’ appears to be constant around 4.5±0.2%0. Our calculations show that SST (t, °C) is related to the coral δ¹⁸O (δ_c) and the sea-water δ¹⁸O (δ_w) by the equation t = 3.0–4.68 (δ_c - δ_w), which is in good agreement with such relationships for corals from the Pacific and Atlantic.     

The middle pleistocene deposits of Waverley Wood Pit,Warwickshire, England

A complex of channels underlying the Baginton-Lillington Gravel (Baginton Formation) at Waverley Wood Quarry, Warwickshire is described. Fossil pollen and plant macrofossils, Coleoptera, Ostracoda, Mollusca and Mammalia are described from the channel-fill deposits. Consideration of all the evidence allows the identification of four separate stages of channel fill which largely occurred under a cool temperate climate. At the top of Channel 2 evidence for a cold, continental climatic episode can be recognised, suggesting that the whole complex was deposited under a fluctuating climate at the end of a temperate stage. At two levels in the channels human artefacts were recovered confirming the presence of Palaeolithic people in Warwickshire during the deposition of the sediments. Amino-acid geochronology suggests an age within the ‘Cromerian Complex’ Stage for the channels. The small vertebrate and molluscan faunas indicate that the deposits are no older than the latter part of the ‘Cromerian Complex’ Stage of East Anglia. The regional stratigraphic significance of the Waverley Wood succession is outlined.     

Fluid infiltration during contact metamorphism of interbedded marble and calc-silicate hornfels, Twin Lakes area, central Sierra Nevada, California

Abstract In the Twin Lakes area, central Sierra Nevada, California, most contact metamorphosed marbles contain calcite + dolomite + forsterite ± diopside ± phlogopite ± tremolite, and most calc-silicate hornfelses contain calcite + diopside + wollastonite + quartz ± anorthite ± K-feldspar ± grossular ± titanite. Mineral-fluid equilibria involving calcite + dolomite + tremolite + diopside + forsterite in two marble samples and wollastonite + anorthite + quartz + grossular in three hornfels samples record P± 3 kbar and T± 630° C. Various isobaric univariant assemblages record CO₂-H₂O fluid compositions of χCO₂= 0.61–0.74 in the marbles and χCO₂= 0.11 in the hornfelses. Assuming a siliceous dolomitic limestone protolith consisting of dolomite + quartz ° Calcite ± K-feldspar ± muscovite ± rutile, all plausible prograde reaction pathways were deduced for marble and hornfels on isobaric T-XCO₂ diagrams in the model system K₂O-CaO-MgO-Al₂O₃-SiO₂-H₂O-CO₂. Progress of the prograde reactions was estimated from measured modes and mass-balance calculations. Time-integrated fluxes of reactive fluid which infiltrated samples were computed for a temperature gradient of 150 °C/km along the fluid flow path, calculated fluid compositions, and estimated reaction progress using the mass-continuity equation. Marbles and hornfelses record values in the range 0.1–3.6 × 10⁴ cm³/cm² and 4.8–12.9 × 10⁴ cm³/cm², respectively. For an estimated duration of metamorphism of 10⁵ years, average in situ metamorphic rock permeabilities, calculated from Darcy's Law, are 0.1–8 × 10^?6 D in the marbles and 10–27 × 10^?6 D in the hornfelses. Reactive metamorphic fluids flowed up-temperature, and were preferentially channellized in hornfelses relative to the marbles. These results appear to give a general characterization of hydrothermal activity during contact metamorphism of small pendants and screens (dimensions ± 1 km or less) associated with emplacement of the Sierra Nevada batholith.     

Deglaciation of the younger dryas or Loch Lomond Stadial ice-field in the northern Highlands,Scotland

When viewed from the air, Scottish ‘hummocky moraine’ can be resolved into a series of linear ridges that resemble those found at the margins of actively retreating glaciers today. Recent work has supported the interpretation of these linear ridges as ice-marginal landforms and the authors believe that the majority of ‘hummocky moraine’ deposits can be interpreted in this way. Consequently the pattern of deglaciation can be established fairly precisely from the pattern of linear ridges. This approach is applied to the landforms of the northern part of the Loch Lomond Stadial ice-field in order to reconstruct the regional pattern of deglaciation. This leads to important inferences about the significance of topographic control during deglaciation and more importantly it provides fresh insight into the environment of the British Isles during the Loch Lomond Stadial.