The equivalence of predictions from universal kriging and intrinsic random-function kriging

A proof is provided that the predictions obtained from kriging based on intrinsic random functions of orderk are identical to those obtained from anappropriate universal kriging model. This is a theoretical result based on known variability measures. It does not imply that people performing traditional universal kriging will get the same predictions as those using intrinsic random functions, because traditionally these methods differ in how variability is modeled. For intrinsic random functions, the same proof shows that predictions do not depend on the specific choice of the generalized covariance function. It is argued that the choice between these methods is really one of modeling and estimating the variability in the data.     

Trace-element distribution in an active hydrothermal system, Roosevelt hot springs thermal area, Utah

Chemical interaction of thermal fluids with reservoir rock in the Roosevelt Hot Springs thermal area, Utah, has resulted in the development of characteristic trace-element dispersion patterns. Multielement analyses of surface rock samples, soil samples and drill cuttings from deep exploration wells provide a three-dimensional perspective of chemical redistribution within this structurally-controlled hot-water geothermal system.Five distinctive elemental suites of chemical enrichment are recognized, each characteristic of a particular combination of physical and chemical conditions within the geothermal system. These are: (1) concentrations of As, Sb, Be, and Hg associated with siliceous material at locations of liquid discharge, fluid mixing or boiling; (2) concentrations of Mn, Ba, W, Be, Cu, Co, As, Sb and Hg in manganese and iron oxide deposits; (3) high concentrations of Hg in argillized rock near fumaroles and lower concentrations in a broad diffuse halo surrounding the thermal center; (4) concentrations of As in sulfides and Li in silicate alteration minerals immediately surrounding high-temperature fluid flow-controlling fractures; (5) deposits of CaCO₃ at depth where flashing of brine to steam has occurred due to pressure release. The geochemical enrichments are not, in general, widespread, pervasively developed zones of regular form and dimension as are typical in many ore-forming hydrothermal systems.As the geothermal system develops, changes and eventually declines through time, the chemical deposits are developed, remobilized or superimposed upon each other, thus preserving within the rocks a record of the history of the geothermal system. Recognition of trace-element distribution patterns during the exploration of a geothermal system may aid definition of the present geometry and interpretation of the history of the system.     

Global high resolution versus Limited Area Model climate change projections over Europe: quantifying confidence level from PRUDENCE results

Four high resolution atmospheric general circulation models (GCMs) have been integrated with the standard forcings of the PRUDENCE experiment: IPCC-SRES A2 radiative forcing and Hadley Centre sea surface temperature and sea-ice extent. The response over Europe, calculated as the difference between the 2071–2100 and the 1961–1990 means is compared with the same diagnostic obtained with nine Regional Climate Models (RCM) all driven by the Hadley Centre atmospheric GCM. The seasonal mean response for 2m temperature and precipitation is investigated. For temperature, GCMs and RCMs behave similarly, except that GCMs exhibit a larger spread. However, during summer, the spread of the RCMs—in particular in terms of precipitation—is larger than that of the GCMs. This indicates that the European summer climate is strongly controlled by parameterized physics and/or high-resolution processes. The temperature response is larger than the systematic error. The situation is different for precipitation. The model bias is twice as large as the climate response. The confidence in PRUDENCE results comes from the fact that the models have a similar response to the IPCC-SRES A2 forcing, whereas their systematic errors are more spread. In addition, GCM precipitation response is slightly but significantly different from that of the RCMs.     

Internal variability of regional climate models

 Two regional climate models have been applied to the task of generating an ensemble of realizations of the year 1982 with observed boundary conditions in areas covering parts of the Mediterranean countries. These realizations were generated by applying boundary conditions from the ECMWF ERA reanalysis project consecutively, carrying over the soil variables from the regional models from one iteration to the next. Monthly mean fields for six iterations of each model have been used as statistical ensembles in order to investigate the internal variability of the regional model dynamics. This internal variability is a necessary consequence of the non-linear physical feedback mechanisms of the RCM being active. A small value of internal variability will give better statistics for climate sensitivity signals, but will make these results less credible. The internal variability is important for the quantitative assessment of a climate sensitivity signal. With the present choice of models and integration domains the internal variabilities of surface fields and precipitation do reach levels that are less than, but in summer of comparable order of magnitude to, corresponding atmospheric variabilities of an atmospheric general circulation model. Received: 26 October 1999 / Accepted: 18 December 2000     

A determination of source properties of large intraplate earthquakes in Alaska

Historically, large and potentially hazardous earthquakes have occurred within the interior of Alaska. However, most have not been adequately studied using modern methods of waveform modeling. The 22 July 1937, 16 October 1947, and 7 April 1958 earthquakes are three of the largest events known to have occurred within central Alaska (M _s =7.3,M _s =7.2 andM _s =7.3, respectively). We analyzed teleseismic body waves to gain information about the focal parameters of these events. In order to deconvolve the source time functions from teleseismic records, we first attempted to improve upon the published focal mechanisms for each event. Synthetic seismograms were computed for different source parameters, using the reflectivity method. A search was completed which compared the hand-digitized data with a suite of synthetic traces covering the complete parameter space of strike, dip, and slip direction. In this way, the focal mechanism showing the maximum correlation between the observed and calculated traces was found. Source time functions, i.e., the moment release as a function of time, were then deconvolved from teleseismic records for the three historical earthquakes, using the focal mechanisms which best fit the data. From these deconvolutions, we also recovered the depth of the events and their seismic moments. The earthquakes were all found to have a shallow foci, with depths of less than 10 km.The 1937 earthquake occurred within a northeast-southwest band of seismicity termed the Salcha seismic zone (SSZ). We confirm the previously published focal mechanism, indicating strike-slip faulting, with one focal plane parallel to the SSZ which was interpreted as the fault plane. Assuming a unilateral fault model and a reasonable rupture velocity of between 2 and 3 km/s, the 21 second rupture duration for this event indicates that all of the 65 km long SSZ may have ruptured during this event. The 1947 event, located to the south of the northwest-southeast trending Fairbanks seismic zone, was found to have a duration of about 11 seconds, thus indicating a rupture length of up to 30 km. The rupture duration of the 1958 earthquake, which occurred near the town of Huslia, approximately 400 km ENE of Fairbanks, was found to be about 9 seconds. This gives a rupture length consistent with the observed damage, an area of 16 km by 64 km.     

A Eulerian air pollution model for Europe with nonlinear chemistry

A long-range transport model with nonlinear chemical reactions is described. The model contains 35 pollutants and 70 chemical reactions. This is a Eulerian model defined on a space domain containing the whole of Europe. The spherical space domain (corresponding to the Earth's surface covered by the model) is mapped into a square plane domain and discretized by using a 32×32 grid. The grid increments are equidistant (both along the Ox axis and along the Oy axis). The choice of values of the physical parameters involved in the model and the numerical treatment of the model are shortly discussed. The model is tested with meteorological data for 1985 and 1989. The numerical results are compared with measurements at stations located in different European countries. Extensive comparisons of ozone concentrations for July 1985 with measurements taken at 24 European stations are also carried out. Results concerning three episodes in July 1985 as well as results obtained in the study of the sensitivity of the ozone concentrations to variations of NO _x and/or anthropogenic VOC emissions are presented. The advantages and the limitations of such a model are discussed. The model is continuously improved by adding new modules to it. The plans for improvements in the near future are outlined.     

Eruptive history and petrologic evolution of the Albano multiple maar (Alban Hills, Central Italy)

A comprehensive volcanological study of the Albano multiple maar (Alban Hills, Italy) using (i) ⁴⁰Ar/³⁹Ar geochronology of the most complete stratigraphic section and other proximal and distal outcrops and (ii) petrographic observations, phase analyses of major and trace elements, and Sr and O isotopic analyses of the pyroclastic deposits shows that volcanic activity at Albano was strongly discontinuous, with a first eruptive cycle at 69±1 ka producing at least two eruptions, and a second cycle with two peaks at 39±1 and 36±1 ka producing at least four eruptions. Contrary to previous studies, we did not find evidence of magmatic or hydromagmatic eruptions younger than 36±1 ka. The activity of Albano was fed by a new batch of primary magma compositionally different from that of the older activity of the Alban Hills; moreover, the REE and ⁸⁷Sr/⁸⁶Sr data indicate that the Albano magma originated from an enriched metasomatized mantle. According to the modeled liquid line of descent, this magma differentiated under the influence of magma/limestone wall rock interaction. Our detailed eruptive and petrologic reconstruction of the Albano Maar evolution substantiates the dormant state of the Alban Hills Volcanic District. Electronic Supplementary Material Supplementary material is available for this article at Editorial responsibility: J. Donnelly-Nolan An erratum to this article can be found at     

An intercomparison of regional climate simulations for Europe: assessing uncertainties in model projections

Ten regional climate models (RCM) have been integrated with the standard forcings of the PRUDENCE experiment: IPCC-SRES A2 radiative forcing and Hadley Centre boundary conditions. The response over Europe, calculated as the difference between the 2071–2100 and the 1961–1990 means can be viewed as an expected value about which various uncertainties exist. Uncertainties are measured here by variance in eight sub-European boxes. Four sources of uncertainty can be evaluated with the material provided by the PRUDENCE project. Sampling uncertainty is due to the fact that the model climate is estimated as an average over a finite number of years (30). Model uncertainty is due to the fact that the models use different techniques to discretize the equations and to represent sub-grid effects. Radiative uncertainty is due to the fact that IPCC-SRES A2 is merely one hypothesis. Some RCMs have been run with another scenario of greenhouse gas concentration (IPCC-SRES B2). Boundary uncertainty is due to the fact that the regional models have been run under the constraint of the same global model. Some RCMs have been run with other boundary forcings. The contribution of the different sources varies according to the field, the region and the season, but the role of boundary forcing is generally greater than the role of the RCM, in particular for temperature. Maps of minimum expected 2m temperature and precipitation responses for the IPCC-A2 scenario show that, despite the above mentioned uncertainties, the signal from the PRUDENCE ensemble is significant.