Analysis of trends in annual streamflow to the Arctic Ocean

There is increasing interest in the magnitude of the flow of freshwater to the Arctic Ocean due to its impacts on the biogeophysical and socio‐economic systems in the north and its influence on global climate. This study examines freshwater flow based on a dataset of 72 rivers that either directly or indirectly contribute flow to the Arctic Ocean or reflect the hydrologic regime of areas contributing flow to the Arctic Ocean. Annual streamflow for the 72 rivers is categorized as to the nature and location of the contribution to the Arctic Ocean, and composite series of annual flows are determined for each category for the period 1975 to 2015. A trend analysis is then conducted for the annual discharge series assembled for each category. The results reveal a general increase in freshwater flow to the Arctic Ocean with this increase being more prominent from the Eurasian rivers than from the North American rivers. A comparison with trends obtained from an earlier study ending in 2000 indicates similar trend response from the Eurasian rivers, but dramatic differences from some of the North American rivers. A total annual discharge increase of 8.7 km³/y/y is found, with an annual discharge increase of 5.8 km³/y/y observed for the rivers directly flowing to the Arctic Ocean. The influence of annual or seasonal climate oscillation indices on annual discharge series is also assessed. Several river categories are found to have significant correlations with the Arctic Oscillation, the North Atlantic Oscillation, or the Pacific Decadal Oscillation. However, no significant association with climate indices is found for the river categories leading to the largest freshwater contribution to the Arctic Ocean.     

Competing space and place identity with landscape change analysis using Web GIS through Singapore historical maps (1828–2015)

The investigation of historical maps is fundamental to understanding the dynamics of landscape changes. In this study, we propose Web GIS tools as a way to compare historical maps efficiently for knowledge production. A key impetus of this article is to contribute to the ongoing efforts to broaden the appeal of ‘mixed‐methods research’ by bridging the epistemological and methodological gaps between GIS and qualitative methodologies in knowledge production. This article proposes a new visualization method for historical landscape change analysis by comparing multiple maps simultaneously with the mash‐up of Web GIS. Based on the analysis of exterior facts represented on the maps, four developmental footprints were investigated, namely, surviving place identity, disappeared place identity replaced by new identity, waning place identity overlapped with new identity, and waning identity by disconnected spatial relationships. To this end, the study argues that Web GIS applications have more potential in spatial knowledge production than traditional desktop GIS.     

Regional and species specific sexual reproductive patterns of three zooxanthellate scleractinian corals across the Eastern Tropical Pacific

Sexual reproduction of zooxanthellate scleractinian corals in the Eastern Tropical Pacific (ETP) is influenced by the interactive effects of regional and local oceanographic conditions, as well as a variety of other environmental factors. Differences in spatial and temporal gamete development and reproductive patterns of three widespread scleractinian corals of this region—Pocillopora damicornis (branching colony morphology, characterized as hermaphrodite broadcaster), Pavona gigantea (massive colonies, characterized as gonochoric broadcaster/sequential co‐sexually hermaphroditic) and Porites panamensis (encrusting colonies, characterized as gonochoric brooder)—were evaluated at local and regional scales across the ETP. This area extended from the Gulf of California (24°N) to the southern coast of Ecuador (–1°S), including offshore islands, using existing data pooled from prior studies. Predictive models were employed based on environmental variables: sea surface temperature, daylight hours, diffuse attenuation co‐efficient at 490 nm and photosynthetically active radiation. Datasets were extracted from satellite images ( https://oceancolor.gsfc.nasa.gov ) and analysed using WAM_STATIST ver. 6.33® software to obtain monthly average data from each site. The spatial (region, sub‐region and site) and seasonal (wet, dry) variation in reproductive activity (% colonies with gametes) differed among the three species; significant interactions were season × sub‐region for P. damicornis, season × site (sub‐region) for P. gigantea, and season × site for P. panamensis, for which sub‐region was not considered as a factor. The predictive models also suggest that gamete production/maturation of P. damicornis and P. gigantea is influenced by local differences in sea surface temperature and daylight hours. Porites panamensis was not correlated with any environmental variable examined. Variations in local and regional reproductive developmental patterns are likely an acclimatization response by each species imposed by the timing and duration of local stressor events. This analysis has provided insights into the diverse local and regional physical drivers that affect species responses and acclimatization in sexual reproduction across the ETP.     

Ceres’ spectral link to carbonaceous chondrites—Analysis of the dark background materials

Ceres’ surface has commonly been linked with carbonaceous chondrites (CCs) by ground‐based telescopic observations, because of its low albedo, flat to red‐sloped spectra in the visible and near‐infrared (VIS/NIR) wavelength region, and the absence of distinct absorption bands, though no currently known meteorites provide complete spectral matches to Ceres. Spatially resolved data of the Dawn Framing Camera (FC) reveal a generally dark surface covered with bright spots exhibiting reflectance values several times higher than Ceres’ background. In this work, we investigated FC data from High Altitude Mapping Orbit (HAMO) and Ceres eXtended Juling (CXJ) orbit (~140 m/pixel) for global spectral variations. We found that the cerean surface mainly differs by spectral slope over the whole FC wavelength region (0.4–1.0 μm). Areas exhibiting slopes ?1 constitute only ~3% of the cerean surface and mainly occur in the bright material in and around young craters, whereas slopes ≥?10% μm^?1 occur on more than 90% of the cerean surface; the latter being denoted as Ceres’ background material in this work. FC and Visible and Infrared Spectrometer (VIR) spectra of this background material were compared to the suite of CCs spectrally investigated so far regarding their VIS/NIR region and 2.7 μm absorption, as well as their reflectance at 0.653 μm. This resulted in a good match to heated CI Ivuna (heated to 200–300 °C) and a better match for CM1 meteorites, especially Moapa Valley. This possibly indicates that the alteration of CM2 to CM1 took place on Ceres.     

China and the changing economic geography of coffee value chains

For the past three centuries, the economic geography of the global coffee sector has been characterized by the supply of beans from tropical countries for consumption in North America and Europe, with various modes of value chain coordination enacted by lead firms to ensure reliable and affordable supply. This pattern is now fundamentally changing, with growth in coffee consumption in emerging markets, including China, exceeding that in established markets. But China is not only a growing consumer market, it is less well known that rapidly increasing agricultural production in Yunnan province of southwest China has also inserted the country as an important source region for coffee, and this has been pivotal in facilitating the emergence of Chinese lead firms in the sector. This article presents the emergence of China, and Chinese firms, at a critical juncture for the structure and governance of the global value chain for coffee. The processes through which this is occurring are outlined, and the implications for regional development prospects across Southeast Asia are discussed. We argue that the changing economic geography of coffee value chains, and their increasing driven‐ness by Chinese actors, is starting to reshape the regional coffee industry in profoundly new ways.     

Phase equilibria modelling of retrograde amphibole and clinozoisite in mafic eclogite from the Tso Morari massif,northwest India: constraining the P–T–M(H2O) conditions of exhumation

Phase equilibria modelling of post‐peak metamorphic mineral assemblages in (ultra)high‐P mafic eclogite from the Tso Morari massif, Ladakh Himalaya, northwest India, has provided new insights into the potential behaviour and source of metamorphic fluid during exhumation, and constrained the P–T conditions of hydration. A series of P–M(H₂O) pseudosections constructed in the Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O (NCKFMASHTO) system show that a number of petrographically distinct hydration episodes occurred during exhumation from peak P–T conditions (~640 °C, 27–28 kbar), resulting in the formation of abundant compositionally zoned amphibole and minor clinozoisite poikiloblasts at the expense of a peak assemblage dominated by garnet and omphacite. Initial hydration is interpreted to have occurred as a result of the destabilization of talc following isothermal decompression to ~23 kbar, which led to the formation of barroisite–winchite amphibole core domains. An episode of fluid infiltration from an external source at ~19 kbar, with or without syn‐decompressional cooling to ~560 °C, resulted in further barroisitic–winchitic amphibole growth, followed by the formation of clinozoisite poikiloblasts. Continued buoyancy‐driven exhumation to the base of the lower crust is constrained to have taken place with no additional fluid input. A final hydration event is characterized by the formation of magnesiohornblende rims on the barroisite–winchite cores, with the former interpreted to have formed during later prograde overprinting in the middle crust associated with the final stages of exhumation. Notably, the vast majority of externally sourced H₂O, comprising just over half of the current bulk rock fluid content, was added during this later hydration event. In a middle crustal setting, this is interpreted as the result of devolatilization reactions occurring in migmatitic host orthogneiss and/or metasedimentary units, or following the crystallization of partial melt.     

Simulation of less‐mobile porosity dynamics in contrasting sediment water interface porous media

Considering heterogeneity in porous media pore size and connectivity is essential to predicting reactive solute transport across interfaces. However, exchange with less‐mobile porosity is rarely considered in surface water/groundwater recharge studies. Previous research indicates that a combination of pore‐fluid sampling and geoelectrical measurements can be used to quantify less‐mobile porosity exchange dynamics using the time‐varying relation between fluid and bulk electrical conductivity. For this study, we use macro‐scale (10 s of cm) advection–dispersion solute transport models linked with electrical conduction in COMSOL Multiphysics to explore less‐mobile porosity dynamics in two different types of observed sediment water interface porous media. Modeled sediment textures contrast from strongly layered streambed deposits to poorly sorted lakebed sands and cobbles. During simulated ionic tracer perturbations, a lag between fluid and bulk electrical conductivity, and the resultant hysteresis, is observed for all simulations indicating differential loading of pore spaces with tracer. Less‐mobile exchange parameters are determined graphically from these tracer time series data without the need for inverse numerical model simulation. In both sediment types, effective less‐mobile porosity exchange parameters are variable in response to changes in flow direction and fluid flux. These observed flow‐dependent effects directly impact local less‐mobile residence times and associated contact time for biogeochemical reaction. The simulations indicate that for the sediment textures explored here, less‐mobile porosity exchange is dominated by variable rates of advection through the domain, rather than diffusion of solute, for typical low‐to‐moderate rate (approximately 3–40 cm/day) hyporheic fluid fluxes. Overall, our model‐based results show that less‐mobile porosity may be expected in a range of natural hyporheic sediments and that changes in flowpath orientation and magnitude will impact less‐mobile exchange parameters. These temporal dynamics can be assessed with the geoelectrical experimental tracer method applied at laboratory and field scales.     

Delimiting the eastern extent of the Altyn Tagh Fault: Insights from structural analyses of seismic reflection profiles

The Altyn Tagh Fault (ATF) serves as a key continental‐scale controlling structural element of the Tibetan Plateau. However, its eastward extent remains controversial. Here we use high‐resolution seismic reflection profiles to investigate the subsurface structures of the easternmost ATF and use these to delimit the easternmost extent of the fault. The structural analyses show an eastward geometric change from transpressional positive flower structures to compressional thrusts, with transpression‐induced shortening magnitudes decreasing eastwards from a maximum of ~5.3 km to being absent. Stratigraphic controls indicate that the deformation took place over the last ~<1.2 Ma. Our wider findings lead us to: (a) reject the suggestion that the ATF previously extended beyond the Kuantan Shan‐Hei Shan to link with the Alxa‐East Mongolia Fault; and (b) propose that the rigid block model used to describe the Tibetan Plateau crust is not consistent with the extent and structural details of the easternmost ATF.     

A luminescence‐based chronology for the Harletz loess sequence,Bulgaria

The Harletz loess‐palaeosol sequence is located in northwestern Bulgaria and represents an important link between well‐studied loess sequences in eastern Romania and further sites to the west of the Carpathians (e.g. Serbia and Hungary). The aim of this study was to establish a chronostratigraphy of the deposits, using various methods of luminescence dating, together with basic stratigraphical field observations as well as magnetic properties. Luminescence dating was carried out using the quartz fine grain fraction and a SAR protocol, and the feldspar coarse grain fraction, applying the MET‐pIRIR protocol. Due to underestimation of the quartz fine grain fraction in the lower parts of the sequence, the resulting chronology is mainly based on the feldspar ages, which are derived from the stimulation temperature at 150 °C. A comparison with nearby sequences from Serbia, Hungary and Romania, and interpretations obtained through the stratigraphical and sedimentological signature of the sequence, supports the established chronology. Our data suggest that the prominent palaeosol (soil complex) in the upper quarter of the sequence was formed during MIS 5. It would follow that large parts of the Last Glacial loess overlying this palaeosol were probably eroded, and that the thick loess accumulation underlying this soil complex can be allocated to the penultimate glacial (MIS 6). A prominent MIS 6 tephra, which has been reported from other sequences in the area, is also present at Harletz.     

Subduction metamorphism of serpentinite‐hosted carbonates beyond antigorite‐serpentinite dehydration (Nevado‐Filábride Complex,Spain)

At sub‐arc depths, the release of carbon from subducting slab lithologies is mostly controlled by fluid released by devolatilization reactions such as dehydration of antigorite (Atg‐) serpentinite to prograde peridotite. Here we investigate carbonate–silicate rocks hosted in Atg‐serpentinite and prograde chlorite (Chl‐) harzburgite in the Milagrosa and Almirez ultramafic massifs of the palaeo‐subducted Nevado‐Filábride Complex (NFC, Betic Cordillera, S. Spain). These massifs provide a unique opportunity to study the stability of carbonate during subduction metamorphism at P–T conditions before and after the dehydration of Atg‐serpentinite in a warm subduction setting. In the Milagrosa massif, carbonate–silicate rocks occur as lenses of Ti‐clinohumite–diopside–calcite marbles, diopside–dolomite marbles and antigorite–diopside–dolomite rocks hosted in clinopyroxene‐bearing Atg‐serpentinite. In Almirez, carbonate–silicate rocks are hosted in Chl‐harzburgite and show a high‐grade assemblage composed of olivine, Ti‐clinohumite, diopside, chlorite, dolomite, calcite, Cr‐bearing magnetite, pentlandite and rare aragonite inclusions. These NFC carbonate–silicate rocks have variable CaO and CO₂ contents at nearly constant Mg/Si ratio and high Ni and Cr contents, indicating that their protoliths were variable mixtures of serpentine and Ca‐carbonate (i.e., ophicarbonates). Thermodynamic modelling shows that the carbonate–silicate rocks attained peak metamorphic conditions similar to those of their host serpentinite (Milagrosa massif; 550–600°C and 1.0–1.4 GPa) and Chl‐harzburgite (Almirez massif; 1.7–1.9 GPa and 680°C). Microstructures, mineral chemistry and phase relations indicate that the hybrid carbonate–silicate bulk rock compositions formed before prograde metamorphism, likely during seawater hydrothermal alteration, and subsequently underwent subduction metamorphism. In the CaO–MgO–SiO₂ ternary, these processes resulted in a compositional variability of NFC serpentinite‐hosted carbonate–silicate rocks along the serpentine‐calcite mixing trend, similar to that observed in serpentinite‐hosted carbonate‐rocks in other palaeo‐subducted metamorphic terranes. Thermodynamic modelling using classical models of binary H₂O–CO₂ fluids shows that the compositional variability along this binary determines the temperature of the main devolatilization reactions, the fluid composition and the mineral assemblages of reaction products during prograde subduction metamorphism. Thermodynamic modelling considering electrolytic fluids reveals that H₂O and molecular CO₂ are the main fluid species and charged carbon‐bearing species occur only in minor amounts in equilibrium with carbonate–silicate rocks in warm subduction settings. Consequently, accounting for electrolytic fluids at these conditions slightly increases the solubility of carbon in the fluids compared with predictions by classical binary H₂O–CO₂ fluids, but does not affect the topology of phase relations in serpentinite‐hosted carbonate‐rocks. Phase relations, mineral composition and assemblages of Milagrosa and Almirez (meta)‐serpentinite‐hosted carbonate–silicate rocks are consistent with local equilibrium between an infiltrating fluid and the bulk rock composition and indicate a limited role of infiltration‐driven decarbonation. Our study shows natural evidence for the preservation of carbonates in serpentinite‐hosted carbonate–silicate rocks beyond the Atg‐serpentinite breakdown at sub‐arc depths, demonstrating that carbon can be recycled into the deep mantle.