Land use change under conditions of high population pressure: the case of Java

A long history of increases in population pressure in Java has caused agricultural land use to expand and intensify. More recent land use changes caused the conversion of prime agricultural land into residential and industrial area. Results of a dynamic, regional-scale, land use change model are presented, defining the spatial distribution of these land use changes. The model is based on multi-scale modelling of the relations between land use and socio-economic and biophysical determinants. Historical validation showed that the model can adequately simulate the pattern of land use change. Future patterns of land use change between 1994 and 2010 are simulated assuming further urbanization. The results suggest that most intensive land use changes will occur in Java's lowland areas.     

Application of Hilbert‐like transforms for enhanced processing of full tensor magnetic gradient data

Commonly, geomagnetic prospection is performed via scalar magnetometers that measure values of the total magnetic intensity. Recent developments of superconducting quantum interference devices have led to their integration in full tensor magnetic gradiometry systems consisting of planar‐type first‐order gradiometers and magnetometers fabricated in thin‐film technology. With these systems measuring directly the magnetic gradient tensor and field vector, a significantly higher magnetic and spatial resolution of the magnetic maps is yield than those produced via conventional magnetometers. In order to preserve the high data quality in this work, we develop a workflow containing all the necessary steps for generating the gradient tensor and field vector quantities from the raw measurement data up to their integration into high­resolution, low­noise, and artefactless two‐dimensional maps of the magnetic field vector. The gradient tensor components are processed by superposition of the balanced gradiometer signals and rotation into an Earth‐centred Earth‐fixed coordinate frame. As the magnetometers have sensitivity lower than that of gradiometers and the total magnetic intensity is not directly recorded, we employ Hilbert‐like transforms, e.g., integration of the gradient tensor components or the conversion of the total magnetic intensity derived by calibrated magnetometer readings to obtain these values. This can lead to a better interpretation of the measured magnetic anomalies of the Earth's magnetic field that is possible from scalar total magnetic intensity measurements. Our conclusions are drawn from the application of these algorithms on a survey acquired in South Africa containing full tensor magnetic gradiometry data.     

Physico-chemical phenomena in comets—III: The continuum of comet Burnham (1960 II)

The new model of the cometary head proposed in papers I and II is developed and applied to comet Burnham. It takes into account the likely existence of a halo of large icy particles surrounding the nucleus. These particles are steadily stripped from the nucleus by evaporating gases. Their terminal velocity and their rate of evaporation set the size of the halo. The existence of the icy halo influences in two ways the photometric characteristics of the coma. This paper establishes the photometric shape of the continuum as reflected by the icy grains, and compares it to the observed continuum of comet Burnham. Paper IV will compare the predictions of the model with the photometric profile of the molecular emission bands of C₂, in the same comet.     

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.     

Occurrence and Distribution of Organic Contaminants in the Aquatic System in Berlin.Part II: Substituted Phenols in Berlin Surface Water

In a surface water screening, 30 representative surface water samples collected from rivers, canals, and lakes in Berlin were investigated for the presence of 22 substituted phenols. The phenols selected include the 11 phenols considered as “priority pollutants” by the US Environmental Protection Agency (US-EPA). Surface water samples were extracted applying solid-phase extraction with styrenedivinylbenzene adsorbent. The recoveries, determined in spiking experiments, were between 80 % and 103 %. After derivatization with N-(tert-butyldimethylsilyl)-N-methyl-trifluoroacetamide (MTBSTFA) the samples were analyzed by capillary gas chromatography-mass spectrometry (GC-MS). Phenol, cresols, 2-ethylphenol, 2-chlorophenol, 4-chloro-3-methylphenol, pentachlorophenol, 2-nitrophenol, and 4-nitrophenol were detected in the surface water samples at concentrations between 0.02 μg/L and 7.8 μg/L, respectively. The distribution of these residues in the Berlin surface waters showed that the phenolic residues, with the exception of pentachlorophenol and 2-ethylphenol, do not originate primarily from municipal sewage treatment plants discharges. Some of the phenols are formed naturally or occur as ubiquitous anthropogenic contaminants in the aquatic system.     

Diffusive equilibrium in thin films provides evidence of suppression of hyporheic exchange and large‐scale nitrate transformation in a groundwater‐fed river

The hyporheic zone of riverbed sediments has the potential to attenuate nitrate from upwelling, polluted groundwater. However, the coarse‐scale (5–10 cm) measurement of nitrogen biogeochemistry in the hyporheic zone can often mask fine‐scale (<1 cm) biogeochemical patterns, especially in near‐surface sediments, leading to incomplete or inaccurate representation of the capacity of the hyporheic zone to transform upwelling NO₃^?. In this study, we utilised diffusive equilibrium in thin‐films samplers to capture high resolution (cm‐scale) vertical concentration profiles of NO₃^?, SO₄^2?, Fe and Mn in the upper 15 cm of armoured and permeable riverbed sediments. The goal was to test whether nitrate attenuation was occurring in a sub‐reach characterised by strong vertical (upwelling) water fluxes. The vertical concentration profiles obtained from diffusive equilibrium in thin‐films samplers indicate considerable cm‐scale variability in NO₃^? (4.4 ± 2.9 mg N/L), SO₄^2? (9.9 ± 3.1 mg/l) and dissolved Fe (1.6 ± 2.1 mg/l) and Mn (0.2 ± 0.2 mg/l). However, the overall trend suggests the absence of substantial net chemical transformations and surface‐subsurface water mixing in the shallow sediments of our sub‐reach under baseflow conditions. The significance of this is that upwelling NO₃^?‐rich groundwater does not appear to be attenuated in the riverbed sediments at <15 cm depth as might occur where hyporheic exchange flows deliver organic matter to the sediments for metabolic processes. It would appear that the chemical patterns observed in the shallow sediments of our sub‐reach are not controlled exclusively by redox processes and/or hyporheic exchange flows. Deeper‐seated groundwater fluxes and hydro‐stratigraphy may be additional important drivers of chemical patterns in the shallow sediments of our study sub‐reach. © 2015 The Authors. Hydrological Processes Published by John Wiley & Sons Ltd.     

Thermodynamic controls on element partitioning between titanomagnetite and andesitic–dacitic silicate melts

Titanomagnetite–melt partitioning of Mg, Mn, Al, Ti, Sc, V, Co, Ni, Cu, Zn, Ga, Zr, Nb, Mo, Hf and Ta was investigated experimentally as a function of oxygen fugacity (fO₂) and temperature (T) in an andesitic–dacitic bulk-chemical compositional range. In these bulk systems, at constant T, there are strong increases in the titanomagnetite–melt partitioning of the divalent cations (Mg²⁺, Mn²⁺, Co²⁺, Ni²⁺, Zn²⁺) and Cu²⁺/Cu⁺ with increasing fO₂ between 0.2 and 3.7 log units above the fayalite–magnetite–quartz buffer. This is attributed to a coupling between magnetite crystallisation and melt composition. Although melt structure has been invoked to explain the patterns of mineral–melt partitioning of divalent cations, a more rigorous justification of magnetite–melt partitioning can be derived from thermodynamic principles, which accounts for much of the supposed influence ascribed to melt structure. The presence of magnetite-rich spinel in equilibrium with melt over a range of fO₂ implies a reciprocal relationship between a(Fe²⁺O) and a(Fe³⁺O_1.5) in the melt. We show that this relationship accounts for the observed dependence of titanomagnetite–melt partitioning of divalent cations with fO₂ in magnetite-rich spinel. As a result of this, titanomagnetite–melt partitioning of divalent cations is indirectly sensitive to changes in fO₂ in silicic, but less so in mafic bulk systems.     

Fractional Crystallization and Liquid Immiscibility Processes in the Alkaline-Carbonatite Complex of Juqui (So Paulo, Brazil)

The Juqui circular intrusion, which is Cretaceous in age (130–135Ma), crops out in the Precambrian gneissic basement in Brazilover an area of 14 km2. It consists of olivine clinopyroxen-itecumulates (with minor olivine gabbros) in the northeastern sector(74 vol.%), whereas ijolites-melteigites-urtites (4%) and nephelinesyenites with minor essexites and syenodiorites (21%) outlinesubannular concentric patterns with an Mg-carbonatite core (1%), in the southwestern part of the complex. Petrographical, bulk rock, and mineral compositional trendsindicate that the origin of the complex can be largely accountedfor by shallow-level fractional crystallization of a carbonatedbasanitic parental magma. Such a magma was generated deep inthe subcontinental lithosphere by low-degree partial meltingof a garnet-phlogopite peridotite source. Mass-balance calculations in agreement with field volume estimatespermit definition of several fractionation stages of the magmaticevolution under nearly closed-system conditions, with inwarddevelopment of zonally arranged side-wall cumulates. These stagesinvolved: (1) fractionation from basanite to essexite magma(liquid fraction F = 33–5%) by crystallization of olivineclinopyroxenite plus minor olivine alkali gabbro cumulates;(2) derivation of the least differentiated mafic nepheline syenite(F = 5–5 %) from essexitic magma by subtraction of a syenodioriteassemblage; (3) exsolution of a carbonatite liquid (5%) froma CO₂-enriched mafic nepheline syenite magma, which also underwentcontinuous fractionation giving rise to ijolite-melteigite-urtitecumulates. The proportion of cumulus clinopyroxene and biotiteand intercumulus nepheline and alkali feldspar in these lastrocks, as well as the absence of alkalis in carbonatite, maybe attributed, at least in part, to loss of alkali-rich hydrousfluids released during and after the unmixing formation of thetwo conjugate liquids. The KD values determined for Mg-carbonatite/nepheline syeniteare lower (1–4–2–9) for light rare earth elements(LREE) than for REE from Eu to Yb (4–6–7–8),in contrast to recent experimental results (Hamilton et al.,1989). A possible explanation is that Juquia Mg-carbonatiterepresents an alreadydifferentiated magma, which underwent extensivefractionation of LREE-enriched calcite. In this way, the highvariability of K₀ REE patterns observed in several alkaline-carbonatitecomplexes can also be accounted for. The remarkably constant initial ⁸⁷Sr/⁸⁶Sr ratios (mostly between0–7052 and 0–7057) support the interpretation ofthe intrusion as having been generated by fractrional crystallizationand liquid immiscibility from a common parental magma. Iligherisotopic ratios (0–7060–0–7078), found mainlyin dykes and in the border facies of the intrusion, may be dueto contamination by the gecissic basement.     

Application of time- and magnitude-predictable model for long-term earthquake prediction in Iran

The regional time- and magnitude-predictable model has been applied successfully in diverse regions of the world to describe the occurrence of main shocks. In the current study, the model has been calibrated against the historical and instrumental catalog of Iranian earthquakes. The Iranian plateau is divided into 15 seismogenic provinces; then, the interevent times for strong main shocks have been determined for each one. The empirical relations reported by Papazachos et al. (Tectonophysics 271:295–323, 1997a) for the Alpine–Himalayan belt (including Iran) were adopted except for the constant terms that were calculated separately for every seismotectonic area. By using the calibrated equations developed for the study area and taking into account the occurrence time and magnitude of the last main shocks in each seismogenic source, the time-dependent conditional probabilities of occurrence P(?t) of the next main shocks during next 10, 20, 30, 40 and 50 years as well as the magnitude of the expected main shocks (M _f) have been estimated. The immediate probability (within next 10 years) of a large main shock is estimated to be high and moderate (>35 %) in all regions except zones 9 (M _f = 5.8) and 15 (M _f = 6.1). However, it should be noted that the probabilities have been estimated for different M _f values in 15 regions. Comparing the model predictions with the actual earthquake occurrence rates shows the good performance of the model for Iranian plateau.