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     

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 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.     

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.     

Contribution of storms to groundwater recharge in the semi‐arid region of Karnataka,India

This paper describes the application of environmental isotopes and injected tracer techniques in estimating the contribution of storms as well as annual precipitation to groundwater recharge and its circulation, in the semi‐arid region of Bagepalli, Kolar district, Karnataka. Environmental isotopes ²H, ¹⁸O and ³H were used to study the effect of storms on the hydrological system, and an isotope balance was used to compute the contribution of a storm component to the groundwater. Some of the groundwater samples collected during the post‐storm periods were highly depleted in stable isotope content with higher deuterium excess relative to groundwater from the pre‐storm periods. Significant variation in deuterium excess in groundwater from the same area, collected in two different periods, indicates the different origin of air masses. The estimated recharge component of a storm event of 600 mm to the groundwater was found to be in the range of 117–165 mm. There was no significant variation in environmental tritium content of post‐storm and pre‐storm groundwater, indicating the fast circulation of groundwater in the system. After completion of the environmental isotope work, an injected radiotracer ³H technique was applied to estimate the direct recharge of total precipitation to the groundwater. The estimated recharge to the groundwater is 33 mm of the 550 mm annual precipitation during 1992. 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.     

Himalayan magmatism and porphyry copper–molybdenum mineralization in the Yulong ore belt, East Tibet

Summary ?The NW–SE-trending Yulong porphyry Cu–Mo ore belt, situated in the Sanjiang0 area of eastern Tibet, is approximately 400 km long and 35 to 70 km wide. Complex tectonic and magmatic processes during the Himalayan epoch have given rise to favorable conditions for porphyry-type Cu–Mo mineralization. Porphyry masses of the Himalayan epoch in the Yulong ore belt are distributed in groups along regional NW–SE striking tectonic lineaments. They were emplaced mainly into Triassic and Lower Permian sedimentary-volcanic rocks. K–Ar und U–Pb isotopic datings give an intrusion age range of 57–26 Ma. The porphyries are mainly of biotite monzogranitic and biotite syenogranitic compositions. Geological and geochemical data indicate that the various porphyritic intrusions in the belt had a common or similar magma source, are metaluminous to peraluminous, Nb–Y–Ba-depleted, I-type granitoids, and belong to the high-K calc-alkaline series. Within the Yulong subvolcanic belt a number of porphyry stocks bear typical porphyry type Cu–Mo alteration and mineralization. The most prominent porphyry Co–Mo deposits include Yulong, Malasongduo, Duoxiasongduo, Mangzong and Zhanaga, of which Yulong is one of the largest porphyry Cu (Mo) deposits in China with approximately 8 × 10⁶ tons of contained Cu metal. Hydrothermal alteration at Yulong developed around a biotite–monzogranitic porphyry stock that was emplaced within Upper Triassic limestone, siltstone and mudstone. The earliest alteration was due to the effects of contact metamorphism of the country rocks and alkali metasomatism (potassic alteration) within and around the porphyry body. The alteration of this stage was accompanied by a small amount of disseminated and veinlet Cu–Mo sulfide mineralization. Later alteration–mineralization zones form more or less concentric shells around the potassic zone, around which are distributed a phyllic or quartz–sericite–pyrite zone, a silicification and argillic zone, and a propylitic zone. Fluid inclusion data indicate that three types of fluids were involved in the alteration–mineralization processes: (1) early high temperature (660–420 °C) and high salinity (30–51 wt% NaCl equiv) fluids responsible for the potassic alteration and the earliest disseminated and/or veinlet Cu–Mo sulfide mineralization; (2) intermediate unmixed fluids corresponding to phyllic alteration and most Cu–Mo sulfide mineralization, with salinities of 30–50 wt% NaCl equiv and homogenization temperatures of 460–280 °C; and (3) late low to moderate temperature (300–160 °C) and low salinity (6–13 wt% NaCl equiv) fluids responsible for argillic and propylitic alteration. Hydrogen and oxygen isotopic studies show that the early hydrothermal fluids are of magmatic origin and were succeeded by increasing amounts of meteoric-derived convective waters. Sulfur isotopes also indicate a magmatic source for the sulfur in the early sulfide mineralization, with the increasing addition of sedimentary sulfur outward from the porphyry stock. Received August 29, 2001; revised version accepted May 1, 2002 Published online: November 29, 2002     

An Instantaneous-Frequency Filtering Technique for High-Frequency Radar Oceanography

High-frequency (HF) radar systems are remote sensing tools that can be used to measure oceanographic parameters. Problems can occur when using the conventional periodogram (PG) method for computing power spectral estimates from backscattered radar signals. Temporal and spatial inhomogeneities within the radar measurement region can cause distortion in the spectra. This paper describes an instantaneous-frequency (IF) filtering technique that has been developed to measure the first-order modulation contained within the radar signal. Successful removal of this modulation is shown to yield an increased quality and quantity of ocean measurements