The use of secondary information in geostatistical target area identification

The identification and characterization of target areas at former bombing ranges is the first step in investigating these sites for residual unexploded ordnance. Traditionally, magnetometer surveys along transects are used in identifying areas with high densities of magnetic anomalies, which are likely former target areas. Combining magnetometer survey data with other data sources may reduce the level of survey data required for site characterization, increasing characterization efficiency. Here, several techniques for incorporating secondary information into kriging estimates of magnetic anomaly density are investigated for a former bombing range located near Pueblo, Colorado. In particular, kriging with external drift, collocated ordinary cokriging, and simple kriging with local means (SKLM) are used to incorporate information from a secondary variable. The secondary variable consists of a grid of crater density values derived from a topographic light detection and ranging (LIDAR) analysis. The craters, which are clearly identifiable in the LIDAR data, were generated through munitions use at the site and are therefore related to the target locations. The results from this study indicate that the inclusion of the secondary information in the kriging estimates does benefit target area characterization and provides a means of elucidating target area details from only limited magnetometer transect data. For the Pueblo site, the use of SKLM with the crater density as a secondary variable and only limited magnetometer transect data, provided results comparable to those obtained from using much larger magnetometer transect data sets.     

Two-phase hydrodynamic and sediment transport modeling of wave-generated sheet flow

This numerical investigation was carried out to advance mechanistic understanding of sediment transport under sheet flow conditions. An Euler–Euler coupled two-phase flow model was developed to simulate fluid–sediment oscillatory sheet flow. Since the concentration of sediment particles is high in such flows, the kinematics of the fluid and sediment phases are strongly coupled. This model includes interaction forces, intergranular stresses and turbulent stress closure. Each phase was modeled via the Reynolds-Averaged Navier–Stokes equations, with interphase momentum conservation accounting for the interaction between the phases. The generation and transformation of turbulence was modeled using the two-equation k–ε$k''–''\epsilon$ turbulence model. Concentration and sediment flux profiles were compared with experimental data for sheet flow conditions considering both symmetric and asymmetric oscillatory flows. Sediment and fluid velocity variations, concentration profiles, sediment flux and turbulence parameters of wave-generated sheet flow were studied numerically with a focus on sediment transport characteristics. In all applications, the model predictions compared well with the experimental data. Unlike previous investigations in which the flow is driven by a horizontal pressure gradient, the present model solves the Navier–Stokes equations under propagating waves. The model’s ability to predict sediment transport under oscillatory sheet flow conditions underscores its potential for understanding the evolution of beach morphology.     

A study of the geoelectrical properties of peatlands and their influence on ground-penetrating radar surveying1

Geophysical surveys and chemical analyses on cores were carried out in three Ontario peatlands, from which we have gained a better understanding of the peat properties that control the geophysical responses. The electrical conductivity depends linearly on the concentration of total dissolved solids in the peat pore waters and the pore waters in turn bear the ionic signatures of the underlying mineral sediments. The ionic concentration, and thus the electrical conductivity, increase linearly from the surface to basement. The average bulk electrical conductivity of peatlands at Ellice Marsh, near Stratford, and at Wally Creek Area Forest Drainage Project, near Cochrane, are of the order of 25 mS/m. The Mer Bleue peatland, near Ottawa, has extremely high electrical conductivity, reaching levels of up to 380 mS/m near the base of the peat. The Mer Bleue peatland water has correspondingly high values of total dissolved solids, which originate from the underlying Champlain Sea glaciomarine clays. The dielectric permittivity in peats is largely controlled by the bulk water content. Ground penetrating radar can detect changes in water content greater than 3%, occurring within a depth interval less than 15 cm. The principal peatland interfaces detected are the near-surface aerobic to anaerobic transition and the peat to mineral basement contact. The potential for the successful detection of the basement contact using the radar can be predicted using the radar instrument specifications, estimates of the peatland depth, and either the bulk peat or the peat pore water electrical conductivities. Predicted depths of penetration of up to 10 m for Ellice Marsh and Wally Creek exceed the observed depths of 1 to 2 m. At Mer Bleue, on the other hand, we observe that, as predicted, a 100 MHz signal will penetrate to the base of a 2 m thick peat but a 200 MHz signal will not.     

A two-dimensional analytical solution of groundwater responses to tidal loading in an estuary and ocean

Previous studies on tidal dynamics of coastal aquifers have focussed on the inland propagation of oceanic tides in the cross-shore direction, a configuration that is essentially one-dimensional. Aquifers at natural coasts can also be influenced by tidal waves in nearby estuaries, resulting in a more complex behaviour of head fluctuations in the aquifers. We present an analytical solution to the two-dimensional depth-averaged groundwater flow equation for a semi-infinite aquifer subject to oscillating head conditions at the boundaries. The solution describes the tidal dynamics of a coastal aquifer that is adjacent to a cross-shore estuary. Both the effects of oceanic and estuarine tides on the aquifer are included in the solution. The analytical prediction of the head fluctuations is verified by comparison with numerical solutions computed using a standard finite-difference method. An essential feature of the present analytical solution is the interaction between the cross- and along-shore tidal waves in the aquifer area near the estuary’s entry. As the distance from the estuary or coastline increases, the wave interaction is weakened and the aquifer response is reduced, respectively, to the one-dimensional solution for oceanic tides or the solution of Sun (Sun H. A two-dimensional analytical solution of groundwater response to tidal loading in an estuary, Water Resour Res 1997;33:1429–35) for two-dimensional non-interacting tidal waves.     

Deep-sea faunal communities associated with a lost intermodal shipping container in the Monterey Bay National Marine Sanctuary,CA

Carrying assorted cargo and covered with paints of varying toxicity, lost intermodal containers may take centuries to degrade on the deep seafloor. In June 2004, scientists from Monterey Bay Aquarium Research Institute (MBARI) discovered a recently lost container during a Remotely Operated Vehicle (ROV) dive on a sediment-covered seabed at 1281 m depth in Monterey Bay National Marine Sanctuary (MBNMS). The site was revisited by ROV in March 2011. Analyses of sediment samples and high-definition video indicate that faunal assemblages on the container’s exterior and the seabed within 10 m of the container differed significantly from those up to 500 m. The container surface provides hard substratum for colonization by taxa typically found in rocky habitats. However, some key taxa that dominate rocky areas were absent or rare on the container, perhaps related to its potential toxicity or limited time for colonization and growth. Ecological effects appear to be restricted to the container surface and the benthos within ∼10 m.     

Atlas of Japan (East) Sea hydrographic properties in summer, 1999

Hydrographic properties from CTD and discrete bottle sample profiles covering the Japan (East) Sea in summer, 1999, are presented in vertical sections, maps at standard depths, maps on isopycnal surfaces, and as property–property distributions. This data set covers most of the Sea with the exception of the western boundary region and northern Tatar Strait, and includes nutrients, pH, alkalinity, and chlorofluorocarbons, as well as the usual temperature, salinity, and oxygen observations.     

Petrophysics and mineral exploration: a workflow for data analysis and a new interpretation framework

As mineral exploration seeks deeper targets, there will be a greater reliance on geophysical data and a better understanding of the geological meaning of the responses will be required, and this must be achieved with less geological control from drilling. Also, exploring based on the mineral system concept requires particular understanding of geophysical responses associated with altered rocks. Where petrophysical datasets of adequate sample size and measurement quality are available, physical properties show complex variations, reflecting the combined effects of various geological processes. Large datasets, analysed as populations, are required to understand the variations. We recommend the display of petrophysical data as frequency histograms because the nature of the data distribution is easily seen with this form of display. A petrophysical dataset commonly contains a combination of overlapping sub-populations, influenced by different geological factors. To understand the geological controls on physical properties in hard rock environments, it is necessary to analyse the petrophysical data not only in terms of the properties of different rock types. It is also necessary to consider the effects of processes such as alteration, weathering, metamorphism and strain, and variables such as porosity and stratigraphy. To address this complexity requires that much more supporting geological information be acquired than in current practice. The widespread availability of field portable instruments means quantitative geochemical and mineralogical data can now be readily acquired, making it unnecessary to rely primarily on categorical rock classification schemes. The petrophysical data can be combined with geochemical, petrological and mineralogical data to derive explanations for observed physical property variations based not only on rigorous rock classification methods, but also in combination with quantitative estimates of alteration and weathering. To understand how geological processes will affect different physical properties, it is useful to define three end-member forms of behaviour. Bulk behaviour depends on the physical properties of the dominant mineral components. Density and, to a lesser extent, seismic velocity show such behaviour. Grain and texture behaviour occur when minor components of the rock are the dominate controls on its physical properties. Grain size and shape control grain properties, and for texture properties the relative positions of these grains are also important. Magnetic and electrical properties behave in this fashion. Thinking in terms of how geological processes change the key characteristics of the major and minor mineralogical components allows the resulting changes in physical properties to be understood and anticipated.     

Application of the SHETRAN basin‐scale,landslide sediment yield model to the Llobregat basin,Spanish Pyrenees

The first application of the SHETRAN basin‐scale, landslide erosion and sediment yield model is carried out for a major landsliding event in the upper 505 km² of the Llobregat basin, in the eastern Spanish Pyrenees, in November 1982. The model simulates the spatial distribution of shallow landslides and their sediment yield. Acknowledging uncertainty in the model parameter evaluation, the aim of the application was not to reproduce the observed occurrence of landslides as accurately as possible with one simulation, but to bracket the observed pattern with several simulations representing uncertainty in the key input conditions. Bounds on the landslide simulations were thus determined as a function of uncertainty in the vegetation root cohesion (used in the model factor of safety calculations). The resulting upper bound considerably overestimates the observed pattern (17 000 landslides compared with an observation of around 700), but it reproduces several of the principal clusters in the observed pattern. The lower bound contains around 500 landslides. The sediment yield estimates (2670–14 630 t km^?2) are comparable to measurements elsewhere in the Pyrenees for extreme events. The results demonstrate an ability to simulate the basin‐scale landslide response to a rainfall event and the resulting sediment yield. They also highlight the need for further research in setting the uncertainty bounds and in avoiding large overestimates of landslide occurrence arising in part from a current inability to model small‐scale controls for a basin of the given size. Copyright © 2006 John Wiley & Sons, Ltd.     

Energy market failure in road transport: Is there scope for ‘no regrets’ greenhouse gas reduction?

The Australian Government policy on reduction of greenhouse gas emissions announced in 1990 includes exploring the scope for immediate, low cost reductions. Such measures can be taken as including ‘no regrets’ policies: those that, in addition to mitigating potential climate change, confer economic gains (including other environmental benefits) which exceed their costs. Some possible ‘no regrets’ opportunities and policies are identified relevant to energy use by the road transport sector over the period to 2020. The MARKALMENSA multi-period linear programming model of the Australian energy sector is used to investigate the cost-effectiveness of these policies.