Gravity field determination using boundary element methods

The Boundary Element Method (BEM), a numerical technique for solving boundary integral equations, is introduced to determine the earth's gravity field. After a short survey on its main principles, we apply this method to the fixed gravimetric boundary value problem (BVP), i.e. the determination of the earth's gravitational potential from measurements of the intensity of the gravity field in points on the earth's surface. We show how to linearize this nonlinear BVP using an implicit function theorem and how to transform the linearized BVP into a boundary integral equation using the single layer representation. A Galerkin method is used to transform the boundary integral equation using the single layer representation. A Galerkin method is used to transform the boundary integral equation into a linear system of equations. We discuss the major problems of this approach for setting up and solving the linear system. The BVP is numerically solved for a bounded part of the earth's surface using a high resolution reference gravity model, measured gravity values of high density, and a 50 50 m² digital terrain model to describe the earth's surface. We obtain a gravity field resolution of 1 1 km² with an accuracy of the order 10^–3 to 10^–4 in about 1 CPU-hour on a Siemens/Fujitsu SIMD vector pipeline machine using highly sophisticated numerical integration techniques and fast equation solvers. We conclude that BEM is a powerful numerical tool for solving boundary value problems and may be an alternative to classical geodetic techniques.     

The Record of an Extreme Flood in the Sediments of Montane Lillooet Lake, British Columbia: Implications for Paleoenvironmental Assessment

Severe rainfall in mid October, 2003 produced the largest floods in almost a century of record on rivers in the Cordillera of southwestern British Columbia. Sediment deposited in Lillooet Lake as a result of this event is clearly distinguished by stratigraphy, colour, texture, magnetic properties, and organic content. Each of these physical properties is related to the lacustrine processes, especially turbid underflow, that distributed the sediment through the lake. The flood, which lasted less than a week, delivered 8–12 times the amount of sediment that accumulates in most entire years in the deepest, central parts of the lake. Recognition of events of this type in the stratigraphic record offers a means of assessing the changing nature of extreme hydroclimatic events, and their relation to more ubiquitous, lower-energy processes.     

Lake Sediment Core Records of Sulphur Accumulation and Sulphur Isotopic Composition in Central Ontario, Canada Lakes

Stable isotopic compositions and concentrations of total sedimentary sulphur (S) were determined in cores from 6 lakes in the acid-sensitive Muskoka-Haliburton region of south-central Ontario. The isotopic composition of S in deep sediment (> ~ 20 cm) was approximately constant in all lakes, and indicated a pre-industrial δ ³⁴S value between +4.0 and +5.3‰, which is similar to current bulk deposition. Similarly, total S concentrations in deep sediment were relatively low (1.9–5 mg S g⁻¹ dwt) and approximately constant with depth within cores. All lakes exhibited up-core increases in total S and decreases in δ ³⁴S at a depth corresponding to the beginning of industrialization in the Great Lakes region ( ~ 1900), resulting in a generally reciprocal depth pattern between total S concentration and δ ³⁴S ratios. While initial shifts in total S and δ ³⁴S were likely due to enhanced SO₄ reduction of newly available anthropogenic SO₄, both the magnitude and pattern of up-core S enrichment and shifts in δ ³⁴S varied greatly among lakes, and did not match changes in S deposition post 1900. Differences between lakes in total S and δ ³⁴S were not related to any single hydrologic (e.g., residence time) or physical (e.g., catchment-area-to-lake area ratio) lake characteristic. This work indicates that sediment cores do not provide consistent records of changes in post-industrial S deposition in this region, likely due to redox-related mobility of S in upper sediment.     

Phase-array decomposition of a seismic section

A seismic re fraction/wide-angle reflection profile is analysed for the presence of correlated events ('phases'). The correlation problem is formulated in terms of temporally, spatially and frequency-local complex covariances. For robustness, the method concentrates on phase rather than amplitude information. This allows a computationally efficient algorithm that can make allowance for signal correlation length and can model curved wavefronts. A statistical test based on residual phase misfit across the analysed subarray is used to assess the probability that a detected event represents a real correlated signal.
With our chosen analysis parameters and confidence level (over 99.9 per cent). 1222 events were detected in the data. Using simple techniques based on 1-D earth models, detected events are associated with a small number of particular wave types. In this way, we have succeeded in classifying almost 95 per cent of the detected events. Those that remain describe those components of the data that are inconsistent with our simple ray paths in the 1-D assumption and with our prescribed tolerance. These include reverberations, near-surface guided waves and reflected waves from strongly laterally inhomogeneous structures. According to our modelling, about 25 per cent of the detected events are consistent with simple P -wave reflected energy, and these are to a very large extent (over 85 per cent) distinct from all the other wave-type models we have used. A direct mapping of the detected events into the offset-depth domain reveals dear internal and external consistencies among the detections for the various wave types. Estimated earth structure is consistent with models from previous analyses based on much larger data sets.
We have thus succeeded in extracting correlated events from the data and decomposing these, approximately but meaningfully, into distinct classes (ray paths)     

Discrimination between climate and human-induced dryland degradation

In this study we present a technique to discriminate between climate or human-induced dryland degradation, based on evaluations of AVHRR NDVI data and rainfall data. Since dryland areas typically have high inter-annual rainfall variations and rainfall has a dominant role in determining vegetation growth, minor biomass trends imposed by human influences are difficult to verify. By performing many linear regression calculations between different periods of accumulated precipitation and the annual NDVI_max, we identify the rainfall period that is best related to the NDVI_max and by this the proportion of biomass triggered by rainfall. Positive or negative deviations in biomass from this relationship, expressed in the residuals, are interpreted as human-induced. We discuss several approaches that use either a temporally fixed NDVI peaking time or an absolute one, a best mean rainfall period for the entire drylands or the best rainfall period for each individual pixel. Advantages and disadvantages of either approach or one of its combinations for discriminating between climate and human-induced degradation are discussed. Depending on the particular land-use either method has advantages. To locate areas with a high likelihood of human-induced degradation we therefore recommend combining results from each approach.     

U/Pb and Pb/Pb zircon ages for arc-related intrusions of the Bolu Massif (W Pontides, NW Turkey): evidence for Late Precambrian (Cadomian) age

We present new U/Pb and Pb/Pb radiometric age data from two tectono-stratigraphic units of the regionally extensive Bolu Massif, in the W Pontides (İstanbul Fragment), N Turkey. A structurally lower unit (Sünnice Group) is cut by small meta-granitic intrusions, whereas the structurally higher unit comprises meta-volcanic rocks (Çaşurtepe Fm) cut by meta-granitic plutons (Tüllükiriş and Kapıkaya plutons). U/Pb single-crystal dating of zircons from the Kapıkaya Pluton yielded a concordant cluster, with a mean ²³⁸U/²⁰⁶Pb age of 565.3 ± 1.9 Ma. Zircons from the Tüllükiriş Pluton (affected by Pb loss) gave a ²⁰⁷Pb/²⁰⁶Pb age of 576 ± 6 Ma age (Late Precambrian). Small meta-granitic intrusions cutting the Sünnice Group yielded a less precise ²⁰⁷Pb/²⁰⁶Pb age of 262 ± 19 Ma (Early Permian). The older ages from the Bolu Massif confirm the existence of latest Precambrian arc magmatism related to subduction of a Cadomian ocean. We infer that the Bolu Massif represents a fragment of a Cadomian active margin. Cadomian orogenic units were dispersed as exotic terranes throughout the Variscan and Tethyan orogens, and the Bolu Massif probably reached its present position prior to latest Palaeozoic time. Our dating results also confirm that NW Turkey was affected by Hercynian magmatism related to subduction of Palaeotethys, as inferred for other areas of the Pontides.     

Hydrothermal Hgbomite Associated with Vanadiferous-Titaniferous(V-Ti) Bearing Magnetite Bands in Bhakatarhalli Chromite Mine,Nuggihalli Greenstone Belt,Western Dharwar Craton,Karnataka, India

Hgbomite, a rare exotic mineral, is found to be associated with the vanadiferous–titaniferous(V-Ti) bearing magnetite bands at Bhakatarhalli, Nuggihalli greenstone belt, western Dharwar Craton, India. We report on a second occurrence of hgbomite from the Dharwar craton in Karnataka, which is the sixth documented occurrence of this mineral from India. We evaluate the chemical characteristics of hgbomite and associated Fe-Ti-minerals in an attempt to identify its formation as a primary hydrothermal mineral in a metamorphosed magnetite layer.We report here the presence of hgbomite as a complex oxide of Fe, Mg, Al and Ti with accessory of Zn, V and Sn. Petrographic studies suggest the(V–Ti) bearing magnetite(Mt) contain spinel, hgbomite, chlorite, martite, ilmenite(Il) and minor amounts of diaspore. The hgbomite displays euhedral to subhedral textures, and is up to 250 μm along the grain boundaries of magnetite and ilmenite. In the samples studied, hgbomite is prismatic, irregular and elongated in shape. The genesis of hgbomite in veins between magnetite and ilmenite implies its precipitation from fluids without involving complicated reactions. Several models were proposed for the formation of hgbomite; however, the subject is still debatable.     

For the deep biosphere,the present is not always the key to the past: what we can learn from the geological record

Microbial life below the Earth's surface (the deep biosphere) has probably varied significantly since the Archaean. Reconstructing changes in deep biosphere activity over geological timescales is necessary to understand its role in biogeochemical cycling. Even for the last few million years, such changes are often not captured by studying the distribution of present activity. However, several studies using samples from scientific drilling have revealed mineralogical, geochemical, isotopic and fossil organic molecule imprints in the sedimentary record that document rather different past deep biosphere conditions. Changing deep biosphere conditions can also be simulated using geochemical models. While some processes occurring in the past can be understood by comparing them with the present deep biosphere, others lack any modern analogue – they are defined as non‐actualistic. A non‐actualistic consideration of the deep biosphere is therefore essential for a better understanding of how Earth and life co‐evolved through time.