Explorative data analysis of heavy metal contaminated soil using multidimensional spatial regression

To obtain data on heavy metal contaminated soil requires laborious and time-consuming data sampling and analysis. Not only has the contamination to be measured, but also additional data characterizing the soil and the boundary conditions of the site, such as pH, land use, and soil fertility. For an integrative approach, combining the analysis of spatial distribution, and of factors influencing the contamination, and its treatment, the Mollifier interpolation was used, which is a non-parametric kernel density regression. The Mollifier was capable of including additional independent variables (beyond the spatial dimensions x and y) in the spatial interpolation and hence explored the combined influence of spatial and other variables, such as land use, on the heavy metal distribution. The Mollifier could also represent the interdependence between different heavy metal concentrations and additional site characteristics. Although the uncertainty measure supplied by the Mollifier at first seems somewhat unusual, it is a valuable feature and supplements the geostatistical uncertainty assessment.     

Ellenbergerite,a new high-pressure Mg-Al-(Ti,Zr)-silicate with a novel structure based on face-sharing octahedra

Ellenbergerite occurs as purple millimetre-size grains associated with talc, kyanite, clinochlore, rutile, and zircon in composite inclusions within decimetre-large pyrope crystals (90–98 mole percent end-member) in the quartzite layer of the Dora Maira massif, Western Alps, from which coesite has been recently reported (Chopin 1984). It is hexagonal, a=12.255(8), c=4.932(4) Å, Z=1, space group P6₃. Mohs hardness 6.5; D_mes 3.15, D_cal 3.10; no cleavage. Uniaxial negative and vividly pleochroic, colourless, colourless to deep lilac with colour zoning. The intensely coloured variety has 1.6789(5), 1.670(1); microprobe analysis yields SiO₂ 39.1, P₂O₅ 0.45, Al₂O₃ 25.1, TiO₂ 4.0, MgO 22.2, FeO 0.20, sum 99.05 wt.% including H₂O 8.0 (coulometrically). The formula calculated on a O₂₈(OH)₁₀ basis (implying 7.5 wt.% H₂O) is Mg_6.71 Fe_0.03 Ti_0.61 Al_6.00 Si_7.92 P_0.08 O₂₈(OH)₁₀ The colour zoning is due to nearly complete TiZr substitution. In addition ellenbergerite may contain more than 8 wt.% P₂O₅ with strictly correlated changes of Si, Mg, Al and Ti+Zr contents, over 80% of which represent the SiAlPMg substitution.The structure has been determined from 1049 observed independent reflections and refined to R=0.034, R_w=0.031, including six of ten protons. It consists of single chains of face-sharing octahedra with one third vacancies extending along the six-fold screw axes, and of pairs of fully occupied face-sharing octahedra linked by edge-sharing to form octahedral double chains parallel to the twofold screw axes, all interconnected by SiO₄ tetrahedra. It may be compared with the dumortierite polymorph with space group P6₃mc derived hypothetically by Moore and Araki (1978). The structural formula is (Mg,Ti,Zr,)₂ Mg₆(Al,Mg)₆ (Si,P)₂ Si₆ O₂₈(OH)₁₀ Face-sharing octahedra are an unusual feature in silicates which results in a dense structure and reflects, considering the common bulk composition, the uncommon high-pressure formation conditions (about 25–30 kbar, 700–800° C). Ti⁴⁺-Fe²⁺ charge transfer between face-sharing octahedra on the six-fold screw axes most likely accounts for the absorption scheme.     

Spatial statistical techniques for aggregating point objects extracted from high spatial resolution remotely sensed imagery

 Using a local maximum filter, individual trees were extracted from a 1 m spatial resolution IKONOS image and represented as points. The spatial pattern of individual trees was determined to represent forest age (a surrogate for forest structure). Point attributes, based on the spatial pattern of trees, were generated via nearest neighbour statistics and used as the basis for aggregating points into forest structure units. The forest structure units allowed for the mapping of a forested area into one of three age categories: young (1–20 years), intermediate (21–120 years), and mature (>120 years). This research indicates a new approach to image processing, where objects generated from the processing of image data (rather than pixels or spectral values) are subjected to spatial statistical analysis to estimate an attribute relating an aspect of forest structure. Received: 22 April 2002 / Accepted: 23 November 2002     

Non-conservative behaviour of molybdenum in coastal waters: Coupling geochemical, biological, and sedimentological processes

Non-conservative behaviour of dissolved Mo was observed during specific time periods in the water column of the Wadden Sea of NW Germany. In July 2005 dissolved Mo declined within 36 h from a level only slightly below seawater (82 nM) to a minimum value of 30 nM, whereas in August 2002 dissolved Mo revealed a tidal cyclicity with maximum values up to 158 nM at low tide. In contrast, cruises in August 2003 and 2004 displayed an almost conservative behaviour of Mo. The decrease in dissolved Mo during July 2005 and elevated values in August 2002 were accompanied by Mo enrichments on aggregates in the water column of the Wadden Sea. Along with Mo, dissolved Mn showed unusual concentration patterns in July 2005, with values distinctly below the common summer level (by a factor of five). A direct relation between the loss of Mo and scavenging by freshly formed MnO_x phases could not be inferred from our data because both metals revealed inverse patterns. Parallel to decreasing dissolved Mo concentrations dissolved Mn showed an increasing trend while particulate Mn decreased. Such finding is compatible with the formation of oxygen-depleted zones in aggregates, which provide suitable conditions for the rapid fixation of Mo and parallel release of Mn by chemically and/or microbially mediated processes. Our assumption is supported by biological (e.g. number of aggregate-associated bacteria) and sedimentological (e.g. aggregate abundance and size) parameters. The production of organic components (e.g. TEP) during breakdown of an algae bloom in July 2005 led to the formation of larger Mo-enriched aggregates, thus depleting the water column in dissolved Mo. After deposition on and incorporation into sandy tidal flats these aggregates are rapidly decomposed by microbial activity. Pore water profiles document that during microbial decomposition of these aggregates, substantial amounts of Mo are released and may replenish and even enrich Mo in the open water column. We postulate a conceptual model for the observed non-conservative behaviour of Mo in coastal waters, which is based on the tight coupling of geochemical, biological, and sedimentological processes.     

Geological indications for Palaeogene uplift in the eastern North Sea Basin

The timing and effect of the Cenozoic uplift of Scandinavia has been investigated using a multi-disciplinary approach involving sedimentological, seismic and biostratigraphic data from the Danish and the adjacent Norwegian parts of the North Sea Basin. It is concluded that significant uplift took place periodically throughout the Palaeogene possibly marking an earlier onset of the so-called “Neogene uplift” of Scandinavia. This conclusion is based on a number of sedimentological observations, including smectite content, grain-size variations, kaolinite thermal stabilities and T_max values supported by seismic reflection geometries and biostratigraphic data. These data indicate several phases of re-working of Palaeogene and older sediments situated further to the east and northeast during the middle to late Eocene and during the middle to late Oligocene. The tectonic patterns were similar during the late Paleocene and the Oligocene with some inversion taking place, whereas no inversion has been observed during the Eocene. Main provenance areas were to the north and northeast during the Paleocene and Oligocene, whereas the Eocene sediments originate mainly from the British Isles to the west. It is proposed that Palaeogene uplift of Scandinavia was associated with regional tectonic movements along crustal zones of weakness, which were reactivated as they accommodated strain induced by the Alpine Orogeny and the opening of the North Atlantic.     

Post‐1850 changes in Glacier Benito,North Patagonian Icefield,Chile

Glacier Benito is a temperate outlet glacier on the west side of the North Patagonian Icefield. Rates of thinning and ablation were obtained using data collected by the British Joint Services Expedition in 1972/73 and subsequent data collected in 2007 and 2011. Ice‐front recession rates were based on dendrochronological dating for the terminal moraines and aerial and satellite imagery of the ice front in 1944, 1998 and 2002. Between the first Benito survey in 1973 and 2007, the lower glacier thinned by nearly 150 m at an average rate of 4.3 m yr^?1 with the rate increasing to 6.1 m yr^?1 between 2007 and 2011, a 28.7% increase during the latter period. Increases in ice movement and ablation were negligible: ice movement for 1973 and 2007 averaged 0.45 m day^?1 and ablation averaged 0.05 m day^?1. Ice front recession along the glacier's centre line from 1886 to 2002 was approximately 1860 m. Retreat rates between 1886 and 1944 averaged 8.9 m yr^?1. Thereafter glacier asymmetry makes measurement along the glacier centre line unrepresentative of areal change until 1998 when symmetry was restored; retreat between 1944 and 1998 was 15.4 m yr^?1. From 1998 to 2002 the rate increased dramatically to 127.2 m yr^?1. Recession from the southern end of Benito's terminal moraine in the 1850s supports an early date for initial retreat of the Icefield's glaciers.     

Propagation of seasonal temperature signals into an aquifer upon bank infiltration

Infiltrating river water carries the temperature signal of the river into the adjacent aquifer. While the diurnal temperature fluctuations are strongly dampened, the seasonal fluctuations are much less attenuated and can be followed into the aquifer over longer distances. In one-dimensional model with uniform properties, this signal is propagated with a retarded velocity, and its amplitude decreases exponentially with distance. Therefore, time shifts in seasonal temperature signals between rivers and groundwater observation points may be used to estimate infiltration rates and near-river groundwater velocities. As demonstrated in this study, however, the interpretation is nonunique under realistic conditions. We analyze a synthetic test case of a two-dimensional cross section perpendicular to a losing stream, accounting for multi-dimensional flow due to a partially penetrating channel, convective-conductive heat transport within the aquifer, and heat exchange with the underlying aquitard and the land surface. We compare different conceptual simplifications of the domain in order to elaborate on the importance of different system elements. We find that temperature propagation within the shallow aquifer can be highly influenced by conduction through the unsaturated zone and into the underlying aquitard. In contrast, regional groundwater recharge has no major effect on the simulated results. In our setup, multi-dimensionality of the flow field is important only close to the river. We conclude that over-simplistic analytical models can introduce substantial errors if vertical heat exchange at the aquifer boundaries is not accounted for. This has to be considered when using seasonal temperature fluctuations as a natural tracer for bank infiltration.