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.     

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.     

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     

Quaternary stratigraphy of the buried valleys of the lower Red Deer River,Alberta, Canada

The Quaternary fills of the buried valleys of southern Alberta and Saskatchewan have provided a wealth of information for the reconstruction of the glacial-interglacial record of the western plains of Canada, and this paper reports on the previously unstudied stratigraphy of the buried Calgary Valley and its former tributaries in the lower Red Deer River area. We attempt to differentiate Empress Group sediments, which potentially relate to pre-glacial, interglacial/ interstadial and post-glacial lake and river deposition, using sedimentology, stratigraphy and palaeoecology. Twenty-nine stratigraphical logs indicate that Empress Group sediments have infilled a considerably large area of badlands and tributary coulees that once drained into the Calgary Valley, located 15 km to the north of Dinosaur Provincial Park. Radiocarbon dates of 52.4 ka, 27.4 ka and > 42.4 ka and glacially modified quartz grains suggest that at least some of the valley fills date to interglacial or interstadial periods and may be mid-Wisconsinan in age. However, outcrops of an older till overlying other valley fills suggest that the buried valleys were only partially excavated during interglacials/interstadials and that older (even pre-glacial) sediments could have survived. Subglacial channels, recognisable on air photographs, largely coincide with buried valley positions due to the preferential excavation of the Quaternary sediment by meltwater and are filled with post-glacial lake sediment from which a radiocarbon date of 16 ka BP was obtained. Pre-glacial and glacial/post-glacial Empress Group sediments are lithologically indistinct but cover a large time span in southern Alberta.     

Timing and tectonic setting of Grenvillian metamorphism—constraints from a transect along Georgian Bay,Ontario*

Systematic mapping of a transect along the well-exposed shores of Georgian Bay, Ontario, combined with the preliminary results of structural analysis, geochronology and metamorphic petrology, places some constraints on the geological setting of high-grade metamorphism in this part of the Central Gneiss Belt. Correlations within and between map units (gneiss associations) have allowed us to recognize five tectonic units that differ in various aspects of their lithology, metamorphic and plutonic history, and structural style. The lowest unit, which forms the footwall to a regional decollement, locally preserves relic pre-Grenvillian granulite facies assemblages reworked under amphibolite facies conditions during the Grenvillian orogeny. Tectonic units above the decollement apparently lack the early granulite facies metamorphism; out-of-sequence thrusting in the south produced a duplex-like structure. Two distinct stages of Grenvillian metamorphism are apparent. The earlier stage (c. 1160–1120 Ma) produced granulite facies assemblages in the Parry Sound domain and upper amphibolite facies assemblages in the Parry Island thrust sheet. The later stage (c. 1040–1020 Ma) involved widespread, dominantly upper amphibolite facies metamorphism within and beneath the duplex. Deformation and metamorphism recently reported from south and east of the Parry Sound domain at c. 1100–1040 Ma have not yet been documented along the Georgian Bay transect. The data suggest that early convergence was followed by a period of crustal thickening in the orogenic core south-east of the transect area, with further advance to the north-west during and after the waning stages of this deformation.     

Modelling the April 15 spectrum of Supernova 1993 J

We present the modeling of the ultraviolet and optical spectra obtained simultaneously on 1993 April 15 with the HST and at Lick Observatory. A Monte Carlo code is employed in the modeling and a comparison is made between models reported by different groups. With an atmosphere similar to the Sun in chemical composition, the observed spectral lines are well reproduced by a power law density structure of index around 20 except the strong H and HeI λ5876 lines which have peculiar absorption profiles. The photospheric velocity is found to be 9500 km/s and the blackbody temperature of the spectrum is 7990 K. For H and HeI λ5876, we suggest a two-component density structure which has a smoother layer located immediately outside the steeply decreasing inner envelope. The power law indices are most probably 20 and 3, respectively, with the transition point at about 13 000 km/s. In addition, this outer smooth layer serves to flatten the far UV spectrum as observed.