Calibrating Environment Canada's MESH Modelling System over the Great Lakes Basin

This paper reports on recent progress towards improved predictions of a land surface-hydrological modelling system, Modélisation Environmentale–Surface et Hydrologie (MESH), via its calibration over the Laurentian Great Lakes Basin. Accordingly, a “global” calibration strategy is utilized in which parameters for all land class types are calibrated simultaneously to a number of sub-basins and then validated in time and space. Model performance was evaluated based on four performance metrics, including the Nash-Sutcliffe (NS) coefficient and simulated compared with observed hydrographs. Results from two calibration approaches indicate that in the model validation mode, the global strategy generates better results than an alternative calibration strategy, referred to as the “individual” strategy, in which parameters are calibrated individually to a single sub-basin with a dominant land type and then validated in a different sub-basin with the same dominant land type. The global calibration strategy was relatively successful despite the large number of calibration parameters (51) and relatively small number of model evaluations (1000) used in the automatic calibration procedure. The NS values for spatial validation range from 0.10 to 0.72 with a median of 0.41 for the 15 sub-basins considered. Results also confirm that a careful model calibration and validation is needed before any application of the model.     

The O’okiep Copper District, Namaqualand, South Africa: a review of the geology with emphasis on the petrogenesis of the cupriferous Koperberg Suite

The O’okiep Copper District is the oldest formal mining area in South Africa. Between 1852 and 2002, the 2,500 km² area yielded two million tons of copper from 32 mines ranging in ore tonnages from 140,000 to 37 million tons. This paper summarizes the calendar of events from the formation of the first primitive crust 1,700–2,000?Ma ago to early Cambrian times ~500?Ma ago, with particular emphasis on the Namaquan (Grenville) Orogeny, notably: the O’okiepian Episode (1,180–1,210?Ma ago) of alpine-type folding, regional granite plutonism, and granulite facies metamorphism and the Klondikean Episode (1,020–1,040?Ma ago) of open and tight folding and the intrusion of the Rietberg Granite and the Koperberg Suite. Almost all of the copper in the O’okiep District occurs in the Koperberg Suite, of which there are 1,700 small bodies that constitute 0.7% of the outcrop area. The suite comprises jotunite, anorthosite, biotite diorite, and hypersthenic rocks ranging from leuconorite to hypersthenite, and it is one of only two world examples of economic copper mineralization in rocks of the anorthosite–charnockite kindred; the second example is Caraiba, Brazil. High I _Sr and low ε _Nd (for a 1,030 Ma intrusion age), and high μ ₂ of 10.1, for Koperberg rock-types indicate a crustal progenitor for the suite, and the presence of jotunite suggests a (subducted) crustal source at ca. 40–50 km depth. The magmatic sulphide paragenesis in the Koperberg Suite is chalcopyrite?+?pyrrhotite (Narrap-type ore) that, in a number of ore-bodies, has been inverted under upper amphibolite facies conditions to bornite?+?Ti-free magnetite (Carolusberg-type ore). Meteoric fluids resulted in supergene Cu enrichment in Koperberg bodies to ~500?m below the pre-Nama peneplane, and lower greenschist facies metamorphism 500–570?Ma ago is reflected by inter alia Hoits-type ore bearing second-generation bornite?+?chalcopyrite(±?covellite?±?chalcocite).     

Fission track thermochronology of the Beni Bousera peridotite massif (Internal Rif,Morocco) and the exhumation of ultramafic rocks in the Gibraltar Arc

The Beni Bousera peridotite massif and its metamorphic surrounding rocks have been analyzed by the fission track (FT) method. The aim was to determine the cooling and uplift history of these mantle and associated crustal rocks after the last major metamorphic event that dates back to the Lower Miocene–Upper Oligocene time (～22–24 Ma). The zircon FT analyses give an average cooling—i.e., below 320 °C—age of ～19.5 Ma. In addition, the apatite FT data give an average cooling—i.e., below 110 °C—age of ～15.5 Ma. Taking into account the thermal properties of the different thermochronological systems used in this work, we have estimated a rate of cooling close to 50 °C/Ma. This cooling rate constrains a denudation rate of about ～2 mm year^?1 from 20 to 15 Ma. These results are similar to those determined in the Ronda peridotite massif of the Betic Cordilleras documenting that some ultrabasic massifs of the internal zones of the two segments of the Gibraltar Arc have a similar evolution. However, Burdigalian sediments occur along the Betic segment (Alozaina area, western Betic segment) unconformably overlying peridotite. At this site, ultramafic rock was exposed to weathering at ages ranging from 20.43 to 15.97 Ma. Since the Beni Bousera peridotite was still at depth until 15.5 Ma, we infer that no simple age projection from massif to massif is possible along the Gibraltar Arc. Moreover, the confined fission track lengths data reveal that a light warming (～100 °C) has reheated the massif during the Late Miocene before the Pliocene–Quaternary tectonic uplift.     

Stable isotope profiles (Ca,O, C) through modern brachiopod shells of T. septentrionalis and G. vitreus: Implications for calcium isotope paleo-ocean chemistry

Brachiopod shells are widely used as an archive to reconstruct elemental and isotopic composition of seawater. Studies, focused on oxygen and carbon isotopes over the last decades, are increasingly extending to the emerging calcium isotope system. To date, only little attention has been paid to test the reliability of fossil brachiopods on their modern counterparts.In this context, the present study investigates two modern brachiopods, Terebratulina septentrionalis (eastern Canada, 5–30 m depth, 7.1 °C seasonal temperature variation, two-layer shell) and Gryphus vitreus (northern Mediterranean, 200 m depth, constant all-year round temperature, three-layer shell). Both species were sampled along the ontogenetic growth direction and calcium, oxygen, and carbon isotopes as well as elemental concentration were measured. Calcium isotopes were analyzed on TIMS. The elemental composition was analyzed by LA-ICP-MS and ICP-AES.The results indicate an intra-specimen δ^44/40Ca variation ranging from 0.16 to 0.33‰, pointing to a fairly homogenous distribution of calcium isotopes in brachiopod shells. However, in the light of the suggested 0.7‰ increase in calcium isotopes over the Phanerozoic such intra-specimen variations constrain ocean reconstruction. δ^44/40Ca values of T. septentrionalis do not seem to be affected by growth rate. Calcium isotopic values of G. vitreus are heavy in the central part of the shell and trend towards lighter values in peripheral areas approaching the maximum isotopic composition of T. septentrionalis. The maximum inter-species δ^44/40Ca difference of 0.62‰ between T. septentrionalis and G. vitreus indicates that care should be taken when using different taxa, species with different strontium content or brachiopods with specialized shell structure, such as G. vitreus, for ocean water reconstruction in terms of Ca isotopic composition. T. septentrionalis may record Ca isotopic fractionation related to seasonal seawater temperature variations in its shell but this is difficult to resolve at the current analytical precision. Average δ¹⁸O-derived temperatures of the two investigated species are close to on-site measured temperatures.     

Dryland Salinity on the Dundas Tableland: A historical appraisal

The groundwater recharge-discharge model, or hydro-equilibrium model, has driven Victorian dryland salinity research for the last 30 years. This paper challenges its application to the Dundas Tableland, south-western Victoria, on historical grounds. In exposing the weak historical foundation of the model for one geographical area so well represented in mainstream investigations, the critique is designed to question the validity of applying the model to other parts of Australia without corresponding historical research. More generally, it is hoped that this contribution from another discipline gives a respectful but significant nudge to the hegemony of physical science in defining and finding solutions for environmental problems.     

Sharp boundary inversion of 2D magnetotelluric data

We consider 2D earth models consisting of laterally variable layers. Boundaries between layers are described by their depths at a set of nodes and interpolated laterally between nodes. Conductivity within each layer is described by values at a set of nodes fixed within each layer, and is interpolated laterally within each layer. Within the set of possible models of this sort, we iteratively invert magnetotelluric data for models minimizing the lateral roughness of the layer boundaries, and the lateral roughness of conductivities within layers, for a given level of data misfit. This stabilizes the inverse problem and avoids superfluous detail. This approach allows the determination of boundary positions between geological units with sharp discontinuities in properties across boundaries, while sharing the stability features of recent smooth conductivity distribution inversions.
We compare sharp boundary inversion results with smooth conductivity distribution inversion results on a numerical example, and on inversion of field data from the Columbia River flood basalts of Washington State. In the synthetic example, where true positions and resistivities are known, sharp boundary inversion results determine both layer boundary locations and layer resistivities accurately. In inversion of Columbia flood basalt data, sharp boundary inversion recovers a model with substantially less internal variation within units, and less ambiguity in both the depth to base of the basalts and depth to resistive basement.     

Anomalous archaeomagnetic directions and site formation processes at archaeological sites in Israel

As part of a program to investigate archaeomagnetic secular variation in Israel and implications for archaeomagnetic dating, we have identified “anomalous” results that have yielded information about site formation processes. Stereonets, vector demagnetization diagrams, and sample location maps are most useful in examining stability of magnetization and consistency of archaeomagnetic directions with each other and with what would be expected from normal secular variation. A collection of examples is presented, including: strong and hard remanence due to vitrification, unstable magnetization due to lightning strikes, superimposed thermal magnetizations, chemical remagnetization, magnetization of a hearth fragment before falling, magnetization of a hot brick after falling, and mechanical deformation of an archaeomagnetic feature. © 1999 John Wiley & Sons, Inc.