Statistical evaluation of methods for the calculation of static formation temperatures in geothermal and oil wells using an extension of the error propagation theory

The feasibility of using the more sophisticated weighted least-squares (WLS) model, as opposed to the traditional ordinary least-squares (OLS), in linear regressions of BHT data to estimate the static formation temperatures (SFT) was investigated. The most commonly used analytical methods (line-source; spherical and radial heat flow; and cylindrical heat source) were applied. Error propagation equations were derived to calculate errors in the time function of each method. These errors were combined with the BHT measurement errors to compute weighting factors for applying the WLS. Intercept uncertainties were estimated for all regressions using sets of synthetic and actual borehole logs taken from geothermal and oil applications. SFT computed with the spherical and radial heat flow method were generally greater than those from the other two methods.     

Detrital zircon reference for the North China block

U–Pb analyses of 250 single detrital zircons from Upper Proterozoic to Ordovician strata collected from the Zhuozi Shan in north-central China provide a detrital zircon reference for the North China block, a major crustal entity in the Asian tectonic collage. The results, which range in age from 1.72 to 2.97 Ga, shed new light on the age of the crystalline basement in North China, much of which is covered by younger sedimentary units. In addition, this detrital zircon reference can be used to help determine the provenance of clastic sedimentary units and for assessing validity of paleogeographic and regional tectonic models that include the complex history of Asian continental amalgamation, terrane accretion, and subsequent translation that is ongoing today.     

Spatial and temporal variability of Holocene temperature in the North Atlantic region

The early Holocene climate of the North Atlantic region was influenced by two boundary conditions that were fundamentally different from the present: the presence of the decaying Laurentide Ice Sheet (LIS) and higher than present summer solar insolation. In order to assess spatial and temporal patterns of Holocene climate evolution across this region, we collated quantitative paleotemperature records at sub-millennial resolution and synthesized their temporal variability using principal components analysis (PCA). The analysis reveals considerable spatial variability, most notably in the time-transgressive expression of the Holocene thermal maximum (HTM). Most of the region, but especially areas peripheral to the Labrador Sea and hence closest to the locus of LIS disintegration, experienced maximum Holocene temperatures that lagged peak summer insolation by 1000-3000 years. Many sites from the northeastern North Atlantic sector, including the Nordic Seas and Scandinavia, either warmed in phase with maximum summer insolation (11,000-9000 years ago) or were less strongly lagged than the Baffin Bay-Labrador Sea region. These spatially complex patterns of Holocene climate development, which are defined by the PCA, resulted from the interplay between final decay of the LIS and solar insolation forcing.     

Ammonium regeneration in intertidal sediments of the Yangtze estuary during air-exposure

The influences of exposure to the atmosphere on ammonium cycle in the intertidal surface sediments were in situ studied with a geochemical approach at a typical station in the Yangtze Estuary during three tidal cycles in September 2003. During an about 8-h emersion period of each diurnal tide, six high-resolution vertical profiles of adsorbed and dissolved ammonium were measured. It was observed that both adsorbed and dissolved ammonium generally had an increasing trend in sediment cores during the exposure. The rate of ammonium regeneration in sediments was estimated using the accumulation amount of ammonium including adsorbed and soluble fractions during the daytime emersion. The calculation result showed that there was relatively high ammonification rate (˜500 nmol N cm⁻³ day⁻¹), which reflected that organic nitrogen in sediments was quickly decomposed with a residence time of ˜52.7 days. Due to the dramatic temperature difference observed in sediment profiles, free convection was considered an important mechanism of regulating the efflux of produced ammonium into overlying waters. The total estimated amount of regenerated ammonium was ˜1.35×10⁵ t N year⁻¹ in the intertidal flat of the Yangtze Estuary, which occupied 7.6% of the total inorganic nitrogen annually transported to the estuarine ecosystem.     

Outlier detection algorithms and their performance in GOCE gravity field processing

The satellite missions CHAMP, GRACE, and GOCE mark the beginning of a new era in gravity field determination and modeling. They provide unique models of the global stationary gravity field and its variation in time. Due to inevitable measurement errors, sophisticated pre-processing steps have to be applied before further use of the satellite measurements. In the framework of the GOCE mission, this includes outlier detection, absolute calibration and validation of the SGG (satellite gravity gradiometry) measurements, and removal of temporal effects. In general, outliers are defined as observations that appear to be inconsistent with the remainder of the data set. One goal is to evaluate the effect of additive, innovative and bulk outliers on the estimates of the spherical harmonic coefficients. It can be shown that even a small number of undetected outliers (<0.2 of all data points) can have an adverse effect on the coefficient estimates. Consequently, concepts for the identification and removal of outliers have to be developed. Novel outlier detection algorithms are derived and statistical methods are presented that may be used for this purpose. The methods aim at high outlier identification rates as well as small failure rates. A combined algorithm, based on wavelets and a statistical method, shows best performance with an identification rate of about 99%. To further reduce the influence of undetected outliers, an outlier detection algorithm is implemented inside the gravity field solver (the Quick-Look Gravity Field Analysis tool was used). This results in spherical harmonic coefficient estimates that are of similar quality to those obtained without outliers in the input data.     

Physical properties of ultrahigh-pressure metamorphic rocks from the Sulu terrain, eastern central China: implications for the seismic structure at the Donghai (CCSD) drilling site

About 30 samples representing major lithologies of Sulu ultrahigh-pressure (UHP) metamorphic rocks were collected from surface exposures and exploration wells, and compressional (Vp) and shear wave (Vs) velocities and their directional dependence (anisotropy) were determined over a range of constant confining pressures up to 600 MPa and temperatures ranging from 20 to 600 °C. Samples range in composition from acidic to ultramafic. P- and S-wave velocities measured at 600 MPa vary from 5.08 to 8.64 km/s and 2.34 to 4.93 km/s, respectively. Densities are in the range from 2.60 to 3.68 g/cm³. To make a direct tie between seismic measurements (refraction and reflection) and subsurface lithologies, the experimental velocity data (corresponding to shallow depths) were used to calculate velocity profiles for the different lithologies and profiles of reflection coefficients at possible lithologic interfaces across the projected 5000-m Chinese Continental Scientific Drilling Program (CCSD) crustal segment. Comparison of calculated in situ velocities with respective intrinsic velocities suggests that the in situ velocities at shallow depths are lowered by an increased abundance of open microcracks. The strongly reflective zone beneath the Donghai drill site can be explained by the impedance contrasts between the different lithologies. Contacts between eclogite/peridotite and felsic rocks (gt-gneiss, granitic gneiss), in particular, may give rise to strong seismic reflections. In addition, shear-induced (lattice preferred orientation (LPO)-related) seismic anisotropy can increase reflectivity. For the explanation of the high velocity bodies (>6.4 km/s) around 1000 m and below 3200-m depth, large proportions of eclogite/peridotite (about 40 and 30 vol.%, respectively) are needed.     

Low-temperature thermodynamic properties of disordered zeolites of the natrolite group

 The heat capacity of paranatrolite and tetranatrolite with a disordered distribution of Al and Si atoms has been measured in the temperature range of 6–309 K using the adiabatic calorimetry technique. The composition of the samples is represented with the formula (Na_1.90K_0.22Ca_0.06)[Al_2.24Si_2.76O₁₀]·nH₂O, where n=3.10 for paranatrolite and n=2.31 for tetranatrolite. For both zeolites, thermodynamic functions (vibrational entropy, enthalpy, and free energy function) have been calculated. At T=298.15 K, the values of the heat capacity and entropy are 425.1 ± 0.8 and 419.1 ±0.8 J K⁻¹ mol⁻¹ for paranatrolite and 381.0 ± 0.7 and 383.2 ± 0.7 J K⁻¹ mol⁻¹ for tetranatrolite. Thermodynamic functions for tetranatrolite and paranatrolite with compositions corrected for the amount of extraframework cations and water molecules have also been calculated. The calculation for tetranatrolite with two water molecules and two extraframework cations per formula yields: C _p (298.15)=359.1 J K⁻¹ mol⁻¹, S(298.15) −S(0)=362.8 J K⁻¹ mol⁻¹. Comparing these values with the literature data for the (Al,Si)-ordered natrolite, we can conclude that the order in tetrahedral atoms does not affect the heat capacity. The analysis of derivatives dC/dT for natrolite, paranatrolite, and tetranatrolite has indicated that the water- cations subsystem within the highly hydrated zeolite may become unstable at temperatures above 200 K. Received: 30 July 2001 / Accepted: 15 November 2001     

Continental crust and lithospheric mantle interaction beneath North China: isotopic evidence from granulite xenoliths in Hannuoba, Sino-Korean craton

Mafic granulite and pyroxenite xenoliths from Cenozoic alkaline basalts at Hannuoba, Hebei Province, North China have been selected for a systematic geochemical and Sr–Nd–Pb isotopic study, which provides a unique opportunity to explore nature of the lower crust and the interaction between the continental crust and lithospheric mantle beneath an Archean craton. The major, compatible and incompatible elements and radiogenic isotopes of these xenoliths suggest great chemical heterogeneity of the lower crust beneath the Hannuoba region. Petrological and geochemical evidences indicate a clear cumulate origin, and most likely, they are related to basaltic underplating in different geological episodes. However, the Sr–Nd–Pb isotopic compositions of the xenoliths reveal a profound enriched source signature (EM I) with some influence of EM II, which implies that some portion of pre-existing, old metasomatized subcontinental lithospheric mantle could have played an important role in their genesis. It is suggested that the interaction between continental crust and subcontinental mantle as manifested by basaltic underplating would be closely related to regional tectonic episodes and geodynamic processes in the deep part of subcontinental lithospheric mantle.