High seas marine protected areas: Benthic environmental conservation priorities from a GIS analysis of global ocean biophysical data

Designing a representative network of high seas marine protected areas (MPAs) requires an acceptable scheme to classify the benthic (as well as the pelagic) bioregions of the oceans. Given the lack of sufficient biological information to accomplish this task, we used a multivariate statistical method with 6 biophysical variables (depth, seabed slope, sediment thickness, primary production, bottom water dissolved oxygen and bottom temperature) to objectively classify the ocean floor into 53,713 separate polygons comprising 11 different categories, that we have termed “seascapes”. A cross-check of the seascape classification was carried out by comparing the seascapes with existing maps of seafloor geomorphology and seabed sediment type and by GIS analysis of the number of separate polygons, polygon area and perimeter/area ratio. We conclude that seascapes, derived using a multivariate statistical approach, are biophysically meaningful subdivisions of the ocean floor and can be expected to contain different biological associations, in as much as different geomorphological units do the same. Less than 20% of some seascapes occur in the high seas while other seascapes are largely confined to the high seas, indicating specific types of environment whose protection and conservation will require international cooperation. Our study illustrates how the identification of potential sites for high seas marine protected areas can be accomplished by a simple GIS analysis of seafloor geomorphic and seascape classification maps. Using this approach, maps of seascape and geomorphic heterogeneity were generated in which heterogeneity hotspots identify themselves as MPA candidates. The use of computer-aided mapping tools removes subjectivity in the MPA design process and provides greater confidence to stakeholders that an unbiased result has been achieved.     

Evolution of rifted ocean ridges

A major question in seafloor tectonics has been, how does the 2-km-deep rift valley characteristic of slow-spreading ridges evolve into the relatively horizontal undulating relief of the rift mountains? Deep-tow studies of the Mid-Atlantic Ridge suggest that the primary mechanism for transformation of the rift valley topography is normal faulting along fault planes which dip away from the valley axis. The faulting occurs in a narrow zone just beyond the outer walls of the rift valley. This model allows for a steady-state evolution of the rift valley into the rift mountains in which the state of stress in the oceanic lithosphere continues to be in horizontal deviatoric tension throughout the entire process. Alternate mechanisms involving reverse faulting or regional tilt may be active but are found to be of less importance. Implications for various dynamic models of the rift valley are discussed.     

Interarc spreading and Cordilleran tectonics as alternates related to the age of subducted oceanic lithosphere

Old, cold oceanic lithosphere is denser and therefore gravitationally more unstable than younger, hotter oceanic lithosphere. Hence, whereas old lithosphere will sink under its own weight, subduction of young lithosphere may require an additional force. Interarc spreading occurs or occurred recently in the western Pacific, in the southern Atlantic, and possibly in the Mediterranean, where the subducted sea floor appears to be more than 50 m.y. old, and in many cases, is more than 100 m.y. old. In most of these regions, the ease with which the old dense lithosphere sinks may have contributed to a seaward migration of the trenches, which led to interarc spreading. Cordilleran tectonics, including high mountains and broad zones of deformation, are present on the margins of the eastern Pacific where the subducted oceanic lithosphere is younger than about 50 m.y. An extra force, which we presume to be necessary to cause subduction of the young lithosphere, may be responsible for the deformation and mountains just as an extra force seems necessary to drive continental collision in Asia. The extensive early Tertiary deformation across a broad zone of western North America may be related to the long-term, continuous subduction of young lithosphere of the Farallon and Kula plates.     

The mineralogic transformation of ferrihydrite induced by heterogeneous reaction with bioreduced anthraquinone disulfonate (AQDS) and the role of phosphate

Bioreduced anthraquinone-2,6-disulfonate (AH₂DS; dihydro-anthraquinone) was reacted with a 2-line, Si-substituted ferrihydrite under anoxic conditions at neutral pH in PIPES buffer. Phosphate (P) and bicarbonate (C); common adsorptive oxyanions and media/buffer components known to effect ferrihydrite mineralization; and Fe(II)_aq (as a catalytic mineralization agent) were used in comparative experiments. Heterogeneous AH₂DS oxidation coupled with Fe(III) reduction occurred within 0.13-1 day, with mineralogic transformation occurring thereafter. The product suite included lepidocrocite, goethite, and/or magnetite, with proportions varing with reductant:oxidant ratio (r:o) and the presence of P or C. Lepidocrocite was the primary product at low r:o in the absence of P or C, with evidence for multiple formation pathways. Phosphate inhibited reductive recrystallization, while C promoted goethite formation. Stoichiometric magnetite was the sole product at higher r:o in the absence and presence of P. Lepidocrocite was the primary mineralization product in the Fe(II)_aq system, with magnetite observed at near equal amounts when Fe(II) was high [Fe(II)/Fe(III)] = 0.5 and P was absent. P had a greater effect on reductive mineralization in the Fe(II)_aq system, while AQDS was more effective than Fe(II)_aq in promoting magnetite formation. The mineral products of the direct AH₂DS-driven reductive reaction are different from those observed in AH₂DS-ferrihydite systems with metal reducing bacteria, particularly in presence of P.     

Source and fractionation controls on subduction-related plutons and dike swarms in southern Patagonia (Torres del Paine area) and the low Nb/Ta of upper crustal igneous rocks

The subduction system in southern Patagonia provides direct evidence for the variability of the position of an active continental arc with respect to the subducting plate through time, but the consequences on the arc magmatic record are less well studied. Here we present a geochemical and geochronological study on small plutons and dykes from the upper crust of the southern Patagonian Andes at ~?51°S, which formed as a result of the subduction of the Nazca and Antarctic plates beneath the South American continent. In situ U–Pb geochronology on zircons and bulk rock geochemical data of plutonic and dyke rocks are used to constrain the magmatic evolution of the retro-arc over the last 30 Ma. We demonstrate that these combined U–Pb and geochemical data for magmatic rocks track the temporal and spatial migration of the active arc, and associated retro-arc magmatism. Our dataset indicates that the rear-arc area is characterized by small volumes of alkaline basaltic magmas at 29–30 Ma that are characterized by low La/Nb and Th/Nb ratios with negligible arc signatures. Subsequent progressive eastward migration of the active arc culminated with the emplacement of calc-alkaline plutons and dikes?~?17–16 Ma with elevated La/Nb and Th/Nb ratios and typical subduction signatures constraining the easternmost position of the southern Patagonian batholith at that time. Geochemical data on the post-16 Ma igneous rocks including the Torres del Paine laccolith indicate an evolution to transitional K-rich calc-alkaline magmatism at 12.5?±?0.2 Ma. We show that trace element ratios such as Nb/Ta and Dy/Yb systematically decrease with increasing SiO₂, for both the 17–16 Ma calc-alkaline and the 12–13 Ma K-rich transitional magmatism. In contrast, Th/Nb and La/Nb monitor the changes in the source composition of these magmas. We suggest that the transition from the common calc-alkaline to K-rich transitional magmatism involves a change in the source component, while the trace element ratios, such as Nb/Ta and Dy/Yb, of derivative higher silica content liquids are controlled by similar fractionating mineral assemblages. Analysis of a global compilation of Nb/Ta ratios of arc magmatic rocks and simple geochemical models indicate that amphibole and variable amounts of biotite exert a major control on the low Dy/Yb and Nb/Ta of derivative granitic liquids. Lastly, we suggest that the low Nb/Ta ratio of silica-rich magmas is a natural consequence of biotite fractionation and that alternative models such as amphibolite melting in subduction zones and diffusive fractionation are not required to explain the Nb/Ta ratio of the upper continental crust.     

Movement of water infiltrated from a recharge basin to wells

Local surface water and stormflow were infiltrated intermittently from a 40-ha basin between September 2003 and September 2007 to determine the feasibility of recharging alluvial aquifers pumped for public supply, near Stockton, California. Infiltration of water produced a pressure response that propagated through unconsolidated alluvial-fan deposits to 125 m below land surface (bls) in 5 d and through deeper, more consolidated alluvial deposits to 194 m bls in 25 d, resulting in increased water levels in nearby monitoring wells. The top of the saturated zone near the basin fluctuates seasonally from depths of about 15 to 20 m. Since the start of recharge, water infiltrated from the basin has reached depths as great as 165 m bls. On the basis of sulfur hexafluoride tracer test data, basin water moved downward through the saturated alluvial deposits until reaching more permeable zones about 110 m bls. Once reaching these permeable zones, water moved rapidly to nearby pumping wells at rates as high as 13 m/d. Flow to wells through highly permeable material was confirmed on the basis of flowmeter logging, and simulated numerically using a two-dimensional radial groundwater flow model. Arsenic concentrations increased slightly as a result of recharge from 2 to 6 μg/L immediately below the basin. Although few water-quality issues were identified during sample collection, high groundwater velocities and short travel times to nearby wells may have implications for groundwater management at this and at other sites in heterogeneous alluvial aquifers.     

Chemical and structural characterization of As immobilization by nanoparticles of mackinawite (FeSm)

The mobility and availability of arsenite, As(III), in anoxic environments is largely controlled by adsorption onto iron sulfides and/or precipitation of arsenic in solid phases. The interaction of As(III) with synthetic mackinawite (FeS_m) in pH 5 and 9 suspensions was investigated using high-angle annular dark field (HAADF) scanning transmission electron microscopy (STEM), STEM elemental mapping, high resolution TEM, and X-ray photoelectron spectroscopy (XPS). At pH 5, arsenic sulfide phases precipitate among the FeS_m particles as discrete particles that are an amorphous hydrous phase of arsenic sulfide. The oxidation state of As in the surface layers of the arsenic sulfide precipitates is ‘realgar-like’ based on XPS results showing that > 75% of the As 3d peak area is due to As with oxidation states between 0 and 2+. Discrete, arsenic sulfide precipitates are absent at pH 9, but elemental mapping in STEM-EDX mode shows that arsenic is uniformly distributed on the FeS_m, suggesting that uptake is caused by the sorption of As(III) oxyanions and/or the precipitation of highly dispersed arsenic sulfides on FeS_m. XPS also revealed that the FeS_m that equilibrated without As(III) has a more oxidized surface composition than the sample at pH 9, as indicated by the higher concentration of O ( three times greater than that at pH 9) and the larger fraction of Fe(III) species making up the total Fe (2p_3/2) peak. These findings provide a better understanding of redox processes and phase transitions upon As(III) adsorption on iron sulfide substrates.     

Erosion of lithospheric mantle beneath the East African Rift system: geochemical evidence from the Kivu volcanic province

This study presents new major and trace element and Sr–Nd isotopic results for a suite of Miocene–Recent mafic lavas from the Kivu volcanic province in the western branch of the East African Rift. These lavas exhibit a very wide range in chemical and isotopic characteristics, due to a lithospheric mantle source region that is heterogeneous on a small scale, probably <1 km. The chemical and isotopic variations are mostly geographically controlled: lavas from Tshibinda volcano, which lies on a rift border fault on the northwestern margin of the province, have higher values of ⁸⁷Sr/⁸⁶Sr, (La/Sm)_n, Ba/Nb, and Zr/Hf than the majority of Kivu (Bukavu) samples. The range of ⁸⁷Sr/⁸⁶Sr at Tshibinda (0.70511–0.70514) overlaps some compositions found in the neighboring Virunga province, while Bukavu group lavas include the lowest ⁸⁷Sr/⁸⁶Sr (0.70314) and highest _Nd (+7.6) yet measured in western rift lavas. The Tshibinda compositions trend towards a convergence for Sr–Nd–Pb isotopic values among western rift lavas. Among Kivu lavas, variations in ¹⁴³Nd/¹⁴⁴Nd correlate with those for certain incompatible trace element ratios (e.g., Th/Nb, Zr/Hf, La/Nb, Ba/Rb), with Tshibinda samples defining one compositional extreme. There are covariations of isotopic and trace element ratios in mafic lavas of the East African Rift system that vary systematically with geographic location. The lavas represent a magmatic sampling of variations in the underlying continental lithospheric mantle, and it appears that a common lithospheric mantle (CLM) source is present beneath much of the East African Rift system. This source contains minor amphibole and phlogopite, probably due to widespread metasomatic events between 500 and 1000 Ma. Lava suites which do not show a strong component of the CLM source, and for which the chemical constraints also suggest the shallowest magma formation depths, are the Bukavu group lavas from Kivu and basanites from Huri Hills, Kenya. The inferred extent of lithospheric erosion therefore appears to be significant only beneath these two areas, which is generally consistent with lithospheric thickness variations estimated from gravity and seismic studies.     

Geochemistry of post-collisional mafic lavas from the North Anatolian Fault zone, Northwestern Turkey

Extensive magmatic activity developed at the northwestern part of the Anatolian block and produced basaltic lavas that are situated along and between the two segments of the North Anatolian Fault zone. This region is a composite tectonic unit formed by collision of continental fragments after consumption of Neotethyan ocean floor during the late Cretaceous. Northwestern Anatolian basalts and evolved lavas exhibit both tholeiitic and calc-alkaline characteristics. Mafic lavas are moderately enriched in LILE (except depleted part of Yuvacık and İznik samples) and depleted in HFSE (but not Zr, Hf) relative to primitive mantle values, suggesting derivation from a MORB-like mantle source that is unexpected in this subduction environment. Sr and Nd isotopes are close to the mantle array and vary beyond analytical error (⁸⁷Sr/⁸⁶Sr 0.70404–0.70546, ¹⁴³Nd/¹⁴⁴Nd 0.51270–0.51289). These geochemical features may result from two possible processes: (1) melting of a MORB-like mantle source that was modified by subduction-released fluids and melts or (2) modification of mafic liquids derived from a dominantly MORB-like source by crustal or lithospheric mantle material. Geochemical characteristics of the lavas (e.g., Ba/Rb, Rb/Sr, Ba/Zr, ⁸⁷Sr/⁸⁶Sr, Sr/P) vary systematically along the fault zone from east to west, consistent with a decrease in the degree of melting from east to west or a change in the nature of the source composition itself. Thus, the difference in incompatible elements and Sr–Nd isotopic ratios seems to result from small-scale mantle heterogeneity in a post-collisional tectonic environment.