Acid mine drainage at Cerro Rico de Potosí I: unabated high-strength discharges reflect a five century legacy of mining

Intensive mining and processing of Ag, Sn, Pb and Zn ores have occurred in various locations within and around the city of Potosí, Bolivia since 1545. Surface and subsurface waters, stream sediments and soils are contaminated with various ecotoxic metals in the headwaters of the economically vital, yet highly impacted, upper Rio Pilcomayo watershed. Previous studies have documented downstream trace metal contamination, however, not addressed specific sources. The AMD discharges identified in this study help link downstream pollution to primary origins. The majority of AMD would be considered high-strength due to metal concentrations and acidity orders of magnitude greater than typical AMD. Discharges from both operating and abandoned portals as well as tailings-related deposits displayed a high degree of heterogeneity with total metal concentrations ranging from 0.11 to 7480, <0.022 to 889, <0.0006 to 65.3, <0.001 to 310, 0.12 to 72,100, 0.3 to 402, <0.012 to 34.8, and 0.24 to 19,600 mg/L of Al, As, Cd, Cu, Fe, Mn, Pb and Zn, respectively. Net acidity and pH ranged from −10 to 246,000 mg/L as CaCO₃ equivalent and 0.90–6.94 standard units, respectively. Data were gathered during two sampling events centered around the most extreme periods of the dry and wet seasons of one water-year. Loadings to local streams were marginally greater for most metals in the wet season. If observed loadings are historically representative, Cerro Rico AMD has contributed thousands of tonnes of ecotoxic metals to the upper Rio Pilcomayo over the last five centuries. Metal and hydrogen ion concentrations in the majority of AMD sampled were several orders of magnitude above discharge limits set by the Bolivian government, yet no action has historically or contemporarily been taken.     

Chemical composition of Russian Arctic precipitation in 2007-2015

The results of observations of the chemical composition of precipitation in the Russian Arctic in 2007-2015 are summarized including the data from NP-35 drifting ice station obtained in 2007-2008 in the framework of the joint program of Arctic and Antarctic Research Institute and Voeikov Main Geophysical Observatory. The qualitative and quantitative difference in the chemical composition of precipitation in the Atlantic, Siberian, and Pacific sectors of the Russian Arctic is revealed. It is found that the concentration of microelements (heavy metals) in precipitation at NP-35 ice station did not exceed 12% of total ions at the lowest mineralization. The comparative analysis is presented of the concentration of sulfate ions in precipitation in Norilsk and at NP-35 ice stations.     

Non-Linear Theory and Power-Law Models for Information Integration and Mineral Resources Quantitative Assessments

Singular physical or chemical processes may result in anomalous amounts of energy release or mass accumulation that, generally, are confined to narrow intervals in space or time. Singularity is a property of different types of non-linear natural processes including cloud formation, rainfall, hurricanes, flooding, landslides, earthquakes, wildfires, and mineralization. The end products of these non-linear processes can be modeled as fractals or multifractals. Hydrothermal processes in the Earth’s crust can result in ore deposits characterized by high concentrations of metals with fractal or multifractal properties. Here we show that the non-linear properties of the end products of singular mineralization processes can be applied for prediction of undiscovered mineral deposits and for quantitative mineral resource assessment, whether for mineral exploration or for regional, national and global planning for mineral resource utilization. In addition to the general theory and framework for the non-linear mineral resources assessment, this paper focuses on several power-law models proposed for characterizing non-linear properties of mineralization and for geoinformation extraction and integration. The theories, methods, and computer system discussed in this paper were validated using a case study dealing with hydrothermal Au mineral potential in southern Nova Scotia, Canada.     

Multi-decadal synthesis of benthic–pelagic coupling in the western arctic: Role of cross-shelf advective processes

Using geographic information systems (GIS) software and geostatistical techniques, we utilized three decades of water-column chlorophyll a data to examine the relative importance of autochthonous versus allochthonous sources of reduced carbon to benthic communities that occur from the northern Bering to the eastern Beaufort Sea shelf. Spatial trend analyses revealed areas of high benthic biomass (>300 g m⁻²) and chlorophyll (>150 mg m⁻²) on both the southern and northern Chukchi shelf; both areas are known as depositional centers for reduced organic matter that originates on the Bering Sea shelf and is advected northward in Anadyr and Bering shelf water masses. We found a significant correlation between biomass and chlorophyll a in the Chukchi Sea, reflective of the strong benthic–pelagic coupling in a system that is utilized heavily by benthic-feeding marine mammals. In contrast, there was no significant correlation between biomass and chlorophyll in the Beaufort Sea, which by comparison, is considerably less productive (biomass and chlorophyll, <75 g m⁻² and <50 mg m⁻², respectively). One notable exception is an area of relatively high biomass (50–100 g m⁻²) and chlorophyll (80 mg m⁻²) near Barter Island in the eastern Beaufort Sea. Compared to other adjacent areas in the Beaufort Sea, the chlorophyll values in the vicinity of Barter Island were considerably higher and likely reflect a long-hypothesized upwelling in that area and close coupling between the benthos and autochthonous production. In the Bering Sea, a drop in benthic biomass in 1994 compared with previous measurements (1974–1993) may support earlier observations that document a decline in biomass that began between the 1980s and 1990s in the Chirikov Basin and south of St. Lawrence Island. The results of this study indicate that the benthos is an excellent long-term indicator of both local and physical advective processes. In addition, this work provides further evidence that secondary production on arctic shelves can be significantly augmented by reduced carbon advected from highly productive adjacent shelves.     

Geometry and kinematics of dunes during steady and unsteady flows in the Calamus River, Nebraska, USA

The geometry and kinematics of river dunes were studied in a reach of the Calamus River, Nebraska. During day-long surveys, dune height, length, steepness, migration rate, creation and destruction were measured concurrently with bedload transport rate, flow depth, flow velocity and bed shear stress. Within a survey, individual dune heights, lengths and migration rates were highly variable, associated with their three-dimensional geometry and changes in their shape through time. Notwithstanding this variability, there were discernible changes in mean dune height, length and migration rate in response to changing discharge over several days. Changes in mean dune height and length lagged only slightly behind changes in discharge. Therefore, during periods of both steady and unsteady flow, mean dune lengths were quite close to equilibrium values predicted by theoretical models. Mean dune steepnesses were also similar to predicted equilibrium values, except during high, falling flows when the steepness was above that predicted. Variations in mean dune height and length with discharge are similar to those predicted by theory under conditions of low mean dune excursion and discharge variation with a short high water period and long low water period. However, the calculated rates of change of height of individual dunes vary considerably from those measured. Rates of dune creation and destruction were unrelated to discharge variations, contrary to previous results. Instead, creations and destructions were apparently the result of local variations in bed shear stress and sediment transport rate. Observed changes in dune height during unsteady flows agree with theory fairly well at low bed shear stresses, but not at higher bed shear stresses when suspended sediment transport is significant.     

Structure Function Analysis of Two-Scale Scalar Ramps. Part II: Ramp Characteristics and Surface Renewal Flux Estimation

Ramp features in the turbulent scalar field are associated with turbulent coherent structures, which dominate energy and mass fluxes in the atmospheric surface layer. Although finer scale ramp-like shapes embedded within larger scale ramp-like shapes can readily be perceived in turbulent scalar traces, their presence has largely been overlooked in the literature. We demonstrate the signature of more than one ramp scale in structure functions of the turbulent scalar field measured from above bare ground and two types of short plant canopies, using structure-function time lags ranging in scale from isotropic to larger than the characteristic coherent structures. Spectral analysis of structure functions was used to characterize different scales of turbulent structures. By expanding structure function analysis to include two ramp scales, we characterized the intermittency, duration, and surface renewal flux contribution of the smallest (i.e., Scale One) and the dominant (i.e., Scale Two) coherent structure scales. The frequencies of the coherent structure scales increase with mean wind shear, implying that both Scale One and Scale Two are shear-driven. The embedded Scale One turbulent structure scale is ineffectual in the surface-layer energy and mass transport process. The new method reported here for obtaining surface renewal-based scalar exchange works well over bare ground and short canopies under unstable conditions, effectively eliminating the α calibration for these conditions and forming the foundation for analysis over taller and more complex surfaces.     

Simulating future Caspian sea level changes using regional climate model outputs

We report on simulations of present-day climate (1961–1990) and future climate conditions (2071–2100, Special Report on Emissions Scenario A2) over the Caspian sea basin with a regional climate model (RCM) nested in time-slice general circulation model (GCM) simulations. We also calculate changes (A2 scenario minus present-day) in Caspian sea level (CSL) in response to changes in the simulated hydrologic budget of the basin. For the present-day run, both the GCM and RCM show a good performance in reproducing the water budget of the basin and the magnitude of multi-decadal changes in CSL. Compared to present-day climate, in the A2 scenario experiment we find an increase in cold season precipitation and an increase in temperature and evaporation, both over land and over the Caspian sea. We also find a large decrease of CSL in the A2 scenario run compared to the present-day run. This is due to increased evaporation loss from the basin (particularly over the sea) exceeding increased cold season precipitation over the basin. Our results suggest that the CSL might undergo large changes under future climate change, leading to potentially devastating consequences for the economy and environment of the region.     

Ocean tide loading (OTL) displacements from global and local grids: comparisons to GPS estimates over the shelf of Brittany, France

In this paper we examine OTL displacements detected by GPS stations of a dedicated campaign and validate ocean tide models. Our area of study is the continental shelf of Brittany and Cotentin in France. Brittany is one of the few places in the world where tides provoke loading displacements of ∼10–12 cm vertically and a few cm horizontally. Ocean tide models suffer from important discrepancies in this region. Seven global and regional ocean tide models were tested: FES2004 corrected for K2, TPXO.7.0, TPXO.6.2, GOT00.2, CSR4.0, NAO.99b and the most recent regional grids of the North East Atlantic (NEA2004). These gridded amplitudes and phases of ocean tides were convolved in order to get the predicted OTL displacements using two different algorithms. Data over a period of 3.5 months of 8 GPS campaign stations located on the north coast of Brittany are used, in order to evaluate the geographical distribution of the OTL effect. We have modified and implemented new algorithms in our GPS software, GINS 7.1. GPS OTL constituents are estimated based on 1-day batch solutions. We compare the observed GPS OTL constituents of M₂, S₂, N₂ and K₁ waves with the selected ocean tide models on global and regional grids. Large phase-lag and amplitude discrepancies over 20° and 1.5 cm in the vertical direction in the semi-diurnal band of M₂ between predictions and GPS/models are detected in the Bay of Mont St-Michel. From a least squares spectral analysis of the GPS time-series, significant harmonic peaks in the integer multiples of the orbital periods of the GPS satellites are observed, indicating the existence of multipath effects in the GPS OTL constituents. The GPS OTL observations agree best with FES2004, NEA2004, GOT00.2 and CSR4.0 tide models.