Variations of trace gases over the Bay of Bengal during the summer monsoon

In situ measurements of near-surface ozone (\(\hbox {O}_{3})\), carbon monoxide (CO), and methane (\(\hbox {CH}_{4})\) were carried out over the Bay of Bengal (BoB) as a part of the Continental Tropical Convergence Zone (CTCZ) campaign during the summer monsoon season of 2009. \(\hbox {O}_{3}\), CO and \(\hbox {CH}_{4}\) mixing ratios varied in the ranges of 8–54 ppbv, 50–200 ppbv and 1.57–2.15 ppmv, respectively during 16 July–17 August 2009. The spatial distribution of mean tropospheric \(\hbox {O}_{3}\) from satellite retrievals is found to be similar to that in surface \(\hbox {O}_{3}\) observations, with higher levels over coastal and northern BoB as compared to central BoB. The comparison of in situ measurements with the Monitoring Atmospheric Composition & Climate (MACC) global reanalysis shows that MACC simulations reproduce the observations with small mean biases of 1.6 ppbv, –2.6 ppbv and 0.07 ppmv for \(\hbox {O}_{3}\), CO and \(\hbox {CH}_{4}\), respectively. The analysis of diurnal variation of \(\hbox {O}_{3}\) based on observations and the simulations from Weather Research and Forecasting coupled with Chemistry (WRF-Chem) at a stationary point over the BoB did not show a net photochemical build up during daytime. Satellite retrievals show limitations in capturing \(\hbox {CH}_{4}\) variations as measured by in situ sample analysis highlighting the need of more shipborne in situ measurements of trace gases over this region during monsoon.     

Identification of buried sinkholes using refraction tomography at Ft. Campbell Army Airfield,Kentucky

Karst aquifers are highly susceptible to contamination, with numerous points of entry for contaminants through recharge features such as sinkholes, swallow holes and solutionally enlarged fractures. These recharge features may be filled or obscured at the surface, requiring the use of geophysical or remote sensing techniques for their identification. This study uses seismic refraction data collected at the Ft. Campbell Army Airfield (CAAF), Kentucky, USA, to test the hypothesis that refraction tomography is a useful tool for imaging bedrock depressions beneath thick overburden (greater than 20 m of unconsolidated sediment). Southeast of the main taxiway of CAAF seismic velocity tomograms imaged a bedrock low, possibly a closed depression, at a depth of 25 m that had been earlier identified through delay-time analysis of the same refraction data. Tomography suggests the bedrock low is about 250-m wide by 10-m deep at its widest point. High rates of contaminant vapor extraction over the western extension of this feature suggest a high concentration of contaminants above, and within, this filled bedrock low, the base of which may contain solutionally enlarged fractures (i.e. karst conduits) that could funnel these contaminants to the upper or lower bedrock aquifers. This study thus demonstrates the viability of seismic refraction tomography as a tool for identification of filled sinkholes and bedrock depressions in karst areas.     

Application of stress mapping in cross‐section to understanding ore geometry,predicting ore zones and development of drilling strategies

Stress mapping is a numerical modelling technique used to determine the distribution and relative magnitude of stress during deformation in a mineralised terrane. It is based on the general principle that fluid flow in the Earth's crust is primarily related to pressure gradients. It is best applied to epigenetic hydrothermal mineral deposits, where fluid flow and fluid flux are enhanced in dilational sections of structures and in sites of enhanced rock permeability due to high fracture density. These are defined by sites of low minimum principal stress (σ₃). Most stress mapping is carried out in two dimensions in plan view using geological maps. This is suitable for terranes with steeply dipping lithostratigraphy and structures in which the distribution of mineral deposits is largely controlled by fault structures portrayed on the maps. However, for terranes with gently dipping sequences and structures, and for situations where deposits are sited in and near the hinges of complex fold structures, stress mapping in cross‐section is preferable. The effectiveness of stress mapping is maximised if mineralisation was late in the evolutionary history of the host terrane, and hence the structural geometry of the terrane and contained deposits were essentially that expressed today. The orientation of syn‐mineralisation far‐field stresses must also be inferred. Two examples of orogenic gold deposits, which meet the above criteria, are used to illustrate the potential of stress mapping in cross‐section. Sunrise Dam, located in the Archaean Yilgarn Craton, is a lode‐gold deposit sited in a thrust‐fold belt. Stress mapping illustrates the heterogeneity of stress distribution in the complex structural geometry of the deposit, and predicts the preferential siting of ore zones around the intersections of more steeply dipping, linking thrusts and banded iron‐formation units, and below the controlling more gently dipping basal thrust, the Sunrise Shear. The Howley Anticline in the Pine Creek block hosts several Palaeoproterozoic gold deposits, sited in complex anticlinal structures in greywacke sequences. Stress mapping indicates that gold ores should develop in the hinge zones of symmetrical anticlines, in the hinge zones and more steeply dipping to overturned limbs of asymmetric anticlines, and in and around thrusts in both anticlines and parasitic synclines. The strong correlation between the predictions of the stress mapping, based on the distribution of low σ₃, and the location of gold ores emphasises the potential of stress mapping in cross‐section, not only as an exploration tool for the discovery of additional resources or deposits, but also as a test of geological models. Knowledge of the potential siting of gold ores and their probable orientations also provides a guide to drilling strategies in both mine‐ and regional‐scale exploration.     

TWO-DIMENSIONAL MASS DISTRIBUTIONS FROM GRAVITY ANOMALIES: A COMPUTER METHOD*

A gravity anomaly profile, considered to be associated with a two-dimensional mass distribution of constant density contrast, may be interpreted with the help of a digital computer. Using the infinite slab formula and anomaly values at selected points of the profile, an initial model of the anomalous mass is built. The gravity effect of the model is computed with a simplified version of the formula derived by Talwani, Worzel and Landisman (1959) for line-integration. After the effect is computed at a point, the model is adjusted depending upon the residual anomaly. The effect of the adjusted model is computed at the following point. Within one cycle of computation the model is adjusted as many times as the number of the selected points. Five such cycles are sufficient to minimize residuals in most cases. A final cycle is computed without adjustment. Local anomaly profiles have been interpreted to obtain the shape of a sedimentary basin and of a faulted block. A regional anomaly profile has been interpreted to obtain crustal structure.     

Specific features of variations in the characteristics of VLF signals when the lunar shadow propagated along the path during the solar eclipse of March 29, 2006

The effects of the solar eclipse of March 29, 2006, in the signals of ULF radio stations, in the intensity of regular radio noise at frequencies of 0.3–10 kHz, and in the number of atmospherics received in Yakutsk mostly from the west have been considered. The observations were performed using a multichannel parallel analyzer-recorder (11 channels in the frequency band 0.47–8.7 kHz), one-point lightning direction and range finder (0.3–100 kHz), narrow-sector radio noise direction-finder (0.3–10 kHz), recorder of signals from VLF radio stations, and broadband radio noise recorder (0.3–100 kHz). A GPS clock was used to synchronize a recorder of signals from VLF radio stations. The effect was observed in radio signals, radio noise, and number of atmospherics from the direction 270° ± 20° counted off clockwise from the north during the last stage of the eclipse (～ 1100–1200 UT), when the lunar shadow approached the line of the nighttime terminator and obscured part of the signal propagation path. The effect was observed as an enhancement of the received signals by a factor of ～1.2, a factor of ～1.4 increase in the number of atmospherics, and a change in the radio station phase values.     

Subglacial bedforms and conditions associated with the 1991 surge of Skeiðarárjökull, Iceland

Much previous research at surge-type glaciers has sought to identify features diagnostic of surge-type behaviour. However, in comparatively little work have subglacial landform–sediment characteristics been used to reconstruct changing basal processes and conditions during surge events. Subglacial bedforms described in this article are associated with the 1991 surge of Skeiðarárjökull, Iceland, and include a series of drumlins with superimposed flutes and basal crevasse-fill ridges. The drumlins were formed by the subglacial erosion of ice-contact fans. Sedimentary evidence indicates a shift from rigid-bed to soft-bed conditions during the surge. The presence of eroded but undeformed fan sediments suggests that they acted as a rigid bed when initially overridden. Subsequent deposition of a layer of deformation till resulted in a change to soft-bed conditions and the generation of flutes and subglacial crevasse-fill ridges. The lack of mixing between this till and the underlying stratified sediments indicates that subglacial sediment deformation was restricted to a thin layer and that its deposition resulted in a cessation of subglacial erosion. The drumlin is therefore a composite of both rigid-bed and soft-bed processes that illustrates changes in basal conditions and processes during the course of the event. The limited time frame in which the drumlin formed and the presence of kettleholes across its surface are distinctive features that may warrant further investigation in the search for features diagnostic of past surge events.     

Coupled HM analysis using zero‐thickness interface elements with double nodes—Part II: Verification and application

In a companion Part I of this paper (Int. J. Numer. Anal. Meth. Geomech. 2008; DOI: 10.1002/nag.735 ), a coupled hydro‐mechanical (HM) formulation for geomaterials with discontinuities based on the finite element method (FEM) with double‐node, zero‐thickness interface elements was developed and presented. This Part II paper includes the numerical solution of basic practical problems using both the staggered and the fully coupled approaches. A first group of simulations, based on the classical consolidation problem with an added vertical discontinuity, is used to compare both the approaches in terms of accuracy and convergence. The monolithic or fully coupled scheme is also used in an application example studying the influence of a horizontal joint in the performance of a reservoir subject to fluid extraction. Results include a comparison with other numerical solutions from the literature and a sensitivity analysis of the mechanical parameters of the discontinuity. Some simulations are also run using both a full non‐symmetric and a simplified symmetric Jacobian matrix. On top of verifying the model developed and its capability to reflect the conductivity changes of the interface with aperture changes, the results presented also lead to interesting observations of the numerical performance of the methods implemented. Copyright © 2008 John Wiley & Sons, Ltd.     

Modelling the non‐linear hydrological behaviour of a small Mediterranean forested catchment

A progressive perceptual understanding approach was used to identify a model structure able to represent the non‐linear behaviour of the hydrological cycle in a small intermittent Mediterranean stream. The initial lumped model structure consisting of a series of four connected water tanks (LU3) progressed to a model with five tanks (LU4), and finally to a semidistributed model structure (SD4) in which spatial variability of the evapotranspiration according to the vegetation cover and to the local aspect was considered. In the final model structure, which gave the best fit (Nash–Sutcliffe efficiency index = 0·78), an additional tank representing the riparian zone was included (SD4‐R). Results showed that the abrupt changes of the riparian water table during summer and the formation of a perched water table during the transition from dry to wet conditions were the main mechanisms leading to the non‐linear hydrological behaviour. The transpiration process from the saturated zone and the spatial variability of evapotranspiration resulted in key factors successfully representing the annual water balance. The spatial and temporal validations carried out for each of the four model structures considered in this study supported the hypothesis adopted during the calibration process. Copyright © 2008 John Wiley & Sons, Ltd.     

Influences of climate,terrain and land cover on stream salinity in southeastern Australia,and implications for management through reforestation

Reforestation of cleared land has the potential to reduce groundwater recharge, salt mobilization and streamflow. Stream salinity change is the net result of changes in stream salt load and streamflow. The net effect of these changes varies spatially as a function of climate, terrain and land cover. Successful natural resource management requires methods to map the spatial variability of reforestation impacts. We investigated salinity data from 2000 bores and streamflow and salinity measurements from 27 catchments in the Goulburn–Broken region in southeast Australia to assess the main factors determining stream salinity and opportunities for management through reforestation. For groundwater systems of similar geology, relationships were found between average annual rainfall and groundwater salinity and between groundwater salinity and low‐flow salinity. Despite its simplicity, we found that the steady‐state component of a simple conceptual coupled water–salt mass balance model (BC2C) adequately explained the spatial variation in streamflow and salinity. The model results suggest the efficiency of afforestation to reduce stream salinity could be increased by more than an order of magnitude through spatial planning. However, appreciable reductions in stream salinity in large rivers through land cover change alone would still require reforestation on an unprecedented scale. Copyright © 2008 John Wiley & Sons, Ltd.     

A Comparative Study of Landsat‐7 and Landsat‐8 Data Using Image Processing Methods for Hydrothermal Alteration Mapping

Comparing spaceborne satellite images of Landsat‐8 Operational Land Imager (OLI) and Landsat‐7 Enhanced Thematic Mapper plus (ETM+) was undertaken to investigate the relative accuracy of mapping hydrothermal alteration minerals. The study investigated the northern part of Rabor, which contains copper mineralization occurrences, and is located in the Kerman Cenozoic magmatic assemblage (KCMA), Iran. Image processing methods of band ratio, principal component analysis (PCA), and spectral angle mapper (SAM) were used to map the distribution of hydrothermally altered rocks associated with the porphyry copper mineralization. The band ratio combination of both sensors for mapping altered areas showed similar outcomes. PCA exposed variations in the spatial distribution of hydroxyl‐bearing minerals. The representation of hydrothermal areas using OLI data was more satisfactory than when using ETM+ data. SAM analysis found similar results for mapping hydroxyl‐bearing zones. Verification of the results came through ground investigation and laboratory studies. Rock samples (n = 56) were collected to validate results using thin sections, X‐ray diffraction (XRD) and spectral analyses. Field observations and laboratory analysis revealed that phyllic and propylitic alterations dominate the alteration zones in the study area. Argillic and iron oxides/hydroxides alterations were observed to a lesser degree. The results indicate that alteration maps prepared by OLI data using PCA for visual interpretation are more suitable than those of ETM+ due to a higher radiometric resolution and lower interference between vegetation and altered areas. As the spectral bandwidth of ETM+ band 7 covers absorption feature of propylitic alteration, better mapping of propylitic alterations is achieved using ETM+ data.