Simulating the diurnal cycle of rainfall in global climate models: resolution versus parameterization

The effects of horizontal resolution and the treatment of convection on simulation of the diurnal cycle of precipitation during boreal summer are analyzed in several innovative weather and climate model integrations. The simulations include: season-long integrations of the Non-hydrostatic Icosahedral Atmospheric Model (NICAM) with explicit clouds and convection; year-long integrations of the operational Integrated Forecast System (IFS) from the European Centre for Medium-range Weather Forecasts at three resolutions (125, 39 and 16 km); seasonal simulations of the same model at 10 km resolution; and seasonal simulations of the National Center for Atmospheric Research (NCAR) low-resolution climate model with and without an embedded two-dimensional cloud-resolving model in each grid box. NICAM with explicit convection simulates best the phase of the diurnal cycle, as well as many regional features such as rainfall triggered by advancing sea breezes or high topography. However, NICAM greatly overestimates mean rainfall and the magnitude of the diurnal cycle. Introduction of an embedded cloud model within the NCAR model significantly improves global statistics of the seasonal mean and diurnal cycle of rainfall, as well as many regional features. However, errors often remain larger than for the other higher-resolution models. Increasing resolution alone has little impact on the timing of daily rainfall in IFS with parameterized convection, yet the amplitude of the diurnal cycle does improve along with the representation of mean rainfall. Variations during the day in atmospheric prognostic fields appear quite similar among models, suggesting that the distinctive treatments of model physics account for the differences in representing the diurnal cycle of precipitation.     

Correction to: Development and Evaluation of Geostatistical Methods for Non-Euclidean-Based Spatial Covariance Matrices

Mathematical Geosciences - The original version of this article unfortunately contained a mistake in equation 9.     

Irrigation produces elevated arsenic in the underlying groundwater of a semi-arid basin in Southwestern Idaho

The shallow aquifer beneath the Western Snake River Plain (Idaho, USA) exhibits widespread elevated arsenic concentrations (up to 120 μg L⁻¹). While semi-arid, crop irrigation has increased annual recharge to the aquifer from approximately 1 cm prior to a current rate of >50 cm year⁻¹. The highest aqueous arsenic concentrations are found in proximity to the water table (all values >50 μg L⁻¹ within 50 m) and concentrations decline with depth. Despite strong vertical redox stratification within the aquifer, spatial distribution of aqueous species indicates that redox processes are not primary drivers of arsenic mobilization. Arsenic release and transport occur under oxidizing conditions; groundwater wells containing dissolved arsenic at >50 μg L⁻¹ exhibit elevated concentrations of O₂ (average 4 mg L⁻¹) and NO₃ (average 8 mg L⁻¹) and low concentrations of dissolved Fe (<20 μg L⁻¹). Sequential extractions and spectroscopic analysis of surficial soils and sediments indicate solid phase arsenic is primarily arsenate and is present at elevated concentrations (4–45 mg kg⁻¹, average: 17 mg kg⁻¹) relative to global sedimentary abundances. The highest concentrations of easily mobilized arsenic (up to 7 mg kg⁻¹) are associated with surficial soils and sediments visibly stained with iron oxides. Batch leaching experiments on these materials using irrigation waters produce pore water arsenic concentrations approximating those observed in the shallow aquifer (up to 152 μg L⁻¹). While As:Cl aqueous phase relationships suggest minor evaporative enrichment, this appears to be a relic of the pre-irrigation environment. Collectively, these data indicate that infiltrating irrigation waters leach arsenic from surficial sediments to the underlying aquifer.     

Plesiosaur remains from the Jurassic-Cretaceous Purbeck Limestone Group of southern England

Plesiosaur fossils from the uppermost Tithonian (Upper Jurassic) to lower Valanginian (Lower Cretaceous) Purbeck Limestone Group of Dorset, England are important because they constitute a rare record from around the Jurassic-Cretaceous boundary. Where diagnostic, the remains are referred to indeterminate plesiosauroids (Cryptocleidoidea) and represent relatively small animals including a previously undescribed juvenile individual. This observation accords with the lagoonal-lacustrine (with periodic marine influx) depositional setting of the Purbeck Limestone Group strata, a palaeoenvironment that is often associated with juvenile plesiosaur occurrences.     

Climatic and lithologic controls on the temporal and spatial variability of CO2 consumption via chemical weathering: An example from the Australian Victorian Alps

A detailed geochemical study on river waters of the Australian Victorian Alps was carried out to determine: (i) the relative significance of silicate, carbonate, evaporite and sulfide weathering in controlling the major ion composition and; (ii) the factors regulating seasonal and spatial variations of CO₂ consumption via silicate weathering in the catchments. Major ion chemistry implies that solutes are largely derived from evaporation of precipitation and chemical weathering of carbonate and silicate lithologies. The input of solutes from rock weathering was determined by calculating the contribution of halite dissolution and atmospheric inputs using local rain and snow samples. Despite the lack of carbonate outcrops in the study area and waters being undersaturated with respect to calcite, the dissolution of vein calcite accounts for up to 67% of the total dissolved cations, generating up to 90% of dissolved Ca and 97% of Mg. Dissolved sulfate has δ³⁴S values of 16 to 20‰_CDT, indicating that it is derived predominantly from atmospheric deposition and minor gypsum weathering and not from bacterial reduction of FeS₂. This militates against sulphuric acid weathering in Victorian rivers. Ratios of Si vs. the atmospheric corrected Na and K concentrations range from ~ 1.1 to ~ 4.3, suggesting incongruent weathering from plagioclase to smectite, kaolinite and gibbsite.Estimated long-term average CO₂ fluxes from silicate weathering range from ~ 0.012 × 10⁶ to 0.039 × 10⁶ mol/km²/yr with the highest values in rivers draining the basement outcrops rather than sedimentary rocks. This is about one order of magnitude below the global average which is due to low relief, and the arid climate in that region. Time series measurements show that exposure to lithology, high physical erosion and long water–rock contact times dominate CO₂ consumption fluxes via silicate weathering, while variations in water temperature are not overriding parameters controlling chemical weathering. Because the atmospheric corrected concentrations of Na, K and Mg act non-conservative in Victorian rivers the parameterizations of weathering processes, and net CO₂ consumption rates in particular, based on major ion abundances, should be treated with skepticism.     

Trace element mapping by LA-ICP-MS: assessing geochemical mobility in garnet

A persistent problem in the study of garnet geochemistry is that the consideration of major elements alone excludes a wealth of information preserved by trace elements, particularly the rare-earth elements (REEs). This is despite the fact that trace elements are generally less vulnerable to diffusive resetting, and are sensitive to a broader spectrum of geochemical interactions involving the entire mineral assemblage, including the growth and/or dissolution of accessory minerals. We outline a technique for the routine acquisition of high-resolution 2D trace element maps by LA-ICP-MS, and introduce an extension of the software package XMapTools for rapid processing of LA-ICP-MS data to visualise and interpret compositional zoning patterns. These methods form the basis for investigating the mechanisms controlling geochemical mobility in garnet, which are argued to be largely dependent on the interplay between element fractionation, mineral reactions and partitioning, and the length scales of intergranular transport. Samples from the Peaked Hill shear zone, Reynolds Range, central Australia, exhibit contrasting trace element distributions that can be linked to a detailed sequence of growth and dissolution events. Trace element mapping is thus employed to place garnet evolution in a specific paragenetic context and derive absolute age information by integration with existing U–Pb monazite and Sm–Nd garnet geochronology. Ultimately, the remarkable preservation of original growth zoning and its subtle modification by subsequent re-equilibration is used to ‘see through’ multiple superimposed events, thereby revealing a previously obscure petrological and temporal record of metamorphism, metasomatism, and deformation.     

Granitoids of the Magba shear zone,West Cameroon,Central Africa: Evidences for emplacement under transpressive tectonic regime

The Magba Shear Zone is made up of granites, migmatites, orthogneiss, metagabbro, mafic dyke and mylonites with coarse grained texture, porphyroblastic, granoblastic, cataclastic and mylonitic texture respectively. Structural features and kinematic indicators testify the syntectonic emplacement of Magba granitoids and also provide detailed information on the relative timing of deformation as follows: (1) D1 of tangential movement immediately followed by (2) the D2 phase which is heterogeneous simple shear in dextral transpressive context with a NW-SE direction (3) D3 tectonic phase is marked by sinistral transpressive tectonic and superposed folding with a NE-SW kinematic direction. Combined ductile NE-SW shear movements and NWSE compressional movements defined a transpressional tectonic regime during the D3 deformation (4) A brittle stage D4 is controlled by transcurrent tectonics and responsible for the emplacement of faults, and joints. The Magba granites would have intruded along sub-vertical mid-crustal feeder channels and were emplaced as a sheet or sheets along the shear zone during the early stage of the C3 shearing, followed by gabbro and mafic dyke at the late stage. Strike-slip dilatancy pumping under transpressive tectonic is suggested as a possible mechanism for the emplacement of the Magba granites.     

Preferential flow,diffuse flow,and perching in an interbedded fractured-rock unsaturated zone

Layers of strong geologic contrast within the unsaturated zone can control recharge and contaminant transport to underlying aquifers. Slow diffuse flow in certain geologic layers, and rapid preferential flow in others, complicates the prediction of vertical and lateral fluxes. A simple model is presented, designed to use limited geological site information to predict these critical subsurface processes in response to a sustained infiltration source. The model is developed and tested using site-specific information from the Idaho National Laboratory in the Eastern Snake River Plain (ESRP), USA, where there are natural and anthropogenic sources of high-volume infiltration from floods, spills, leaks, wastewater disposal, retention ponds, and hydrologic field experiments. The thick unsaturated zone overlying the ESRP aquifer is a good example of a sharply stratified unsaturated zone. Sedimentary interbeds are interspersed between massive and fractured basalt units. The combination of surficial sediments, basalts, and interbeds determines the water fluxes through the variably saturated subsurface. Interbeds are generally less conductive, sometimes causing perched water to collect above them. The model successfully predicts the volume and extent of perching and approximates vertical travel times during events that generate high fluxes from the land surface. These developments are applicable to sites having a thick, geologically complex unsaturated zone of substantial thickness in which preferential and diffuse flow, and perching of percolated water, are important to contaminant transport or aquifer recharge.     

ISCO for Groundwater Remediation: Analysis of Field Applications and Performance

A critical analysis of in situ chemical oxidation (ISCO) projects was performed to characterize situations in which ISCO is being implemented, how design and operating parameters are typically employed, and to determine the performance results being achieved. This research involved design of a database, acquisition and review of ISCO project information, population of the database, and analyses of the database using statistical methods. Based on 242 ISCO projects included in the database, ISCO has been used to treat a variety of contaminants; however, chlorinated solvents are by far the most common. ISCO has been implemented at sites with varied subsurface conditions with vertical injection wells and direct push probes being the most common delivery methods. ISCO has met and maintained concentrations below maximum contaminant levels (MCLs), although not at any sites where dense nonaqueous phase liquids (DNAPL) were presumed to be present. Alternative cleanup levels and mass reduction goals have also been attempted, and these less stringent goals are met with greater frequency than MCLs. The use of pilot testing is beneficial in heterogeneous geologic media, but not so in homogeneous media. ISCO projects cost $220,000 on average, and cost on average $94/yd³ of target treatment zone. ISCO costs vary widely based on the size of the treatment zone, the presence of DNAPL, and the oxidant delivery method. No case studies were encountered in which ISCO resulted in permanent reductions to microbial populations or sustained increases in metal concentrations in groundwater at the ISCO-treated site.