Water-rock interactions drive chemostasis

The western U.S. is experiencing shifts in recharge due to climate change, and it is currently unclear how hydrologic shifts will impact geochemical weathering and stream concentration–discharge (C–Q) patterns. Hydrologists often use C–Q analyses to assess feedbacks between stream discharge and geochemistry, given abundant stream discharge and chemistry data. Chemostasis is commonly observed, indicating that geochemical controls, rather than changes in discharge, are shaping stream C–Q patterns. However, few C–Q studies investigate how geochemical reactions evolve along groundwater flowpaths before groundwater contributes to streamflow, resulting in potential omission of important C–Q controls such as coupled mineral dissolution and clay precipitation and subsequent cation exchange. Here, we use field observations—including groundwater age, stream discharge, and stream and groundwater chemistry—to analyse C–Q relations in the Manitou Experimental Forest in the Colorado Front Range, USA, a site where chemostasis is observed. We combine field data with laboratory analyses of whole rock and clay x-ray diffraction and soil cation-extraction experiments to investigate the role that clays play in influencing stream chemistry. We use Geochemist's Workbench to identify geochemical reactions driving stream chemistry and subsequently suggest how climate change will impact stream C–Q trends. We show that as groundwater age increases, C–Q slope and stream solute response are not impacted. Instead, primary mineral dissolution and subsequent clay precipitation drive strong chemostasis for silica and aluminium and enable cation exchange that buffers calcium and magnesium concentrations, leading to weak chemostatic behaviour for divalent cations. The influence of clays on stream C–Q highlights the importance of delineating geochemical controls along flowpaths, as upgradient mineral dissolution and clay precipitation enable downgradient cation exchange. Our results suggest that geochemical reactions will not be impacted by future decreasing flows, and thus where chemostasis currently exists, it will continue to persist despite changes in recharge.     

Evaluation of redox-sensitive metals in marine surface sediments influenced by the oxygen minimum zone of the Humboldt Current System,Northern Chile

Upwelling coastal systems can be used to understand how dissolved oxygen and biological productivity control the accumulation of redox-sensitive metals in marine sediments. The aluminium (Al), cadmium (Cd), iron (Fe), nickel (Ni), molybdenum (Mo), vanadium (V), total organic carbon (TOC), total nitrogen (TN) and total sulfur (TS) contents in surficial sediment collected from different water depths (30, 70, and 120?m) in three northern Chilean bays influenced by coastal upwelling and oxygen minimum zones (OMZs) were measured. Principal component analysis (PCA), cluster analysis, and Spearman?s rank correlation were used to identify the mechanisms responsible for the redox-sensitive metal accumulation. The content of redox-sensitive metals and organic components in sediment increased with increasing water column depth, whereas lithogenic metals decreased. In the Mejillones del Sur and Caldera bays, the enrichment factors of the redox-sensitive metals showed enrichment for all metals with depth. The Cd and V enrichments are mainly the product of biogenic flow to the seabed, and the Mo and Ni enrichments are due to preservation under low subsurface oxygen conditions. Sulfate reduction is not an important mechanism in the accumulation of redox-sensitive metals in the sediment of the three bays. The PCA showed that the behaviors of the redox-sensitive metals and organic components reflect differences in the effects of the OMZ in sediment along the coast of northern Chile, with a more intense OMZ in Mejillones del Sur bay and weaker OMZs in Caldera and Inglesa bays. However, the high degree of enrichment in redox-sensitive metals in Caldera Bay can be attributed to the intense activity of the mining industry near the bay, a situation that produces geochemical behavior similar to that observed in Mejillones del Sur Bay.     

Ancient subsurface structure beneath crater Clavius:constraint by recent high-precision gravity and topography data

With the increasing precision of the GRAIL gravity field models and topography from LOLA, it is possible to investigate the substructure beneath crater Clavius. An admittance between gravity and topography data is commonly used to estimate selenophysical parameters, including load ratio, crustal thickness and density, and elastic thickness. Not only a surface load, but also a subsurface load is considered in estimation. The algorithm of particle swarm optimization(PSO) with a swarm size of 400 is employed as well.Results indicate that the observed admittance is best-fitted by the modeled admittance based on a spherical shell model, which was proved to be unsatisfactory in the previous study. The best-fitted load ratio f is around-0.194. Such a small load ratio conforms to the direct proportion between the nearly uncompensated topography and its corresponding negative gravity anomaly. It also indicates that a surface load dominates all the loads. Constrained within 2σSTD, a small crustal thickness(～30 km) and a crustal density of ～2587 kg m-3are found, quite close to the results from previous GRAIL research. Considering the well constrained crustal thickness and density, the best-fitted elastic thickness(～7 km) is rational. This result is slightly smaller than the previous study(～12 km). Such difference can be attributed to the difference in crustal density used and the precision of gravity and topography data. Considering that the small difference between the modeled gravity anomaly and observations is quite small, a parameter inversed here could be an indicator of the subsurface structure beneath Clavius.     

Mesoscale eddies off Peru in altimeter records: Identification algorithms and eddy spatio-temporal patterns

Relatively little is known about coherent vortices in the eastern South-Pacific along the Peruvian coast, even with regard to basic facts about their frequency of occurrence, longevity and structure. This study addresses these issues with nearly 15 years of relatively high-resolution satellite altimetry measurements.We first compare two distinct automated methods for eddy identification. The objective validation protocol shows that the rarely-used geometrical or “winding-angle method”, based on the curvature of the streamline functions, is more accurate than the commonly-used “Okubo–Weiss algorithm”, which defines a vortex as a simple connected region with values of Okubo–Weiss parameter weaker than a given threshold.We then investigate vortices off Peru using more than 20,000 mesoscale eddies identified by the winding-angle method. Coherent eddies, characterized by a high ratio of vorticity to deformation rate, are typically formed along the coast and propagate westward at 3–6 cm s⁻¹. The vortices have a mean radius of 80 km, increasing northward, and are most frequently observed off of Chimbote (9°S) and south of San Juan (15°S). The mean eddy lifetime is about 1 month, but if eddies survive at least 2 months, the probability for surviving an additional week (or month) is constant at 90% (or 67%). Anticyclonic eddies tend to propagate northwestward whereas cyclonic vortices migrate southwestward. In general, cyclones and anticyclones are similar, except for eddies surviving at least 6 months. In this case, after a similar 3–4 months of radius and amplitude growth, amplitudes (or sizes) decay particularly rapidly for anticyclonic (or cyclonic) eddies. In terms of intensity, cyclonic eddies show a rapid decay during the first 3 months before arriving at a quasi-constant value, whereas anticyclones exhibit steady decline. Finally, eddy temporal variations were examined at seasonal and interannual scales in the “coastal” region favorable to the formation of energetic mesoscale structures. On seasonal scales, eddy activity is maximal in fall and minimum in spring. At interannual scales, the eddy activity index was maximal during the strong El Niño of 1997–1998 but another strong maximum of eddy activity also occurred late in 2004. These temporal variations are probably associated with the intensification of the upwelling thermal front and with the passage of coastal-trapped waves which generate baroclinic instabilities. Further investigation of the mechanisms involved on the eddy genesis is needed.     

Paleolakes, paleofloods, and depressions in Aurorae and Ophir Plana, Mars: Connectivity of surface and subsurface hydrological systems

The plains of Aurorae and Ophir in the equatorial region of Mars display geomorphic evidence indicative of extensive but generally short-lived paleohydrological processes. Elaver Vallis in Aurorae Planum south of Ganges Chasma is an outflow channel system >180 km long, and here inferred to have formed by cataclysmic spillover flooding from a paleolake(s) contained in the Morella crater basin. Ganges Cavus is an enormous 5-km-deep depression of probable collapse origin located in the Morella basin. The fluid responsible for the infilling of the Morella basin likely emerged at least partially through Ganges Cavus or its incipient depression, and it may have been supplied also from small-scale springs in the basin. Similar paleohydrological processes are inferred also in Ophir Planum. It is reasonable to assume that water, sometimes sediment-laden and/or mixed with gases, was the responsible fluid for these phenomena although some of the observed features could be explained by non-aqueous processes such as volcanism. Water emergence may have occurred as consequences of ground ice melting or breaching of cryosphere to release water from the underlying hydrosphere. Dike intrusion is considered to be an important cause of formation for the cavi and smaller depressions in Aurorae and Ophir Plana, explaining also melting of ground ice or breaching of cryosphere. Alternatively, the depressions and crater basins may have been filled by regional groundwater table rising during the period(s) when cryosphere was absent or considerably thin. The large quantities of water necessary for explaining the paleohydrological processes in Aurorae and Ophir Plana could have been derived through crustal migration from the crust of higher plains in western Ophir Planum where water existed in confined aquifers or was produced by melting of ground ice due to magmatic heating or climatic shift, or from a paleolake in Candor Chasma further west.     

Early Cretaceous (late Berriasian to early Aptian) palaeoceanographic change along the northwestern Tethyan margin (Vocontian Trough,southeastern France): δ13C, δ18O and Sr-isotope belemnite and whole-rock records

Stable carbon, oxygen, and strontium isotope records were obtained from uppermost Hauterivian to lowermost Aptian belemnite rostra, which were collected in well-dated sections from the Vocontian Trough (southeastern France). This data set complements previously published belemnite-isotope records from the uppermost Berriasian-Hauterivian interval from the same basin. The belemnite carbon and oxygen isotope record is compared to the carbonate bulk-rock isotope record from the same sections, and from additional Italian sections. With regards to their long-term trends, both belemnite and whole-rock δ¹⁸O records are well correlated, except for the uppermost Hauterivian-lower Barremian interval, within which they deviate. This discrepancy is interpreted to be linked to the latest Hauterivian Faraoni oceanic anoxic event and its early Barremian aftermath. The Faraoni level is characterized by enhanced sea-water stratification, probably induced by the onset of a warmer and more humid climate along the northern Tethyan margin. The early Barremian was characterized by stronger vertical sea-water mixing reflected by a decrease in density contrast between sea-surface and deeper waters. The belemnite oxygen isotope record shows a more stable evolution with smaller fluctuations than its bulk-rock counterpart, which indicates that deeper water masses were not as much subjected to density fluctuations as sea-surface water. The comparison of belemnite and bulk-rock carbon isotope records allows observing the impact of regional influence exerted by platform carbonate ooze shedding on the carbon cycle. Discrepancies in the two records are observed during time of photozoan carbonate platform growth. The strontium isotopic record shows a gradual increase from the uppermost Berriasian to the uppermost lower Barremian followed by a rapid decrease until the uppermost Barremian and a renewed small increase within the lowermost Aptian. The major inflection point in the uppermost lower Barremian appears to predate the onset in the formation of the Ontong-Java volcanic plateau.     

A Comparison of Fuzzy AHP and Ideal Point Methods for Evaluating Land Suitability

This article compares two fuzzy approaches to land suitability evaluations, Analytical Hierarchy Process (AHP) and Ideal Point. The methods were evaluated using a case study which models the opportunities for wheat production under irrigation conditions in the north‐western region of Jeffara Plain, Libya. A number of relevant soil and landscape criteria were identified through a review of the literature and their weights specified as a result of discussions with local experts. The results of the Fuzzy AHP showed that the majority of the study area has membership values to the set of suitability between 0.40 and 0.50, while the results of the Ideal Point approach revealed most of the study area to have membership values between 0.30 and 0.40. While the Fuzzy AHP and Ideal Point approaches accommodate the continuous nature of many soil properties and produce more intuitive distributions of land suitabilities values, the Fuzzy AHP approach was found to be better than Fuzzy Ideal Point. This was due to the latter's tendency to be biased towards positive and negative ideal values.