Fiber Optic Measurements of Soil Moisture in a Waste Rock Pile

The design and construction of a waste rock pile influences water infiltration and may promote the production of contaminated mine drainage. The objective of this project is to evaluate the use of an active fiber optic distributed temperature sensing (aFO-DTS) protocol to measure infiltration and soil moisture within a flow control layer capping an experimental waste rock pile. Five hundred meters of fiber optic cable were installed in a waste rock pile that is 70 m long, 10 m wide, and was covered with 0.60 m of fine compacted sand and 0.25 m of non-reactive crushed waste rock. Volumetric water content was assessed by heating the fiber optic cable with 15-min heat pulses at 15 W/m every 30 min. To test the aFO-DTS system 14 mm of recharge was applied to the top surface of the waste rock pile over 4 h, simulating a major rain event. The average volumetric water content in the FCL increased from 0.10 to 0.24 over the duration of the test. The volumetric water content measured with aFO-DTS in the FCL and waste rock was within ±0.06 and ±0.03, respectively, compared with values measured using 96 dielectric soil moisture probes over the same time period. Additional results illustrate how water can be confined within the FCL and monitored through an aFO-DTS protocol serving as a practical means to measure soil moisture at an industrial capacity.     

Observation of snowfall with a low-power FM-CW K-band radar (Micro Rain Radar)

Quantifying snowfall intensity especially under arctic conditions is a challenge because wind and snow drift deteriorate estimates obtained from both ground-based gauges and disdrometers. Ground-based remote sensing with active instruments might be a solution because they can measure well above drifting snow and do not suffer from flow distortions by the instrument. Clear disadvantages are, however, the dependency of e.g. radar returns on snow habit which might lead to similar large uncertainties. Moreover, high sensitivity radars are still far too costly to operate in a network and under harsh conditions. In this paper we compare returns from a low-cost, low-power vertically pointing FM-CW radar (Micro Rain Radar, MRR) operating at 24.1?GHz with returns from a 35.5?GHz cloud radar (MIRA36) for dry snowfall during a 6-month observation period at an Alpine station (Environmental Research Station Schneefernerhaus, UFS) at 2,650?m height above sea level. The goal was to quantify the potential and limitations of the MRR in relation to what is achievable by a cloud radar. The operational MRR procedures to derive standard radar variables like effective reflectivity factor (Z _e) or the mean Doppler velocity (W) had to be modified for snowfall since the MRR was originally designed for rain observations. Since the radar returns from snowfall are weaker than from comparable rainfall, the behavior of the MRR close to its detection threshold has been analyzed and a method is proposed to quantify the noise level of the MRR based on clear sky observations. By converting the resulting MRR-Z _e into 35.5?GHz equivalent Z _e values, a remaining difference below 1?dBz with slightly higher values close to the noise threshold could be obtained. Due to the much higher sensitivity of MIRA36, the transition of the MRR from the true signal to noise can be observed, which agrees well with the independent clear sky noise estimate. The mean Doppler velocity differences between both radars are below 0.3?ms^?1. The distribution of Z _e values from MIRA36 are finally used to estimate the uncertainty of retrieved snowfall and snow accumulation with the MRR. At UFS low snowfall rates missed by the MRR are negligible when comparing snow accumulation, which were mainly caused by intensities between 0.1 and 0.8 mm?h^?1. The MRR overestimates the total snow accumulation by about 7%. This error is much smaller than the error caused by uncertain Z _e?Csnowfall rate relations, which would affect the MIRA36 estimated to a similar degree.     

Circumglobal wave train and the summer monsoon over northwestern India and Pakistan: the explicit role of the surface heat low

This study examines the influence of the mid-latitude circulation on the surface heat low (HL) and associated monsoon rainfall over northwestern India and Pakistan using the ERA40 data and high resolution (T106L31) climate model ECHAM5 simulation. Special emphasis is given to the surface HL which forms over Pakistan and adjoining areas of India, Iran and Afghanistan during the summer season. A heat low index (HLI) is defined to depict the surface HL. The HLI displays significant correlations with the upper level mid-latitude circulation over western central Asia and low level monsoon circulation over Arabian Sea and acts as a bridge connecting the mid-latitude wave train to the Indian summer monsoon. A time-lagged singular value decomposition analysis reveals that the eastward propagation of the mid-latitude circumglobal wave train (CGT) influences the surface pressure anomalies over the Indian domain. The largest low (negative) pressure anomalies over the western parts of the HL region (i.e., Iran and Afghanistan) occur in conjunction with the upper level anomalous high that develops over western-central Asia during the positive phase of the CGT. The composite analysis also reveals a significant increase in the low pressure anomalies over Iran and Afghanistan during the positive phase of CGT. The westward increasing low pressure anomalies with its north?Csouth orientation provokes enormous north?Csouth pressure gradient (lower pressure over land than over sea). This in turn enables the moist southerly flow from the Arabian Sea to penetrate farther northward over northwestern India and Pakistan. A monsoon trough like conditions develops over northwestern India and Pakistan where the moist southwesterly flow from the Arabian Sea and the Persian Gulf converge. The convergence in association with the orographic uplifting expedites convection and associated precipitation over northwestern India and Pakistan. The high resolution climate model ECHAM5 simulation also underlines the proposed findings and mechanism.     

A scale-selective multilevel method for long-wave linear acoustics

A new method for the numerical integration of the equations for one-dimensional linear acoustics with large time steps is presented. While it is capable of computing the “slaved” dynamics of short-wave solution components induced by slow forcing, it eliminates freely propagating compressible short-wave modes, which are under-resolved in time. Scale-wise decomposition of the data based on geometric multigrid ideas enables a scale-dependent blending of time integrators with different principal features. To guide the selection of these integrators, the discrete-dispersion relations of some standard second-order schemes are analyzed, and their response to high wave number low frequency source terms are discussed. The performance of the new method is illustrated on a test case with “multiscale” initial data and a problem with a slowly varying high wave number source term.     

Convectively coupled Kelvin and easterly waves in a regional climate simulation of the tropics

This study evaluates the performance of a regional climate model in simulating two types of synoptic tropical weather disturbances: convectively-coupled Kelvin and easterly waves. Interest in these two wave modes stems from their potential predictability out to several weeks in advance, as well as a strong observed linkage between easterly waves and tropical cyclogenesis. The model is a recent version of the weather research and forecast (WRF) system with 36-km horizontal grid spacing and convection parameterized using a scheme that accounts for key convective triggering and inhibition processes. The domain spans the entire tropical belt between 45°S and 45°N with periodic boundary conditions in the east–west direction, and conditions at the meridional/lower boundaries specified based on observations. The simulation covers 6 years from 2000 to 2005, which is long enough to establish a statistical depiction of the waves through space-time spectral filtering of rainfall data, together with simple lagged-linear regression. Results show that both the horizontal phase speeds and three-dimensional structures of the waves are qualitatively well captured by the model in comparison to observations. However, significant biases in wave activity are seen, with generally overactive easterly waves and underactive Kelvin waves. Evidence is presented to suggest that these biases in wave activity (which are also correlated with biases in time–mean rainfall, as well as biases in the model’s tropical cyclone climatology) stem in part from convection in the model coupling too strongly to rotational circulation anomalies. Nevertheless, the model is seen to do a reasonable job at capturing the genesis of tropical cyclones from easterly waves, with evidence for both wave accumulation and critical layer processes being importantly involved.     

Kinetics of H2-O2-H2O redox equilibria and formation of metastable H2O2 under low temperature hydrothermal conditions

Hydrothermal experiments were conducted to evaluate the kinetics of H_2(aq) oxidation in the homogeneous H₂-O₂-H₂O system at conditions reflecting subsurface/near-seafloor hydrothermal environments (55-250 °C and 242-497 bar). The kinetics of the water-forming reaction that controls the fundamental equilibrium between dissolved H_2(aq) and O_2(aq), are expected to impose significant constraints on the redox gradients that develop when mixing occurs between oxygenated seawater and high-temperature anoxic vent fluid at near-seafloor conditions. Experimental data indicate that, indeed, the kinetics of H_2(aq)-O_2(aq) equilibrium become slower with decreasing temperature, allowing excess H_2(aq) to remain in solution. Sluggish reaction rates of H_2(aq) oxidation suggest that active microbial populations in near-seafloor and subsurface environments could potentially utilize both H_2(aq) and O_2(aq), even at temperatures lower than 40 °C due to H_2(aq) persistence in the seawater/vent fluid mixtures. For these H₂-O₂ disequilibrium conditions, redox gradients along the seawater/hydrothermal fluid mixing interface are not sharp and microbially-mediated H_2(aq) oxidation coupled with a lack of other electron acceptors (e.g. nitrate) could provide an important energy source available at low-temperature diffuse flow vent sites.More importantly, when H_2(aq)-O_2(aq) disequilibrium conditions apply, formation of metastable hydrogen peroxide is observed. The yield of H₂O_2(aq) synthesis appears to be enhanced under conditions of elevated H_2(aq)/O_2(aq) molar ratios that correspond to abundant H_2(aq) concentrations. Formation of metastable H₂O₂ is expected to affect the distribution of dissolved organic carbon (DOC) owing to the existence of an additional strong oxidizing agent. Oxidation of magnetite and/or Fe⁺⁺ by hydrogen peroxide could also induce formation of metastable hydroxyl radicals (•OH) through Fenton-type reactions, further broadening the implications of hydrogen peroxide in hydrothermal environments.     

Stable Cu isotope fractionation in soils during oxic weathering and podzolization

Copper stable isotope ratios are fractionated during various biogeochemical processes and may trace the fate of Cu during long-term pedogenetic processes. We assessed the effects of oxic weathering (formation of Cambisols) and podzolization on Cu isotope ratios (δ⁶⁵Cu). Two Cambisols (oxic weathered soils without strong vertical translocations of soil constituents) and two Podzols (soils showing vertical translocation of organic matter, Fe and Al) were analyzed for Cu concentrations, partitioning of Cu in seven fractions of a sequential extraction and δ⁶⁵Cu values in bulk soil. Cu concentrations in the studied soils were low (1.4-27.6 μg g⁻¹) and Cu was mainly associated with strongly bound Fe oxide- and silicate-associated forms. Bulk δ⁶⁵Cu values varied between −0.57‰ and 0.44‰ in all studied horizons. The O horizons had on average significantly lighter Cu isotope compositions (−0.21‰) than the A horizons (0.13‰) which can either be explained by Cu isotope fractionation during cycling through the plants or deposition of isotopically light Cu from the atmosphere. Oxic weathering without pronounced podzolization in both Cambisols and a weakly developed Podzol (Haplic Podzol 2) caused no significant isotope fractionation in the single profiles, while a slight tendency to lower δ⁶⁵Cu values with depth was visible in all four profiles. This is the opposite depth distribution of δ⁶⁵Cu values to that we observed in hydromorphic soils (soils which show indication of redox changes because of the influence of water saturation) in a previous study. In a more pronounced Podzol (Haplic Podzol 1), δ⁶⁵Cu values and Cu concentrations decreased from Ah to E horizons and increased again deeper in the soil. Humus-rich sections of the Bhs horizon had higher Cu concentrations (2.8 μg g⁻¹) and a higher δ⁶⁵Cu value (−0.18‰) than oxide-rich sections (1.9 μg g⁻¹, −0.35‰) suggesting Cu translocation between E and B horizons as organo-Cu complexes. The different depth distributions in oxic weathered and hydromorphic soils and the pronounced vertical differences in δ⁶⁵Cu values in Haplic Podzol 1 indicate a promising potential of δ⁶⁵Cu values to improve our knowledge of the fate of Cu during long-term pedogenetic processes.     

Summer temperature in the eastern part of southern South America: its variability in the twentieth century and a teleconnection with Oceania

The 1907–2001 summer-to-summer surface air temperature variability in the eastern part of southern South America (SSA, partly including Patagonia) is analysed. Based on records from instruments located next to the Atlantic Ocean (36°S–55°S), we define indices for the interannual and interdecadal timescales. The main interdecadal mode reflects the late-1970s cold-to-warm climate shift in the region and a warm-to-cold transition during early 1930s. Although it has been in phase with the Pacific Decadal Oscillation (PDO) index since the 1960s, they diverged in the preceding decades. The main interannual variability index exhibits high spectral power at ~3.4 years and is representative of temperature variability in a broad area in the southern half of the continent. Eleven-years running correlation coefficients between this index and December-to-February (DJF) Niño3.4 show significant decadal fluctuations, out-of-phase with the running correlation with a DJF index of the Southern Annular Mode. The main interannual variability index is associated with a barotropic wavetrain-like pattern extending over the South Pacific from Oceania to SSA. During warm (cold) summers in SSA, significant anticyclonic (cyclonic) anomalies tend to predominate over eastern Australia, to the north of the Ross Sea, and to the east of SSA, whereas anomalous cyclonic (anticyclonic) circulation is observed over New Zealand and west of SSA. This teleconnection links warm (cold) SSA anomalies with dry (wet) summers in eastern Australia. The covariability seems to be influenced by the characteristics of tropical forcing; indeed, a disruption has been observed since late 1970s, presumably due to the PDO warm phase.     

Modeling Ambiguity in Census Microdata Allocations to Improve Demographic Small Area Estimates

This article describes a methodology for allocating demographic microdata to small enumeration areas such as census tracts, in the presence of underlying ambiguities. Maximum Entropy methods impute population weights that are constrained to match a set of census tract‐level summary statistics. Once allocated, the household characteristics are summarized to revise estimates of tract‐level demographic summary statistics, and to derive measures of ambiguity. The revised summary statistics are compared with original tract summaries within a context of expected variation. Allocation ambiguity is quantified for each household as a function of the distribution of imputed sample weights over all census tracts, and by computed metrics of confusion and variety of allocation to any census tract. The process reported here allows differentiation of households with regard to inherent ambiguity in the allocation decision. Ambiguity assessment represents an important component that has been neglected in spatial allocation work to date but can be seen as important additional knowledge for demographers and users of small area estimates. For the majority of tested variables, the revised tract level summaries correlate highly with original tract summary statistics. In addition to assessments for individual households, it is also possible to compute average allocation ambiguity for individual tracts, and to associate this with demographic characteristics not utilized in the allocation process.     

Road network selection for medium scales using an extended stroke-mesh combination algorithm

The road network is an essential feature class in topographic maps and databases. Road network selection for smaller scales forms a prerequisite for all other generalization operators and is thus a fundamental operation in the overall process of topographic map and database production. The objective of this paper was to develop an algorithm for automated road network selection from a large-scale (1:10,000) to a medium-scale database (1:50,000). The project was pursued in collaboration with swisstopo, the national mapping agency of Switzerland. Three algorithms (a stroke-based, a mesh-based, and a combined stroke-mesh algorithm) were implemented from the literature and analyzed using swisstopo’s large-scale TLM3D spatial database, with requirements set forth by expert cartographers. Initial experiments showed that the combination algorithm performed best, yet still it could not meet all requirements. Therefore, extensions to the basic stroke-mesh algorithm were developed, significantly improving the selection result with real-world, large test databases. Three extensions introduce modifications to the stroke-mesh combination algorithm. Furthermore, two extensions include external feature classes, ensuring accessibility of points of interest and appropriate network density representation in settlement areas, respectively. The results were evaluated by expert cartographers, who concluded that the proposed approach is ready to be deployed in production at swisstopo.