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.     

Unusual biomarker maturation ratio changes through the oil window, a consequence of varied thermal history

Saturated and aromatic biomarker ratios continue to change systematically through the oil window and into the gas-condensate window to high vitrinite reflectances (R_o = 1.16%) in mature marine and lacustrine Mesozoic clastic samples from a South African basin. Two of the ratios reverse above R_o = 0.9%. These unusual maturation effects result from isolated periods of high rates of maturation increase. The basin cooled regionally after the break-up of Gondwana but high heating rates prevailed during the late Cretaceous-early Tertiary, as Africa moved across a hotspot, and again in the late Tertiary as a result of a possible hotspot and hydrothermal event.     

Kinetic O effects in the hydrologic cycle: Indirect evidence and implications

The abundances of ¹⁸O and deuterium in the present and past hydrologic cycle have proven to be an important tool in Earth systems science. In contrast, the abundance of ¹⁷O in precipitation has thus far been assumed to carry no additional information to that of ¹⁸O. Here, we demonstrate, using known constraints on oxygen isotope abundances from the O₂ cycle and existing data about the natural abundance of ¹⁷O in water, that the relationship between the discrimination against ¹⁷O and ¹⁸O in water may vary. This relationship, presented here as θ = ln (¹⁷α)/ln (¹⁸α), is found to be 0.511 ± 0.005 for kinetic transport effects and 0.526 ± 0.001 for equilibrium effects, with very low temperature sensitivity. As a result, the ¹⁷Δ of precipitation is controlled primarily by kinetic effects during evaporation of the initial vapor and, in contrast to the deuterium excess, is independent of the temperature at the evaporation (and condensation) site. This makes ¹⁷Δ a unique tracer that complements ¹⁸O and deuterium, and may allow for a decoupling of changes in the temperature of the ocean, that serves as the vapor source, from changes in the relative humidity above it. In addition, the ¹⁷Δ of ice caps is influenced by the kinetic effects in ice formation, and therefore measurement of ice ¹⁷Δ can be used as an additional constraint for better understanding and parameterization of these effects.     

Hydrologic testing during drilling: application of the flowing fluid electrical conductivity (FFEC) logging method to drilling of a deep borehole

Drilling of a deep borehole does not normally allow for hydrologic testing during the drilling period. It is only done when drilling experiences a large loss (or high return) of drilling fluid due to penetration of a large-transmissivity zone. The paper proposes the possibility of conducting flowing fluid electrical conductivity (FFEC) logging during the drilling period, with negligible impact on the drilling schedule, yet providing important information on depth locations of both high- and low-transmissivity zones and their hydraulic properties. The information can be used to guide downhole fluid sampling and post-drilling detailed testing of the borehole. The method has been applied to the drilling of a 2,500-m borehole at Åre, central Sweden, firstly when the drilling reached 1,600 m, and then when the drilling reached the target depth of 2,500 m. Results unveil eight hydraulically active zones from 300 m down to borehole bottom, with depths determined to within the order of a meter. Further, the first set of data allows the estimation of hydraulic transmissivity values of the six hydraulically conductive zones found from 300 to 1,600 m, which are very low and range over one order of magnitude.     

Evaluating Tidal Wetland Restoration Performance Using National Estuarine Research Reserve System Reference Sites and the Restoration Performance Index (RPI)

Evaluations of tidal wetland restoration efforts suffer from a lack of appropriate reference sites and standardized methods among projects. To help address these issues, the National Estuarine Research Reserve System (NERRS) and the NOAA Restoration Center engaged in a partnership to monitor ecological responses and evaluate 17 tidal wetland restoration projects associated with five reserves. The goals of this study were to (1) determine the level of restoration achieved at each project using the restoration performance index (RPI), which compares change in parameters over time between reference and restoration sites, (2) compare hydrologic and excavation restoration projects using the RPI, (3) identify key indicator parameters for assessing restoration effectiveness, and (4) evaluate the value of the NERRS as reference sites for local restoration projects. We found that the RPI, modified for this study, was an effective tool for evaluating relative differences in restoration performance; most projects achieved an intermediate level of restoration from 2008 to 2010, and two sites became very similar to their paired reference sites, indicating that the restoration efforts were highly effective. There were no differences in RPI scores between hydrologic and excavation restoration project types. Two abiotic parameters (marsh platform elevation and groundwater level) were significantly correlated with vegetation community structure and thus can potentially influence restoration performance. Our results highlight the value of the NERRS as reference sites for assessing tidal wetland restoration projects and provide improved guidance for scientists and restoration practitioners by highlighting the RPI as a trajectory analysis tool and identifying key monitoring parameters.     

Sedimentology of the continental end‐Permian extinction event in the Sydney Basin,eastern Australia

Upper Permian to Lower Triassic coastal plain successions of the Sydney Basin in eastern Australia have been investigated in outcrop and continuous drillcores. The purpose of the investigation is to provide an assessment of palaeoenvironmental change at high southern palaeolatitudes in a continental margin context for the late Permian (Lopingian), across the end‐Permian Extinction interval, and into the Early Triassic. These basins were affected by explosive volcanic eruptions during the late Permian and, to a much lesser extent, during the Early Triassic, allowing high‐resolution age determination on the numerous tuff horizons. Palaeobotanical and radiogenic isotope data indicate that the end‐Permian Extinction occurs at the top of the uppermost coal bed, and the Permo‐Triassic boundary either within an immediately overlying mudrock succession or within a succeeding channel sandstone body, depending on locality due to lateral variation. Late Permian depositional environments were initially (during the Wuchiapingian) shallow marine and deltaic, but coastal plain fluvial environments with extensive coal‐forming mires became progressively established during the early late Permian, reflected in numerous preserved coal seams. The fluvial style of coastal plain channel deposits varies geographically. However, apart from the loss of peat‐forming mires, no significant long‐term change in depositional style (grain size, sediment‐body architecture, or sediment dispersal direction) was noted across the end‐Permian Extinction (pinpointed by turnover of the palaeoflora). There is no evidence for immediate aridification across the boundary despite a loss of coal from these successions. Rather, the end‐Permian Extinction marks the base of a long‐term, progressive trend towards better‐drained alluvial conditions into the Early Triassic. Indeed, the floral turnover was immediately followed by a flooding event in basinal depocentres, following which fluvial systems similar to those active prior to the end‐Permian Extinction were re‐established. The age of the floral extinction is constrained to 252.54 ± 0.08 to 252.10 ± 0.06 Ma by a suite of new Chemical Abrasion Isotope Dilution Thermal Ionization Mass Spectrometry U‐Pb ages on zircon grains. Another new age indicates that the return to fluvial sedimentation similar to that before the end‐Permian Extinction occurred in the basal Triassic (prior to 251.51 ± 0.14 Ma). The character of the surface separating coal‐bearing pre‐end‐Permian Extinction from coal‐barren post‐end‐Permian Extinction strata varies across the basins. In basin‐central locations, the contact varies from disconformable, where a fluvial channel body has cut down to the level of the top coal, to conformable where the top coal is overlain by mudrocks and interbedded sandstone–siltstone facies. In basin‐marginal locations, however, the contact is a pronounced erosional disconformity with coarse‐grained alluvial facies overlying older Permian rocks. There is no evidence that the contact is everywhere a disconformity or unconformity.     

Methodological challenges to confirmatory latent variable models of social vulnerability

Natural Hazards - Socially vulnerable communities experience disproportionately negative outcomes following natural disasters and underscoring a need for well-validated measures to identify those...     

LA‐ICP‐MS Pb isotope test of meteorite provenance: A terrestrial origin for Lovina

Lovina, classified as an ungrouped ataxite, is controversial and its identity as a meteorite has been questioned. In this work, we use Pb isotopes on targeted troilite nodules in Lovina as a test of its antiquity and provenance. Although precise ages cannot be obtained, LA‐ICP‐MS offers a rapid, straightforward procedure to establish the source of lead, whether ancient (meteoritic) or modern (terrestrial). For nine pristine, unweathered nodules in Lovina, we find a lead isotopic composition of: ²⁰⁶Pb/²⁰⁸Pb = 0.492 ± 0.003 (2σ, MSWD 0.79; 95%) and ²⁰⁷Pb/²⁰⁶Pb = 0.852 ± 0.003 (2σ, MSWD 1.09; 95%) with no detectable uranium. All lead compositions of the troilite fall in the range expected for modern environmental and mantle lead and are distinctly different from the primordial Canyon Diablo Troilite (CDT) composition of ancient meteoritic troilite. Although the origin of Lovina remains unknown, we conclude that lead in the Lovina troilite is unsupported by U decay and originated from a terrestrial source.