Fate and transport of pollutants through a municipal solid waste landfill leachate in Sri Lanka

This study focuses on the characterization of leachate generated from Gohagoda dumpsite in Kandy, Sri Lanka, assessment of its spatial and temporal variations, and identification of subsurface canals and perched water bodies in the wetland system affected by the leachate flow. Leachate samples were collected monthly throughout dry and rainy seasons from different points of the leachate drainage channel over a period of 1 year and they were tested for quality parameters: pH, temperature, electrical conductivity, total dissolved soils, alkalinity, hardness, total solids, volatile solids, total suspended solids, volatile suspended solids, biochemical oxygen demand (BOD₅), chemical oxygen demand, nitrate-nitrogen, nitrite-nitrogen, phosphates, ammonium-nitrogen, chloride, dissolved organic carbon, total organic carbon and heavy metals. Sequential soil extraction procedures were performed for the characterization of leachate-affected local soil. A geophysical survey using direct current resistivity technique was conducted at locations downstream of the dumpsite. Leachate characteristics indicated that the leachate is in the methanogenic phase and the results strongly suggest that the leachate may be polluting the river where the leachate is discharged directly. Leachate exceeds the allowable limits of Sri Lankan wastewater discharge standards for many of the parameters. Significant difference (P < 0.05) was observed for most of organic and inorganic parameters among all sampling locations. Many parameters showed a negative correlation with pH. The affected soils showed high heavy metal concentrations. Resistivity study confirmed a confined leachate flow at the near surface with few subsurface canals. However, no separate subsurface plume movement was observed. The results of this research can effectively be used for the establishment of an efficient and effective treatment method for the Gohagoda landfill leachate.     

Can a Training Image Be a Substitute for a Random Field Model?

In most multiple-point simulation algorithms, all statistical features are provided by one or several training images (TI) that serve as a substitute for a random field model. However, because in practice the TI is always of finite size, the stochastic nature of multiple-point simulation is questionable. This issue is addressed by considering the case of a sequential simulation algorithm applied to a binary TI that is a genuine realization of an underlying random field. At each step, the algorithm uses templates containing the current target point as well as all previously simulated points. The simulation is validated by checking that all statistical features of the random field (supported by the simulation domain) are retrieved as an average over a large number of outcomes. The results are as follows. It is demonstrated that multiple-point simulation performs well whenever the TI is a complete (infinitely large) realization of a stationary, ergodic random field. As soon as the TI is restricted to a limited domain, the statistical features cannot be obtained exactly, but integral range techniques make it possible to predict how much the TI should be extended to approximate them up to a prespecified precision. Moreover, one can take advantage of extending the TI to reduce the number of disruptions in the execution of the algorithm, which arise when no conditioning template can be found in the TI.     

Analysis of coprolites from the extinct mountain goat Myotragus balearicus

Humans colonized the Balearic Islands 5–4 ka ago. They arrived in a uniquely adapted ecosystem with the Balearic mountain goat Myotragus balearicus (Bovidae, Antilopinae, Caprini) as the only large mammal. This mammal went extinct rapidly after human arrival. Several hypotheses have been proposed to explain the extinction of M. balearicus. For the present study ancient DNA analysis (Sanger sequencing, Roche-454, Ion Torrent), and pollen and macrofossil analyses were performed on preserved coprolites from M. balearicus, providing information on its diet and paleo-environment. The information retrieved shows that M. balearicus was heavily dependent on the Balearic box species Buxus balearica during at least part of the year, and that it was most probably a browser. Hindcast ecological niche modelling of B. balearica shows that local distribution of this plant species was affected by climate changes. This suggests that the extinction of M. balearicus can be related to the decline and regional extinction of a plant species that formed a major component of its diet. The vegetation change is thought to be caused by increased aridity occurring throughout the Mediterranean. Previous hypotheses relating the extinction of M. balearicus directly to the arrival of humans on the islands must therefore be adjusted.     

Age and impacts of the caldera-forming Aniakchak II eruption in western Alaska

The mid-Holocene eruption of Aniakchak volcano (Aniakchak II) in southwest Alaska was among the largest eruptions globally in the last 10,000 years (VEI-6). Despite evidence for possible impacts on global climate, the precise age of the eruption is not well-constrained and little is known about regional environmental impacts. A closely spaced sequence of radiocarbon dates at a peatland site over 1000 km from the volcano show that peat accumulation was greatly reduced with a hiatus of approximately 90–120 yr following tephra deposition. During this inferred hiatus no paleoenvironmental data are available but once vegetation returned the flora changed from a Cyperaceae-dominated assemblage to a Poaceae-dominated vegetation cover, suggesting a drier and/or more nutrient-rich ecosystem. Oribatid mites are extremely abundant in the peat at the depth of the ash, and show a longer-term, increasingly wet peat surface across the tephra layer. The radiocarbon sample immediately below the tephra gave a date of 1636–1446 cal yr BC suggesting that the eruption might be younger than previously thought. Our findings suggest that the eruption may have led to a widespread reduction in peatland carbon sequestration and that the impacts on ecosystem functioning were profound and long-lasting.     

Geometry and Properties of the Excavation-Induced Fractures at the Meuse/Haute-Marne URL Drifts

The performance and safety assessment and technology demonstration are the main objectives of research programs for feasibility studies for deep geological repository of radioactive waste. In this context, the French national radioactive waste management agency (ANDRA) started to develop the Meuse/Haute-Marne underground research laboratory (URL) at Bure, nearly 300 km East of Paris. The host formation consists of a Callovo-Oxfordian claystone found between 420 and 550 m below ground, overlain and underlain by poorly permeable carbonate formations. One of the major concerns related to performance assessment is the excavation-induced fractures which may provide groundwater preferential pathway for radionuclide migration. The extent of the fractures possibly acting significantly in the radionuclide migration is known as the excavation damaged zone (EDZ). A scientific study on the EDZ characterization is performed at the main level of the URL (?490 m). Observations such as structural analysis on core, overcored resin-filled samples, geological survey of the drift face and sidewalls, were made to better understand the fracture network characteristics, extent and its generation. Pulse and constant head test hydraulic conductivity measurements were performed with multi packer system to estimate the extension of the EDZ hydraulic conductivity. Fractures exhibited high transmissivity near the excavation walls, but farther from the exaction walls, shear fractures showed hydraulic conductivity values reflecting values of undisturbed or slightly disturbed rock mass condition. The major findings in terms of geometry and properties of excavation-induced fractures are discussed in detail in this paper. For example, it is observed that the shape of the fracture network depends on the orientation of the drift in relation to the orientation of the in situ stress field.     

Multifrequency algorithms for precise point positioning: MAP3

We present the new MAP3 algorithms to perform static precise point positioning (PPP) from multifrequency and multisystem GNSS observations. MAP3 represents a two-step strategy in which the least squares theory is applied twice to estimate smoothed pseudo-distances, initial phase ambiguities, and slant ionospheric delay first, and the absolute receiver position and its clock offset in a second adjustment. Unlike the classic PPP technique, in our new approach, the ionospheric-free linear combination is not used. The combination of signals from different satellite systems is accomplished by taking into account the receiver inter-system bias. MAP3 has been implemented in MATLAB and integrated within a complete PPP software developed on site and named PCube. We test the MAP3 performance numerically and contrast it with other external PPP programs. In general, MAP3 positioning accuracy with low-noise GPS dual-frequency observations is about 2.5 cm in 2-h observation periods, 1 cm in 10 h, and 7 mm after 1 day. This means an improvement in the accuracy in short observation periods of at least 7 mm with respect to the other PPP programs. The MAP3 convergence time is also analyzed and some results obtained from real triple-frequency GPS and GIOVE observations are presented.     

Charge distribution on the atmospheric aerosol over the North Atlantic and the Indian Ocean

The distribution of electric charge on the marine aerosol was determined near the sea surface of the Indian Ocean and the North Atlantic during the final stage of the Snellius II-expedition. Mean values for small ion concentrationsn ⁺=455 cm^–3 andn ^–=340 cm^–3 were found over the Atlantic, whilen ⁺=310 andn ^–=250 cm^–3 were the mean values over the Indian Ocean. The ration ⁺/n^– increased from 1.2 to 2.0 with decreasing wind velocity. At wind velocities below 5 m/s 75% of the net space charge near the ocean surface was found to be carried by small ions.     

Thresholds of intrabed flow and other interactions of turbidity currents with soft muddy substrates

Controlled laboratory experiments reveal that the lower part of turbidity currents has the ability to enter fluid mud substrates, if the bed shear stress is higher than the yield stress of the fluid mud and the density of the turbidity current is higher than the density of the substrate. Upon entering the substrate, the turbidity current either induces mixing between flow‐derived sediment and substrate sediment, or it forms a stable horizontal flow front inside the fluid mud. Such ‘intrabed’ flow is surrounded by plastically deformed mud; otherwise it resembles the front of a ‘bottom‐hugging’ turbidity current. The ‘suprabed’ portion of the turbidity current, i.e. the upper part of the flow that does not enter the substrate, is typically separated from the intrabed flow by a long horizontal layer of mud which originates from the mud that is swept over the top of the intrabed flow and then incorporated into the flow. The intrabed flow and the mixing mechanism are specific types of interaction between turbidity currents and muddy substrates that are part of a larger group of interactions, which also include bypass, deposition, erosion and soft sediment deformation. A classification scheme for these types of interactions is proposed, based on an excess bed shear stress parameter, which includes the difference in the bed shear stress imposed by the flow and the yield stress of the substrate and an excess density parameter, which relies on the density difference between the flow and the substrate. Based on this classification scheme, as well as on the sedimentological properties of the laboratory deposits, an existing facies model for intrabed turbidites is extended to the other types of interaction involving soft muddy substrates. The physical threshold of flow‐substrate mixing versus stable intrabed flow is defined using the gradient Richardson number, and this method is validated successfully with the laboratory data. The gradient Richardson number is also used to verify that stable intrabed flow is possible in natural turbidity currents, and to determine under which conditions intrabed flow is likely to be unstable. It appears that intrabed flow is likely only in natural turbidity currents with flow velocities well below ca 3·5 m s^?1, although a wider range of flows is capable of entering fluid muds. Below this threshold velocity, intrabed flow is stable only at high‐density gradients and low‐velocity gradients across the upper boundary of the turbidity current. Finally, the gradient Richardson number is used as a scaling parameter to set the flow velocity limits of a natural turbidity current that formed an inferred intrabed turbidite in the deep‐marine Aberystwyth Grits Group, West Wales, United Kingdom.     

Change of Support for Estimating Local Block Grade Distributions

An important aspect in mineral resource evaluation is the reduction of variance when post-processing the grade distributions defined on the support (volume) of the available data into distributions defined on the support of the proposed selective mining units. Although the volume-variance relationship is well understood for the estimation of global grade distributions, it is still an unsolved issue for local estimation studies based on non-parametric geostatistical methods, such as indicator kriging, for which the support correction is not inherent to the method. To clarify this relationship, the local change of support problem is examined in the scope of two parametric models (multi-Gaussian and discrete Gaussian models). It is shown that the variance reduction factor between point and block-support local distributions depends on the block being considered and is less than the global variance reduction factor. As a consequence, post-processing the local point-support grade distributions on the basis of the latter systematically understates the importance of the change of support at the local scale and makes selective mining appear more economically attractive than it really is. In the light of these results, a methodology is proposed to post-process the local point-support distributions obtained via non-parametric (indicator) methods into block-support distributions. An application to simulated data indicates that this methodology provides an accurate estimation at the block support when dealing with diffusion-type random fields.