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.     

Long-term seismicity of Istanbul and of the Marmara Sea region

We examine the long-term seismicity of the Marmara Sea region in Turkey over a period of twenty centuries to attempt to answer two questions: (1) how well recorded in history is the earthquake activity of this region? and (2) does seismicity over this long period of time differ from that over the present century? We study this densely-populated and fast-developing region because it is one of the most tectonically active regions on the continents that appears today to be relatively quiescent.     

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.     

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.     

River sediments,great floods and centennial‐scale Holocene climate change

A new analysis of all 346 published ¹⁴C dated Holocene alluvial units in Britain offers a unique insight into the regional impacts of global change and shows how surprisingly sensitive British rivers have been to relatively modest but repeated changes in climate. Fourteen major but probably brief periods of flooding are identified bracketed within the periods 400–1070, 1940–3940, 7520–8100 and at ca. 10 420 cal. yr BP. There is a strong correspondence between climatic deteriorations inferred from mire wet shifts and major periods of flooding, especially at ca. 8000 cal. yr BP and since ca. 4000 cal. yr BP. The unusually long and complete British record also demonstrates that alterations in land cover have resulted in a step change in river basin sensitivity to variations in climate. This has very important implications for assessing and mitigating the impact of increasing severe flooding. In small and medium‐sized river basins land use is likely to play a key role in either moderating or amplifying the climatic signal. Copyright © 2003 John Wiley & Sons, Ltd.     

Late Quaternary vegetation changes around Lake Rutundu,Mount Kenya,East Africa: evidence from grass cuticles,pollen and stable carbon isotopes

Woody, subalpine shrubs and grasses currently surround Lake Rutundu, Mount Kenya. Multiple proxies, including carbon isotopes, pollen and grass cuticles, from a 755‐cm‐long core were used to reconstruct the vegetation over the past 38 300 calendar years. Stable carbon‐isotope ratios of total organic carbon and terrestrial biomarkers from the lake sediments imply that the proportion of terrestrial plants using the C₄ photosynthetic pathway was greater during the Late Pleistocene than in the Holocene. Pollen data show that grasses were a major constituent of the vegetation throughout the Late Pleistocene and Holocene. The proportion of grass pollen relative to the pollen from other plants was greatest at the last glacial maximum (LGM). Grass cuticles confirm evidence that C₄ grass taxa were present at the LGM and that the majority followed the cold‐tolerant NADP‐MEC₄ subpathway. Copyright © 2003 John Wiley & Sons, Ltd.     

Palaeotemperature reconstructions of the European permafrost zone during marine oxygen isotope Stage 3 compared with climate model results

A palaeotemperature reconstruction based on periglacial phenomena in Europe north of approximately 51 °N, is compared with high‐resolution regional climate model simulations of the marine oxygen isotope Stage 3 (Stage 3) palaeoclimate. The experiments represent Stage 3 warm (interstadial), Stage 3 cold (stadial) and Last Glacial Maximum climatic conditions. The palaeotemperature reconstruction deviates considerably for the Stage 3 cold climate experiments, with mismatches up to 11 °C for the mean annual air temperature and up to 15 °C for the winter temperature. However, in this reconstruction various factors linking climate and permafrost have not been taken into account. In particular a relatively thin snow cover and high climatic variability of the glacial climate could have influenced temperature limits for ice‐wedge growth. Based on modelling the 0 °C mean annual ground temperature proves to be an appropriate upper temperature limit. Using this limit, mismatches with the Stage 3 cold climate experiments have been reduced but still remain. We therefore assume that the Stage 3 ice wedges were generated during short (decadal time‐scale) intervals of extreme cold climate, below the mean temperatures indicated by the Stage 3 cold climate model simulations. Copyright © 2003 John Wiley & Sons, Ltd.     

Estimation of snow water equivalent using microwave radiometry over Arctic first‐year sea ice

The magnitude and spatial distribution of snow on sea ice are both integral components of the ocean–sea‐ice–atmosphere system. Although there exists a number of algorithms to estimate the snow water equivalent (SWE) on terrestrial surfaces, to date there is no precise method to estimate SWE on sea ice. Physical snow properties and in situ microwave radiometry at 19, 37 and 85 GHz, V and H polarization were collected for a 10‐day period over 20 first‐year sea ice sites. We present and compare the in situ physical, electrical and microwave emission properties of snow over smooth Arctic first‐year sea ice for 19 of the 20 sites sampled. Physical processes creating the observed vertical patterns in the physical and electrical properties are discussed. An algorithm is then developed from the relationship between the SWE and the brightness temperature measured at 37 GHz (55°) H polarization and the air temperature. The multiple regression between these variables is able to account for over 90% of the variability in the measured SWE. This algorithm is validated with a small in situ data set collected during the 1999 field experiment. We then compare our data against the NASA snow thickness algorithm, designed as part of the NASA Earth Enterprise Program. The results indicated a lack of agreement between the NASA algorithm and the algorithm developed here. This lack of agreement is attributed to differences in scale between the Special Sensor Microwave/Imager and surface radiometers and to differences in the Antarctic versus Arctic snow physical and electrical properties. Copyright © 2003 John Wiley & Sons, Ltd.