A paleoclimatic context for the drought of 1999–2001 in the northern Great Plains of North America

The least annual precipitation in the western interior of North America occurs in the northern Great Plains, including an area that encompasses parts of south-eastern Alberta, south-western Saskatchewan and eastern Montana. During 1999–2001, most climate stations in this region had record low precipitation. This paper examines this three-year drought in the context of historical climate records from Medicine Hat, Alberta and Havre, Montana and summer (June–July) and annual (August–July) precipitation reconstructed from standardized tree-ring widths (residual chronologies) from Pinus contorta (lodgepole pine) sampled in the Cypress Hills of Alberta and Saskatchewan and the Bears Paw Mountains of north-central Montana. Drought is operationally defined as precipitation in the lower 10th and 20th percentiles. Plots of reconstructed precipitation and cumulative departure from median values indicate a shift in climate variability prior to the twentieth century, when EuroCanadians settled in this region. The eighteenth and nineteenth centuries are characterized by sustained periods of progressively wetter and drier conditions, including prolonged drought. Various archival sources document the significant impacts of these prolonged droughts. While drought was frequent in the twentieth century, it tended to be of short duration and the impacts also were ameliorated by intervening periods of relatively high precipitation. Increasing aridity in response to global warming could expose a larger area of the northern Great Plains to the impacts of drought.     

Flow Cross K-function: a bivariate flow analytical method

Spatial flow data represent meaningful interaction activities between pairs of corresponding locations, such as daily commuting, animal migration, and merchandise shipping. Despite recent advances in flow data analytics, there is a lack of literature on detecting bivariate or multivariate spatial flow patterns. In this paper we introduce a new spatial statistical method called Flow Cross K-function, which combines the Cross K-function that detects marked point patterns and the Flow K-function that detects univariate flow clustering patterns. Flow Cross K-function specifically assesses spatial dependence of two types of flow events, in other words, whether one type of flows is spatially associated with the other, and if so, whether this is according to a clustering or dispersion trend. Both a global version and a local version of Flow Cross K-function are developed. The former measures the overall bivariate flow patterns in the study area, while the latter can identify anomalies at local scales that may not follow the global trend. We test our method with carefully designed synthetic data that simulate the extreme situations. We exemplify the usefulness of this method with an empirical study that examines the distributions of taxi trip flows in New York City.     

Influence of microclimate and geomorphological factors on alpine vegetation in the Western Swiss Alps

Among the numerous environmental factors affecting plant communities in alpine ecosystems, the influence of geomorphic processes and landforms has been minimally investigated. Subjected to persistent climate warming, it is vital to understand how these factors affect vegetation properties. Here, we studied 72 vegetation plots across three sites located in the Western Swiss Alps, characterized by high geomorphological variability and plant diversity. For each plot, vascular plant species were inventoried and ground surface temperature, soil moisture, topographic variables, earth surface processes (ESPs) and landform morphodynamics were assessed. The relationships between plant communities and environmental variables were analysed using non-metric multi-dimensional scaling (NMDS) and multivariate regression techniques (generalized linear model, GLM, and generalized additive model, GAM). Landform morphodynamics, growing degree days (sum of degree days above 5°C) and mean ground surface temperature were the most important explanatory variables of plant community composition. Furthermore, the regression models for species cover and species richness were significantly improved by adding a morphodynamics variable. This study provides complementary support that landform morphodynamics is a key factor, combined with growing degree days, to explain alpine plant distribution and community composition. © 2019 John Wiley & Sons, Ltd. © 2019 John Wiley & Sons, Ltd.     

A pattern-based approach for multiple removal applied to a 3D Gulf of Mexico data set

Surface‐related multiples are attenuated for one sail line and one streamer of a 3D data set (courtesy of Compagnie Générale de Géophysique). The survey was carried out in the Gulf of Mexico in the Green Canyon area where salt intrusions close to the water‐bottom are present. Because of the complexity of the subsurface, a wavefield method incorporating the full 3D volume of the data for multiple removal is necessary. This method comprises modelling of the multiples, where the data are used as a prediction operator, and a subtraction step, where the model of the multiples is adaptively removed from the data with matching filters. The accuracy of the multiple model depends on the source/receiver coverage at the surface. When this coverage is not dense enough, the multiple model contains errors that make successful subtraction more difficult. In these circumstances, one can either (1) improve the modelling step by interpolating the missing traces, (2) improve the subtraction step by designing methods that are less sensitive to modelling errors, or (3) both. For this data set, the second option is investigated by predicting the multiples in a 2D sense (as opposed to 3D) and performing the subtraction with a pattern‐based approach. Because some traces and shots are missing for the 2D prediction, the data are interpolated in the in‐line direction using a hyperbolic Radon transform with and without sparseness constraints. The interpolation with a sparseness constraint yields the best multiple model. For the subtraction, the pattern‐based technique is compared with a more standard, adaptive‐subtraction scheme. The pattern‐based approach is based on the estimation of 3D prediction‐error filters for the primaries and the multiples, followed by a least‐squares estimation of the primaries. Both methods are compared before and after prestack depth migration. These results suggest that, when the multiple model is not accurate, the pattern‐based method is more effective than adaptive subtraction at removing surface‐related multiples while preserving the primaries.     

An improved process-based representation of stream solute transport in the soil and water assessment tools

Hydrological models have long been used to study the interactions between land, surface and groundwater systems, and to predict and manage water quantity and quality. The soil and water assessment tool (SWAT), a widely used hydrological model, can simulate various ecohydrological processes on land and subsequently route the water quality constituents through surface and subsurface waters. So far, in-stream solute transport algorithms of the SWAT model have only been minimally revised, even though it has been acknowledged that an improvement of in-stream process representation can contribute to better model performance with respect to water quality. In this study, we aim to incorporate a new and improved solute transport model into the SWAT model framework. The new process-based model was developed using in-stream process equations from two well established models—the One-dimensional Transport with Inflow and Storage model and the Enhanced Stream Water Quality Model. The modified SWAT model (Mir-SWAT) was tested for water quality predictions in a study watershed in Germany. Compared to the standard SWAT model, Mir-SWAT improved dissolved oxygen (DO) predictions by removing extreme low values of DO (<6 mg/L) simulated by SWAT. Phosphate concentration peaks were reduced during high flows and a better match of daily predicted and measured values was attained using the Mir-SWAT model (R² = 0.17, NSE = −0.65, RSR = 1.29 with SWAT; R² = 0.28, NSE = −0.04, RSR = 1.02 with Mir-SWAT). In addition, Mir-SWAT performed better than the SWAT model in terms of Chlorophyll-a content particularly during winter months, improving the NSE and RSR for monthly average Chl-a by 74 and 42%, respectively. With the new model improvements, we aim to increase confidence in the stream solute transport component of the model, improve the understanding of nutrient dynamics in the stream, and to extend the applicability of SWAT for reach-scale analysis and management.     

Computation of the potential of a simple inhomogeneous equipotential layer with a given density on a given surface

Summary The potential of a simple equipotential layer, defined by the equation of a surface, on which the layer is spread, and by the density function in the form of harmonic expansions, is expressed with the aid of the parameters of the surface and of the density function. The integral of the product of three Legendre associated functions is treated with the aid of Gaunt's integral [3, 4].     

Sinusoidal oscillations radiating from a cylindrical source in thermal conduction or groundwater flow: Closed‐form solutions

The diffusion equation governs thermal conduction and groundwater flow phenomena. In this paper, we study the two‐dimensional radial propagation of a sinusoidal perturbation radiating from a cylindrical source within an infinite slab of homogeneous material. The solution of this problem has several applications. For instance, it can be used to determine the hydraulic diffusivity of the subsurface based on measurements of the hydraulic head around a vertical well during its development. For thermal problems, it can be used to determine the thermal diffusivity based on measurements of the temperature distribution around a cylindrical heat source generating a sinusoidal power per unit length. In this paper, we present a comprehensive analytical solution of this problem and we compare these solutions with numerical solutions. Two approximate analytical solutions, which can be relevant in practice, are also presented. Finally, we give an upper bound for the survival time of the transient part of the solution and we provide an estimate of the radius of influence of the sinusoidal solicitation. Copyright © 2010 John Wiley & Sons, Ltd.     

Analysis and Estimation of Multi-unit Deposits: Application to a Porphyry Copper Deposit

This work concerns mineral deposits consisting of geological bodies whose metal grades have different characteristics in terms of distribution and variogram, which means that estimating grades by ordinary kriging may produce unrealistic spatial continuity. This paper proposes a method based on the indicators of the geological objects (hereafter called units) and their product with the metal grade. This is illustrated by an application to a porphyry copper deposit. The aim of this paper is essentially to promote the use of variogram ratios to analyze and characterize deposits.     

A Leaky Aquifer below Champlain Sea Clay: Closed‐Form Solutions for Natural Seepage

Closed‐form solutions are proposed for natural seepage in semiconfined (leaky) aquifers such as those existing below the massive Champlain Sea clay layers in the Saint‐Lawrence River Valley. The solutions are for an ideal horizontal leaky aquifer below an ideal aquitard that may have either a constant thickness and a constant hydraulic head at its surface, or a variable thickness and a variable hydraulic head at its surface. A few simplifying assumptions were needed to obtain the closed‐form solutions. These have been verified using a finite element method, which did not make any of the assumptions but gave an excellent agreement for hydraulic heads and groundwater velocities. For example, the difference between the two solutions was smaller than 1 mm for variations in the 5 to 8 m range for the hydraulic head in the semiconfined aquifer. Note that fitting the hydraulic head data of monitoring wells to the theoretical solutions gives only the ratio of the aquifer and aquitard hydraulic conductivities, a clear case of multiple solutions for an inverse problem. Consequently, field permeability tests in the aquitard and the aquifer, and pumping tests in the aquifer, are still needed to determine the hydraulic conductivity values.     

Dynamics and fluxes of organic carbon and nitrogen in two Guiana Shield river basins impacted by deforestation and mining activities

Deforestation and mining activities have proven to be very damaging to rivers because these activities disturb the environmental characteristics of rivers. Thus, the concentrations of dissolved organic carbon (DOC), particulate organic carbon (POC), particulate nitrogen (PN), and Chlorophyll‐a (Chl‐a) were measured monthly during 2 hydrological years in the Maroni and Oyapock Rivers to assess the dynamics and fluxes of organic carbon and nitrogen in these 2 Guiana Shield basins, which have been strongly (Maroni) and weakly (Oyapock) impacted by deforestation and mining activities. The 2‐year time series show that DOC, POC, PN, and Chl‐a concentrations vary seasonally with discharge in both rivers, indicating a hydrologically dominated control. Temporal patterns of DOC, POC, and PN indicate that these variables show maximum concentrations in rising waters due to the yield of organic matter and nitrogen accumulated in soils, which are incorporated into the rivers during rainfall. However, the Chl‐a concentrations were at a maximum during low‐water stages. The C/N and C/Chl‐a ratios also showed a seasonal trend, with lower values during the low water periods due to an increase in algal biomass. During high water, the POC in both rivers is the result of terrestrial organic matter, whereas during low‐water autochthonous organic matter can reach up to 34% of the POC. The mean annual fluxes of TOC and PN were higher (4.56 × 10⁵ tonC year^?1 and 1.77 × 10⁴ tonN year^?1, respectively) in the Maroni River than those (1.84 × 10⁵ tonC year^?1 and 0.54 × 10⁴ tonN year^?1, respectively) in the Oyapock River. However, the specific fluxes of DOC, POC, and PN from both basins were nearly the same. Although gold mining activities are performed in both basins, there is no conclusive evidence regarding the impact of these activities on the dynamics of organic matter and particulate nitrogen in the Maroni and Oyapock Rivers.