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.     

Refining accretionary orogen models for the Tasmanides of eastern Australia

The well-known southwest-to-northeast younging of stratigraphy over a present-day cross strike distance of >1500 km in the southern Tasmanides of eastern Australia has been used to argue for models of accretionary orogenesis behind a continually eastwards-rolling paleo-Pacific plate. However, these accretionary models need modification, since the oldest (ca 530 Ma) outcrops of Cambrian supra-subduction zone rocks occur in the outboard New England Orogen, now ～900 km east of the next oldest (520–510 Ma) supra-subduction zone rocks. This is not consistent with simple, continuous easterly rollback. Instead, the southern Tasmanides contain an early history characterised by a westwards-migrating margin between ca 530 and ca 520 Ma, followed by rapid eastwards rollback of the paleo-Pacific plate from 520 to 502 Ma that opened a vast backarc basin ～2000 km across that has never been closed. From the Ordovician through to the end of the Carboniferous, the almost vertical stacking of continental margin arcs (within a hundred kilometres of each other) in the New England Orogen indicates a constant west-dipping plate boundary in a Gondwana reference frame. Although the actual position of the boundary is inferred to have undergone contraction-related advances and extension-related retreats, these movements are estimated to be ～250 km or less. Rollback in the early Permian was never completely reversed, so that late Permian–Triassic to Cretaceous arcs lie farther east, in the very eastern part of eastern Australia, with rifted fragments occurring in the Lord Howe Rise and in New Zealand. The northern Tasmanides are even more anomalous, since they missed out on the middle Cambrian plate boundary retreat seen in the south. As a result, their Cambrian-to-Devonian history is concentrated in a ～300 km wide strip immediately west of Precambrian cratonic Australia and above Precambrian basement. The presence in this narrow region of Ordovician to Carboniferous continental margin arcs and backarc basins also implies a virtually stationary plate boundary in a Gondwana frame of reference. This bipolar character of the Tasmanides suggests the presence of a segmented paleo-Pacific Plate, with major transform faults propagating into the Tasmanides as tear faults that were favourably oriented for the formation of local supra-subduction zone systems and for subsequent intraplate north–south shortening. In this interpretation of the Tasmanides, Lower–Middle Ordovician quartz-rich turbidites accumulated as submarine fan sequences, and do not represent multiple subduction complexes developed above subduction zones lying behind the plate boundary. Indeed, the Tasmanides are characterised by the general absence of material accreted from the paleo-Pacific plate and by the dominance of craton-derived, recycled sedimentary rocks.     

Internal gravity waves in a slowly varying,dissipative medium

Nonlinear internal gravity waves in a slightly dissipative, slightly compressible fluid are discussed for the case when the properties of the medium vary slowly on a scale determined by the local wave structure. A two‐timing technique is used to obtain transport equations which describe the changes in amplitude, phase and mean flow of a wave packet. Various solutions of these transport equations are discussed, with relevance to critical layer absorption.     

Crustal complexity in the Lachlan Orogen revealed from teleseismic receiver functions

There is an ongoing debate about the tectonic evolution of southeast Australia, particularly about the causes and nature of its accretion to a much older Precambrian core to the west. Seismic imaging of the crust can provide useful clues to address this issue. Seismic tomography imaging is a powerful tool often employed to map elastic properties of the Earth's lithosphere, but in most cases does not constrain well the depth of discontinuities such as the Mohorovi?i? (Moho). In this study, an alternative imaging technique known as receiver function (RF) has been employed for seismic stations near Canberra in the Lachlan Orogen to investigate: (i) the shear-wave-velocity profile in the crust and uppermost mantle, (ii) variations in the Moho depth beneath the Lachlan Orogen, and (iii) the nature of the transition between the crust and mantle. A number of styles of RF analyses were conducted: H-K stacking to obtain the best compressional–shear velocity (V _P /V _S) ratio and crustal thickness; nonlinear inversion for the shear-wave-velocity structure and inversion of the observed variations in RFs with back-azimuth to investigate potential dipping of the crustal layers and anisotropy. The thick crust (up to 48 km) and the mostly intermediate nature of the crust?mantle transition in the Lachlan Orogen could be due to the presence of underplating at the base of the crust, and possibly to the existing thick piles of Ordovician mafic rocks present in the mid and lower crust. Results from numerical modelling of RFs at three seismic stations (CAN, CNB and YNG) suggest that the observed variations with back-azimuth could be related to a complex structure beneath these stations with the likelihood of both a dipping Moho and crustal anisotropy. Our analysis reveals crustal thickening to the west beneath CAN station which could be due to slab convergence. The crustal thickening may also be related to the broad Macquarie volcanic arc, which is rooted to the Moho. The crustal anisotropy may arise from a strong N–S structural trend in the eastern Lachlan Orogen and to the preferred crystallographic orientation of seismically anisotropic minerals in the lower and middle crust related to the paleo-Pacific plate convergence.     

An answer to an important controversy and the need for caution when using simple models to predict inelastic earthquake response of buildings with torsion

This paper presents evidence that the extension of conclusions based on the widely used simplified, one story, eccentric systems of the shear‐beam type, to actual, nonsymmetric buildings and consequent assessments of the pertinent code provisions, can be quite erroneous, unless special care is taken to match the basic properties of the simplified models to those of the real buildings. The evidence comes from comparisons of results obtained using three variants of simplified models, with results from the inelastic dynamic response of three‐ and five‐story eccentric buildings computed with detailed MDOF systems, where the members are idealized with the well‐known plastic hinge model. In addition, a convincing answer is provided on a pertinent hanging controversy: For frame‐type buildings, designed in accordance with the dynamic provisions of modern codes (such as EC8 or IBC2000), which allow reduced shears at the stiff edge due to torsion, the frames at the flexible sides are the critical elements in terms of ductility demands. Copyright © 2009 John Wiley & Sons, Ltd.     

Manganese in runoff from upland catchments: temporal patterns and controls on mobilization / Teneur en manganèse de l'écoulement en bassins versants d'altitude: variations temporelles et contrôles de la mobilisation

Abstract

Knowledge of the hydrochemical dynamics of the trace metal manganese (Mn) in upland catchments is required for water quality management. Stream water Mn and other solutes and flow were monitored in two upland catchments in northern England with different soils: one dominated by peat (HS7), the other by mineral soils (HS4). Maximum Mn concentrations occurred at different times in the two catchments: in summer baseflow at HS4 and during late summer storm events at HS7. A two-component chemical mixing model was used to identify the hydrological processes controlling Mn concentrations in stream water. This approach was more successful for HS4 than HS7, probably because of different processes of Mn release in the two catchments and also difficulties in selecting conservative solutes. Factor analysis of the stream water chemistry data set for each catchment was more useful in identifying the controls on Mn release into runoff. The factors indicate that the main source of Mn at HS4 is the hydrological pathway supplying summer baseflow, whereas at HS7 Mn is released during the rewetting of dried peat soils. Manganese concentrations in stream water in upland catchments appear to depend on soil type and antecedent moisture conditions. This has implications for the design of sampling strategies in upland catchments and also for managing the quality of water supplies from such areas.     

Pulse generation models for daily intermittent streamflows / Modèles de génération d’impulsions pour des débits journaliers intermittents

Abstract

Abstract Generating pulses and then converting them into flow are two main steps of daily streamflow generation. Three pulse generation models have been proposed on the basis of Markov chains for the purpose of generating daily intermittent streamflow time series in this study. The first one is based on two two-state Markov chains, whereas the second uses a three-state Markov chain. The third model uses harmonic analysis and fits Fourier series to the three-state Markov chain. Results for a daily intermittent streamflow data series show a good performance of the proposed models.     

Assessing uncertainties in a conceptual water balance model using Bayesian methodology / Estimation bayésienne des incertitudes au sein d’une modélisation conceptuelle de bilan hydrologique

Abstract

Abstract The aim of this study was to estimate the uncertainties in the streamflow simulated by a rainfall–runoff model. Two sources of uncertainties in hydrological modelling were considered: the uncertainties in model parameters and those in model structure. The uncertainties were calculated by Bayesian statistics, and the Metropolis-Hastings algorithm was used to simulate the posterior parameter distribution. The parameter uncertainty calculated by the Metropolis-Hastings algorithm was compared to maximum likelihood estimates which assume that both the parameters and model residuals are normally distributed. The study was performed using the model WASMOD on 25 basins in central Sweden. Confidence intervals in the simulated discharge due to the parameter uncertainty and the total uncertainty were calculated. The results indicate that (a) the Metropolis-Hastings algorithm and the maximum likelihood method give almost identical estimates concerning the parameter uncertainty, and (b) the uncertainties in the simulated streamflow due to the parameter uncertainty are less important than uncertainties originating from other sources for this simple model with fewer parameters.     

Integration of remote sensing data into hydrological models for reservoir management

Abstract

The purpose of this paper is to present the methodology set up to derive catchment soil moisture from Earth Observation (EO) data using microwave spaceborne Synthetic Aperture Radar (SAR) images from ERS satellites and to study the improvements brought about by an assimilation of this information into hydrological models. The methodology used to derive EO data is based on the appropriate selection of land cover types for which the radar signal is mainly sensitive to soil moisture variations. Then a hydrological model is chosen, which can take advantage of the new information brought by remote sensing. The assimilation of soil moisture deduced from EO data into hydrological models is based principally on model parameter updating. The main assumption of this method is that the better the model simulates the current hydrological system, the better the following forecast will be. Another methodology used is a sequential one based on Kalman filtering. These methods have been put forward for use in the European AIMWATER project on the Seine catchment upstream of Paris (France) where dams are operated to alleviate floods in the Paris area.     

Seismic stability analysis of piled embankments: A multiphase approach

The seismic stability analysis of an embankment lying over a soft foundation soil reinforced by a group of vertical piles is performed within the framework of the upper bound kinematic approach of yield design. The analysis is based on a previously developed ‘multiphase’ model of the reinforced ground, which explicitly accounts for the shear and bending resistances of the piles. Making use of appropriate failure mechanisms involving shear zones across which the reinforcements are continuously deforming, along with ‘plastic hinge’ surfaces, upper bound estimates to the critical seismic coefficient of the structure are derived. The results, which are confirmed by the simulations obtained from a finite element elastoplastic code, give clear evidence of the key role played by the bending strength capacities of the piles in ensuring the stability of the pile reinforced embankment. Copyright © 2009 John Wiley & Sons, Ltd.