A Framework and Software Tool to Support Collaborative Landscape Analysis: Fitting Square Pegs into Square Holes

For landscape models to be applied successfully in management situations, models must address appropriate questions, include relevant processes and interactions, be perceived as credible and involve people affected by decisions. We propose a framework for collaborative model building that can address these issues, and has its roots in adaptive management, computer‐supported collaborative work and landscape ecology. Models built through this framework integrate a variety of information sources, address relevant questions, and are customized for the particular landscape and policy environment under study. Participants are involved in the process from the start, and because their input is incorporated, they feel ownership of the resulting models, increasing the chance of model acceptance and application. There are two requirements for success: a tool that supports rapid model prototyping and modification, that makes a clear link between a conceptual and implemented model, and that has the ability to implement a wide range of model types; and a core team with skills in communication, research and analysis, and knowledge of ecology and forestry in addition to modelling. SELES (Spatially Explicit Landscape Event Simulator) is a tool for building and running models of landscape dynamics. It combines discrete event simulation with a spatial database and a relatively simple modelling language to allow rapid development of landscape simulations, and provides a high‐level means of specifying complex model behaviours ranging from management actions to natural disturbance and succession. We have applied our framework in several forest modelling projects in British Columbia, Canada. We have found that this framework increases the interest by local experts and decision‐makers to participate actively in the model building process. The workshop process and resulting models have efficiently provided insight into the dynamics of large landscapes over long time frames. The use of SELES has facilitated this process by providing a flexible, transparent environment in which models can be rapidly implemented and refined. As a result, model findings may be more readily incorporated into decision‐support systems designed to assist resource managers in making informed decisions.     

The frequency response of the horizontal magnetic field for a conductive channel

Summary. Babour & Mosnier's results showing no frequency dependence between the anomalous horizontal magnetic field above a conductor and the difference between the horizontal magnetic fields above and below the conductor over a wide range of frequencies led them to conclude that this effect is due to current channelling. A two-dimensional numerical model of a conductive channel, with a uniform horizontal source field, shows the same effect over a wide range of frequencies. Thus local induction can show the same effect.     

Lithogeochemical localisation of disseminated gold in the White River area, Yukon, Canada

Gold mineralisation in the White River area, 80 km south of the highly productive Klondike alluvial goldfield, is hosted in amphibolite facies gneisses in the same Permian metamorphic pile as the basement for the Klondike goldfield. Hydrothermal fluid which introduced the gold was controlled by fracture systems associated with middle Cretaceous to early Tertiary extensional faults. Gold deposition occurred where highly fractured and chemically reactive rocks allowed intense water–rock interaction and hydrothermal alteration, with only minor development of quartz veins. Felsic gneisses were sericitised with recrystallisation of hematite and minor arsenic mobility, and extensively pyritised zones contain gold and minor arsenic (ca 10 ppm). Graphitic quartzites (up to 5 wt.% carbon) caused chemical reduction of mineralising fluids, with associated recrystallisation of metamorphic minerals (graphite, pyrrhotite, pyrite, chalcopyrite) in host rocks and veins, and introduction of arsenic (up to 1 wt.%) to form arsenopyrite in veins and disseminated through host rock. Veins have little or no hydrothermal quartz, and up to 19 wt.% carbon as graphite. Late-stage oxidation of arsenopyrite in some graphitic veins has formed pharmacosiderite. Gold is closely associated with disseminated and vein sulphides in these two rock types, with grades of up to 3 ppm on the metre scale. Other rock types in the White River basement rocks, including biotite gneiss, hornblende gneiss, pyroxenite, and serpentinite, have not developed through-going fracture systems because of their individual mineralogical and rheological characteristics, and hence have been little hydrothermally altered themselves, have little hydrothermal gold, and have restricted flow of fluids through the rock mass. Some small post-metamorphic quartz veins (metre scale) have been intensely fractured and contain abundant gold on fractures (up to 40 ppm), but these are volumetrically minor. The style of gold mineralisation in the White River area is younger than, and distinctly different from, that of the Klondike area. Some of the mineralised zones in the White River area resemble, mineralogically and geochemically, nearby coeval igneous-hosted gold deposits, but this resemblance is superficial only. The White River mineralisation is an entirely new style of Yukon gold deposit, in which host rocks control the mineralogy and geochemistry of disseminated gold, without quartz veins.     

Application of numerical modelling to a case of compaction-driven dolomitization: a Jurassic palaeohigh in the Po Plain, Italy

Dolomitization of a carbonate platform can occur at different times and in different diagenetic environments, from synsedimentary to deep burial settings. Numerical simulations are valuable tools to test and select the model that, among different hypotheses compatible with field and geochemical data, best honour mass balance, kinetic and thermodynamic constraints. Moreover, the simulation can predict the distribution of the dolomitized bodies in the subsurface and evaluate porosity changes; valuable information for the oil industry. This study is the first attempt to reproduce and investigate the compaction dolomitization model. The diagenetic study of the Jurassic carbonate basin and palaeohigh system of the Po Plain indicates that the carbonates of the palaeohighs were dolomitized by basin compaction fluids. The main goal of the simulations is to evaluate the origin and evolution of the dolomitizing fluids and to provide insights regarding the distribution of the potential reservoir‐dolomitized bodies in the Po Plain. The modelling process is subdivided into two steps: basin modelling and reactive transport modelling. The SEBE3 basin simulator (Eni proprietary) was used to create a three‐dimensional model of the compacting system. The results include compaction fluid flow rate from the basin to the palaeohigh, compaction duration and a determination of the total amount of fluid introduced into the palaeohigh. These data are then used to perform reactive transport modelling with the TOUGHREACT code. Sensitivities on dolomite kinetic parameters suggest that dolomitization was an efficient process even at low temperatures, with differences mainly related to the dynamics of the process. Fluid composition is one of the main constraints, the sea water derived compaction fluid is proven to be efficient for dolomitization due to its relatively high Mg content. Simulations also confirmed that permeability is the most important factor influencing fluid flow and, consequently, the dolomite distribution in the formation. Permeable fractured zones have a strong influence, diverting the dolomitizing fluids from their normal path towards overlying or lateral zones. Moreover, the simulations showed that, after dolomite replacement is complete, the dolomitizing fluids can precipitate dolomite cement, causing over‐dolomitization, with related localized plugging effects in the zone of influx. Mass balance calculations indicate that in the dolomitization compaction model, the amount of compaction water fluxed from the basin to the carbonate is the main constraint on dolomitization efficiency. This observation implies that the ratio between the volume of the basin undergoing compaction and the volume of the palaeohigh is a limiting factor on the final size of the dolomitized bodies. An isolated palaeohigh could be an ideal site for pervasive replacement dolomitization due to the large volume of compaction fluids available compared with the carbonate rock volume. In the case of large platforms, the more permeable margin lithofacies are the most likely sites for compaction model dolomitization. The combined use of a basin simulator and reactive transport modelling has proved to be a successful method to verify model reliability and it provides insights into the volumetric distribution of diagenetic products.     

Geoecology of ecosystem recovery at an inactive coal mine site,New Zealand

Ecological restoration of the Wangaloa coal mine in southern New Zealand is hindered by a range of geoecological factors. The site has some substrate acidification (down to pH 1) and acid rock drainage with discharge waters initially down to pH 4, although this has since risen to ca. pH 6. Surface and ground waters develop elevated sulfate (up to 700 mg/kg) during oxidation of pyrite in coal and quartz in waste rock. Coal has elevated boron content (up to 450 mg/kg) and surface waters on coal-rich waste rock have up to 6 mg/L dissolved boron. Evaporation causes formation of salt encrustations dominated by gypsum with minor boron salts. Boron is bioavailable and may be at toxic levels (>200 mg/kg) in some plants. Quartz-rich waste rock is readily eroded, and develops a cm-scale low-nutrient quartz pebble armouring layer with low water retention properties. All waste rocks including loess siltstone have low nutrient contents, and low moisture retention properties, that are barely sufficient for plant establishment. Native plants introduced to the site during rehabilitation have grown on loess substrate (up to fivefold increase in height over 3 years), with poor or no growth on coal-rich and quartz-rich substrates. In contrast, natural colonisation of manuka (Leptospermum scoparium) has been most effective at revegetation on even the most hostile substrates. This natural revegetation has been facilitated by islands of manuka established accidentally during 60 years of mining history. Manuka from local genetic stock is most viable for this revegetation, and introduced manuka seedlings have had a 70% mortality rate. Natural plant colonisation is the key step in overall ecosystem recovery, and invertebrates rapidly colonise beneath new shrubs irrespective of the nature of the substrate from vegetation islands that have high invertebrate numbers and species richness.     

Exploring the effects of data quality,data worth,and redundancy of CO2 gas pressure and saturation data on reservoir characterization through PEST inversion

This study examined the impacts of reservoir properties on carbon dioxide (CO₂) migration after subsurface injection and evaluated the possibility of characterizing reservoir properties using CO₂ monitoring data such as spatial–temporal distributions of gas pressure, which can be reasonably monitored in practice. The injection reservoir was assumed to be located 1,400–1,500 m below the ground surface such that CO₂ remained in the supercritical state. The reservoir was assumed to contain layers with alternating conductive and resistive properties, which is analogous to actual geological formations such as the Mount Simon Sandstone unit. The CO₂ injection simulation used a cylindrical grid setting in which the injection well was situated at the center of the domain, which extended out 8,000 m from the injection well. The CO₂ migration was simulated using the latest version of the simulator, subsurface transport over multiple phases (the water–salt–CO₂–energy module), developed by Pacific Northwest National Laboratory. A nonlinear parameter estimation and optimization modeling software package, Parameter ESTimation (PEST), is adopted for automated reservoir parameter estimation. The effects of data quality, data worth, and data redundancy were explored regarding the detectability of reservoir parameters using gas pressure monitoring data, by comparing PEST inversion results using data with different levels of noises, various numbers of monitoring wells and locations, and different data collection spacing and temporal sampling intervals. This study yielded insight into the use of CO₂ monitoring data for reservoir characterization and how to design the monitoring system to optimize data worth and reduce data redundancy. The feasibility of using CO₂ saturation data for improving reservoir characterization was also discussed.