Magnetic unmixing of first-order reversal curve diagrams using principal component analysis

We describe a quantitative magnetic unmixing method based on principal component analysis (PCA) of first-order reversal curve (FORC) diagrams. For PCA we resample FORC distributions on grids that capture diagnostic signatures of single-domain (SD), pseudo-single-domain (PSD), and multidomain (MD) magnetite, as well as of minerals such as hematite. Individual FORC diagrams are recast as linear combinations of end-member (EM) FORC diagrams, located at user-defined positions in PCA space. The EM selection is guided by constraints derived from physical modeling and imposed by data scatter. We investigate temporal variations of two EMs in bulk North Atlantic sediment cores collected from the Rockall Trough and the Iberian Continental Margin. Sediments from each site contain a mixture of magnetosomes and granulometrically distinct detrital magnetite. We also quantify the spatial variation of three EM components (a coarse silt-sized MD component, a fine silt-sized PSD component, and a mixed clay-sized component containing both SD magnetite and hematite) in surficial sediments along the flow path of the North Atlantic Deep Water (NADW). These samples were separated into granulometric fractions, which helped constrain EM definition. PCA-based unmixing reveals systematic variations in EM relative abundance as a function of distance along NADW flow. Finally, we apply PCA to the combined dataset of Rockall Trough and NADW sediments, which can be recast as a four-EM mixture, providing enhanced discrimination between components. Our method forms the foundation of a general solution to the problem of unmixing multi-component magnetic mixtures, a fundamental task of rock magnetic studies.The research leading to these results has received funding from the European Research Council under the European Union's Seventh Framework Programme (FP/2007-2013)/ERC grant agreement 320750.This is the author accepted manuscript. The final version is available from Wiley via http://dx.doi.org/10.1002/2015GC00590

A geochemical study of the impact of irrigation and aquifer lithology on groundwater in the Upper Yakima River Basin, Washington, USA

The Yakima River, a major tributary of the Columbia River, is currently overallocated in its surface water usage in part because of large agricultural water use. As a result, groundwater availability and surface water/groundwater interactions have become an important issue in this area. In several sub-basins, the Yakima River water is diverted and applied liberally to fields in the summer creating artificial recharge of shallow groundwater. Major ion, trace element, and stable isotope geochemistry of samples from 26 groundwater wells from a transect across the Yakima River and 24 surface waters in the Kittitas sub-basin were used to delineate waters with similar geochemical signatures and to identify surface water influence on groundwater. Major ion chemistry and stable isotope signatures combined with principal component analysis revealed four major hydrochemical groups. One of these groups, collected from shallow wells within the sedimentary basin fill, displays temporal variations in NO3 and SO4 along with high δ18O and δD values, indicating significant contribution from Yakima River and/or irrigation water. Two other major hydrochemical groups reflect interaction with the main aquifer lithologies in the basin: the Columbia River basalts (high-Na groundwaters), and the volcaniclastic rocks of the Ellensburg Formation (Ca–Mg–HCO3 type waters). The fourth major group has interacted with the volcaniclastic rocks and is influenced to a lesser degree by surface waters. The geochemical groupings constrain a conceptual model for groundwater flow that includes movement of water between underlying Columbia River basalt and deeper sedimentary basin fill and seasonal input of irrigation water