Using constrained inversion of gravity and magnetic field to produce a 3D litho‐prediction model

Geologically constrained inversion of gravity and magnetic field data of the Victoria property (located in Sudbury, Canada) was undertaken in order to update the present three‐dimensional geological model. The initial and reference model was constructed based on geological information from over 950 drillholes to constrain the inversion. In addition, downhole density and magnetic susceptibility measured in six holes were statistically analysed to derive lower and upper bounds on the physical properties attributed to the lithological units in the reference model. Constrained inversion of the ground gravity and the airborne magnetic data collected at the Victoria property were performed using GRAV3D and MAG3D, respectively. A neural network was trained to predict lithological units from the physical properties measured in six holes. Then, the trained network was applied on the three‐dimensional distribution of physical properties derived from the inversion models to produce a three‐dimensional litho‐prediction model. Some of the features evident in the lithological model are remnants of the constraints, where the data did not demand a significant change in the model from the initial constraining model (e.g., the thin pair of diabase dykes). However, some important changes away from the initial model are evident; for example, a larger body was predicted for quartz diorite, which may be related to the prospective offset dykes; a new zone was predicted as sulfide, which may represent potential mineralisation; and a geophysical subcategory of metabasalt was identified with high magnetic susceptibility and high density. The litho‐prediction model agrees with the geological expectation for the three‐dimensional structure at Victoria and is consistent with the geophysical data, which results in a more holistic understanding of the subsurface lithology.     

Potential of legacy 2D seismic data for deep targeting and structural imaging at the Neves–Corvo massive sulphide-bearing deposit,Portugal

Seismic methods are becoming an established choice for deep mineral exploration after being extensively tested and employed for the past two decades. To investigate whether the early European mineral-exploration datasets had potential for seismic imaging that was overlooked, we recovered a low-fold legacy seismic dataset from the Neves–Corvo mine site in the Iberian Pyrite Belt in southern Portugal. This dataset comprises six 4–6 km long profiles acquired in 1996 for deep targeting. Using today's industry-scale processing algorithms, the world-class, ca. 150 Mt, Lombador massive sulphide and other smaller deposits were better imaged. Additionally, we also reveal a number of shallow but steeply dipping reflections that were absent in the original processing results. This study highlights that legacy seismic data are valuable and should be revisited regularly to take advantage of new processing algorithms and the experiences gained from processing such data in hard-rock environments elsewhere. Remembering that an initial processing job in hard rock should always aim to first obtain an overall image of the subsurface and make reflections visible, and then subsequent goals of the workflow could be set to, for example understanding relative amplitude ratios. The imaging of the known mineralization implies that this survey could likely have been among one of the pioneer studies in the world that demonstrated the capability of directly imaging massive sulphide deposits using the seismic method.     

Lithostratigraphy, Sedimentology, and Provenance of the Balfour Formation (Beaufort Group) in the Fort Beaufort–Alice Area, Eastern Cape Province, South Africa

Abstract: The Balfour Formation has a pronounced lithological variation that is characterized by alternating sandstone- and mudstone-dominated members. The sandstone-dominated Oudeberg and Barberskrans Members are composed of lithofacies that range from intraformational conglomerates to fine-grained sediments, whereas the mudstone-dominated members (Daggaboersnek, Elandsberg, and Palingkloof) are dominated by the facies Fm and Fl. Petrography, geochemistry, and a paleocurrent analysis indicated that the source rock of the Balfour Formation was to south east and the rocks had a transitional/dissected magmatic arc signature. The sandstones-rich members were deposited by seasonal and ephemeral high-energy, low-sinuous streams, and the fine-grained-rich members were formed by ephemeral meandering streams. The paleoclimates have been equated to present temperate climates; they were semiarid becoming arid towards the top of the Balfour Formation. This has been determined by reconstructing the paleolatitude of the Karoo Basin, geochemistry, paleontology, sedimentary structures, and other rock properties, like color.     

Experimental and modeling study on Cr(VI) transfer from soil into surface runoff

With the development of modern agriculture, large amount of fertilizer and pesticide outflow from farming land causes great waste and serious pollution to surface water and groundwater, and threatens ecological environment and even human life. In this paper, laboratory experiments are conducted to simulate adsorbed Cr(VI) transfer from soil into runoff. A two-layer in-mixing analytical model is applied to analyze laboratory experimental results. A data assimilation (DA) method via the ensemble Kalman filter (EnKF) is used to update parameters and improve modeling results. In comparison with the experimental data, DA updated modeling results are much better than those without the updating. To make predictions better, the inflation method with a constant inflation factor via DA method was used to compensate the fast decrease of ensemble spread partially related to filter inbreeding. Based on the used rainfall and relevant physical principles, the updated value of the incomplete mixing coefficient γ is about 14.0 times of the value of the incomplete mixing coefficient α in experiment 1 and about 7.4 times in experiment 2, while the difference between the flow rate of runoff and infiltration is not so large even after reaching stable infiltration condition. The results indicate the loss of Cr(VI) in soil solute is mainly due to infiltration, rather than surface runoff. With the increase of mixing layer depth, soil adsorption capacity will increase and the loss of soil solute will decrease. The study results provide information for reducing and even preventing the agricultural nonpoint source pollution.     

On the use of IPCC-class models to assess the impact of climate on Living Marine Resources

The study of climate impacts on Living Marine Resources (LMRs) has increased rapidly in recent years with the availability of climate model simulations contributed to the assessment reports of the Intergovernmental Panel on Climate Change (IPCC). Collaboration between climate and LMR scientists and shared understanding of critical challenges for such applications are essential for developing robust projections of climate impacts on LMRs. This paper assesses present approaches for generating projections of climate impacts on LMRs using IPCC-class climate models, recommends practices that should be followed for these applications, and identifies priority developments that could improve current projections. Understanding of the climate system and its representation within climate models has progressed to a point where many climate model outputs can now be used effectively to make LMR projections. However, uncertainty in climate model projections (particularly biases and inter-model spread at regional to local scales), coarse climate model resolution, and the uncertainty and potential complexity of the mechanisms underlying the response of LMRs to climate limit the robustness and precision of LMR projections. A variety of techniques including the analysis of multi-model ensembles, bias corrections, and statistical and dynamical downscaling can ameliorate some limitations, though the assumptions underlying these approaches and the sensitivity of results to their application must be assessed for each application. Developments in LMR science that could improve current projections of climate impacts on LMRs include improved understanding of the multi-scale mechanisms that link climate and LMRs and better representations of these mechanisms within more holistic LMR models. These developments require a strong baseline of field and laboratory observations including long time series and measurements over the broad range of spatial and temporal scales over which LMRs and climate interact. Priority developments for IPCC-class climate models include improved model accuracy (particularly at regional and local scales), inter-annual to decadal-scale predictions, and the continued development of earth system models capable of simulating the evolution of both the physical climate system and biosphere. Efforts to address these issues should occur in parallel and be informed by the continued application of existing climate and LMR models.     

Simulating flow in karst aquifers at laboratory and sub-regional scales using MODFLOW-CFP

Groundwater flow in a well-developed karst aquifer dominantly occurs through bedding planes, fractures, conduits, and caves created by and/or enlarged by dissolution. Conventional groundwater modeling methods assume that groundwater flow is described by Darcian principles where primary porosity (i.e. matrix porosity) and laminar flow are dominant. However, in well-developed karst aquifers, the assumption of Darcian flow can be questionable. While Darcian flow generally occurs in the matrix portion of the karst aquifer, flow through conduits can be non-laminar where the relation between specific discharge and hydraulic gradient is non-linear. MODFLOW-CFP is a relatively new modeling program that accounts for non-laminar and laminar flow in pipes, like karst caves, within an aquifer. In this study, results from MODFLOW-CFP are compared to those from MODFLOW-2000/2005, a numerical code based on Darcy’s law, to evaluate the accuracy that CFP can achieve when modeling flows in karst aquifers at laboratory and sub-regional (Woodville Karst Plain, Florida, USA) scales. In comparison with laboratory experiments, simulation results by MODFLOW-CFP are more accurate than MODFLOW 2005. At the sub-regional scale, MODFLOW-CFP was more accurate than MODFLOW-2000 for simulating field measurements of peak flow at one spring and total discharges at two springs for an observed storm event.     

Miocene silicic volcanism in southwestern Idaho: geochronology,geochemistry, and evolution of the central Snake River Plain

New ⁴⁰Ar-³⁹Ar geochronology, bulk rock geochemical data, and physical characteristics for representative stratigraphic sections of rhyolite ignimbrites and lavas from the west-central Snake River Plain (SRP) are combined to develop a coherent stratigraphic framework for Miocene silicic magmatism in this part of the Yellowstone ‘hotspot track’. The magmatic record differs from that in areas to the west and east with regard to its unusually large extrusive volume, broad lateral scale, and extended duration. We infer that the magmatic systems developed in response to large-scale and repeated injections of basaltic magma into the crust, resulting in significant reconstitution of large volumes of the crust, wide distribution of crustal melt zones, and complex feeder systems for individual eruptive events. Some eruptive episodes or ‘events’ appear to be contemporaneous with major normal faulting, and perhaps catastrophic crustal foundering, that may have triggered concurrent evacuations of separate silicic magma reservoirs. This behavior and cumulative time-composition relations are difficult to relate to simple caldera-style single-source feeder systems and imply complex temporal-spatial development of the silicic magma systems. Inferred volumes and timing of mafic magma inputs, as the driving energy source, require a significant component of lithospheric extension on NNW-trending Basin and Range style faults (i.e., roughly parallel to the SW–NE orientation of the eastern SRP). This is needed to accommodate basaltic inputs at crustal levels, and is likely to play a role in generation of those magmas. Anomalously high magma production in the SRP compared to that in adjacent areas (e.g., northern Basin and Range Province) may require additional sub-lithospheric processes. Electronic supplementary material The online version of this article (doi: ) contains supplementary material, which is available to authorized users. This paper constitutes part of a special issue dedicated to Bill Bonnichsen on the petrogenesis and volcanology of anorogenic rhyolites.