Using the potassium-argon laser experiment (KArLE) to date ancient,low-K chondritic meteorites

Several laboratories have been investigating the feasibility of in situ K-Ar dating for use in future landing planetary missions. One drawback of these laboratory demonstrations is the insufficient analogy of the analyzed analog samples with expected future targets. We present the results obtained using the K-Ar laser experiment (KArLE) on two old and K-poor chondritic samples, Pułtusk and Hvittis, as better lunar analogs. The KArLE instrument uses laser ablation to vaporize rock samples and quantifies K content by laser-induced breakdown spectroscopy (LIBS), Ar by quadrupole mass spectrometry (QMS), and ablated mass by laser profilometry. We performed 64 laser ablations on the chondrites to measure spots with a range of K₂O and Ar content and used the data to construct isochrons to determine the chondrite formation age. The KArLE isochron ages on Pułtusk and Hvittis are 5059 ± 892 Ma and 4721 ± 793 Ma, respectively, which is within the uncertainty of published reference ages, and interpreted as the age of their formation. The uncertainty (2σ) on the KArLE ages obtained in this study is better than 20% (18% for Pułtusk and 17% for Hvittis). The precision, which compares our obtained ages to the reference ages, is also better than 20% (11% for Pułtusk and 4% for Hvittis). These results are encouraging for understanding the limits of this technique to measure ancient planetary samples and for guiding future improvements to the instrument.     

Exploring full-waveform LiDAR parameters for tree species classification

Precise tree species classification with high density full-waveform LiDAR data is a key research topic for automated forest inventory. Most approaches constrain to geometric features and only a few consider intensity values. Since full-waveform data offers a much larger amount of deducible information this study explores a high number of parameter and feature combinations. Those variables having the highest impact on species differentiation are determined. To handle the large amount of airborne full-waveform data and to extract a comprehensive number of variable combinations an improved algorithm was developed. The full-waveform point parameters amplitude, width, range corrected intensity and total number of targets within a beam are transferred into raster covering a test site of 10 km². It was possible to isolate the three most important variables based on the intensity, the width and the total number of targets. Up to six tree species were classified with an overall accuracy of 57%, limiting to the four main species accuracy was improved to 78% and constraining just to conifers and broadleaved trees even 91% could be classified correctly.     

Uniformly mantle-like δ<Superscript>18</Superscript>O in zircons from oceanic plagiogranites and gabbros

Lower ocean crust is primarily gabbroic, although 1–2% felsic igneous rocks that are referred to collectively as plagiogranites occur locally. Recent experimental evidence suggests that plagiogranite magmas can form by hydrous partial melting of gabbro triggered by seawater-derived fluids, and thus they may indicate early, high-temperature hydrothermal fluid circulation. To explore seawater–rock interaction prior to and during the genesis of plagiogranite and other late-stage magmas, oxygen-isotope ratios preserved in igneous zircon have been measured by ion microprobe. A total of 197 zircons from 43 plagiogranite, evolved gabbro, and hydrothermally altered fault rock samples have been analyzed. Samples originate primarily from drill core acquired during Ocean Drilling Program and Integrated Ocean Drilling Program operations near the Mid-Atlantic and Southwest Indian Ridges. With the exception of rare, distinctively luminescent rims, all zircons from ocean crust record remarkably uniform δ¹⁸O with an average value of 5.2 ± 0.5‰ (2SD). The average δ¹⁸O(Zrc) would be in magmatic equilibrium with unaltered MORB [δ¹⁸O(WR) ~ 5.6–5.7‰], and is consistent with the previously determined value for equilibrium with the mantle. The narrow range of measured δ¹⁸O values is predicted for zircon crystallization from variable parent melt compositions and temperatures in a closed system, and provides no indication of any interactions between altered rocks or seawater and the evolved parent melts. If plagiogranite forms by hydrous partial melting, the uniform mantle-like δ¹⁸O(Zrc) requires melting and zircon crystallization prior to significant amounts of water–rock interactions that alter the protolith δ¹⁸O. Zircons from ocean crust have been proposed as a tectonic analog for >3.9 Ga detrital zircons from the earliest (Hadean) Earth by multiple workers. However, zircons from ocean crust are readily distinguished geochemically from zircons formed in continental crustal environments. Many of the >3.9 Ga zircons have mildly elevated δ¹⁸O (6.0–7.5‰), but such values have not been identified in any zircons from the large sample suite examined here. The difference in δ¹⁸O, in combination with newly acquired lithium concentrations and published trace element data, clearly shows that the >3.9 Ga detrital zircons did not originate by processes analogous to those in modern mid-ocean ridge settings.     

Speciation and thermodynamic properties for cobalt chloride complexes in hydrothermal fluids at 35-440 °C and 600 bar: An in-situ XAS study

Aqueous Co(II) chloride complexes play a crucial role in cobalt transport and deposition in ore-forming hydrothermal systems, ore processing plants, and in the corrosion of special Co-bearing alloys. Reactive transport modelling of cobalt in hydrothermal fluids relies on the availability of thermodynamic properties for Co complexes over a wide range of temperature, pressure and salinity. Synchrotron X-ray absorption spectroscopy was used to determine the speciation of cobalt(II) in 0-6 m chloride solutions at temperatures between 35 and 440 °C at a constant pressure of 600 bar. Qualitative analysis of XANES spectra shows that octahedral species predominate in solution at 35 °C, while tetrahedral species become increasingly important with increasing temperature. Ab initio XANES calculations and EXAFS analyses suggest that in high temperature solutions the main species at high salinity (Cl:Co >> 2) is CoCl₄²⁻, while a lower order tetrahedral complex, most likely CoCl₂(H₂O)_2(aq), predominates at low salinity (Cl:Co ratios ∼2). EXAFS analyses further revealed the bonding distances for the octahedral Co(H₂O)₆²⁺ (^octCo-O = 2.075(19) Å), tetrahedral CoCl₄²⁻ (^tetCo-Cl = 2.252(19) Å) and tetrahedral CoCl₂(H₂O)_2(aq) (^tetCo-O = 2.038(54) Å and ^tetCo-Cl = 2.210(56) Å). An analysis of the Co(II) speciation in sodium bromide solutions shows a similar trend, with tetrahedral bromide complexes becoming predominant at higher temperature/salinity than in the chloride system. EXAFS analysis confirms that the limiting complex at high bromide concentration at high temperature is CoBr₄²⁻. Finally, XANES spectra were used to derive the thermodynamic properties for the CoCl₄²⁻ and CoCl₂(H₂O)_2(aq) complexes, enabling thermodynamic modelling of cobalt transport in hydrothermal fluids. Solubility calculations show that tetrahedral CoCl₄²⁻ is responsible for transport of cobalt in hydrothermal solutions with moderate chloride concentration (∼2 m NaCl) at temperatures of 250 °C and higher, and both cooling and dilution processes can cause deposition of cobalt from hydrothermal fluids.     

Source, attenuation and potential mobility of arsenic at New Britannia Mine, Snow Lake, Manitoba

This study was to investigate the source, mobility and attenuation of As at the New Britannia Mine, Snow Lake, Manitoba. One major source of As contamination was determined to be an arsenopyrite residue stockpile (ARS) containing refractory Au in a waste rock impoundment. It appears that As is still moving through glacial clay at the base of the ARS into a confined aquifer even though the pile was capped in the year 2000. Arsenic is also being mobilized from a deposit of tailings, which formed following spills by previous owners, Nor Acme. Arsenic from the tailings is being mobilized by oxidation of arsenopyrite and reduction of arsenate to the more mobile arsenite by arsenate-reducing bacteria. This contamination is affecting a shallow unconfined aquifer and surface water flowing from the tailings through wetlands towards Snow Lake. Arsenic is being attenuated by adsorption to hydrated ferric oxides (HFO) in the tailings, wetland soils and aquatic plants. Although As in surface water, soils and plants along the flow path from the mine to Snow Lake are above Canadian drinking water guidelines, efficient natural attenuation by HFO in soils and plants of the wetlands have limited the concentration in Snow Lake to below drinking water standards.     

Paleo-climate of the central European uplands during the last glacial maximum based on glacier mass-balance modeling

During the last glacial maximum (LGM), glaciers existed in scattered mountainous locations in central Europe between the major ice masses of Fennoscandia and the Alps. A positive degree-day glacier mass-balance model is used to constrain paleo-climate conditions associated with reconstructed LGM glacier extents of four central European upland regions: the Vosges Mountains, the Black Forest, the Bavarian Forest, and the Giant Mountains. With reduced precipitation (25–75%), reflecting a drier LGM climate, the modeling yields temperature depressions of 8–15°C. To reproduce past glaciers more severe cooling is required in the west than in the east, indicating a strong west–east temperature anomaly gradient.     

Long-term spatio-temporal hydrochemical and 222Rn tracing to investigate groundwater flow and water–rock interaction in the Gran Sasso (central Italy) carbonate aquifer

In the Gran Sasso fissured carbonate aquifer (central Italy), a long-term (2001–2007) spatio-temporal hydrochemical and ²²²Rn tracing survey was performed with the goal to investigate groundwater flow and water–rock interaction. Analyses of the physico-chemical parameters, and comparisons of multichemical and characteristic ratios in space and time, and subsequent statistical analyses, permitted a characterisation of the hydrogeology. At the regional scale, groundwater flows from recharge areas to the springs located at the aquifer boundaries, with a gradual increase of mineralisation and temperature along its flowpaths. However, the parameters of each group of springs may significantly deviate from the regional trend owing to fast flows and to the geological setting of the discharge spring areas, as corroborated by statistical data. Along regional flowpaths, the effects of seasonal recharge and lowering of the water table clearly cause changes in ion concentrations over time. This conceptual model was validated by an analysis of the ²²²Rn content in groundwater. ²²²Rn content, for which temporal variability depends on seasonal fluctuations of the water table, local lithology and the fracture network at the spring discharge areas, was considered as a tracer of the final stages of groundwater flowpaths.     

The Vestan cataclysm: Impact‐melt clasts in howardites and the bombardment history of 4 Vesta

Crystalline impact‐melt samples were created in high‐temperature environments by relatively large craters and, as such, give additional constraints on the nature of the impacts that created them. This article provides new ⁴⁰Ar‐³⁹Ar ages of impact‐melt clasts in howardites and shows that these clasts formed on the HED parent body, 4 Vesta, within the time period 3.3–3.8 Ga. Rather than resulting from an increased number of impacts, however, impact‐melted material in howardites may result from unusually high‐velocity impacts occurring in the asteroid belt during this period. This scenario is similar to the late heavy bombardment of the Moon, pointing to an unusual dynamical event at this time across the inner solar system. Therefore, impact‐melt rocks in howardites uniquely record a Vestan cataclysm.     

The interplay between basement fabric,rifting, syn-rift folding,and inversion in the Rio do Peixe Basin,NE Brazil

The Rio do Peixe Basin is part of a series of aborted Cretaceous rifts formed within the Proterozoic Borborema Province, onshore NE Brazil in response to rifting between Africa and South America. The basin is remarkably well-imaged and comprises fault-bounded depocentres, the main ones being the NE-oriented Brejo das Freiras and the E–W-oriented Sousa half-grabens. These grabens and their bounding faults are influenced by Neoproterozoic basement shear zones and present a complex framework of secondary normal faults and folds. Recent workers also interpret large reverse faults and regional post-rift shortening driven by far-field stresses from the Andes. For those reasons, the basin represents an ideal setting to investigate the multiphase history of rift basins. We thus combine borehole-calibrated 2D and 3D seismic and magnetic data with section restoration and numerical modelling to investigate the architecture and evolution of this basin. We aim to understand: (i) the controls of the basement fabric in 3D fault architecture and kinematics and (ii) how syn-rift faults controlled the geometry and development of fault-related folds. By doing this, we also investigate the timing, kinematics, and magnitude of inversion in the basin to explore its multiphase history. We demonstrate that (i) the basement fabric controlled not only the strike of faults but also their geometry and polarity at depth, (ii) folds in the syn-rift sequence are attributed simply to syn-rift extension along stepped and/or curved faults, and (iii) inverted and/or reverse faults occur within the basin, but these are minor and appear to have formed during rifting. We explain this minor inversion by a change in plate kinematics related to the onset of the nearby transform margin to the north. These results have implications for understanding the 3D evolution of oblique grabens, the role of structural inheritance, and the recognition of inversion- versus extension-related folds in rift basins worldwide.     

Automated Extraction of Natural Drainage Density Patterns for the Conterminous United States through High-Performance Computing

Hydrographic networks form an important data foundation for cartographic base mapping and for hydrologic analysis. Drainage density patterns for these networks can be derived to characterize local landscape, bedrock and climate conditions, and further inform hydrological and geomorphological analysis by indicating areas where too few headwater channels are represented. Natural drainage density patterns are not consistently available in existing hydrographical data bases for the United States because compilation and capture criteria historically varied, along with climate, during the period of data collection over the various terrain types throughout the country. This paper demonstrates an automated workflow that is being tested in a high-performance computing environment by the U.S. Geological Survey (USGS) to map natural drainage density patterns at the 1:24,000-scale (24K) for the conterminous United States. Hydrographic network drainage patterns may be extracted from elevation data to guide corrections for existing hydrographic network data. The paper describes three stages in this workflow including data pre-processing, natural channel extraction, and generation of drainage density patterns from extracted channels. The workflow is implemented in parallel fashion by simultaneously executing procedures on multiple subbasin watersheds within the U.S. National Hydrography Dataset (NHD). Pre-processing defines parameters needed for the extraction process. Extraction proceeds in standard fashion: filling sinks, developing flow direction and weighted flow accumulation rasters. Drainage channels with assigned Strahler stream order are extracted within a subbasin and simplified. Drainage density patterns are then estimated with 100-m resolution and subsequently smoothed with a low-pass filter. The extraction process is found to be of better quality in higher slope terrains. Concurrent processing through the high-performance computing environment is shown to facilitate and refine the choice of drainage density extraction parameters and more readily improve extraction procedures than conventional processing.