Formation waters discharge to rivers near oil sands projects

Groundwater discharges in the western Canadian oil sands region impact river water quality. Mapping groundwater discharges to rivers in the oil sands region is important to target water quality monitoring efforts and to ensure injected wastewater and steam remain sequestered rather than eventually resurfacing. Saline springs composed of Pleistocene‐aged glacial meltwater exist in the region, but their spatial distribution has not been mapped comprehensively. Here we show that formation waters discharge into 3 major rivers as they flow through the Athabasca Oil Sands Region adjacent to many active oil sands projects. These discharges increase river chloride concentrations from river headwaters to downstream reaches by factors of ~23 in the Christina River, ~4 in the Clearwater River, and ~5 in the Athabasca River. Our survey provides further evidence for the substantial impact of formation water discharges on river water quality, even though they comprise less than ~2% of total streamflow. Geochemical evidence supporting formation water discharges as the leading control on river salinity include increases in river chloride concentrations, Na/(Na + Ca) ratios, Cl/(Cl + SO₄) ratios and decreases in ⁸⁷Sr/⁸⁶Sr ratios; each mixing trend is consistent with saline groundwater discharges sourced from Cretaceous or Devonian aquifers. These regional subsurface‐to‐surface connections signify that injected wastewater or steam may potentially resurface in the future, emphasizing the critical importance of mapping groundwater flow paths to understand present‐day streamflow quality and to predict the potential for injected fluids to resurface.     

Post-Marquette discharge from Glacial Lake Agassiz into the Superior basin

Glacial Lake Agassiz, the largest of the North American glacial lakes, discharged through several different outlets during its history, although the timing and location of discharge remain controversial. However, one discharge event is well established based on extensive onshore observations: drainage through the Nipigon Lake area into the Superior basin about 10,700 years ago, following retreat of ice of the Marquette advance from the basin. High-resolution, single-channel seismic-reflection data collected with a small airgun were acquired to test the hypothesis that the Post-Marquette drainage event left diagnostic stratigraphic and geomorphic signatures beneath Lake Superior. The unique bathymetry of northwestern Lake Superior, where water depth plunges off Nipigon and Black Bays, makes this location ideal for the characterization of the post-Marquette depositional and erosional features. According to our hypothesis, the initial, sudden discharge of high-velocity water would have eroded channels through the bays. The steep and sudden drop-off into the Superior basin would have caused the flow to slow and drop much of the sediment it was carrying. Our results confirm the existence of both erosional features and depositional sediment packages related to Lake Agassiz discharge at this time. The erosional features include deeply incised bedrock channels in the bays. The depositional features comprise subaqueous fans that are thickest in the deep water areas adjacent to the outlets and thin lakeward and laterally away from the channels. The seismic character of the basal units of the fans, proximal to the channels, is chaotic and only very weakly stratified, suggesting that these deposits represent coarse sediment laid down during the initial, high-energy stages of the flood. These sediments are overlain by a stratified package which is interpreted as the fine grained sediment associated with the later, low-energy stages of the flood. The combination of channels and subaqueous fans is inferred to be diagnostic of high-energy Lake Agassiz discharge into the Superior basin, and they serve as analogs for hypothesized discharge at other times.     

Geography in interdisciplinarity: Towards a third conversation

The ‘era of interdisciplinarity’ heralds collaborative inquiry as effective for addressing complex issues at the nexus of disciplinary interests. Geographers have long argued that they are particularly well-suited to contribute to interdisciplinary endeavors because of the breadth and depth that the discipline enfolds. However, within the literature about geography and interdisciplinarity, we find only two rather limiting conversations. The first conversation is concerned with the role(s) and position of geography within academia and focuses on what geographers can do to distinguish themselves while also improving their interactions with scholars from other disciplines. The second conversation largely revolves around how best to conduct interdisciplinarity ‘in the lab and classroom’ and focuses on practical issues associated with making collaborative research operate smoothly for multiple, disciplinary participants. We propose opening up intellectual space for a third conversation about the benefits, challenges and contributions of individual disciplines in interdisciplinary environments. Using survey data, we analyze how geography was perceived by collaborators from various disciplines in an interdisciplinary Urban Ecology program at the University of Washington. We offer this pilot study as a heuristic for others wishing to perform similar small-scale reflexive exercises and advance this “third conversation”.     

The internal structure of Jupiter family cometary nuclei from Deep Impact observations: The “talps” or “layered pile” model

We consider the hypothesis that the layering observed on the surface of Comet 9P/Tempel 1 from the Deep Impact spacecraft and identified on other comet nuclei imaged by spacecraft (i.e., 19P/Borrelly and 81P/Wild 2) is ubiquitous on Jupiter family cometary nuclei and is an essential element of their internal structure. The observational characteristics of the layers on 9P/Tempel 1 are detailed and considered in the context of current theories of the accumulation and dynamical evolution of cometary nuclei. The works of Donn [Donn, B.D., 1990. Astron. Astrophys. 235, 441-446], Sirono and Greenberg [Sirono, S.-I., Greenberg, J.M., 2000. Icarus 145, 230-238] and the experiments of Wurm et al. [Wurm, G., Paraskov, G., Krauss, O., 2005. Icarus 178, 253-263] on the collision physics of porous aggregate bodies are used as basis for a conceptual model of the formation of layers. Our hypothesis is found to have implications for the place of origin of the JFCs and their subsequent dynamical history. Models of fragmentation and rubble pile building in the Kuiper belt in a period of collisional activity (e.g., [Kenyon, S.J., Luu, J.X., 1998. Astron. J. 115, 2136-2160; 1999a. Astron. J. 118, 1101-1119; 1999b. Astrophys. J. 526, 465-470; Farinella, P., Davis, D.R., Stern, S.A., 2000. In: Mannings, V., Boss, A.P., Russell, S.S. (Eds.), Protostars and Planets IV. Univ. of Arizona Press, Tucson, pp. 1255-1282; Durda, D.D., Stern, S.J., 2000. Icarus 145, 220-229]) following the formation of Neptune appear to be in conflict with the observed properties of the layers and irreconcilable with the hypothesis. Long-term residence in the scattered disk [Duncan, M.J., Levison, H.F., 1997. Science 276, 1670-1672; Duncan, M., Levison, H., Dones, L., 2004. In: Festou, M., Keller, H.U., Weaver, H.A. (Eds.), Comets II. Univ. of Arizona Press, Tucson, pp. 193-204] and/or a change in fragmentation outcome modeling may explain the long-term persistence of primordial layers. In any event, the existence of layers places constraints on the environment seen by the population of objects from which the Jupiter family comets originated. If correct, our hypothesis implies that the nuclei of Jupiter family comets are primordial remnants of the early agglomeration phase and that the physical structure of their interiors, except for the possible effects of compositional phase changes, is largely as it was when they were formed. We propose a new model for the interiors of Jupiter family cometary nuclei, called the talps or “layered pile” model, in which the interior consists of a core overlain by a pile of randomly stacked layers. We discuss how several cometary characteristics—layers, surface texture, indications of flow, compositional inhomogeneity, low bulk density low strength, propensity to split, etc., might be explained in terms of this model. Finally, we make some observational predictions and suggest goals for future space observations of these objects.     

Metallurgy at the Crossroads: New Analyses of Copper-based Objects at Tianshanbeilu, Eastern Xinjiang, China

Tianshanbeilu is the largest Bronze Age site in eastern Xinjiang, China. Stretching across the entire second millennium BC, it performed a prominent role in connecting the Hexi corridor, Central China and the steppe. A further insight into the metallurgical tradition and the metal supply network is of vital importance to improve our understanding of its multi-connected nature. This paper offers a new set of chemical and isotopic data on the copper-based objects at Tianshanbeilu, including alloying elements, trace elements (impurities) and lead isotopes. Combining the concentrations of arsenic and antimony reveals that arsenic was introduced to copper partially due to the use of specific minerals tethrahedrite-tennantite. Lead isotopes demonstrate that multiple sources of copper were employed at Tianshanbeilu and a majority of them are characterized by common lead, which appears rather different from those of the Central Plains and the Hexi corridor, but highly consistent with local ores. Surprisingly, one object at Tianshanbeilu contains the well-known highly radiogenic lead. This object undoubtedly marks the westernmost boundary of the distribution of the highly radiogenic lead. We also anticipate that more lead isotopic analyses in NW China will further contribute to the study of the highly radiogenic lead in Central China.     

Evidence from olivine/melt element partitioning that nonbridging oxygen in silicate melts are not equivalent

Partitioning of Ca, Mn, Mg, and Fe²⁺ between olivine and melt has been used to examine the influence of energetically nonequivalent nonbridging oxygen in silicate melts. Partitioning experiments were conducted at ambient pressure in air and 1400°C with melts in equilibrium with forsterite-rich olivine (Fo >95 mol%). The main compositional variables of the melts were NBO/T and Na/(Na+Ca). In all melts, the main structural units were of Q⁴, Q³, and Q² type with nonbridging oxygen, therefore, in the Q³ and Q² units.For melts with high Q³/Q²-abundance ratio (corresponding to NBO/T near 1), increasing Na/(Na+Ca) [and Na/(Na+Ca+Mn+Mg+Fe²⁺)] results in a systematic decrease of the partition coefficients, K_Ca^ol/melt, K_Mn^ol/melt, K_Mg^ol/melt, and K_Fe2+^ol/melt, because of ordering of the network-modifying Ca, Mn, Mg, and Fe²⁺ among nonbridging oxygen in Q³ and Q² structural units. This decrease is more pronounced the smaller the ionic radius of the cation. With decreasing Q³/Q² abundance ratio (less-polymerized melts) this effect becomes less pronounced.Activity-composition relations among network-modifying cations in silicate melts are, therefore, governed by availability of energetically nonequivalent nonbridging oxygen in individual Qⁿ-species in the melt. As a result, any composition change that enhances abundance of highly depolymerized Qⁿ-species will cause partition coefficients to decrease.     

Effect of monochloramine on isolated fulvic acid

Monochloramine is interesting both as a selective oxidant of fulvic acid and as a drinking water disinfectant. In this study, the controlled reaction of aquatic fulvic acid with monochloramine did not result in products detectable by ether extraction-gas chromatography-flame ionization detection techniques. Evidence of a reaction was shown by bleaching of the fulvic acid solution (decreased absorbance at 465 nm) and chlorine substitution. Chlorine-containing products were quantified by the total carbon adsorbable organic halide (TOX) parameter. Bleaching and organic chlorine formation were much less extensive than in the reaction of free chlorine with fulvic acid. Monochloramine was shown to produce an organic chlorine fraction more hydrophilic and with higher molecular size than that produced by chlorine and fulvic acid. Results suggest that monochloramine may be useful tool for the investigation of certain humic functionalities because it reacts rather selectively and to a small extent with aquatic fulvic acid.     

Multidimensional Investigation of Bedrock Heterogeneity/Unconformities at a DNAPL‐Impacted Site

Organic solvent (i.e., dense nonaqueous phase liquid, DNAPL) migration in the subsurface is known to be extremely sensitive to geologic heterogeneity. There is often a focus on heterogeneity that results from changing depositional conditions over short spatial scales. Similar or even more extreme spatial heterogeneity can result postdeposition due to erosional processes. This study applies a synergistic approach based on a combination of high‐resolution lithologic logs of continuous cores, borehole geophysical logs, surface electrical resistivity, and seismic refraction tomography models to assess spatial heterogeneity in a shallow bedrock sequence subject to multiple unconformities and contaminated with a mixture of organic chemicals. The persistence of DNAPL in the source zone and an associated dissolved‐phase plume led to variable impacts on formation resistivity across the study site. Seismic refraction in combination with electrical resistivity tomography improved interpretation of highly irregular erosional boundaries by delineating sharp lateral transitions in lithologic composition near the source zone and across the dissolved‐phase plume. Electrical resistivity was effective at differentiating between clean and mud‐rich sandstones and their unconformable contact with an underlying dolostone. Geophysical measurements revealed eroded dolostone mounds encased by a network of younger mud‐rich sandstones channelized by clean semi‐lithified sand, all of which was buried beneath variable glacial drift. Our synergistic multidimensional approach resulted in the development of a detailed three‐dimensional shallow bedrock geospatial model, which has led to an improved understanding of DNAPL migration and contaminant plume heterogeneity.     

Global distribution of coastal cliffs

Previous studies have estimated that coastal cliffs exist on about 80% of the global shoreline, but have not been validated on a global scale. This study uses two approaches to capture information on the worldwide existence and erosion of coastal cliffs: a detailed literature survey and imagery search, and a GIS-based global mapping analysis. The literature and imagery review show coastal cliffs exist in 93% of the combined recognized independent coastal states and non-independent coastal regions worldwide (total of 213 geographic units). Additionally, cliff retreat rates have been quantified in at least one location within 33% of independent coastal states and 15% of non-independent regions. The GIS-based mapping used the near-global Shuttle Radar Topography Mission 3 arc second digital elevation model and Arctic Coastal Dynamics Database to obtain near-global backshore coastal elevations at 1 km alongshore intervals comprising about 1,340,000 locations (81% of the world vector shoreline). Backshore coastal elevations were compared with the mapped distribution of European coastal cliffs to produce a model training set, and this relationship was extended globally to map the likelihood of coastal cliff locations. About 21% of the transects (17% of the world vector shoreline) were identified as mangroves and eliminated as potential cliff locations. The results were combined with estimates of cliff percentages for Greenland and Antarctica from the literature, extending the global coverage to estimate cliff occurrence across 89% of the world vector shoreline. The results suggest coastal cliffs likely exist on about 52% of the global shoreline. © 2018 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.