Comparison of methods to model the gravitational gradients from topographic data bases

A number of methods have been developed over the last few decades to model the gravitational gradients using digital elevation data. All methods are based on second-order derivatives of the Newtonian mass integral for the gravitational potential. Foremost are algorithms that divide the topographic masses into prisms or more general polyhedra and sum the corresponding gradient contributions. Other methods are designed for computational speed and make use of the fast Fourier transform (FFT), require a regular rectangular grid of data, and yield gradients on the entire grid, but only at constant altitude. We add to these the ordinary numerical integration (in horizontal coordinates) of the gradient integrals. In total we compare two prism, two FFT and two ordinary numerical integration methods using 1" elevation data in two topographic regimes (rough and moderate terrain). Prism methods depend on the type of finite elements that are generated with the elevation data; in particular, alternative triangulations can yield significant differences in the gradients (up to tens of Eötvös). The FFT methods depend on a series development of the topographic heights, requiring terms up to 14th order in rough terrain; and, one popular method has significant bias errors (e.g. 13 Eötvös in the vertical–vertical gradient) embedded in its practical realization. The straightforward numerical integrations, whether on a rectangular or triangulated grid, yield sub-Eötvös differences in the gradients when compared to the other methods (except near the edges of the integration area) and they are as efficient computationally as the finite element methods.     

Subsidence in the super-deep Pattani and Malay basins of Southeast Asia: a coupled model incorporating lower-crustal flow in response to post-rift sediment loading

Two Early Cenozoic rifts in Southeast Asia (beneath the Pattani and Malay basins) experienced only limited upper-crustal extension (β≤1.5); yet very thick post-rift sequences are present, with 6–12 km of Late Cenozoic terrestrial and shallow-marine sediment derived from adjacent sources. Conventional post-rift backstripping requires depth-dependent lithospheric thinning by β=2–4 to explain these tremendous thicknesses. We assess an alternative explanation for this post-rift subsidence, involving lower-crustal flow from beneath these basins in response to lateral pressure-gradients induced by the sediment loads and the negative loads arising from the erosion of their sediment sources. We calculate that increased rates of erosion in western Thailand in the Early Miocene placed the crust in a non-steady thermal state, such that the depth (and thus, the pressure) at the base of the brittle upper crust subsequently varied over time. Following such a perturbation, thermal and mass-flux steady-state conditions took millions of years to re-establish. In the meantime, the lateral pressure-gradient caused net outflow of lower crust, thinning the crust beneath the depocentre by several kilometres (mimicking the isostatic effect of greater crustal extension having occurred beforehand) and thickening it beneath the sediment source region. The local combination of hot crust and high rates of surface processes, causing lower-crustal flow to be particularly vigorous and thus making its effects more readily identifiable, means that the Pattani and Malay basins represent a set of conditions different from basins in many other regions. However, lower-crustal flow induced by surface processes will also occur to some extent, but less recognisably, in many other continental crustal provinces, but its effects may be mistaken for those of other processes, such as larger-magnitude stretching and/or depth-dependent stretching.     

Carbon budget for a subtropical seagrass dominated coastal lagoon: How important are seagrasses to total ecosystem net primary production?

It has been assumed that because seagrasses dominate macrophyte biomass in many estuaries they also dominate primary production. We tested this assumption by developing three carbon budgets to examine the contribution of autotrophic components to the total ecosystem net primary production (TENPP) of Lower Laguna Madre, Texas. The first budget coupled average photosynthetic parameters with average daily irradiance to calculate daily production. The second budget used average photosynthetic parameters and hourly in situ irradiance to estimate productivity. The third budget integrated temperature-adjusted photosynthetic parameters (using Q₁₀=2) and hourly in situ irradiance to estimate productivity. For each budget TENPP was calculated by integrating production from each autotroph based on the producers’ areal distribution within the entire Lower Laguna Madre. All budgets indicated that macroalgae account for 33–42% of TENPP and seagrasses consistently accounted for about 33–38%. The contribution by phytoplankton was consistently about 15–20%, and the contribution from the benthic microalgae varied between 8% and 36% of TENPP, although this may have been underestimated due to our exclusion of the within bed microphytobenthos component. The water column over the seagrass beds was net heterotrophic and consequently was a carbon sink consuming between 5% and 22% of TENPP, TENPP ranged between 5.41×10¹⁰ and 2.53×10¹¹ g C yr⁻¹, depending on which budget was used. The simplest, most idealized budget predicted the highest TENPP, while the more realistic budgets predicted lower values. Annual production rates estimated using the third budget forHalodule urightii andThalassia testudinum compare well with field data. Macroalgae and microalgae contribute 50–60% of TENPP, and seagrass may be more important as three-dimensional habitat (i.e., structure) than as a source of organic carbon to the water column in Lower Laguna Madre.     

Oxygen fugacity and geochemical variations in the martian basalts: implications for martian basalt petrogenesis and the oxidation state of the upper mantle of Mars

The oxygen fugacity of the Dar al Gani 476 martian basalt is determined to be quartz-fayalite-magnetite (QFM) −2.3 ± 0.4 through analysis of olivine, low-Ca pyroxene, and Cr-spinel and is in good agreement with revised results from Fe-Ti oxides that yield QFM −2.5 ± 0.7. This estimate falls within the range of oxygen fugacity for the other martian basalts, QFM −3 to QFM −1. Oxygen fugacity in martian basalts correlates with ⁸⁷Sr/⁸⁶Sr, ¹⁴³Nd/¹⁴⁴Nd, and La/Yb ratios, indicating that the mantle source of the basalts is reduced and that assimilation of crust-like material controls the oxygen fugacity. This allows constraints to be placed on the oxidation state of the martian mantle and on the nature of assimilated crustal material. The assimilated material may be the product of early and extensive hydrothermal alteration of the martian crust, or it may be amphibole- or phlogopite-bearing basaltic rock within the crust. In either case, water may play a significant role in the oxidation of basaltic magmas on Mars, although it may be secondary to assimilation of ferric iron-rich material.     

G-Probe-1 - An International Proficiency Test for Microprobe Laboratories - Report on Round 1 : February 2002 (TB-1 Basaltic Glass)

Results are presented for round one of a new international proficiency testing programme designed for microprobe laboratories involved in the routine analysis of silicate minerals. The sample used for this round was TB-1, a basaltic glass fused and prepared by the USGS. Thirty nine laboratories contributed data to this round, the majority of major element results being undertaken by EPMA and the majority of trace elements by LA-ICP-MS. Assigned values were derived from the median of results produced by nine selected laboratories that analysed powdered material by conventional ICP-MS, INAA and XRF techniques using bulk powders of the sample. Submitted microprobe results were evaluated using a target precision calculated using the Horwitz function, adopting the same criteria as those used for "applied" geochemistry laboratories in the companion GeoPT proficiency testing programme for laboratories involved in the routine bulk analysis of silicate rocks. An evaluation of results from participating microprobe laboratories indicated that overall, data were compatible with this precision function. A comparison between the performance of bulk and microprobe techniques used in the analysis of the basaltic glass showed remarkably good agreement, with significant bias only observed for the major oxide MgO.     

Changes in the mode of Southern Ocean circulation over the last glacial cycle revealed by foraminiferal stable isotopic variability

Benthic foraminiferal oxygen and carbon isotopic records from Southern Ocean sediment cores show that during the last glacial period, the South Atlantic sector of the deep Southern Ocean filled to roughly 2500 m with water uniformly low in δ¹³C, resulting in the appearance of a strong mid-depth nutricline similar to those observed in glacial northern oceans. Concomitantly, deep water isotopic gradients developed between the Pacific and Atlantic sectors of the Southern Ocean; the δ¹³C of benthic foraminifera in Pacific sediments remained significantly higher than those in the Atlantic during the glacial episode. These two observations help to define the extent of what has become known as the ‘Southern Ocean low δ¹³C problem’. One explanation for this glacial distribution of δ¹³C calls upon surface productivity overprints or changes in the microhabitat of benthic foraminifera to lower glacial age δ¹³C values. We show here, however, that glacial-interglacial δ¹³C shifts are similarly large everywhere in the deep South Atlantic, regardless of productivity regime or sedimentary environment. Furthermore, the degree of isotopic decoupling between the Atlantic and Pacific basins is proportional to the magnitude of δ¹³C change in the Atlantic on all time scales. Thus, we conclude that the profoundly altered distribution of δ¹³C in the glacial Southern Ocean is most likely the result of deep ocean circulation changes. While the characteristics of the Southern Ocean δ¹³C records clearly point to reduced North Atlantic Deep Water input during glacial periods, the basinal differences suggest that the mode of Southern Ocean deep water formation must have been altered as well.     

Geo–ecology and management of sensitive montane landscapes

Montane (alpine) areas are generally of high value for nature conservation. Such environments and the habitats they support are dynamic and often fragile. They are vulnerable to disturbance from a range of human activities and are responsive to climate changes over short and long timescales. Biodiversity and conservation values are closely linked to geological history, geomorphological processes and soils, and it is crucial that management systems are based on understanding these links.
There are many similarities between the Cairngorm Mountains (Scotland), the Giant Mountains (Czech Republic) and Abisko Mountains (Sweden) in terms of geology, geomorphology, ecology, links with biodiversity and high conservation importance. Comparable pressures and management issues involve, to varying degrees, a history of human use and impacts from deforestation, pasturing, grazing, recreation and atmospheric pollution. Landscape change therefore involves a complex interplay between natural and anthropogenic factors. Managing such change requires better understanding of the geo–ecological processes involved and the factors that determine landscape sensitivity. This is illustrated through a simple framework and examples from the three areas. Comparison of landscape sensitivity between similar montane areas, but in different geographic locations and climatic environments, should allow more informed management planning and a precautionary approach in advance of further changes in human activity and from predicted global warming scenarios.     

Fluvial sediment transport and deposition following the 1991 eruption of Mount Pinatubo

The 1991 eruption of Mount Pinatubo generated extreme sediment yields from watersheds heavily impacted by pyroclastic flows. Bedload sampling in the Pasig–Potrero River, one of the most heavily impacted rivers, revealed negligible critical shear stress and very high transport rates that reflected an essentially unlimited sediment supply and the enhanced mobility of particles moving over a smooth, fine-grained bed. Dimensionless bedload transport rates in the Pasig–Potrero River differed substantially from those previously reported for rivers in temperate regions for the same dimensionless shear stress, but were similar to rates identified in rivers on other volcanoes and ephemeral streams in arid environments. The similarity between volcanically disturbed and arid rivers appears to arise from the lack of an armored bed surface due to very high relative sediment supply; in arid rivers, this is attributed to a flashy hydrograph, whereas volcanically disturbed rivers lack armoring due to sustained high rates of sediment delivery. This work suggests that the increases in sediment supply accompanying massive disturbance induce morphologic and hydrologic changes that temporarily enhance transport efficiency until the watershed recovers and sediment supply is reduced.     

Flow and fracture maps for basaltic rock deformation at high temperatures

Understanding how the strength of basaltic rock varies with the extrinsic conditions of stress state, pressure and temperature, and the intrinsic rock physical properties is fundamental to understanding the dynamics of volcanic systems. In particular it is essential to understand how rock strength at high temperatures is limited by fracture. We have collated and analysed laboratory data for basaltic rocks from over 500 rock deformation experiments and plotted these on principal stress failure maps. We have fitted an empirical flow law (Norton’s law) and a theoretical fracture criterion to these data. The principal stress failure map is a graphical representation of ductile and brittle experimental data together with flow and fracture envelopes under varying strain rate, temperature and pressure. We have used these maps to re-interpret the ductile–brittle transition in basaltic rocks at high temperatures and show, conceptually, how these failure maps can be applied to volcanic systems, using lava flows as an example.     

Short‐term erosion rates from a 7Be inventory balance

Detailed soil erosion studies bene?t from the ability to quantify the magnitude of erosion over time scales appropriate to the process. An inventory balance for ⁷Be was used to calculate sediment erosion in a 30·73 m² plot during a series of runoff‐producing thunderstorms occurring over three days at the Deep Loess Research Station in Treynor, Iowa, USA. The inventory balance included determination of the pre‐ and post‐storm ⁷Be inventories in the soil, the atmospheric in?ux of ⁷Be during the event, and pro?les of the ⁷Be activity in the soil following the atmospheric deposition. The erosion calculated in the plot using the ⁷Be inventory balance was 0·058 g cm^?2, which is 23 per cent of the annual average erosion determined using ¹³⁷Cs inventories. The calculated erosion from the mass balance is similar to the 0·059 g cm^?2 of erosion estimated from the amount of sediment collected at the outlet of the 6 ha ?eld during the study period and the delivery ratio (0·64). The inventory balance of ⁷Be provides a new means for evaluating soil erosion over the time period most relevant to quantifying the prediction of erosion from runoff. Copyright © 2003 John Wiley & Sons, Ltd.