A statistical test on the local effects of spatially structured variance

Spatial variance is an important characteristic of spatial random variables. It describes local deviations from average global conditions and is thus a proxy for spatial heterogeneity. Investigating instability in spatial variance is a useful way of detecting spatial boundaries, analysing the internal structure of spatial clusters and revealing simultaneously acting geographic phenomena. Recently, a corresponding test statistic called ‘Local Spatial Heteroscedasticity’ (LOSH) has been proposed. This test allows locally heterogeneous regions to be mapped and investigated by comparing them with the global average mean deviation in a data set. While this test is useful in stationary conditions, its value is limited in a global heterogeneous state. There is a risk that local structures might be overlooked and wrong inferences drawn. In this paper, we introduce a test that takes account of global spatial heterogeneity in assessing local spatial effects. The proposed measure, which we call ‘Local Spatial Dispersion’ (LSD), adapts LOSH to local conditions by omitting global information beyond the range of the local neighbourhood and by keeping the related inferential procedure at a local level. Thereby, the local neighbourhoods might be small and cause small-sample issues. In the view of this, we recommend an empirical Bayesian technique to increase the data that is available for resampling by employing empirical prior knowledge. The usefulness of this approach is demonstrated by applying it to a Light Detection and Ranging-derived data set with height differences and by making a comparison with LOSH. Our results show that LSD is uncorrelated with non-spatial variance as well as local spatial autocorrelation. It thus discloses patterns that would be missed by LOSH or indicators of spatial autocorrelation. Furthermore, the empirical outcomes suggest that interpreting LOSH and LSD together is of greater value than interpreting each of the measures individually. In the given example, local interactions can be statistically detected between variance and spatial patterns in the presence of global structuring, and thus reveal details that might otherwise be overlooked.     

Slab melt as metasomatic agent in island arc magma mantle sources, Negros and Batan (Philippines)

Abstract Two new cases of association of adakites with ‘normal’ island arc lavas and transitional adakites are recognized in the islands of Batan and Negros in northern and central Philippines, respectively. The Batan lavas are related to the subduction of the middle Miocene portion of the South China Sea basin along the Manila trench; those of Negros come from the almost aseismic subduction of the middle Miocene Sulu Sea crust along the Negros trench. The occurrence of the Batan adakites is consistent with previous findings showing adakitic glass inclusions within minerals of mantle xenoliths associated with Batan arc lavas. The similarity of adakite ages (1.09 Ma) and that of the metasomatized xenoliths (1 Ma) suggests that both are linked to the same slab‐melting and metasomatic event. Earlier Sr, Pb and Nd‐isotopic studies, however, also reveal the presence of an important sediment contribution to the Batan lava geochemistry. Thus, the role played by slab melts, assumed to have mid‐ocean ridge basalts‐like (MORB) isotopic characteristics, in enriching the Batan subarc mantle is largely masked by the sediment input. The Negros adakites are present only in Mount Cuernos, the volcanic center nearest to the Negros trench. Batch partial melting calculations show that the Negros adakites could be derived from a garnet amphibolitic source with normal‐MORB (N‐MORB) geochemistry. This is supported by the MORB‐like isotopic characteristics of the Mount Cuernos lavas. The volcanic rocks from the other volcanoes consist of normal arc and transitional adakitic lavas that have slightly higher Sr‐ and Pb‐isotopic ratios, probably due to slight sediment input. Mixing of adakites and normal arc lavas to produce transitional adakites is only partly supported by trace element geochemistry and not by field evidence. The transitional adakites can be modeled as partial melts of an adakite‐enriched mantle. Trace element enrichment of non‐adakitic lavas could reflect the interaction of their mantle source with uprising slab melts, as metasomatic mantle minerals scavenge certain trace elements from the adakitic fluids. Therefore, in arcs beneath which thick (up to 2 km) continent‐derived detrital sediments are involved in subduction, like in Batan, the sediment signature can overwhelm the slab melt input. In arcs like Negros where slow subduction could cause a more efficient scraping of thinner (approximately 1 km) detrital sediments, the contribution of slab melts is easier to detect.     

Collision, subduction and accretion events in the Philippines: A synthesis

Abstract The Philippines preserves evidence of the superimposition of tectonic processes in ancient and present‐day collision and subduction zone complexes. The Baguio District in northern Luzon, the Palawan–Central Philippine region and the Mati–Pujada area in southeastern Mindanao resulted from events related to subduction polarity reversal leading to trench initiation, continent‐arc collision and autochthonous oceanic lithosphere emplacement, respectively. Geological data on the Baguio District in Northern Luzon reveal an Early Miocene trench initiation for the east‐dipping Manila Trench. This followed the Late Oligocene cessation of subduction along the west‐dipping proto‐East Luzon Trough. The Manila Trench initiation, which is modeled as a consequence of the counter‐clockwise rotation of Luzon, is attributed to the collision of the Palawan microcontinental block with the Philippine Mobile Belt. In the course of rotation, Luzon onramped the South China Sea crust, effectively converting the shear zone that bounded them into a subduction zone. Several collision‐related accretionary complexes (e.g. Romblon, Mindoro) are present in the Palawan–Central Philippine region. The easternmost collision zone boundary is located east of the Romblon group of islands. The Early Miocene southwestward shift of the collision boundary from Romblon to Mindoro started to end by the Pliocene. Continuous interaction between the Palawan microcontinental block and the Philippine Mobile Belt is presently taken up again along the collisional boundary east of the Romblon group of islands. The Mati–Pujada Peninsula area, on the other hand, is underlain by the Upper Cretaceous Pujada Ophiolite. This supra‐subduction zone ophiolite is capped by chert and pelagic limestones which suggests its derivation from a relatively deep marginal basin. The Pujada Ophiolite could be a part of a proto‐Molucca Sea plate. The re‐interpretation of the geology and tectonic settings of the three areas reaffirm the complex geodynamic evolution of the Philippine archipelago and addresses some of its perceived geological enigmas.     

Investigations on Aral Sea Regressions from Mirabilite Deposits and Remote Sensing

Remote sensing techniques including radar (Topex/Poseidon, Jason-1 and Envisat) and laser altimetry (Icesat), and moderate resolution spectro-radiometer (MODIS) images, are used to estimated current level and surface extent time variations of the Aral Sea. During the Holocene several phases of regression occurred, leading to desiccation of the Aral Sea. During the last 50 years, Aral Sea has drastically shrunk due to intense use of river’s water for irrigation purposes. It is currently separated into four distinct water bodies, namely, the Small Aral in the North, the Tchebas Bay in the North West, and the South West and the South East basins. The Kulandy strait connected the SW and SE basins until very recent times. These basins are now almost separated and salinity becomes very high (140–180 g/l) in the Eastern part. Rubanov discovered past deposits of mirabilite in the years 1970–1980. We investigate the significance of these deposits in the light of current evolution of the four water bodies that constitute the heritage of Aral Sea contemporary desiccation. Using remote sensing techniques, we have attempted to calculate the water balance of south Aral Sea during the last 3 years. We conclude in strong probability that the Kulandy strait carries water most of the time from the Eastern Basin to the Western Basin. We have demonstrated that it should have been the same process in the past to explain the Mirabilite deposit, but unfortunately, due to recent artificial water monitoring of the Aral Sea (dam in the Berg’s strait, new reservoirs in the Amu Darya’s delta), it is impossible to make definitive conclusion from actual Aral Sea water balance.     

Upward continuation of Dome-C airborne gravity and comparison with GOCE gradients at orbit altitude in east Antarctica

An airborne gravity campaign was carried out at the Dome-C survey area in East Antarctica between the 17th and 22nd of January 2013, in order to provide data for an experiment to validate GOCE satellite gravity gradients. After typical filtering for airborne gravity data, the cross-over error statistics for the few crossing points are 11.3 mGal root mean square (rms) error, corresponding to an rms line error of 8.0 mGal. This number is relatively large due to the rough flight conditions, short lines and field handling procedures used. Comparison of the airborne gravity data with GOCE RL4 spherical harmonic models confirmed the quality of the airborne data and that they contain more high-frequency signal than the global models. First, the airborne gravity data were upward continued to GOCE altitude to predict gravity gradients in the local North-East-Up reference frame. In this step, the least squares collocation using the ITGGRACE2010S field to degree and order 90 as reference field, which is subtracted from both the airborne gravity and GOCE gravity gradients, was applied. Then, the predicted gradients were rotated to the gradiometer reference frame using level 1 attitude quaternion data. The validation with the airborne gravity data was limited to the accurate gradient anomalies (T_XX, T_YY, T_ZZ and T_XZ) where the long-wavelength information of the GOCE gradients has been replaced with GOCO03s signal to avoid contamination with GOCE gradient errors at these wavelengths. The comparison shows standard deviations between the predicted and GOCE gradient anomalies T_XX, T_YY, T_ZZ and T_XZ of 9.9, 11.5, 11.6 and 10.4 mE, respectively. A more precise airborne gravity survey of the southern polar gap which is not observed by GOCE would thus provide gradient predictions at a better accuracy, complementing the GOCE coverage in this region.     

Aqueous and petroleum fluid flow associated with sand injectites

Field, petrographic and fluid inclusion characteristics of sand injectites from five outcrop localities and from the subsurface of the Tertiary of the south Viking Graben are described. Although the case studies are from a wide variety of sedimentological, stratigraphic and tectonic settings, and hence their diagenetic evolutions differ significantly, it is possible and useful to assign diagenetic events to three distinct phases of fluid flow associated with sand injectites in sedimentary basins. Firstly, there is fluid flow associated with the injection of the fluid–sediment mix during shallow burial. Early diagenetic imprints in sand injectites reveal that basinal fluids, which may be released during movement along deeper‐seated faults, can be associated with this process and thus the injection process may reveal information on the timing of basin‐scale movement of fluids. Secondly, following the injection process, basinal fluids continue to migrate through uncemented injectites and mix with the ambient meteoric and/or marine pore fluids that invade injectites from the overlying and surrounding host sediments. Early, often pervasive, carbonate cementation is common within sand injectites and rapidly turns sand injectites into flow barriers during shallow (<1 km) burial. If early carbonate cementation is not pervasive, fluid inclusions in late quartz cement (～>2 km of burial) reveal additional information on fluid flow associated with sand injectites during deeper burial. The latest phase of fluid flow occurs when sand injectites are reactivated as preferential fluid conduits during phases of deformation, when well‐cemented subvertical sand injectites become sites of focussed brittle deformation (fracturing). This study shows that sand injectites are a common and volumetrically important type of structural heterogeneity in sedimentary basins and that long‐lived fluid flow associated with sand injectites in very different settings can be assessed and compared systematically using a combination of petrography and fluid inclusion studies.     

Copper Biosorption by Spent Coffee Ground: Equilibrium,Kinetics, and Mechanism

The potential use of spent coffee ground (SCG) for the removal of copper has been investigated as a low‐cost adsorbent for the biosorption of heavy metals. Adsorption batch experiments were conducted to determine isotherms and kinetics. The biosorption equilibrium data were found to fit well the Freundlich model and an experimental maximum biosorption capacity of copper ions 0.214 mmol/g was achieved. The biosorption kinetics of SCG was studied at different adsorbate concentrations (0.1–1.0 mM) and stirring speeds (100–400/min). The results showed an increase in the copper ion uptake with raising the initial metal concentration and the kinetic data followed the pseudo‐second order rate expression. The effect of stirring speed was a significant factor for the external mass transfer resistance at 100/min and coefficients were estimated by the Mathews and Weber model. Biosorption of copper ions onto SCG was observed to be related mainly with the release of calcium and hydrogen ions suggesting that biosorption performance by SCG can be attributed to ion‐exchange mechanism with calcium and hydrogen ions neutralizing the carboxyl and hydroxyl groups of the biomass.     

Interpreting characteristic drainage timescale variability across Kilombero Valley,Tanzania

We explore seasonal variability and spatiotemporal patterns in characteristic drainage timescale (K) estimated from river discharge records across the Kilombero Valley in central Tanzania. K values were determined using streamflow recession analysis with a Brutsaert–Nieber solution to the linearized Boussinesq equation. Estimated K values were variable, comparing between wet and dry seasons for the relatively small catchments draining upland positions. For the larger catchments draining through valley bottoms, K values were typically longer and more consistent across seasons. Variations in K were compared with long‐term averaged, Moderate‐resolution Imaging Spectroradiometer‐derived monthly evapotranspiration. Although the variations in K were potentially related to evapotranspiration, the influence of data quality and analysis procedure could not be discounted. As such, even though recession analysis offers a potential approach to explore aquifer release timescales and thereby gain insight to a region's hydrology to inform water resources management, care must be taken when interpreting spatiotemporal shifts in K in connection with process representation in regions like the Kilombero Valley. © 2014 The Authors. Hydrological Processes published by John Wiley & Sons Ltd.     

Water contents of clinopyroxenes from sub‐arc mantle peridotites

One poorly constrained reservoir of the Earth's water budget is that of clinopyroxene in metasomatised, mantle peridotites. This study presents reconnaissance Sensitive High‐Resolution, Ion Microprobe–Stable Isotope (SHRIMP–SI) determinations of the H₂O contents of (dominantly) clinopyroxenes in rare mantle xenoliths from four different subduction zones, i.e. Mexico, Kamchatka, Philippines, and New Britain (Tabar‐Feni island chain) as well as one intra‐plate setting (western Victoria). All of the sub‐arc xenoliths have been metasomatised and carry strong arc trace element signatures. Average measured H₂O contents of the pyroxenes range from 70 ppm to 510 ppm whereas calculated bulk H₂O contents range from 88 ppm to 3 737 ppm if the variable presence of amphibole is taken into account. In contrast, the intra‐plate, continental mantle xenolith from western Victoria has higher water contents (3 447 ppm) but was metasomatised by alkali and/or carbonatitic melts and does not carry a subduction‐related signature. Material similar to the sub‐arc peridotites can either be accreted to the base of the lithosphere or potentially be transported by convection deeper into the mantle where it will lose water due to amphibole breakdown.     

Evaluation of the controls affecting the quality of spatial data derived from historical aerial photographs

This paper is concerned with the fundamental controls affecting the quality of data derived from historical aerial photographs typically used in geomorphological studies. A short review is provided of error sources introduced into the photogrammetric workflow. Data‐sets from two case‐studies provided a variety of source data and hence a good opportunity to evaluate the influence of the quality of archival material on the accuracy of coordinated points. Based on the statistical weights assigned to the measurements, precision of the data was estimated a priori, while residuals of independent checkpoints provided an a posteriori measure of data accuracy. Systematic discrepancies between the two values indicated that the routinely used stochastic model was incorrect and overoptimistic. Optimized weighting factors appeared significantly larger than previously used (and accepted) values. A test of repeat measurements explained the large uncertainties associated with the use of natural objects for ground control. This showed that the random errors not only appeared to be much larger than values accepted for appropriately controlled and targeted photogrammetric networks, but also small undetected gross errors were induced through the ‘misidentification’ of points. It is suggested that the effects of such ‘misidentifications’ should be reflected in the stochastic model through selection of more realistic weighting factors of both image and ground measurements. Using the optimized weighting factors, the accuracy of derived data can now be more truly estimated, allowing the suitability of the imagery to be judged before purchase and processing. Copyright © 2010 John Wiley & Sons, Ltd.