Sedimentary basin formation associated with a silicic large igneous province: stratigraphy and provenance of the Mesoproterozoic Roopena Basin,Gawler Range Volcanics

The Mesoproterozoic Gawler Silicic Large Igneous Province (SLIP) in the Gawler Craton and Curnamona Province, southern Australia, comprises extensive felsic and lesser mafic volcanic sequences, with only limited sedimentary successions. The Roopena Basin is a rare example of a synvolcanic sedimentary basin that formed within the Gawler SLIP in the eastern Gawler Craton. It is a north–south-trending basin with a preserved area of 75 km², bound by the Roopena and Wizzo Well faults, and contains three units of the lower Gawler Range Volcanics; the Angle Dam Dacite, Fresh Well Formation and Roopena Basalt. The Angle Dam Dacite is a porphyritic lava and the oldest part of the volcanic succession, directly overlying basement. The Fresh Well Formation overlies the Angle Dam Dacite or basement, comprises three coarsening-upwards volcaniclastic packages of claystone, siltstone, fine-grained to coarse-grained lithic sandstone and conglomerate deposited in a fluvio-lacustrine setting, and contains three tuff members. The Roopena Basalt is interlayered with the Fresh Well Formation, and comprises auto-brecciated lavas that exhibit only local interaction with water or wet sediment. Sharp basal contacts of the prograding packages within the Fresh Well Formation provide evidence of rapid flooding events within the basin. New detrital zircon geochronology of a sandstone within the Fresh Well Formation yielded a maximum depositional age of ca 1580 Ma, with provenance dominated by felsic volcanic units of the 1635–1605 Ma St Peter Suite. Sedimentation in the Gawler SLIP appears to have occurred in isolated basins with limited areal extent. It was largely restricted to the eastern Gawler Craton, and as well as the Roopena Basin, and includes similar basins at the Olympic Dam and Prominent Hill iron oxide–copper–gold ± uranium (IOCG ± U) deposits. The coincidence of sedimentation and mafic volcanism in the eastern Gawler Craton suggests that this region underwent extension at this time, although high-temperature metamorphism and compressional deformation occurred in some parts of the Gawler Craton and Curnamona Province synchronous with the Gawler SLIP. The Roopena Basin sedimentary rocks and underlying basement contain hematite–chlorite–sericite–white mica assemblages, permissive of hematite-style IOCG mineral deposits; however, no significant ore deposit has yet been discovered in the Roopena Basin.     

Relationships between magmatism and deformation in northern Yorke Peninsula and southeastern Proterozoic Australia

The ca 1600–1580 Ma time interval is recognised as a significant period of magmatism, deformation and mineralisation throughout eastern Proterozoic Australia. Within the northern Yorke Peninsula in South Australia, this period was associated with the emplacement of multiple phases of the Tickera Granite, an intensely foliated quartz alkali-feldspar syenite, a leucotonalite and an alkali-feldspar granite. These granites belong to the broader Hiltaba Suite that was emplaced at shallow crustal levels throughout the Gawler Craton. Geochemical and isotopic analysis suggests these granite phases were derived from a heterogeneous source region. The syenite and alkali-feldspar granite were derived from similar source regions, likely the underlying ca 1850 Ma Donington Suite and/or the ca 1750 Ma Wallaroo Group metasediments with some contamination from an Archean basement. The leucotonalite is sourced from a similar but more mafic/lower crustal source. Phases of the Tickera Granite were emplaced synchronously with deformation that resulted in development of a prominent northeast-trending structural grain throughout the Yorke Peninsula region. This fabric is associated with composite events resulting from folding, shearing and faulting within the region. The intense deformation and intrusion of granites within this period resulted in mineralisation throughout the region, as seen in Wheal Hughes and Poona mines. The Yorke Peninsula shares a common geological history with the Curnamona Province, which was deformed during the ca 1600–1585 Ma Olarian Orogeny, and resulted in development of early isoclinal and recumbent folds overprinted by an upright fold generation, a dominant northeast-trending structural grain, mineralisation, and spatially and temporally related intrusions. This suggests correlation of parts of the Gawler Craton with the Curnamona Province, and that the Olarian Orogeny also affected the southeastern Gawler Craton.     

Hydrogeochemical exploration for volcanic-hosted massive sulfide deposits in semi-arid Australia

Research on hydrogeochemistry for mineral exploration for inland Australia includes development of weathering models and extensive mine-scale and regional groundwater data. Mineral saturation indices for groundwater, activity–activity plots and reaction modelling simulate weathering of volcanic-hosted massive sulfide (VHMS) deposits in deeply weathered environments. At 10 m or more below surface, dissolved O₂ is very low and other solutes such as sulfate, carbonate and nitrate are more likely oxidants. Modelling indicates that these processes differ from oxic weathering of highly eroded terrains, and provide the framework to develop robust hydrogeochemical exploration procedures in covered terrains. Sulfide weathering potentially occurs in two or more phases that effect surrounding groundwaters in differing manners. Deeper oxidative alteration of sulfides (e.g. bornite to chalcopyrite), occurring tens to hundreds of metres below surface, uses sulfate and carbonate as oxidants, causing neutral to alkaline conditions. In this zone, only pyritic massive sulfides potentially generate acidic conditions. Thus, deep sulfide-rich rocks are indicated by sulfate-depleted groundwater. Closer to the surface, sulfides are oxidised to soluble sulfates by dissolved nitrate, with much less acid production than if dissolved oxygen was the main oxidant. Thus, in shallow groundwater, sulfides are indicated by sulfate enrichment and nitrate depletion. Elements are released from sulfides and wall rocks by acid or alkaline conditions. The derived FeS (pH–Eh + Fe + Mn) and AcidS (Li + Mo + Ba + Al) indices distinguish sulfide systems through tens of metres of cover. VHMS systems are distinguished from other non-economic sulfide deposits where there is little transported cover, using various dissolved elements, including Zn, Pb and Cu. Elsewhere, ‘patchiness’ and limited aerial extent of metal signals are due to adsorption effects, that intensify with depth. Other elements such as Mn and Co have lesser diminution effects, but are less selective indicators for VHMS. There is exploration potential for elements such as Pt or Ag. These varying sulfide indicators have moderate utility, even for large-scale (～5 km spacing) sampling. Results indicate that hydrogeochemistry can add value to greenfields exploration for VHMS ore deposits in deeply weathered terrains. It is also moderately successful at indicating the presence of sulfide-rich systems (whether magmatic or hydrothermal) under >100 m cover, thus providing a rapid and cost-effective regional prospectivity tool for deeply buried terrains. Such mineral exploration tools will encourage exploration investment for more difficult regions of Australia and in other deeply weathered regions of the world.     

An apatite U-Pb thermal history map for the northern Gawler Craton,South Australia

Apatite U-Pb thermochronology was applied to granitoid basement samples across the northern Gawler Craton to unravel the Proterozoic, post-orogenic, cooling history and to examine the role of major fault zones during cooling. Our observations indicate that cooling following the ~2500 Ma Sleaford Orogeny and ~1700 Ma Kimban Orogeny is restricted to the Christie and Wilgena Domains of the central northern Gawler Craton. The northern Gawler Craton mainly records post-Hiltaba Event(~1590 Ma) U-Pb cooling ages. Cooling following the ~1560 Ma Kararan Orogeny is preserved within the Coober Pedy Ridge,Nawa Domain and along major shear zones within the south-western Fowler Domain. The Nawa Domain samples preserve U-Pb cooling ages that are 150 Ma younger than the samples within the Coober Pedy Ridge and Fowler Domain, indicating that later(~1300 Ma) fault movement within the Nawa Domain facilitated cooling of these samples, caused by arc collision in the Madura Province of eastern Western Australia. When compared to~(40)Ar/~(39) Ar from muscovite, biotite and hornblende, our new apatite U-Pb ages correlate well, particularly in regions of higher data density. Our data also preserve a progressive younging of U-Pb ages from the nucleus of the craton to the periphery with a stark contrast in U-Pb ages across major structures such as the Karari Shear Zone and the Southern Overthrust, which indicates the timing of reactivation of these major crustal structures. Although this interpolation was based solely on thermochronological data and did not take into account structural or other geological data, these maps are consistent with the structural architecture of the Gawler Craton and reveal the thermal footprint of known tectonic and magmatic events in the Gawler Craton.     

Retreat dynamics of the eastern sector of the British–Irish Ice Sheet during the last glaciation

The findings of BRITICE-CHRONO Transect 2 through the North Sea Basin and eastern England are reported. We define ice-sheet marginal oscillation between ~31 and 16 ka, with seven distinctive former ice-sheet limits (L1–7) constrained by Bayesian statistical analysis. The southernmost limit of the North Sea Lobe is recorded by the Bolders Bank Formation (L1; 25.8–24.6 ka). L2 represents ice-sheet oscillation and early retreat to the northern edge of the Dogger Bank (23.5–22.2 ka), with the Garret Hill Moraine in north Norfolk recording a significant regional readvance to L3 at 21.5–20.8 ka. Ice-marginal oscillations at ~26–21 ka resulted in L1, L2 and L3 being partially to totally overprinted. Ice-dammed lakes related to L1–3, including Lake Humber, are dated at 24.1–22.3 ka. Ice-sheet oscillation and retreat from L4 to L5 occurred between 19.7 and 17.3 ka, with grounding zone wedges marking an important transition from terrestrial to marine tidewater conditions, triggered by the opening of the Dogger Lake spillway between 19.9 and 17.5 ka. L6 relates to ice retreat under glacimarine conditions and final ice retreat into the Firth of Forth by 15.8 ka. L7 (~15 ka) represents an ice retreat from Bosies Bank into the Moray Firth.     

Prospective modelling of 3D hyporheic exchange based on high‐resolution topography and stream elevation

River managers and scientists interested in hyporheic processes need adequate tools for characterizing hyporheic exchange flow (HEF) at local sites where only poor information on subsurface properties are available. This study evaluates a three‐dimensional modelling approach, on the basis of detailed surface parameterization and a simplified subsurface structure, for comparison of potential HEF characteristics at three experimental reaches at the channel‐unit scale. First, calibration is conducted to determine the best fit‐of‐heads given the model simplification, then the structure of residuals are used to evaluate the origin of the misfit, and finally, a sensitivity analysis is conducted to identify inter‐site differences in HEF. Results show that such an approach can highlight potential magnitude differences in HEF characteristics between reaches. The sensitivity analysis is successful in delineating the small area of exchange that remains under conditions of high groundwater discharge. In this case, however, the calibrated model performs poorly in representing the exchange pattern at the sediment–water interface, thus suggesting that the approach is less adequate for a deterministic simulation of observed heads. The summary statistics are in the range of similar published models, for which the reported indicator is the root mean square error on heads normalized by the head drop over the reach. We recommend, however, that modellers use a more comparable indicator, such as a measure of the residuals normalized by a measure of observed vertical head differences. Overall, when subsurface data are unavailable or sparse, a three‐dimensional groundwater model based on high‐resolution topographic data combined with a sensitivity analysis appears as a useful tool for a preliminary characterization of HEF. Copyright © 2013 John Wiley & Sons, Ltd.     

Up‐scaling potential impacts on water flows from agricultural water interventions: opportunities and trade‐offs in the Osman Sagar catchment,Musi sub‐basin,India

Agricultural water management (AWM) has been shown to improve and secure yields in the tropics and has been suggested as an important way to combat poverty in the region. In this paper, we describe potential impacts on upstream and downstream flows of extensive AWM interventions, using the watershed development programme of the Osman Sagar catchment of Musi sub‐basin, Andhra Pradesh semi‐arid India, as an example. Various AWM interventions are compared with a non‐intervention state and the current state of the study area, using 31 years of data by application of the calibrated and validated ARCSWAT 2005 (Version 2.1.4a) modelling tool. Different AWM interventions contribute to improved livelihoods of upstream smallholder farmers by increasing soil moisture availability and groundwater recharge, which can subsequently be used for irrigation. The result is higher crop production and hence larger incomes. Moreover, lower flow intensities and sediment losses reduced by 30–50%, reducing the risk of flooding and sediment accumulation in the Osman Sagar drinking water reservoir. On the other hand, AWM interventions are predicted to result in reduced total water inflows to the Osman Sagar reservoir from 11% of the total annual rainfall (754 mm) recorded at present, to 8% if AWM interventions were implemented at large scale throughout the catchment. A cost–benefit analysis of AWM interventions showed that the highest net economic returns were achieved at intermediate intervention levels (only in‐situ AWM). Copyright © 2012 John Wiley & Sons, Ltd.     

Sea ice as an evaluation target for HEM modelling and inversion

We investigate helicopter electromagnetic (HEM) inversion schemes applied to synthetic and measured HEM sea ice profiling data. Direct HEM data-to-ice-thickness inversion is compared to three different formal, least squares layered earth inversion algorithms.By making several approximations, it is possible to directly invert a single channel measurement (i.e., the in-phase or quadrature component of a single frequency measurement) to obtain an estimate of sea ice thickness. Measurements from multiple input channels, however, can be used in a layered earth inversion to simultaneously recover several model parameters such as sea ice thickness, sea ice conductivity and sub-ice bathymetry. Synthetic data sets for a particular two-frequency HEM system showed that simple least squares inversion algorithms produce reliable estimates of sea ice thickness in cases where the ice is thicker than 3 m. These methods could also recover acceptable estimates of sea ice thickness when a thin, conductive, partially melted sea ice layer was present, and could determine shallow, sub-ice bathymetry in brackish water. As expected, 1D transformations and inversions of synthetic data for a three-dimensional pressure ridge keel structure contained artifacts, notably broadening of the apparent width of the keel.Prior to inverting a field data set acquired over rather thin (~ 0.5 m) Antarctic sea ice, we found it necessary to recalibrate the phase angle of the measurements using a phasor diagram-based method. Direct transformation of a single channel from the recalibrated data set produced more accurate estimates of sea-ice thickness than formal inversion of multi-channel data. We suggest that the least squares inversion methods are inferior in this situation because of the particular characteristics of the two-frequency HEM system used in this evaluation; the extreme differences in sensitivity of high and low frequency data components, the overall low sensitivity to sea ice conductivity (especially for thin ice), and the partially low signal-to noise ratios of the measurements. The data sets used in this study will be made available to the public to allow alternate inversion approaches to be applied and evaluated. It is suggested that inclusion of parameter bounds and other forms of regularization could help to improve the inversion results.     

Air quality planning in California’s changing climate

California is home to some of the worst air quality in the nation and ninety percent of the state’s population lives in areas that are out of attainment with at least one of the National Ambient Air Quality Standards. Increasing temperatures associated with climate change will make meeting air quality standards more difficult. Under a changing climate, additional emission reductions will be needed to achieve clean air standards. These additional emission reductions and associated costs are called the “climate penalty.” Air quality planning is the process of assessing the emission reductions needed to meet air quality standards and outlining the programs and policies that will be implemented to achieve these emission reductions. This paper reviews the challenges that a changing climate will pose for air quality planning in California and identifies opportunities for adaptation. While state air quality regulators in California are taking enormous strides to address global warming, less work is happening at the regional, air district level. Air districts are the agencies responsible for developing air quality improvement plans. An important first step for regional air quality regulators will be to quantify the climate penalty and understand their region’s vulnerability to climate change. Limitations in regulatory authority could impede measures to improve preparedness. Regional agencies will likely need to look to state and federal agencies for additional emission reductions.