ALBERT

Description: Modern biotechnology, including the application of transgenic techniques to produce Genetically Modified Organisms (GMOs), can play a significant role in increasing agricultural production in a sustainable way, but its products need to be tailored for the developing world. In sub-Saharan Africa, the capacity to develop GMOs and ensure they meet stringent regulatory requirements is somewhat limited. Most African governments contribute little to science and technology either financially or through strong policies. This leaves the determination of research and development priorities in the hands of international funding agencies. Whereas funding from the United States is generally supportive of GM technology, the opposite is true of funding from European sources. African countries are thus pulled in two different directions. One alternative to this dilemma might be for countries in the sub-Saharan Africa region to develop stronger South-South collaborations, but these need to be supported with adequate funding. African governments as well as external funding agencies are urged to consider the important role that biotechnology, including GM technology, can play in contributing to sustainable development in Africa, and to provide adequate support to the development of capacity to research, develop and commercialize GMOs in the region.

Description: SUMMARY In his comment, De Basabe criticises our paper ignoring the advantages of unstructured element mesh used in the finite element method, and argues that the Gaussian quadrature grid (GQG) approach is limited to a homogenous or constant layered geological model and does not have the spectral accuracy. In this reply, we give our response to his criticism and comment, and further clarify the accuracy and capability of the GQG approach.

Description: SUMMARY High-porosity channels are important pathways for melt migration in the mantle. To better understand the lithology and geometry of high-porosity melt channels, we conduct high-order accurate numerical simulations of reactive dissolution along a solubility gradient in an upwelling and viscously deformable porous column. In contrast to earlier studies, we assume the dissolution reaction to be at equilibrium, consider a finite soluble mineral abundance, and employ a high-order accurate numerical scheme. Using sustained perturbations in porosity and soluble mineral abundance at the inflow boundary, we explore the structure of steady-state high-porosity melt channels and their associated lithologies over a range of parameters, including soluble mineral abundance, solubility gradient, amplitude and lateral variation in inflow melt flux, melt fraction and upwelling rate. In general, high-porosity dunite channels are transient and shallow parts of pathways for melt migration in the mantle. The lower parts of a high-porosity channel are orthopyroxene-bearing dunite, harzburgite and possibly lherzolite. A wide orthopyroxene-free dunite channel may contain two or three high-porosity melt channels. The depth of dunite channel initiation depends on the solubility gradient, soluble mineral abundance, inflow melt flux, melt suction rate and upwelling rate. The amplitude and length scale of lateral variation in porosity and orthopyroxene abundance at the base of the upwelling column are important in determining the size and dimension of dunite channels, the strength of melt focusing and the melt suction rate, the presence of compacting boundary layer, as well as the number of high-porosity melt branches within a dunite channel. The spatial relations among the high-porosity melt channels, dunite and harzburgite channels documented in our numerical simulations may shed new light on a number of field, petrological and geochemical observations related to melt migration in the mantle.

Description: SUMMARY Geothermal energy is becoming an important clean energy source, however, the stimulation of a reservoir for an Enhanced Geothermal System (EGS) is associated with seismic risk due to induced seismicity. Seismicity occurring due to the water injection at depth have to be well recorded and monitored. To mitigate the seismic risk of a damaging event, an appropriate alarm system needs to be in place for each individual experiment. In recent experiments, the so-called traffic-light alarm system, based on public response, local magnitude and peak ground velocity, was used. We aim to improve the pre-defined alarm system by introducing a probability-based approach; we retrospectively model the ongoing seismicity in real time with multiple statistical forecast models and then translate the forecast to seismic hazard in terms of probabilities of exceeding a ground motion intensity level. One class of models accounts for the water injection rate, the main parameter that can be controlled by the operators during an experiment. By translating the models into time-varying probabilities of exceeding various intensity levels, we provide tools which are well understood by the decision makers and can be used to determine thresholds non-exceedance during a reservoir stimulation; this, however, remains an entrepreneurial or political decision of the responsible project coordinators. We introduce forecast models based on the data set of an EGS experiment in the city of Basel. Between 2006 December 2 and 8, approximately 11 500 m 3 of water was injected into a 5-km-deep well at high pressures. A six-sensor borehole array, was installed by the company Geothermal Explorers Limited (GEL) at depths between 300 and 2700 m around the well to monitor the induced seismicity. The network recorded approximately 11 200 events during the injection phase, more than 3500 of which were located. With the traffic-light system, actions where implemented after an M L 2.7 event, the water injection was reduced and then stopped after another M L 2.5 event. A few hours later, an earthquake with M L 3.4, felt within the city, occurred, which led to bleed-off of the well. A risk study was later issued with the outcome that the experiment could not be resumed. We analyse the statistical features of the sequence and show that the sequence is well modelled with the Omori–Utsu law following the termination of water injection. Based on this model, the sequence will last 31+29/−14 years to reach the background level. We introduce statistical models based on Reasenberg and Jones and Epidemic Type Aftershock Sequence (ETAS) models, commonly used to model aftershock sequences. We compare and test different model setups to simulate the sequences, varying the number of fixed and free parameters. For one class of the ETAS models, we account for the flow rate at the injection borehole. We test the models against the observed data with standard likelihood tests and find the ETAS model accounting for the on flow rate to perform best. Such a model may in future serve as a valuable tool for designing probabilistic alarm systems for EGS experiments.

Description: SUMMARY In principle, we can never know the true Green's function, which is a major error source in seismic waveform inversion. So far, many studies have devoted their efforts to obtain a Green's function as precise as possible. In this study, we propose a new strategy to cope with this problem. That is to say, we introduce uncertainty of Green's function into waveform inversion analyses. Due to the propagation law of errors, the uncertainty of Green's function results in a data covariance matrix with significant off-diagonal components, which naturally reduce the weight of observed data in later phases. Because the data covariance matrix depends on the model parameters that express slip distribution, the inverse problem to be solved becomes non-linear. Applying the developed inverse method to the teleseismic P -wave data of the 2006 Java, Indonesia, tsunami earthquake, we obtained a reasonable slip-rate distribution and moment-rate function without the non-negative slip constraint. The solution was independent of the initial values of the model parameters. If we neglect the modelling errors due to Green's function as in the conventional formulation, the total slip distribution is much rougher with significant opposite slip components, whereas the moment-rate function is much smoother. If we use a stronger smoothing constraint, more plausible slip distribution can be obtained, but then the moment-rate function becomes even smoother. By comparing the observed waveforms with the synthetic waveforms, we found that high-frequency components were well reproduced only by the new formulation. The modelling errors are essentially important in waveform inversion analyses, although they have been commonly neglected.

Description: The Peruvian anchovy fishery is the largest worldwide in terms of catches. The fishery started during the mid 1950s, and since then it has been highly dependent on natural stock fluctuations, due to the sensitivity of anchovy stocks to ocean-climate variability. The main driver of anchovy stock variability is the El Niño Southern Oscillation (ENSO), and three extreme ENSO warm events were recorded in 1972–1973, 1983–1984 and 1997–1998. This study investigates the evolution of coping strategies developed by the anchovy fisheries to deal with climate variability and extreme ENSO events. Results showed eight coping strategies to reduce impacts on the fishery. These included: decentralized installation of anchovy processing factories; simultaneous ownership of fishing fleet and processing factories; use of low-cost unloading facilities; opportunistic utilization of invading fish populations; low cost intensive monitoring; rapid flexible management; reduction of fishmeal price uncertainty through controlled production based on market demand; and decoupling of fishmeal prices from those of other protein-rich feed substitutes like soybean. This research shows that there are concrete lessons to be learned from successful adaptations to cope with climate change-related extreme climatic events that impact the supply of natural resources. The lessons can contribute to improved policies for coping with climate change in the commercial fishery sector.

Description: SUMMARY Geological and geodetic records show that motion is distributed over a broad and complex system of faults in southern California, with many large earthquakes occurring off the main San Andreas fault. This unusually complex plate boundary geometry is imposed by the big bend of the San Andreas Fault and the dynamics of the underlying lithosphere, but lithosphere rheology and distribution of applied forces are still poorly constrained. Geodynamic modelling is thus very uncertain, and distribution of lithospheric deformations highly controversial. Here, exploiting the property that deformation orients olivine crystals, we show that deep lithospheric deformations can be mapped directly with a new imaging approach of seismic anisotropy, with a resolution fine enough for detailed geodynamic interpretation. This method relies entirely on finite-frequency effects in the splitting of SKS waves. We build an extensive data set of more than 3400 SKS splitting measurements from the analysis of SKS waves recorded by all the broad-band stations in southern California. Resolution tests demonstrate that, using this data set, we are able to resolve anisotropic structures smaller than the size of the first Fresnel zone of SKS waves. The good vertical resolution in our tomographic images results from the short inter-station spacing in southern California. In our 3-D anisotropy model, we find no evidence for localized lithospheric shear deformation beneath the San Andreas Fault. Instead, the lithospheric plate boundary is localized along a broad shear zone beneath the East California Shear Zone. Therefore, surface and deep deformation patterns are poorly correlated and most likely decoupled. Active N–S convergence in the Transverse Ranges results in strongly coherent E–W alignment of olivine fast axes in the shallow lithosphere. At the top of the asthenosphere, the observed low level of anisotropy suggests weak lithospheric basal tractions and thus a strong decoupling between the lithosphere and the asthenosphere.

Description: SUMMARY Earthquakes often occur on faults separating materials with different elastic properties. On theoretical grounds, it is expected that earthquakes on such bimaterial interfaces might have a preferred rupture propagation direction, that being the direction of motion of the more compliant material. The goal of this paper is to determine whether a large sample of natural earthquakes on a bimaterial interface exhibits this tendency. Since the creeping section of San Andreas Fault has a large across-fault velocity contrast and has produced thousands of microearthquakes over the last few decades, the rupture directions of ∼3100 magnitude 0.5–3.0 earthquakes were studied using spectral ratio analysis. The spectral ratios of all earthquake pairs in spatially defined clusters were fitted with synthetic spectral ratios at qualified stations. The synthetics were computed from a simple moving point source model in which each modelled earthquake has four parameters: two rupture lengths (one to the SE and one to the NW) and their propagation velocities. The resolution of rupture directivity increases with event size, such that nearly 900 events, mostly those larger than ∼70 m, appear reasonably well resolved. The inversion results suggest that ∼40 per cent of the well-resolved events are roughly bilateral, although more than ∼80 per cent of the 144 events classified as "strongly unilateral" rupture to the SE, consistent with the theoretical prediction. For those rupture halves that were large enough for the propagation speed to be somewhat resolved, that speed was greater by roughly 10 per cent for those halves propagating to the SE, qualitatively consistent with numerical and laboratory experiments. We also found that events with nearby foreshocks within several hours tend to rupture away from those foreshocks, whether to the NW or to the SE, indicating that asymmetry of prior stressing history can exert a stronger influence on rupture directivity than the material contrast.

Description: Sludge produced by municipal wastewater treatment plants (WWTPs) amounts to only a few percent by volume of the processed wastewater, but its handling accounts for up to 50% of total operating costs. Moreover, the need to achieve a sustainable sludge management strategy has become of great concern. It follows that as conventional and more traditional options, such as land spreading for agricultural purposes, are progressively restricted, and often legally banned, the development of innovative systems to maximize the recovery of useful materials and/or energy is required. A change toward more sustainable procedures can be promoted through an integrated approach, including the assessment of management routes capable of maximizing the recycle/recovery benefits, through low energy impact systems, and the development of operational systems appropriate to local circumstances. Based on the above considerations, an integrated system is proposed in this paper. It includes Anaerobic digestion, Dewatering/Drying, and Pyrolysis/Gasification processes which are efficiently coupled for the recovery of products for material reuse and/or energy purposes. Such an integrated system should also allow the recovery of one or more materials, depending on the combination of processes which best fit specific local situations.

Description: Based on an analysis of the working principles of the hydraulic variable pitch system of a wind turbine, a novel Petri net model and reliability evaluation method are proposed. First, Petri net theory is adopted to build a model for each discrete state of the operation of the hydraulic pitch system of the wind turbine and at the same time a fault Petri net model is established. Then through qualitative analysis and quantitative calculations based on the fault Petri net, the system reliability indexes are obtained. During the qualitative analysis process, in order to more conveniently find the minimal cut sets of the fault Petri net, a Visual C++ 6.0-based algorithm is compiled and the minimal cut sets are tested correctly with another method. During the quantitative calculation process, the fault probability has been obtained from the equations according to the fault probability of libraries and transitions between different states. Not only does the proposed Petri net describe the structure, function and operation of the hydraulic pitch system with a graphic language, but the fault Petri net model can also clearly express the logical relations among faults. The novel Petri net model offers simple calculations and the prospect of broad applicability and the new reliability evaluation method provides an important reference for the performance evaluation of these systems.

Description: SUMMARY We present results on the composition, structure and particle size distribution (PSD) of pulverized and damaged granitic rocks in a 42-m-deep core adjacent to the San Andreas Fault near Littlerock, CA. The cored section is composed of pulverized granites and granodiorites, and is cut by numerous mesoscopic secondary shears. The analysis employs XRD, XRF, thin sections and semi-automated particle size analyser methods, including a novel calibration method. The mean particle size for the majority of samples falls between 50 and 470 μm. The PSDs can be fitted by a power law, with D -values ranging between 2.5 and 3.1, as well as by a superposition of two Gaussians. Fracture surface energy calculations based on the observed particle distributions provide very low values, indicating that the part of the total earthquake energy budget expended for breaking or shattering rocks is small. Shear deformation is likely to dominate near secondary faults. The most pronounced fault-related alteration occurs along the secondary shears, and is a function of both composition and depth. The alteration to clay appears to be the result of fluid–rock interaction and brittle deformation under low temperature conditions, rather than of surface-related weathering. The particle size reduction noted in the core reflects multiple mechanisms of comminution. The zones of pulverization that lack significant weathering likely result from repeating episodes of dynamic dilation and contraction.

Description: SUMMARY Broad-band seismic data from the southern African seismic experiment and the AfricaArray network are used to investigate the seismic velocity structure of the upper mantle beneath southern Africa, and in particular beneath the Kaapvaal Craton. A two-plane approximation method that includes a finite frequency sensitivity kernel is employed to measure Rayleigh wave phase velocities, which are inverted to obtain a quasi-3-D shear wave velocity model of the upper mantle. We find phase velocities for the Kaapvaal Craton and surrounding mobile belts that are comparable to those reported by previous studies, and we find little evidence for variation from east to west across the Namaqua–Natal Belt, a region not well imaged in previous studies. A high-velocity upper-mantle lid is found beneath the Kaapvaal Craton and most of southern Africa. For the Kaapvaal Craton, the thickness of the lid (∼150–200 km) is consistent with the lid thicknesses reported in many previous studies. The cratonic lid is underlain by a ∼100-km thick low-velocity zone with a 3.9 per cent maximum velocity reduction. By comparing the velocity model to those published for other Archean cratons, we find few differences, and therefore conclude that there is little evidence in the shear wave velocity structure of the mantle to indicate that the southern African plateau is supported by an upper-mantle thermal anomaly.

Description: Water, being a primary element in the diet and a necessary resource for the agriculture, can be considered a basic need for humans. In addition, also industrial practices need a growing amount of water. Since human population is continuously growing at a rate that, in the last two centuries, approximates well the exponential, water demand is increasing. However, the water resources on the Earth are finite. For this reason, even disregarding the potential threats due to the climate change, this situation appears as one of the biggest challenges of the current era. Actually, several small-scale regions already face water sustainability problems, and the scarcity of water resources is expected to spread to wider areas in the near future, if the actual trends of development and population growth do not change. The situation is exacerbated as the climate is already changing, due to the anthropogenic emissions of greenhouse gases in the atmosphere, and its rate is expected to increase by the end of this century. The effects of these changes will increase the natural variability of the climate, exacerbating the extreme climatic phenomena (drought and flood events) and increasing the difficulty of managing water resources, especially in the most vulnerable regions.

Description: SUMMARY To better constrain Baltica's position within Pangea, we conducted a palaeomagnetic study of Permo-Triassic dykes from the Oslo Graben, as a follow-up to an initial, but rather limited, study by Torsvik and colleagues in 1998. The age of these so-called Lunner dykes had previously been determined as ∼240 Ma in that study, but details in their analyses and new 40 Ar/ 39 Ar ages reveal that there may have been some argon loss in the initially dated dyke minerals and that a combined (weighted mean) age of 271 ± 2.7 (2σ) Myr for the dykes is preferable. We find two major components of magnetization in our samples: one carried by an Fe-sulphide (likely pyrrhotite) and the other carried by low-Ti magnetite; these magnetization components may be found together (superposed) in a given sample or they may occur apart. Micronmetre-sized crystals of Ti-Fe oxides, observed with a scanning electron microscope (SEM) show exsolution lamellae, formed upon cooling from intrusion temperatures. Assuming that the submicronmetre-sized (Ti)-magnetite grains that carry a stable remanence are of the same generation as the observed larger grains, we interpret the magnetite remanence in the dykes as of primary, thermoremanent origin. The sulphide remanence appears to be slightly younger, as seen by the SEM observations of pyrite framboids and a Fe-sulphide grain invading a Ti-magnetite grain. Moreover, the sulphide mineralization is likely of region-wide hydrothermal origin. The magnetizations carried by the pyrrhotite and magnetite have nearly identical directions and so, must be nearly of the same age. For this study, we sampled 56 sites including 39 dykes, 10 baked-contact rocks and 7 host rocks removed from the immediate dyke contacts. The dykes and the contact rocks have the same SW and up directions of magnetization, and contain the Fe-sulphide or the magnetite magnetization or both, as diagnosed by their relative unblocking temperatures. However, all the sampled carbonate and igneous host rocks far away from the dykes also have the same directions. Thus, all of the 10 originally planned contact tests are inconclusive. The new palaeopoles of this study are a few degrees apart; the magnetite pole (from dykes only, N = 25) is located at 51°N, 164°E, K = 69, A 95 = 3.5°, whereas the pole calculated from iron sulphide magnetic directions (all rock types, N = 20) is at 54°N, 166°E, K = 112, A 95 = 3.1°. All directions are of reversed polarity, suggesting that the magnetization was acquired during the Kiaman Reversed Superchron. The palaeomagnetic mean result from the magnetite-bearing sites implies a palaeolatitude of Oslo of 23°N, whereas the palaeolatitude calculated from the pyrrhotite magnetizations is 25–27°N, depending on choice of host lithologies. As noted in many previous publications, the palaeomagnetic poles for the late Palaeozoic and Early-Middle Triassic are in conflict with classical Pangea reconstructions. The poles with ages of 250 ± 10 Ma, in particular, previously showed a discrepancy of some 25° or more, when the Gondwana and Laurussia continents are restored to their juxtapositions in the Pangea-A fit, before the opening of the Atlantic Ocean. Proposed solutions to this conundrum have been controversial, involving doubts about (1) the geocentric coaxial dipole field model, (2) the reliability of the palaeomagnetic results or their ages, or (3) the validity of the Pangea-A reconstruction, leading to proposals of a Pangea B reconstruction in which Gondwana is displaced some 3500 km to the east with respect to Laurussia. The significance of our new result for this Pangea controversy resides in its improved age within an early Guadelupian (mid to late Permian) time interval where few results exist from well-dated igneous rocks in either Baltica or Laurentia. There are quite a few results from sedimentary rocks, but these may be suspected to suffer inclination shallowing, and are therefore less suitable to settle a palaeolatitudinal argument. Our new result of the magnetite magnetization, granted it is primary and acquired at about 270 Ma, combined with a new ∼265 Ma result from Argentina and selected other poles from igneous rocks, leaves enough room for the north–south configuration of Pangea A at 270 Ma and avoids the overlap between Baltica and Gondwana that necessitated Pangea B, at least for the Late Permian.

Description: This paper provides data from two years of monitoring of the chemical quality of rivers and streams in North West England from the clean headwaters to polluted rivers just above the tidal reach and covers 26 sites including the Ribble, Wyre and the tributary rivers of the Calder and Douglas. Across the basins that include areas of rural, urban and industrial typologies, data is presented for three of the priority substances in the Water Framework Directive i.e., nickel (Ni), cadmium (Cd), and lead (Pb). Average concentrations are low and well below the Environmental Quality Standards values for all three of these substances. Cadmium and Pb appear in approximately equal proportions in the dissolved (0.45 µm) whilst Ni occurs predominantly in the dissolved form (92%). Regional inputs of these metals arise mostly from diffuse sources as the storm-flow concentrations are generally greater than at base-flow condition. Greater concentrations of Ni are transported at the headwaters and smaller tributary sites under storm flow condition than for the main stream of the Ribble. For Ni, amounts increase as the river proceeds from its headwaters down towards the Ribble and Wyre estuaries, whilst Cd and Pb show consistent values throughout the catchment. There is annual cycling of dissolved concentrations of Cd, Pb and Ni for the clean headwater streams that gives maxima during the latter half of the year when the river flow is greater. For the impacted sites the pattern is less distinct or absent. Our estimates suggest that the Ribble estuary receives 550 t y−1 of dissolved Ni, 16 t y−1 of dissolved Cd and 240 t y−1 of dissolved Pb.

Description: The use of remote sensing for environmental policy development is now quite common and well-documented, as images from remote sensing platforms are often used to focus attention on emerging environmental issues and spur debate on potential policy solutions. However, its use in policy implementation and evaluation has not been examined in much detail. Here we examine the use of remote sensing to support the implementation and enforcement of policies regarding the conservation of forests and wetlands in the USA. Specifically, we focus on the “Roadless Rule” and “Travel Management Rules” as enforced by the US Department of Agriculture Forest Service on national forests, and the “No Net Loss” policy and Clean Water Act for wetlands on public and private lands, as enforced by the US Environmental Protection Agency and the US Army Corps of Engineers. We discuss several national and regional examples of how remote sensing for forest and wetland conservation has been effectively integrated with policy decisions, along with barriers to further integration. Some of these barriers are financial and technical (such as the lack of data at scales appropriate to policy enforcement), while others are political.

Description: The aim of this paper is to analyze two Finnish Joint Regional Authorities for Water Supply—namely the Raisio-Naantali Joint Municipal Authority for Water Supply (established in 1957) and the Tuusula Region Joint Municipal Authority for Water Supply (established in 1967)—for assessing the development of supra-municipal water governance. The above two cases make it possible to analyze and assess water policies in settings where the owners are groups of municipalities. The analysis is based on two separately conducted case studies. The study data consist of several types of materials: Annual reports, local government documents, etc. The conducted interviews were semi-structured with some themes defined beforehand. The studies describe two authorities in the context of historical development and as a part of local development.

Description: SUMMARY We investigate 3-D basin structures and site resonance frequencies in the İzmit Bay area of Turkey by new geophysical surveys based on 239 single station microtremor and 405-point gravity measurements. A fundamental resonance frequency map of the İzmit Bay was produced from the main peak in the horizontal-to-vertical component spectral ratio (HVSR) of microtremors. The HVSR analysis of the microtremor data reveals single, double, broad peaked or no peak type HVSR curves varying in accordance with the surface geology and spatial extent of the three basins present in the İzmit Bay area. In the deepest sections of the İzmit, Gölcük and Derince basins of the İzmit Bay, the fundamental resonance frequencies are dominantly 0.24–0.30 Hz. These resonance frequencies should be taken into consideration along with higher mode frequencies to construct earthquake resistant structures in the İzmit Bay area. The 3-D gravimetric bedrock depth map of the İzmit basin shows that the basin has an asymmetric shape with its deepest section coinciding with the surface trace of the North Anatolian Fault. The northern shoulder of the basin has a gentle dip on the Kocaeli Peneplain side and the southern shoulder is much steeper, and it is bounded by the Samanlı Mountains. We derive a power-law relationship that correlates the fundamental site resonance frequencies with the sedimentary cover thickness obtained from the gravity and shear wave velocity data in the İzmit Bay. We use this relationship to estimate bedrock depths beneath Gölcük and Derince basins. Our estimation of maximum basin depths is 1400 m for the İzmit and is 800 m for the Gölcük and Derince basins. Finally, we have analysed a converted Sp phase from a local earthquake recording made at site CMP to calculate and verify the sediment thickness estimations obtained from our gravimetric and microtremor analyses. This calculation shows close agreement with that of the gravimetric and microtremor results. Our results show that the basins in the İzmit Bay area have a very thick sedimentary cover with very low shear velocities underlined by hard bedrock, forming a sharp impedance contrast. We anticipate that these results will be a key contribution to the quantitative assessment of seismic hazard for the İzmit Bay area before the occurrence of strong earthquakes in the Marmara region.

Description: Unconferencing is a method for organizing social learning which could be suitable to trigger sustainability learning processes. An unconference is defined as participant-driven meeting that tries to avoid one or more aspects of a conventional conference, such as top-down organization, one-way communication and power-relationships based on titles, formal hierarchies and status. This paper presents a case study on the application of unconferencing in a large Swiss university (ETH Zurich) where an unconference was conducted to engage students, academics, staff and external experts in a mutual learning process aimed at the development of project ideas for reducing its CO2 emissions. The study analyzes how the unconferencing format initiated and promoted sustainability oriented group processes during the unconference, and in how far the projects which were developed contributed to a reduction of the university’s CO2 emissions.

Description: SUMMARY To constrain the forecast horizon of geomagnetic data assimilation, it is of interest to quantify the range of predictability of the geodynamo. Following earlier work in the field of dynamic meteorology, we investigate the sensitivity of numerical dynamos to various perturbations applied to the magnetic, velocity and temperature fields. These perturbations result in some errors, which affect all fields in the same relative way, and grow at the same exponential rate  λ=τ −1 e , independent of the type and the amplitude of perturbation. Errors produced by the limited resolution of numerical dynamos are also shown to produce a similar amplification, with the same exponential rate. Exploring various possible scaling laws, we demonstrate that the growth rate is mainly proportional to an advection timescale. To better understand the mechanism responsible for the error amplification, we next compare these growth rates with two other dynamo outputs which display a similar dependence on advection: the inverse  τ −1 SV  of the secular-variation timescale, characterizing the secular variation of the observable field produced by these dynamos; and the inverse  (τ mag diss ) −1  of the magnetic dissipation time, characterizing the rate at which magnetic energy is produced to compensate for Ohmic dissipation in these dynamos. The possible role of viscous dissipation is also discussed via the inverse  (τ kin diss ) −1  of the analogous viscous dissipation time, characterizing the rate at which kinetic energy is produced to compensate for viscous dissipation. We conclude that  τ e  tends to equate  τ mag diss  for dynamos operating in a turbulent regime with low enough Ekman number, and such that  τ mag diss 〈 τ kin diss . As these conditions are met in the Earth's outer core, we suggest that  τ e  is controlled by magnetic dissipation, leading to a value  τ e =τ mag diss ≈ 30  yr . We finally discuss the consequences of our results for the practical limit of predictability of the geodynamo.

Description: Understanding the socio-economic factors that are associated with fishers’ willingness to delay gratification may be useful for designing appropriate fisheries management and conservation policies. We aim to identify the predictors of low discounting behaviour among fishers, which is analogous to having a longer-term outlook. We base our empirical study on two small-scale tropical reef fisheries in Sabah, Malaysia, and Fiji. We use an experimental approach to identify fishers with low discount rates, and then use a logistic regression model to identify predictors of low discount rates. We find that 42% of the respondents have low discount rates, and that site and village level variables are significant predictors of low discount rates. Within Sabah and Fiji, boat ownership and relative catch differentiate low discounting from non-low discounting fishers, but these variables have contradictory effects in Sabah and Fiji. Overall, our results imply that a substantial proportion of reef fishers may be willing to engage in conservation initiatives; however, local socio-cultural, economic, and ecological conditions have to be considered first during the process of designing management interventions.

Description: SUMMARY A closed-form solution is given for a 2-D, transverse electric mode, magnetotelluric (MT) problem. The model system consists of a finite vertical thin conductor with variable integrated conductivity over a perfectly conducting base. A notable property of the solution is that the frequency response possesses a single pole in the complex plane. Systems with finitely many resonances play a central role in the 1-D MT inverse problem based on finite data sets, but until now, no 2-D system of this kind was known. The particular model is shown to be just one of a large class of thin conductors with same the property, and further examples are given. The solutions of the induction problem for members of this family can often be written in compact closed form, making them the simplest known solutions to the 2-D MT problem.

Description: Peelle’s Pertinent Puzzle (PPP) was described in 1987 in the context of estimating fundamental parameters that arise in nuclear interaction experiments. In PPP, generalized least squares (GLS) parameter estimates fell outside the range of the data, which has raised concerns that GLS is somehow flawed and has led to suggested alternatives to GLS estimators. However, there have been no corresponding performance comparisons among methods, and one suggested approach involving simulated data realizations is statistically incomplete. Here we provide performance comparisons among estimators, introduce approximate Bayesian computation (ABC) using density estimation applied to simulated data realizations to produce an alternative to the incomplete approach, complete the incompletely specified approach, and show that estimation error in the assumed covariance matrix cannot always be ignored.

Description: This study used an affordable ground-based portable LiDAR system to provide an understanding of the structural differences between old-growth and secondary-growth Southeastern pine. It provided insight into the strengths and weaknesses in the structural determination of portable systems in contrast to airborne LiDAR systems. Portable LiDAR height profiles and derived metrics and indices (e.g., canopy cover, canopy height) were compared among plots with different fire frequency and fire season treatments within secondary forest and old growth plots. The treatments consisted of transitional season fire with four different return intervals: 1-yr, 2-yr, 3-yr fire return intervals, and fire suppressed plots. The remaining secondary plots were treated using a 2-yr late dormant season fire cycle. The old growth plots were treated using a 2-yr growing season fire cycle. Airborne and portable LiDAR derived canopy cover were consistent throughout the plots, with significantly higher canopy cover values found in 3-yr and fire suppressed plots. Portable LiDAR height profile and metrics presented a higher sensitivity in capturing subcanopy elements than the airborne system, particularly in dense canopy plots. The 3-dimensional structures of the secondary plots with varying fire return intervals were dramatically different to old-growth plots, where a symmetrical distribution with clear recruitment was visible. Portable LiDAR, even though limited to finer spatial scales and specific biases, is a low-cost investment with clear value for the management of forest canopy structure.

Description: SUMMARY The eikonal equation in an attenuating medium has the form of a complex-valued Hamilton–Jacobi equation and must be solved in terms of the complex-valued traveltime (complex-valued action function). A very suitable approximate method for calculating the complex-valued traveltime right in real space is represented by the perturbation from the reference traveltime calculated along real-valued reference rays to the complex-valued traveltime defined by the complex-valued Hamilton–Jacobi equation. The real-valued reference rays are calculated using the reference Hamiltonian function. The perturbation Hamiltonian function is parametrized by one or more perturbation parameters, and smoothly connects the reference Hamiltonian function with the Hamiltonian function corresponding to a given complex-valued Hamilton–Jacobi equation. Both the reference Hamiltonian function and the perturbation Hamiltonian function may be constructed in different ways, yielding differently accurate perturbation expansions of traveltime. All present perturbation methods use reference rays calculated in a reference anisotropic non-attenuating medium. In this paper, the reference Hamiltonian function is constructed directly using the Hamiltonian function corresponding to a given complex-valued Hamilton–Jacobi equation, and the perturbation Hamiltonian function is linear with respect to the perturbation parameter. The direct construction of the reference Hamiltonian function from the given complex-valued Hamilton–Jacobi equation is very general and accurate, especially for homogeneous Hamiltonian functions of degree N =−1 with respect to the slowness vector.

Description: SUMMARY The lithospheric structure in Oman has been determined by analyzing teleseismic P -receiver functions recorded at broad-band and short-period seismic stations of the Oman Seismological Network. Lithospheric structure is obtained by jointly inverting receiver functions and Rayleigh wave group velocities derived from continental-scale surface wave tomography of Pasyanos. We observe relatively thick crust (40–48 km) within the ophiolite formed mountains in northern Oman. The crustal thickness is about 35 km within the passive continental margin of the southern Oman region. Uppermost (〈5 km) crustal shear wave velocities are faster in the northern ophiolite regions compared to the southern Oman region, while shear velocities in the middle crust are faster in the Southern Oman region compared to the ophiolite region. This observation coincides well with cretaceous to Eocene marine platform sequences overlying Precambrian to Cambrian basement of the southern part. Joint inversion analysis shows that the Moho depth of Oman varies from 34 km in the southern region to 48 km in the northern part.

Description: This study integrates photovoltaic (PV) system, building structure, and heat flow mechanism to propose the notion of ventilated Building-Integrated Photovoltaic (BIPV) walls. The energy-saving potential of the ventilated BIPV walls was investigated via engineering considerations and computational fluid dynamics (CFD) simulations. The results show that the heat removal rate and indoor heat gain of the proposed ventilated BIPV walls were dominantly affected by outdoor wind velocity and airflow channel width. Correlations for predicting the heat removal rate and indoor heat gain, the reduction ratio of the indoor heat gain, CO2 reduction, and induced indoor air exchange are introduced.

Description: The goal of the two year Ground European Network for Earth Science Interoperations-Digital Repositories (GENESI-DR) project was to build an open and seamless access service to Earth science digital repositories for European and world-wide science users. In order to showcase GENESI-DR, one of the developed technology demonstrators focused on fast search, discovery, and access to remotely sensed imagery in the context of post-disaster building damage assessment. This paper describes the scenario and implementation details of the technology demonstrator, which was developed to support post-disaster damage assessment analyst activities. Once a disaster alert has been issued, response time is critical to providing relevant damage information to analysts and/or stakeholders. The presented technology demonstrator validates the GENESI-DR project data search, discovery and security infrastructure and integrates the rapid urban area mapping and the near real-time orthorectification web processing services to support a post-disaster damage needs assessment analysis scenario. It also demonstrates how the GENESI-DR SOA can be linked to web processing services that access grid computing resources for fast image processing and use secure communication to ensure confidentiality of information.

Description: In this paper, two main approaches for automatic building detection and localization using high spatial resolution imagery and LiDAR data are compared and evaluated: thresholding-based and object-based classification. The thresholding-based approach is founded on the establishment of two threshold values: one refers to the minimum height to be considered as building, defined using the LiDAR data, and the other refers to the presence of vegetation, which is defined according to the spectral response. The other approach follows the standard scheme of object-based image classification: segmentation, feature extraction and selection, and classification, here performed using decision trees. In addition, the effect of the inclusion in the building detection process of contextual relations with the shadows is evaluated. Quality assessment is performed at two different levels: area and object. Area-level evaluates the building delineation performance, whereas object-level assesses the accuracy in the spatial location of individual buildings. The results obtained show a high efficiency of the evaluated methods for building detection techniques, in particular the thresholding-based approach, when the parameters are properly adjusted and adapted to the type of urban landscape considered.

Description: A review of the literature published from 2008 to 2010 on topics related to chemicals and allied products is presented. The review considered several sections such as waste management, physicochemical treatment, aerobic treatment, anaerobic treatment, air emissions, soils and groundwater, and reuse.

Description: Accurate electric load forecasting has become the most important issue in energy management; however, electric load demonstrates a seasonal/cyclic tendency from economic activities or the cyclic nature of climate. The applications of the support vector regression (SVR) model to deal with seasonal/cyclic electric load forecasting have not been widely explored. The purpose of this paper is to present a SVR model which combines the seasonal adjustment mechanism and a chaotic immune algorithm (namely SSVRCIA) to forecast monthly electric loads. Based on the operation procedure of the immune algorithm (IA), if the population diversity of an initial population cannot be maintained under selective pressure, then IA could only seek for the solutions in the narrow space and the solution is far from the global optimum (premature convergence). The proposed chaotic immune algorithm (CIA) based on the chaos optimization algorithm and IA, which diversifies the initial definition domain in stochastic optimization procedures, is used to overcome the premature local optimum issue in determining three parameters of a SVR model. A numerical example from an existing reference is used to elucidate the forecasting performance of the proposed SSVRCIA model. The forecasting results indicate that the proposed model yields more accurate forecasting results than the ARIMA and TF-ε-SVR-SA models, and therefore the SSVRCIA model is a promising alternative for electric load forecasting.

Description: We present experimental data on the partitioning of Li, Be, B, K, Mg, Sr, Ga, Rb, Cs, Ba, La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu, Th, U, Hf, Zr, Nb and Ta between lawsonite and fluid, and zoisite and fluid at 3·0–3·5 GPa and 650–850°C. The aim is to provide data bearing on the trace element contents of fluids released during dehydration of subducting oceanic crust. Experimental trace element partition coefficients for lawsonite indicate a preference for the light rare earth elements (LREE) over the heavy REE (HREE) and for Be. These characteristics are consistent with the chemical composition of lawsonite in natural rocks. Experimental trace element partition coefficients for zoisite indicate a preference for HREE relative to LREE. This observation, consistent with earlier experimental data, is the reverse of the observed trace element compositions of natural zoisites, indicating the influence of other factors on the trace element contents of this phase. Lattice strain theory explains well the experimentally derived partitioning of divalent cations in the Ca-site between lawsonite and fluid. However, the weak relative fractionation of REE between lawsonite and fluid cannot be explained by lattice strain theory, as previously observed for zoisite–fluid REE partitioning. We combine our experimental data with thermodynamic models of mineral stability to model the compositions of fluids released during subduction of altered normal mid-ocean ridge basalt. The low La/Sm ratio associated with very high Ba/Th in arc magmas can be explained only if allanite is stable in the subducting oceanic crust. This suggests that the crustal fluid component involved in arc magma petrogenesis results from processes occurring in the warm, top part of the subducting slab. Decreasing lawsonite modal proportion with depth is associated with a large release of fluid characterized by low B/Be ratios that could explain the decreasing B/Be ratios in arc magmas with increasing distance from the trench. This implies that an important Be input in arc magma originates from the fluid generated during oceanic crust dehydration.

Description: To contribute to our understanding of the mechanisms and pathways of fluid movement through deeply subducted crust, we investigate high-pressure veins cutting eclogite-facies (~2·0 GPa and ~600°C) metagabbros of the Monviso Ophiolite, Italian Western Alps. The veins consist mainly of omphacite with minor garnet, rutile, talc and accessory zircon. Most of the vein minerals have major and trace element compositions that are comparable with the host-rock minerals, and vein and host-rock zircons have similar Hf isotopic compositions. These observations support the conclusions of previous studies that these veins largely formed from a locally sourced hydrous fluid during prograde or peak metamorphism. However, the bulk-rock Cr and Ni contents of the veins are significantly higher than those of the surrounding host eclogites. We also document distinct Cr-rich (up to weight per cent levels) zones in omphacite, garnet and rutile in some vein samples. Vein garnet and talc also have relatively high MgO and Ni contents. X-ray maps of vein garnet and rutile grains reveal complex internal zoning features, which are largely defined by micrometre-scale variations in Cr content. Some grains have concentric and oscillatory zoning in Cr, whereas others feature a chaotic fracture-like pattern. These Cr-rich zones are associated with high concentrations of Ni, B, As, Sb, Nb, Zr and high ratios of light rare earth elements (LREE) to middle REE (MREE) compared with low-Cr vein and host-rock minerals. Petrological and mass-balance constraints verify that the Cr-rich zones in the veins were not derived from internally sourced fluids, but represent precipitates from an external fluid. The external source that is consistent with the distinctive trace element characteristics of the vein components is antigorite serpentinite, which forms the structural basement of the high-pressure metagabbros. We propose at least two separate growth mechanisms for the Monviso veins. Most vein infillings were formed during progressive prograde metamorphism from locally derived fluid. Influx of the serpentinite-derived or other external fluid was transient and episodic and was probably achieved via brittle fractures, which preferentially formed along the pre-existing vein structures. The dehydration of serpentinite at high pressures in subduction zones may provide crucial volatiles and trace elements for arc magmas. Our results indicate that the movement of these fluids through subducted oceanic crust is likely to be highly channeled and transient so the progressive development of vein systems in mafic rocks may also be crucial for forming channelways for long-distance fluid flow at depth in subduction zones.

Description: Siwi caldera, in the Vanuatu arc (Tanna island), is a rare volcanic complex where both persistent eruptive activity (Yasur volcano) and rapid block resurgence (Yenkahe horst) can be investigated simultaneously during a post-caldera stage. Here we provide new constraints on the feeding system of this volcanic complex, based on a detailed study of the petrology, geochemistry and volatile content of Yasur–Siwi bulk-rocks and melt inclusions, combined with measurements of the chemical composition and mass fluxes of Yasur volcanic gases. Major and trace element analyses of Yasur–Siwi volcanic rocks, together with literature data for other volcanic centers, point to a single magmatic series and possibly long-lived feeding of Tanna volcanism by a homogeneous arc basalt. Olivine-hosted melt inclusions show that the parental basaltic magma, which produces basaltic-trachyandesites to trachyandesites by ~50–70% crystal fractionation, is moderately enriched in volatiles (~1 wt % H 2 O, 0·1 wt % S and 0·055 wt % Cl). The basaltic-trachyandesite magma, emplaced at between 4–5 km depth and the surface, preserves a high temperature (1107 ± 15°C) and constant H 2 O content (~1 wt %) until very shallow depths, where it degasses extensively and crystallizes. These conditions, maintained over the past 1400 years of Yasur activity, require early water loss during basalt differentiation, prevalent open-system degassing, and a relatively high heat flow (~10 9 W). Yasur volcano releases on average ≥ 13·4 x 10 3 tons d –1 of H 2 O and 680 tons d –1 of SO 2 , but moderate amounts of CO 2 (840 tons d –1 ), HCl (165 tons d –1 ), and HF (23 tons d –1 ). Combined with melt inclusion data, these gas outputs constrain a bulk magma degassing rate of ~5 x 10 7 m 3 a –1 , about a half of which is due to degassing of the basaltic-trachyandesite. We compute that 25 km 3 of this magma have degassed without erupting and have accumulated beneath Siwi caldera over the past 1000 years, which is one order of magnitude larger than the accumulated volume uplift of the Yenkahe resurgent block. Hence, basalt supply and gradual storage of unerupted degassed basaltic-trachyandesite could easily account for (or contribute to) the Yenkahe block resurgence.

Description: Primitive basalts are rarely found in arcs. The active NW Rota-1 volcano in the Mariana arc has erupted near-primitive lavas, which we have sampled with ROV Hyper-Dolphin (HPD). Samples from the summit (HPD480) and eastern flank (HPD488) include 17 magnesian basalts (51–52 wt % SiO 2 ) with 7·5–9·5 wt % MgO and Mg-number of 61–67, indicating little fractionation. Olivine phenocrysts are as magnesian as Fo 93 and contain 0·4 wt % NiO; the Cr/(Cr + Al) values of spinels are mostly 0·5–0·8, indicating equilibrium with depleted mantle. There are three petrographic groups, based on phenocryst populations: (1) cpx–olivine basalt (COB); (2) plagioclase–olivine basalt (POB); (3) porphyritic basalt. Zr/Y and Nb/Yb are higher in POB (3·1–3·2 and 1·2–1·5, respectively) than in COB (Zr/Y = 2·8–3·0 and Nb/Yb = 0·7–0·9), suggesting that POB formed from lower degrees of mantle melting, or that the COB mantle source was more depleted. On the other hand, COB have Ba/Nb (70–80) and Th/Nb (0·4–0·5) that are higher than for POB (Ba/Nb = 30–35 and Th/Nb = 0·1–0·2), and also have steeper light rare earth element (LREE)-enriched patterns. Moreover, COB have enriched 87 Sr/ 86 Sr and 143 Nd/ 144 Nd, and higher Pb isotope values, suggesting that COB has a greater subduction component than POB. 176 Hf/ 177 Hf between COB and POB are similar and Hf behavior in COB and POB is similar to that of Zr, Y and HREE, suggesting that Hf is not included in the subduction component, which produced the differences between COB and POB. The calculated primary basaltic magmas of NW Rota-1 volcano (primary COB and POB magmas) indicate segregation pressures of 2–1·5 GPa (equivalent to 65–50 km depth). These magmas formed by 24–18% melting of mantle peridotite having Mg-number ~89·5. Diapiric ascent of hydrous peridotite mixed heterogeneously with sediment melts may be responsible for the NW Rota-1 basalts. These two basalt magma types are similar to those found at Sumisu and Torishima volcanoes in the Izu–Bonin arc, with COB representing wetter and POB representing drier magmas, where subduction zone-derived melt components are coupled with the water contents.

Description: The origin of andesite is an important issue in petrology because andesite is the main eruptive product at convergent margins, corresponds to the average crustal composition and is often associated with major Cu–Au mineralization. In this study we present petrographic, mineralogical, geochemical and isotopic data for basaltic andesites of the latest Pleistocene Pilavo volcano, one of the most frontal volcanoes of the Ecuadorian Quaternary arc, situated upon thick (30–50 km) mafic crust composed of accreted Cretaceous oceanic plateau rocks and overlying mafic to intermediate Late Cretaceous–Late Tertiary magmatic arcs. The Pilavo rocks are basaltic andesites (54–57·5 wt % SiO 2 ) with a tholeiitic affinity as opposed to the typical calc-alkaline high-silica andesites and dacites (SiO 2 59–66 wt %) of other frontal arc volcanoes of Ecuador (e.g. Pichincha, Pululahua). They have much higher incompatible element contents (e.g. Sr 650–1350 ppm, Ba 650–1800 ppm, Zr 100–225 ppm, Th 5–25 ppm, La 15–65 ppm) and Th/La ratios (0·28–0·36) than Pichincha and Pululahua, and more primitive Sr ( 87 Sr/ 86 Sr ~0·7038–0·7039) and Nd ( Nd ~ +5·5 to +6·1) isotopic signatures. Pilavo andesites have geochemical affinities with modern and recent high-MgO andesites (e.g. low-silica adakites, Setouchi sanukites) and, especially, with Archean sanukitoids, for both of which incompatible element enrichments are believed to result from interactions of slab melts with peridotitic mantle. Petrographic, mineral chemistry, bulk-rock geochemical and isotopic data indicate that the Pilavo magmatic rocks have evolved through three main stages: (1) generation of a basaltic magma in the mantle wedge region by flux melting induced by slab-derived fluids (aqueous, supercritical or melts); (2) high-pressure differentiation of the basaltic melt (at the mantle–crust boundary or at lower crustal levels) through sustained fractionation of olivine and clinopyroxene, leading to hydrous, high-alumina basaltic andesite melts with a tholeiitic affinity, enriched in incompatible elements and strongly impoverished in Ni and Cr; (3) establishment of one or more mid-crustal magma storage reservoirs in which the magmas evolved through dominant amphibole and clinopyroxene (but no plagioclase) fractionation accompanied by assimilation of the modified plutonic roots of the arc and recharge by incoming batches of more primitive magma from depth. The latter process has resulted in strongly increasing incompatible element concentrations in the Pilavo basaltic andesites, coupled with slightly increasing crustal isotopic signatures and a shift towards a more calc-alkaline affinity. Our data show that, although ultimately originating from the slab, incompatible element abundances in arc andesites with primitive isotopic signatures can be significantly enhanced by intra-crustal processes within a thick juvenile mafic crust, thus providing an additional process for the generation of enriched andesites.

Description: To meet European Water Framework Directive requirements, data are needed on faecal indicator organism (FIO) concentrations in rivers to enable the more heavily polluted to be targeted for remedial action. Due to the paucity of FIO data for the UK, especially under high-flow hydrograph event conditions, there is an urgent need by the policy community for generic models that can accurately predict FIO concentrations, thus informing integrated catchment management programmes. This paper reports the development of regression models to predict base- and high-flow faecal coliform (FC) and enterococci (EN) concentrations for 153 monitoring points across 14 UK catchments, using land cover, population (human and livestock density) and other variables that may affect FIO source strength, transport and die-off. Statistically significant models were developed for both FC and EN, with greater explained variance achieved in the high-flow models. Both land cover and, in particular, population variables are significant predictors of FIO concentrations, with r2 maxima for EN of 0.571 and 0.624, respectively. It is argued that the resulting models can be applied, with confidence, to other UK catchments, both to predict FIO concentrations in unmonitored watercourses and evaluate the likely impact of different land use/stocking level and human population change scenarios.

Description: SUMMARY We analyse the spatial distribution of the intensity data points surveyed after the M w 6.3, 2009 L’Aquila (central Italy) earthquake, with the aim to recognize and quantify finite-fault and directivity effects. The study is based on the analysis of the residuals, evaluated with respect to attenuation-with-distance models, calibrated for L’Aquila earthquake. We apply a non-parametric approach considering both the epicentral and the rupture distance, which accounts for the finite extension of the source. Then, starting from a simplified kinematic rupture model of the L’Aquila fault, we compute four directivity predictors proposed in literature, and assess their correlation with intensity residuals. We derive a so-called Intensity Directivity Factor by the correlation between theoretical predictors and observed residuals that allows us to identify and quantify the intensity data points affected by forward and backward directivity during L’Aquila earthquake. We find that the effects are more pronounced in the forward directivity direction and increments up to 1 MCS intensity unit are expected. Moreover, the directivity predictor that accounts for radiation pattern poorly correlates with residuals. These results show that the spatial distribution of the L’Aquila macroseismic field is affected by source effects and in particular that directivity-induced amplification effects can be recognized. We show that the quasi-unilateral rupture propagation along the fault can explain the high-intensity patterns observed along specific direction at relatively large distance from the instrumental epicentre, in accordance with the seismological source models derived from the analysis of instrumental observations.

Description: SUMMARY Temporal variations in the elastic behaviour of the Earth’s crust can be monitored through the analysis of the Earth’s seismic response and its evolution with time. This kind of analysis is particularly interesting when combined with the reconstruction of seismic Green’s functions from the cross-correlation of ambient seismic noise, which circumvents the limitations imposed by a dependence on the occurrence of seismic events. In fact, because seismic noise is recorded continuously and does not depend on earthquake sources, these cross-correlation functions can be considered analogously to records from continuously repeating doublet sources placed at each station, and can be used to extract observations of variations in seismic velocities. These variations, however, are typically very small: of the order of 0.1 per cent. Such accuracy can be only achieved through the analysis of the full reconstructed waveforms, including later scattered arrivals. We focus on the method known as Moving-Window Cross-Spectral Analysis that has the advantage of operating in the frequency domain, where the bandwidth of coherent signal in the correlation function can be clearly defined. We investigate the sensitivity of this method by applying it to microseismic noise cross-correlations which have been perturbed by small synthetic velocity variations and which have been randomly contaminated. We propose threshold signal-to-noise ratios above which these perturbations can be reliably observed. Such values are a proxy for cross-correlation convergence, and so can be used as a guideline when determining the length of microseismic noise records that are required before they can be used for monitoring with the moving-window cross-spectral technique.

Description: SUMMARY The young (〈 5 Ma) Corinth Rift is an ideal natural laboratory to investigate rift deformation mechanisms by comparing extension rates determined by various methods spanning different depth and time ranges. Corinth Rift geodetic extension rates averaged over 5–100 yr have been interpreted to increase from ∼5 mm yr –1 or less in the east to 〉10–15 mm yr –1 in the west. We quantify total upper-crust and whole-crust extension on three profiles across the Corinth Rift. Whole-crust extension is greater across the central rift (∼11–21 km) than across the western part of the rift (∼5–13 km). This correlates with the overall rift morphology, which shows maximum basement subsidence, sediment accumulation, rift width and greatest summed Late Quaternary fault displacements in the central basin, but contrasts with the pattern of geodetic extension rates which are greater to the west of the central basin. The E–W increase in strain rates interpreted from geodetic data cannot have persisted over rift history to produce the observed rift morphology. We suggest the discrepancy between short-term and long-term extension patterns is related to shifts in the loci of maximum extension due to fault growth and linkage during Corinth Rift history, and is likely a characteristic of rift development in general. Total upper-crust and whole-crust extension estimates in the western rift, where extension estimates are best constrained, are within error. We propose that uniform pure-shear extension is a viable extension mechanism in the western rift and crustal extension estimates do not require the existence of a major active N–S dipping detachment fault.

Description: A comparison of the community structure of juvenile hermatypic corals of 2 to 37 m depth at the fringing reefs of Curaçao between 1975 and 2005 shows a decline of 54.7% in juvenile coral abundance and a shift in species composition. Agaricia species and Helioseris cucullata, the most common juveniles in 1975, showed the largest decline in juvenile abundance (a 9 and 120 fold decrease in density respectively) with Helioseris cucullata being nearly extirpated locally. In 2005, Porites astreoides contributed most colonies to the juvenile coral community, increasing from 8.2% (in 1975) to 19.9% of the total juvenile community. Between 1975 and 2005, juveniles of brooding species decreased in relative abundance while the abundance of juveniles of broadcast spawning species increased or remained the same. These data illustrate the magnitude of the changes that have occurred in only three decades in the composition of juvenile coral communities.

Description: SUMMARY This paper describes an innovative method of volcano deformation measurements, applied to camera images taken from the 2004–2008 eruption period at Mount St. Helens. Dome growth was thought to be characterized by sustained, near-linear rates of a solid dacite plug. Through spatial digital image correlation (DIC) analysis of the camera images, new evidences arise that the deformation and strain rate of the spine was more complex. DIC yielded cumulative and incremental displacements, strain and shear planes at decimetre resolution. It was found that dome extrusion rates are highly non-linear, decelerating prior to partial collapse, followed by a pronounced dome extrusion increase and direction change. Associated processes have been identified through DIC, such as shallow landslides and reworking of talus apron material. The work highlights the strengths of camera strain monitoring, and illustrates that dome growth and collapse is a very dynamic process complexly interplaying with the surrounding.

Description: SUMMARY Direct detection of Rayleigh-wave azimuthal anisotropy is reported by applying an array analysis to broad-band seismic data in Southern California, USA. Our approach has excellent resolution for frequencies between 30 and 60 mHz and good resolution between 10 and 30 mHz. Limitation from array size limits accuracy below 10 mHz and complicated wave propagation effects lead to difficulty above 60 mHz. Between 30 and 60 mHz, azimuthal anisotropy of Rayleigh-wave phase velocity is detected cleanly. Phase velocity results show that the 2θ components dominate the 4θ components. Strength of anisotropy hovers around 1 per cent for the 2θ components with standard error ( 1 σ ) of about 0.4 per cent. Fast axes orientations show stable orientation in the direction 110° clockwise from north and its opposite direction 290° . Depth sensitivity kernels for P -wave and S -wave anisotropy indicate that anisotropy is confined to depths in the upper 100 km and may even be to shallower depths. We speculate that the fast axes alignment may be associated with a strong shearing process over the thickness of lithosphere due to obliquely collisional motions between the Pacific Plate and the North America Plate in the region.

Description: This paper explores the politics around the role of agency in the UK climate change debate. Government interventions on the demand side of consumption have increasingly involved attempts to obtain greater traction with the values, attitudes and beliefs of citizens in relation to climate change and also in terms of influencing consumer behaviour at an individual level. With figures showing that approximately 40% of the UK’s carbon emissions are attributable to household and transport behaviour, policy initiatives have progressively focused on the facilitation of “sustainable behaviours”. Evidence suggests however, that mobilisation of pro-environmental attitudes in addressing the perceived “value-action gap” has so far had limited success. Research in this field suggests that there is a more significant and nuanced “gap” between context and behaviour; a relationship that perhaps provides a more adroit reflection of reasons why people do not necessarily react in the way that policy-makers anticipate. Tracing the development of the UK Government’s behaviour change agenda over the last decade, we posit that a core reason for the limitations of this programme relates to an excessively narrow focus on the individual. This has served to obscure some of the wider political and economic aspects of the debate in favour of a more simplified discussion. The second part of the paper reports findings from a series of focus groups exploring some of the wider political views that people hold around household energy habits, purchase and use of domestic appliances, and transport behaviour-and discusses these insights in relation to the literature on the agenda’s apparent limitations. The paper concludes by considering whether the aims of the Big Society approach (recently established by the UK’s Coalition Government) hold the potential to engage more directly with some of these issues or whether they merely constitute a “repackaging” of the individualism agenda.

Description: An adsorption microcalorimeter was designed and built in our laboratory and used for the determination of differential adsorption heats in different samples of porous solids: activated carbon granules, activated carbon pellets, an activated carbon monolith and a zeolite sample. This work shows the relationship between adsorption heat and the pore size of different porous solids using adsorption of NH3, CO and N2O. The result shows that the thermal effect can be related with textural properties and superficial chemical groups of the studied porous solids. The values of differential heats of N2O adsorption in the investigated systems have shown that this interaction is weaker than that with CO. Small amounts of N2O are chemisorbed in the investigated systems. For the room temperature adsorption of N2O, the strongest active sites for the interaction with Brönsted acid groups in the ACM structure were identified. The values determined are between −60 kJ/mol and −110 kJ/mol for ZMOR and ACM, respectively, for the adsorption of N2O and −95 kJ/mol and −130 kJ/mol for the adsorption of CO.

Description: SUMMARY Zhou & Greenhalgh have recently presented an application of the Gaussian quadrature grid to seismic modelling in which the authors propose a meshing scheme that partitions the domain independently of the discontinuities in the media parameters. This comment aims to clarify the implications that this strategy has on the accuracy.

Description: The development of genetically engineered (GE) crops has focused predominantly on enhancing conventional pest control approaches. Scientific assessments show that these GE crops generally deliver significant economic and some environmental benefits over their conventional crop alternatives. However, emerging evidence indicates that current GE crops will not foster sustainable cropping systems unless the negative environmental and social feedback effects are properly addressed. Moreover, GE crop innovations that promote more sustainable agricultural systems will receive underinvestment by seed and chemical companies that must understandably focus on private returns for major crops. Opportunities to promote crops that convey multi-faceted benefits for the environment and the poor are foundational to a sustainable food system and should not be neglected because they also represent global public goods. In this paper, we develop a set of criteria that can guide the development of GE crops consistent with contemporary sustainable agriculture theory and practice. Based on those principles, we offer policy options and recommendations for reforming public and private R&D and commercialization processes to further the potential contributions of GE crops to sustainable agriculture. Two strategies that would help achieve this goal would be to restore the centrality of the public sector in agricultural R&D and to open the technology development process to more democratic participation by farmers and other stakeholders.

Description: Landscape scale conservation planning efforts have been in place for the past several decades to maintain biodiversity. Objectives of past efforts have been to identify areas to create reserves based on species diversity, land ownership, and landscape context. Risk analysis has not often been included in these spatial analyses. Datasets such as the Southwest Regional Gap Analysis (SWReGAP) are now available as are processes that allow risk analysis to be viewed in a spatial context in relations to factors that affect habitats over broad scales. We describe a method to include four spatial datasets to provide coarse scale delineation on areas to focus conservation including species numbers, key habitats, land management and factors that influence habitats. We used the SWReGAP management status dataset to identify management categories for long-term intent of management for biodiversity. The New Mexico Department of Game and Fish identified a set of 290 Species of Greatest Conservation Need (SGCN). Species occurrences for these species were associated with hydrologic unit codes from the National Hydrography Dataset (NHD). Key habitats were identified by using the SWReGAP land cover dataset and NHD derivatives. Factors that influence habitats were identified and scored for 89 land cover types and 23 aquatic habitats identified by the NMDGF. Our final model prioritizes landscapes that are within key habitats, have high numbers of terrestrial and aquatic Species of Greatest Conservation Need taxa, may be potentially altered by multiple effects that influence habitats, and lack long-term legally-binding management plans protecting them from anthropogenic degradation. Similar to other efforts, riparian and aquatic habitats were identified as the most important for conservation. This information may be displayed spatially, allowing land managers and decision makers to understand the ecological context where multiple effects of potential factors may influence some habitats greater than others, and repeat process with CWCS revisions.

Description: SUMMARY We obtained the crustal thickness ( H ) and average Vp / Vs ratio (κ) for the central (CNCC) and western North China Craton (WNCC) by H –κ stacking of receiver functions. Our results show that H and κ varies significantly but similarly in the northern and southern CNCC, reflecting possible effects of the Phanerozoic cratonotic reactivation on the region. Of the WNCC, the south is featured by κ values typical for Precambrian shields (∼1.77) and variations of H mirroring surface topography; the north presents large fluctuations of κ (1.68–1.93) and a crustal thickening (up to 50 km) uncorrelated with topography. Such N–S differences were thought to be associated with a ∼1.95-Ga suturing event occurred in the north with little effects on the south. Our observations suggest that distinct structural heterogeneities of the NCC not only reflect diverse regional tectonics in Phanerozoic, but also provide evidence for long-term survival of crustal fabrics in stable continental interior.

Description: SUMMARY The aim of this paper is to present an alternative method for the analysis of resistivity data. The methodology was developed during a study to evaluate if electrical resistivity can be used as a tool for analysing subsurface gas dynamics and gas emissions from landfills. The main assumption of this study was that variations in time of resistivity data correspond to variations in the relative amount of gas and water in the soil pores. Field measurements of electrical resistivity, static chamber gas flux and weather data were collected at a landfill in Helsingborg, Sweden. The resistivity survey arrangement consisted of nine lines each with 21 electrodes in an investigation area of 16 ×20 m. The ABEM Lund Imaging System provided vertical and horizontal resistivity profiles every second hour. The data were inverted in Res3Dinv using L 1 -norm–based optimization method with a standard least-squares formulation. Each horizontal soil layer was then represented as a linear interpolated raster model. Different areas underneath the gas flux measurement points were defined in the resistivity model of the uppermost soil layer, and the vertical extension of the zones could be followed at greater depths in deeper layer models. The average resistivity values of the defined areas were calculated and plotted on a time axis, to provide graphs of the variation in resistivity with time in a specific section of the ground. Residual variation of resistivity was calculated by subtracting the resistivity variations caused by the diurnal temperature variations from the measured resistivity data. The resulting residual resistivity graphs were compared with field data of soil moisture, precipitation, soil temperature and methane flux. The results of the study were qualitative, but promising indications of relationships between electrical resistivity and variations in the relative amount of gas and water in the soil pores were found. Even though more research and better data quality is necessary for verification of the results presented here, we conclude that this alternative methodology of working with resistivity data seems to be a valuable and flexible tool for this application.

Description: SUMMARY We analysed high-frequency body waves of local earthquakes to image the damage zone of the Calico fault in the eastern California shear zone. We used generalized ray theory and finite difference methods to compute synthetic seismograms for a low-velocity fault zone (FZ) to model the direct and FZ-reflected P and S traveltimes of local earthquakes recorded by a temporary array across the fault. The low velocity zone boundaries were determined by apparent traveltime delays across the fault. The velocity contrast between the fault zone and host rock was constrained by the traveltime delays of P and S waves and differential traveltimes between the direct and FZ-reflected waves. The dip and depth extent of the low velocity zone were constrained by a systematic analysis of direct P traveltimes of events on both sides of the fault. We found that the Calico fault has a ∼1.3-km-wide low velocity zone in which the P - and S -wave velocity decreased 40 and 50 per cent, respectively, with respect to the host rock. The low velocity zone dips 70° northeast and extends 3 km in depth.

Description: SUMMARY High-grade metamorphic rocks of the upper amphibolite to granulite facies in the Higher Himalayan crystalline (HHC) were studied for their appropriateness for palaeomagnetic investigations. Sampling and fieldwork was accomplished on a south–north transect at ∼87°E in the Khumbu region south of Everest. Paragenesis of iron-sulphides and oxides was found to be variable due to different source rocks and metamorphic conditions in this area. Unblocking temperature ranges, isothermal remanence acquisition, low-temperature magnetic properties and microscopic observations confirm the presence of pyrrhotite and magnetite. Two remanence components were separated with unblocking temperature intervals of 250–330 °C and 370–550 °C, which can be attributed to characteristic remanent magnetizations (ChRM) of monoclinic pyrrhotite (ChRM pyr ) and magnetite (ChRM mag ), respectively. Mainly the ChRM pyr showed a consistent demagnetization behaviour, probably residing in pyrrhotite particles within the single domain to pseudo-single domain range. Normal and reverse polarities were observed, partly co-existing within the unblocking spectrum of single specimens; this (i) supports a thermoremanent (TRM) origin of the ChRM pyr and (ii) indicates that the time span of remanence acquisition was long enough to average out palaeosecular variation. The age of the ChRM pyr can be related to the last metamorphic cooling event at the Early to Late Miocene boundary (16 Ma). The ChRM mag is probably a secondary chemical or thermochemical remanence acquired by hydrothermal processes during retrogression. ChRM pyr and ChRM mag directions are similar suggesting that pyrrhotite and magnetite were magnetized approximately contemporaneously. Separation of the sampling area into three sectors reveals a south–north trend with (i) variable block rotations in the southernmost part close to the Main Central Thrust (five sites), (ii) no rotation in the area between Lukla and Namche (10 sites) and (iii) a significant clockwise rotation ( N = 5 sites, D / I = 42.3°/29.4°, k = 20.7, α 95 = 17.2°) along the Imja Khola north of Namche. We interpret the clockwise rotation as an apparent rotation due to tilting by crustal doming in the Everest region. This doming obviously occurred north of a fault zone located at around Namche and was not affecting the region south of it. The clockwise sense of the apparent vertical-axis rotation indicates that the area north of Namche is located at the western side of the dome structure where westward tilting around a subhorizontal axis occurred (best fit: tilt direction 277°, tilt angle 51°).

Description: SUMMARY The deep structure of the Upper Rhine Graben (URG), a continental rift in SW Germany and E France, is still poorly known. This deficit impedes a full understanding of the geodynamic evolution of this prominent rift. We study the lithosphere–asthenosphere structure using teleseismic waveforms obtained from the passive broad-band TIMO project across the central URG. The recovered, crust-corrected traveltime residuals relative to the iasp91 earth model are tiny (mostly less than 0.2–0.3 s). The average measured slowness (〈1 s deg –1 ) and backazimuth (〈5°) deviations are also very small and do not show any systematic wave front anomalies. These observed perturbation values are smaller than expected ones from synthetic 3-D ray tracing modelling with anomalies exceeding 2–3 per cent seismic velocity in the mantle. Thus there is no significant hint for any deep-seated anomaly such as a mantle cushion, etc. This result means that the rifting process did not leave behind a lower lithospheric signature, which could be clearly verified with high-resolution teleseismic experiments. The only significant traveltime perturbation at the central URG is located at its western side in the upper crust around a known geothermal anomaly. The upper crustal seismic anomaly with traveltime delays of 0.2–0.3 s cannot be explained with increased temperature alone. It is possibly related to a zone of highly altered granite. In the west of our network a traveltime anomaly (0.6–0.7 s delay) related with the Eifel plume is confirmed by the TIMO data set.

Description: Mid-ocean ridge basalts (MORB) from the Arctic Ocean have been significantly less studied than those from other oceans. The Arctic ridges (Gakkel Ridge and Lena Trough) are ultraslow-spreading ridges with low melt productivity and are thus the best locations to investigate mantle heterogeneity. We report the major and trace element and Sr–Nd–Pb–Hf isotope compositions of basalts generated along the Lena Trough and the westernmost part of the Gakkel Ridge in the Arctic Ocean. Basalts from the northern Lena Trough and westernmost Gakkel Ridge (NLT–WGR) have compositions close to normal MORB. The geochemical composition of the NLT–WGR lavas confirms a binary mixing model involving melts from a depleted MORB mantle source and a Spitsbergen amphibole-bearing subcontinental lithospheric mantle (SCLM) source. In contrast, in the central part of the Lena Trough (CLT), the basalts are alkalic with relatively high Mg-number (60–65), high SiO 2 (51·0–51·6 wt %), Al 2 O 3 (18·1–18·4 wt %), Na 2 O (4·0–4·2 wt %), K 2 O (1·0–1·6 wt %), K 2 O/TiO 2 (0·6–0·9) and (La/Sm) PM (1·4–1·8), and low FeO (6·5–6·8 wt %) contents. These basalts display isotope variations with 87 Sr/ 86 Sr ranging from 0·70361 to 0·70390, 143 Nd/ 144 Nd from 0·51283 to 0·51290 ( Nd + 3·7 to +5·2), 176 Hf/ 177 Hf from 0·28313 to 0·28322 ( Hf + 11·6 to +14·9) and 206 Pb/ 204 Pb from 17·752 to 17·884, 207 Pb/ 204 Pb from 15·410 to 15·423 and 208 Pb/ 204 Pb from 37·544 to 37·670. These isotope compositions clearly distinguish the CLT lavas from those generated along the Gakkel Ridge. For the CLT lavas, involvement of a phlogopite- or amphibole- and (possibly garnet)-bearing SCLM source component is proposed. Owing to SCLM contamination along the entire length of the Lena Trough, we classify the Lena Trough as an ocean–continent transition boundary. Magmatism similar to that observed in the Lena Trough would be expected to occur wherever ocean spreading initiates.

Description: The urban heat island effect is linked to the built environment and threatens human health during extreme heat events. In this study, we analyzed whether characteristic land uses within an urban area are associated with higher or lower surface temperatures, and whether concentrations of “hot” land uses exacerbate this relationship. Zonal statistics on a thermal remote sensing image for the City of Toronto revealed statistically significant differences between high average temperatures for commercial and resource/industrial land use (29.1 °C), and low average temperatures for parks and recreational land (25.1 °C) and water bodies (23.1 °C). Furthermore, higher concentrations of either of these land uses were associated with more extreme surface temperatures. We also present selected neighborhoods to illustrate these results. The paper concludes by recommending that municipal planners and decision-makers formulate policies and regulations that are specific to the problematic land uses, in order to mitigate extreme heat.

Description: SUMMARY Recent revisions to the satellite-derived vertical gravity gradient (VGG) data reveal more detail of the ocean bottom and have allowed us to develop a non-linear inversion method to detect seamounts in VGG data. We approximate VGG anomalies over seamounts as sums of individual, partially overlapping, elliptical polynomial functions, which allows us to form a non-linear inverse problem by fitting the polynomial model to the observations. Model parameters for a potential seamount include geographical location, peak VGG amplitude, major and minor axes of the elliptical base, and the azimuth of the major axis. The non-linear inversion is very sensitive to the initial values for the location and amplitude; hence, they are constrained by the centre and amplitude of the uppermost contours obtained with a 1-Eötvös contour interval. With these initial conditions from contouring, we execute a step-wise and fully automated inversion and obtain optimal model estimates for potential seamounts; these are statistically evaluated for significance using the Akaike Information Criterion and  F  tests. A logarithmic barrier technique is applied to ensure positivity of all seamount amplitudes. After automatic and manual inspections of the model parameters we estimate actual heights and basal ellipses of the inspected potential seamounts directly from the predicted bathymetry grid. In this study, we find globally 24 643 potential seamounts ( h  ≥ 0.1 km) that are located away from continental margins; 8458 potential seamounts are taller than 1 km. Although our global estimate is significantly lower than predictions from previous studies, a first-order reconciliation of the size-frequency statistics obtained from those studies reveals that the previous counts are systematically overestimated. Because of the ambiguity of gravity signals due to small seamounts of  h  〈 1 km and the overlap with abyssal hills, we estimate the global seamount census to lie in the 40 000–55 000 range. The seamount data from this study are accessible from http://www.soest.hawaii.edu/PT/SMTS .

Description: The objective of this study is to validate an approach based on the change detection in multitemporal TerraSAR images (X-band) for mapping soil moisture in the Sahelian area. In situ measurements were carried out simultaneously with TerraSAR-X acquisitions on two study sites in Niger. The results show the need for comparing the difference between the rainy season image and a reference image acquired in the dry season. The use of two images enables a reduction of the roughness effects. The soils of plateaus covered with erosion crusts are dry throughout the year while the fallows show more significant moisture during the rainy season. The accuracy on the estimate of soil moisture is about 2.3% (RMSE) in comparison with in situ moisture contents.

Description: Negligence of former industrial sites (“brownfields”) has adversely impacted urban landscapes throughout the industrialized world. Brownfield redevelopment has recently emerged as a sustainable land use strategy and impetus for urban revitalization. This study presents a system dynamics model of the redevelopment process that illustrates how delays compound before realizing financial benefits from investment in these core urban areas. We construct a dynamic hypothesis, in which brownfield redevelopment activities are dependent upon funding and in turn bolster tax base through job creation. Drawing on previous studies, barriers to brownfield redevelopment are explored, including fear of liability, regulatory concerns, and uncertain cleanup standards and funding mechanisms. We model a case study of redevelopment in the State of Michigan (USA), which is informed by data from the Department of Environmental Quality (MDEQ) and U.S. Conference of Mayors brownfield surveys. Stock-flow structures represent phases of redevelopment, with diverted streams for sites in which no contamination was found (false alarms) and those with excess contamination level. The model is used to examine the point at which cumulative tax revenues from redeveloped areas exceed cumulative expenditures on brownfield redevelopment under different levels of funding availability.

Description: Typically the aim in the construction process is to calculate the energy, space and cost efficiency in the design phase. These factors’ influence on decision making extends to the whole building process. How these decisions affect the use of the building and user satisfaction as well as maintenance is still not that well understood. This study analyses different schools and day care centers and their energy as well as primary energy use. The buildings are located in southern Finland. Each building has had different objectives with respect to energy efficiency in the design phase. Our objective was to find out how those decisions made in the design and construction phase have influenced the overall energy performance of the building compared to existing building stock of similar building type. The results show that the studied buildings had lower thermal energy consumption compared to existing building stock. Thus the special attention in the design phase allowed achieving the desired goal. However, for the electricity consumption such a correlation could not be found. One of the reasons could be also different service level of buildings (more equipment). Also other quality values could not be compared since such data were not available from the existing building stock. As many earlier studies have indicated users have a high influence on the energy consumption. In the future, when feed-back from the users are obtained it will be interesting to analyze the results and compare what kind of influence that user behavior will have on the overall energy consumption of the studied buildings.

Description: SUMMARY The growth of the Earthscope/USArray Transportable Array (TA) has prompted the development of new methods in surface wave tomography that track phase fronts across the array and map the traveltime field for each earthquake or for each station from ambient noise. Directionally dependent phase velocities are determined locally by measuring the gradient of the observed traveltime field without the performance of a formal inversion. This method is based on the eikonal equation and is, therefore, referred to as ‘eikonal tomography’. Eikonal tomography is a bent-ray theoretic method, but does not account for finite frequency effects such as wave interference, wave front healing, or backward scattering. This shortcoming potentially may lead to both systematic bias and random error in the phase velocity measurements, which would be particularly important at the longer periods studied with earthquakes. It is shown here that eikonal tomography can be improved by using amplitude measurements to construct a geographically localized correction via the Helmholtz equation. This procedure should be thought of as a finite-frequency correction that does not require the construction of finite-frequency kernels and is referred to as ‘Helmholtz tomography’. We demonstrate the method with Rayleigh wave measurements following earthquakes between periods of 30 and 100 s in the western US using data from the TA. With Helmholtz tomography at long periods (〉50 s): (1) resolution of small-scale isotropic structures, which correspond to known geological features, is improved, (2) uncertainties in the isotropic phase velocity maps are reduced, (3) the directionally dependent phase velocity measurements are less scattered, (4) spurious 1-psi azimuthal anisotropy near significant isotropic structural contrasts is reduced, and (5) estimates of 2-psi anisotropy are better correlated across periods.

Description: The Jordan River basin is subject to extreme and increasing water scarcity. Management of transboundary water resources in the basin is closely intertwined with political conflicts in the region. We have jointly developed with stakeholders and experts from the riparian countries, a new dynamic consensus database and—supported by hydro-climatological model simulations and participatory scenario exercises in the GLOWA (Global Change and the Hydrological Cycle) Jordan River project—a basin-wide Water Evaluation and Planning (WEAP) tool, which will allow testing of various unilateral and multilateral adaptation options under climate and socio-economic change. We present its validation and initial (climate and socio-economic) scenario analyses with this budget and allocation tool, and invite further adaptation and application of the tool for specific Integrated Water Resources Management (IWRM) problems.

Description: SUMMARY Geomagnetic polarity transitions may be significantly more complex than are currently depicted in many sedimentary and lava-flow records. By splicing together paleomagnetic results from earlier studies at Steens Mountain with those from three newly studied sections of Oregon Plateau flood basalts at Catlow Peak and Poker Jim Ridge 70–90 km to the southeast and west, respectively, we provide support for this interpretation with the most detailed account of a magnetic field reversal yet observed in volcanic rocks. Forty-five new distinguishable transitional (T) directions together with 30 earlier ones reveal a much more complex and detailed record of the 16.7 Ma reversed (R)-to-normal (N) polarity transition that marks the end of Chron C5Cr. Compared to the earlier R-T-N-T-N reversal record, the new record can be described as R-T-N-T-N-T-R-T-N. The composite record confirms earlier features, adds new west and up directions and an entire large N-T-R-T segment to the path, and fills in directions on the path between earlier directional jumps. Persistent virtual geomagnetic pole (VGP) clusters and separate VGPs have a preference for previously described longitudinal bands from transition study compilations, which suggests the presence of features at the core–mantle boundary that influence the flow of core fluid and distribution of magnetic flux. Overall the record is consistent with the generalization that VGP paths vary greatly from reversal to reversal and depend on the location of the observer. Rates of secular variation confirm that the flows comprising these sections were erupted rapidly, with maximum rates estimated to be 85–120 m ka −1 at Catlow and 130–195 m ka −1 at Poker Jim South. Paleomagnetic poles from other studies are combined with 32 non-transitional poles found here to give a clockwise rotation of the Oregon Plateau of 11.4°± 5.6° with respect to the younger Columbia River Basalt Group flows to the north and 14.5°± 4.6° with respect to cratonic North America (95 per cent confidence interval).

Description: SUMMARY The principal aim of performing a survey or experiment is to maximize the desired information within a data set by minimizing the post-survey uncertainty on the ranges of the model parameter values. Using Bayesian, non-linear, statistical experimental design (SED) methods we show how industrial scale amplitude variations with offset (AVO) surveys can be constructed to maximize the information content contained in AVO crossplots, the principal source of petrophysical information from seismic surveys. The design method allows offset dependent errors, previously not allowed in non-linear geoscientific SED methods. The method is applied to a single common-midpoint gather. The results show that the optimal design is highly dependent on the ranges of the model parameter values when a low number of receivers is being used, but that a single optimal design exists for the complete range of parameters once the number of receivers is increased above a threshold value. However, when acquisition and processing costs are considered we find that a design with constant spatial receiver separation survey becomes close to optimal. This explains why regularly-spaced, 2-D seismic surveys have performed so well historically, not only from the point of view of noise attenuation and imaging in which homogeneous data coverage confers distinct advantages, but also to provide data to constrain subsurface petrophysical information.

Description: SUMMARY Tsunami generated by the Late Bronze Age (LBA) eruption of Thera were simulated using synthetic tide records produced for selected nearshore (∼20 m depths) sites of northern Crete, the Cyclades Islands, SW Turkey and Sicily. Inundation distances inland were also calculated along northern Crete. Modelling was performed by incorporating fully non-linear Boussinesq wave theory with two tsunamigenic mechanisms. The first involved the entry of pyroclastic flows into the sea, assuming a thick (55 m; 30 km 3 ) flow entering the sea along the south coast of Thera in three different directions all directed towards northern Crete, then a thin pyroclastic flow (1 m; 1.2 km 3 ) entering the sea along the north coast of Thera directed towards the Cyclades Islands. Flows were modelled as a solid block that slowly decelerates along a horizontal surface. The second mechanism assumed caldera collapse, of 19 km 3 and 34 km 3 modelled as a dynamic landslide producing rapid vertical displacements. Calculated nearshore wave amplitudes varied from a few metres to 28 m along northern Crete from pyroclastic flows, and up to 19 m from caldera collapse (34 km 3 volume). Inundation distances on Crete were 250–450 m. Waves produced by pyroclastic flows were highly focused, however, as a function of sea entry direction. Smaller volume pyroclastic flows produced nearshore wave amplitudes up to 4 m in the Cyclades islands north of Thera. Wave amplitudes in the Cyclades from smaller volume caldera collapse (19 km 3 ) were up to 24 m, whereas in SW Turkey were as low as 2.1 and 0.8 m (Didim and Fethye where LBA tsunami deposits have been found). Wave amplitudes for the larger volume caldera collapse (34 km 3 ) were generally 2.5–3 times larger than those generated by the smaller volume collapse (19 km 3 ). These results provide estimates for understanding possible consequences of tsunami impact in LBA coastal zones, thus providing criteria at archaeological sites for detecting inundation damage, as well as for contemporary hazard assessment; they also provide additional criteria for deciphering homogenite layers in the abyssal stratigraphy of the Ionian and eastern Mediterranean Seas.

Description: SUMMARY The first archaeomagnetic secular variation (SV) curves for the whole Southern Balkan Peninsula are presented. These are based on all data within a 700 km circle centred at Thessaloniki (40.60 o N, 23.00 o E). This data set consists of 325 directional and 625 intensity data mainly from Greece, Bulgaria, Serbia and southern Hungary. Some data from southern Italy are also included. The sliding moving window technique, was used to calculate a continuous SV curve for intensity while the directional SV curves were calculated using the bivariate extension of the Fisher statistics. These curves are well constrained and clearly show the main features of the geomagnetic field variation in this region during the last eight millennia. Comparisons with the predictions of the SCHA.DIF.3K and SCHA.DIF.8K regional and the CALS7K.2 and ARCH3K.1 global geomagnetic field models show a good agreement for the last 3000 years but differences for older times. The Balkan SV curves identify several rapid changes of the geomagnetic field in eastern Europe and can be used as reference curves for archaeomagnetic dating in the Balkan Peninsula.

Description: SUMMARY Receiver function analysis reveals a 7 km step-like change in crustal thickness across the Taranaki–Ruapehu Line (TR-line) of North Island, New Zealand. The TR-line runs east–west between the active andesite volcanoes of Mt Taranaki and Ruapehu and marks the southern-most extent of subduction zone volcanism in New Zealand. North of the TR-line receiver functions show a strong and sharp P -to- S ( Ps ) conversion at 25 ± 1.5 km depth, which is interpreted as a shallow Moho. At the TR-line the Moho Ps conversion deepens across a step to ∼32 km depth and weakens. Further south the Moho deepens to 〉35 km. Most of the 7 km step-change in crustal thickness occurs over a lateral distance of ∼8 km, yet there is little surface or topographic manifestation of the line. Given Fresnel zone considerations, the dip of the Moho offset could vary between 90° and 45°. Gravity, seismic attenuation and electrical data all show that the TR-line is not only a step in crustal thickness but also a profound lithospheric boundary as mantle properties, such as attenuation ( Q p −1 ), implied density, and electrical resistivity change abruptly across the line. An east–west oriented cluster of earthquakes with hypocentral depths of 20–40 km is centred on the Moho step. We propose that the Moho step, the earthquakes and the rapid change in mantle properties across the TR boundary are causally related. Processes that could be responsible for the phenomena described here include the rapid removal of mantle lithosphere and lower crust to the north of the TR-line.

Description: Both organic and conventional farming processes require energy input in the form of diesel fuel for farming equipment, animal feed, and fertilizer compounds. The most significant difference between the two methods is the use in conventional farming of mineral fertilizers and pesticides that are minimally employed in organic management. It is argued that organic farming is more environmentally friendly, given that synthetic fertilizers mainly used at conventional farms are replaced with animal manure and cover crops. Nutrient uptake by plants is additionally enhanced by the effective use of rhizobia and other types of plant growth-promoting bacteria, in combination with arbuscular mycorrhizal fungi. This article aims to compare the amounts and/or types of energy and nutrients required for both farming systems and provide feasible suggestions for the sustainable use of farm resources in combination with good crop yields.

Description: In 2009, just 27% of American teens with mobile phones reported using their devices to access the internet. However, teens from lower income families and minority teens were significantly more likely to use their phones to go online. Together, these surprising trends suggest a potential narrowing of the digital divide, offering internet access to those without other means of going online. This is an important move, as, in today’s society, internet access is central to active citizenship in general and teen citizenship in particular. Yet the cost of this move toward equal access is absorbed by those who can least afford it: Teenagers from low income households. Using survey and focus group data from a national study of “Teens and Mobile Phone Use” (released by Pew and the University of Michigan in 2010), this article helps identify and explain this and other emergent trends for teen use (as well as non-use) of the internet through mobile phones.

Description: Forest management planning in Finland is currently adopting a new-generation forest inventory method, which is based on interpretation of airborne laser scanning data and digital aerial images. The inventory method is based on a systematic grid, where the grid elements serve as inventory units, for which the laser and aerial image data are extracted and the forest variables estimated. As an alternative or a complement to the grid elements, image segments can be used as inventory units. The image segments are particularly useful as the basis for generation of the silvicultural treatment and cutting units since their boundaries should follow the actual stand borders, whereas when using grid elements it is typical that some of them cover parts of several forest stands. The proportion of the so-called mixed cells depends on the size of the grid elements and the average size and shape of the stands. In this study, we carried out automatic segmentation of two study areas on the basis of laser and aerial image data with a view to delineating micro-stands that are homogeneous in relation to their forest attributes. Further, we extracted laser and aerial image features for both systematic grid elements and segments. For both units, the feature set used for estimating the forest attributes was selected by means of a genetic algorithm. Of the features selected, the majority (61–79%) were based on the airborne laser scanning data. Despite the theoretical advantages of the image segments, the laser and aerial features extracted from grid elements seem to work better than features extracted from image segments in estimation of forest attributes. We conclude that estimation should be carried out at grid level with an area-specific combination of features and estimates for image segments to be derived on the basis of the grid-level estimates.

Description: SUMMARY We investigate the slip distribution of the 2008 May 12 Wenchuan  M s 7.9 earthquake using GPS data and InSAR measurements under a linear inversion scheme, with emphasis on the effect of three factors, including constraint on rake, different discretizations, and layered elastic model. Within our inversion parameterization context, we find the most influential factor would be constraint on rake. Without constraint on rake, the slip model seems physically wrong under the depth of 15 km, due to the limited depth resolution of the geodetic data used, especially the one orbit of InSAR measurements. Thus it is necessary to add  a priori  to the slip rake to obtain a reasonable fault source model. Different discretizations of the subfault patches have a notable impact on the slip distribution. Also, the layered elastic model predicts more slip at depth than does the half-space model, by about 15–20 per cent. The characteristics of slip distribution established through our inversions include some points as follows: (1) the most reliable results would be expected at depth range of 0–15 km; (2) three peak-slip asperities are inverted, at Yingxiu county (near epicentre), Yuejiashan county, and Beichuan county, respectively; (3) the inverted rake along the Yingxiu-Beichuan fault changes from dominantly thrusting motion at the southwest segments to dominant right-lateral or even pure right-lateral strike slip at the northeast segments, and (4) only thrust slip is occurred on the Guanxian-Jiangyou fault.

Description: SUMMARY Recordings from a temporary seismic network in the Rwenzori region of the East African Rift show seven earthquakes at a depth range between 53 and 60 km. Evidence from receiver functions indicates a crustal thickness of about 32 km in this region suggesting that the earthquakes are located well within the mantle beneath the rift. The occurrence of faulting at this depth is of relevance for models of continental-rift initiation and related deep magmatic processes.

Description: Narrow channels located in the longwave IR CO2 absorption region between approximately 13.2 and 14.5 μm, the well known CO2 slicing channels, have been proven to be quite effective for the estimates of cloud heights and effective cloud amounts as well as atmospheric temperature profiles. The designs of some of the near-future multi-channel earth observing satellite sensors cannot accommodate these longwave IR channels. Based on the analysis of the multi-channel imaging data collected with the NASA Moderate Resolution Imaging SpectroRadiometer (MODIS) instrument and on theoretical cloud radiative transfer modeling, we have found that narrow channels located at the midwave IR region between approximately 4.2 and 4.55 μm, where the combined CO2 and N2O absorption effects decrease rapidly with increasing wavelength, have similar properties as the longwave IR CO2 slicing channels. The scattering of solar radiation by clouds on the long wavelength side of the 4.3 μm CO2 absorption makes only a small contribution to the upwelling radiances. In order to retain the crucial cloud and temperature sensing capabilities, future satellite sensors should consider including midwave IR CO2 and N2O slicing channels if the longwave IR channels cannot be implemented on the sensors. The hyperspectral data covering the 3.7-15.5 mm wavelength range and measured with the Infrared Atmospheric Sounding Interferometer (IASI) can be used to further assess the utility of midwave IR channels for satellite remote sensing.

Description: SUMMARY We present an analysis on thermal properties, density and porosity of clastic, chemical/biochemical and intrasedimentary volcanic rocks collected from petroleum exploration wells of the Po Basin (Northern Italy). Moreover, we investigate the applicability of the Hashin–Shtrikman's model for a mineral aggregate in combination with the Zimmerman's model that takes pore shape into account, to calculate the bulk thermal conductivity. In case of macroscopically isotropic rocks, deviations between predicted and measured values range from −2.2 per cent to 6.9 per cent, and significantly decrease if a proper pore aspect ratio is chosen. Regarding the volumetric heat capacity, approximate estimates were obtained by means of a weighted average of the volumetric heat capacity of the mineral grains and that of the pore-filling water. The differences between the computed and measured values range from –6.2 per cent to 4.9 per cent and, on average, the computed volumetric heat capacity is lower by 1.6 per cent. The water loss during compaction and the temperature increase with depth are main factors controlling thermal properties. An anisotropy effect occurs in the case of rocks rich in sheet silicates. Due to rotation of these minerals, the vertical thermal conductivity of sheet silicates decreases exponentially with burial depth from 2.13 W m −1 K −1 at the surface to 0.52 W m −1 K −1 at 4.5 km. The laboratory data allow the formulation of compaction curves for the different sedimentary rock types. Examples of estimations of in situ vertical thermal conductivity and heat flow are finally given for two petroleum wells at which lithostratigraphic information is known in good detail.

Description: SUMMARY We present forward and adjoint spectral-element simulations of coupled acoustic and (an)elastic seismic wave propagation on fully unstructured hexahedral meshes. Simulations benefit from recent advances in hexahedral meshing, load balancing and software optimization. Meshing may be accomplished using a mesh generation tool kit such as CUBIT, and load balancing is facilitated by graph partitioning based on the SCOTCH library. Coupling between fluid and solid regions is incorporated in a straightforward fashion using domain decomposition. Topography, bathymetry and Moho undulations may be readily included in the mesh, and physical dispersion and attenuation associated with anelasticity are accounted for using a series of standard linear solids. Finite-frequency Fréchet derivatives are calculated using adjoint methods in both fluid and solid domains. The software is benchmarked for a layercake model. We present various examples of fully unstructured meshes, snapshots of wavefields and finite-frequency kernels generated by Version 2.0 ‘Sesame’ of our widely used open source spectral-element package SPECFEM3D.

Description: This article reviews and compares assessments of three biodiesel fuels: (1) transesterified lipids, (2) hydrotreated vegetable oils (HVO), and (3) woody biomass-to-liquid (BTL) Fischer-Tropsch diesel and selected feedstock options. The article attempts to rank the environmental performance and costs of fuel and feedstock combinations. Due to inter-study differences in goal and study assumptions, the ranking was mostly qualitative and intra-study results are emphasized. Results indicate that HVO made from wastes or by-products such as tall oil, tallow or used cooking oil outperform transesterified lipids and BTL from woody material, both with respect to environmental life cycle impacts and costs. These feedstock options are, however, of limited availability, and to produce larger volumes of biofuels other raw materials must also be used. BTL from woody biomass seems promising with good environmental performance and the ability not to compete with food production. Production of biofuels from agricultural feedstock sources requires much energy and leads to considerable emissions due to agrochemical inputs. Thus, such biodiesel fuels are ranked lowest in this comparison. Production of feedstock is the most important life cycle stage. Avoiding detrimental land use changes and maintaining good agricultural or forestry management practices are the main challenges to ensure that biofuels can be a sustainable option for the future transport sector.

Description: SUMMARY Allègre et al. recently presented new experimental data regarding the dependence of the streaming potential coupling coefficient with the saturation of the water phase. Such experiments are important to model the self-potential response associated with the flow of water in the vadose zone and the electroseismic/seismoelectric conversions in unsaturated porous media. However, the approach used to interpret the data is questionable and the conclusions reached by Allègre et al. likely incorrect

Description: SUMMARY Revil and Linde recently commented our paper concerning streaming potential (SP) measurements in unsaturated sand during drainage experiments. This comment suggests that the approach used to infer SP coefficients was inappropriate for unsaturated conditions, and therefore yielded wrong conclusions and ‘unphysical’ results regarding the behaviour of the relative SP coefficient. This reply argues that even if in Allègre et al . we neglected some secondary electrokinetic sources, the resulting conclusions are still representative of the behaviour of the true SP coefficient, and that the remarks of Revil & Linde arose from a misunderstanding of the drainage experiment conditions. We also find support for our results from a comparison between our observations and previous experimental studies.

Description: Environmental monitoring programs must efficiently describe state shifts. We propose using maximum entropy modeling to select dissimilar sampling sites to capture environmental variability at low cost, and demonstrate a specific application: sample site selection for the Central Plains domain (453,490 km2) of the National Ecological Observatory Network (NEON). We relied on four environmental factors: mean annual temperature and precipitation, elevation, and vegetation type. A “sample site” was defined as a 20 km × 20 km area (equal to NEON’s airborne observation platform [AOP] footprint), within which each 1 km2 cell was evaluated for each environmental factor. After each model run, the most environmentally dissimilar site was selected from all potential sample sites. The iterative selection of eight sites captured approximately 80% of the environmental envelope of the domain, an improvement over stratified random sampling and simple random designs for sample site selection. This approach can be widely used for cost-efficient selection of survey and monitoring sites.

Description: SUMMARY A small-aperture (40 m) short-period seismic array was installed during four months on the Vercors massif (Western French Alps) at the top of a limestone column which collapsed one month and a half later. During this monitoring period, 193 seismic events were recorded by the seven seismometers of the array. Signal analysis yielded three main types of nearby seismic events to be identified from temporal and spectral characteristics: microearthquakes (single or multiple events), individual block falls and rock falls. It turned out that 60 per cent of the 193 events were classified as microearthquakes, exhibiting distinct  P  and  S  waves, while 17 per cent of these events remained unclassified. Out of the microearthquakes, 40 events with a good signal-to-noise ratio were selected and processed.  P-  and  S -wave traveltimes were picked on the records and the inferred hypocentral distances agree with the two zones of the scarp exhibiting fresh ruptures after the fall. Polarization analysis of the 3-C records, along with numerical simulations, allowed discriminating between the two possible rupture mechanisms (toppling and sliding). Shear rupture (sliding) was the predominant mode in the lower part of the column whereas traction rupture (toppling) affected the upper part. Finally, the comparison between the ground motions recorded on the column and on the rock mass showed a systematic amplification on the column. Signal processing and numerical modelling both suggest that this amplification resulted from the excitation of the natural frequencies of the column and is particularly high (〉3) for microearthquakes occurring at the column-to-mass interface.

Description: SUMMARY The kinematic fractal source model presented in this study is able to simulate broad-band accelerograms with spectral amplitudes proportional to a fraction of the directivity coefficient C d in the far-field approximation. This approach is based on a composite source description, where subevents are generated using a fractal distribution of sizes, and, by summation, produces k -square space distribution of the slip. Each elementary source is described as a crack-type slip model growing circularly from a nucleation point when the rupture front reaches it. In order to better control the directivity effect, the location of the nucleation point for an elementary source is assumed to be scale-dependent. For the larger sources, the nucleation point is located near the intercept of the crack with the rupture front, whereas for smaller sources, it is randomly chosen within the crack. For simplicity, a constant rupture velocity is assumed. Each subevent is set up with scale-dependent rise-time, assuming a boxcar source–time function, hence filtering out its own high frequency radiation. The resulting mean slip-velocity functions are very similar to the ones derived from dynamic rupture modelling. Ground motion synthetics are computed by convolving the slip-velocity functions with the Green's functions. It is demonstrated that, in the far-field approximation, accelerogram spectra follow the ω 2 model with amplitudes controlled by the frequency-dependent directivity effects. In particular, spectral amplitudes at high-frequencies are proportional to a fraction of C d . These results were verified for few earthquake magnitudes. In addition, a validation exercise was made in the near-fault region by modelling the complete wavefield of strong ground motion at a few receivers and for several rupture scenarios. The synthetic strong-motion parameters are compared to the ones predicted by empirical attenuation relationships. It is shown that calculated standard deviations are in good agreement with the empirical ones, as well as the ground-motion parameter amplitudes predicted as a function of distance for whole interval of source distance considered in modelling. Minor differences were found in peak ground-accelerations computed at large distance from the fault, a problem related to the simplified response of the medium.

Description: Landsat MSS and ETM+ images were analyzed to detect 25-year land-cover change (1976–2001) in the critical Taguibo Watershed in Mindanao Island, Southern Philippines. This watershed has experienced historical modifications of its land-cover due to the presence of logging industries in the 1950s, and continuous deforestation due to illegal logging and slash-and-burn agriculture in the present time. To estimate the impacts of land-cover change on watershed runoff, land-cover information derived from the Landsat images was utilized to parameterize a GIS-based hydrologic model. The model was then calibrated with field-measured discharge data and used to simulate the responses of the watershed in its year 2001 and year 1976 land-cover conditions. The availability of land-cover information on the most recent state of the watershed from the Landsat ETM+ image made it possible to locate areas for rehabilitation such as barren and logged-over areas. We then created a “rehabilitated” land-cover condition map of the watershed (re-forestation of logged-over areas and agro-forestation of barren areas) and used it to parameterize the model and predict the runoff responses of the watershed. Model results showed that changes in land-cover from 1976 to 2001 were directly related to the significant increase in surface runoff. Runoff predictions showed that a full rehabilitation of the watershed, especially in barren and logged-over areas, will be likely to reduce the generation of a huge volume of runoff during rainfall events. The results of this study have demonstrated the usefulness of multi-temporal Landsat images in detecting land-cover change, in identifying areas for rehabilitation, and in evaluating rehabilitation strategies for management of tropical watersheds through its use in hydrologic modeling.

Description: SUMMARY Geodynamic models incorporating metamorphic phase transformations almost invariably assume the validity of the Boussinesq approximation that violates conservation of mass. In such models metamorphic density changes take place without volumetric effects. We assess the impact of the Boussinesq approximation by comparing models of orogeny accompanied by lower crustal eclogitization both with and without the approximation. Our results demonstrate that the approximation may cause errors approaching 100 per cent in characteristic measures of orogenic shape. Mass conservation errors in Boussinesq models amplify with model time. Mass conservative models of metamorphism are therefore essential to understand long-term tectonic evolution and to assess the importance of the different geodynamic processes.

Description: Community forestry and related small and medium forest enterprises (SMFEs) can contribute towards the achievement of REDD+ goals, since they can promote sustainable use and conservation of forests and, therefore, a reduction in forest-related carbon emissions. Additionally, they can improve the quality of life of forest-dependant people by generating alternative sources of income and employment. However, SMFEs often face a number of challenges, including non-conducive policy environments, inadequate business skills, and moreover, limited access to financial services. In this paper, we propose to direct a portion of REDD+ readiness efforts towards promoting the generation of an enabling environment for SMFEs that includes: the construction of an adequate Business Environment (BE), the provision of Business Development Services (BDS) and better access to Financial Services (FS). With the application of this framework, SMFEs will be more likely to proliferate and succeed, leading to enhanced community resilience and empowerment, in addition to increasing the likelihood of forest carbon stock permanence and the long term achievement of REDD+ goals. Opportunities and challenges of applying this approach in Latin America are discussed.

Description: The Western Alpine Sesia–Lanzo Zone (SLZ) is a sliver of eclogite-facies continental crust exhumed from mantle depths in the hanging wall of a subducted oceanic slab. Eclogite-facies felsic and basic rocks sampled across the internal SLZ show different degrees of retrograde metamorphic overprint associated with fluid influx. The weakly deformed samples preserve relict eclogite-facies mineral assemblages that show partial fluid-induced compositional re-equilibration along grain boundaries, brittle fractures and other fluid pathways. Multiple fluid influx stages are indicated by replacement of primary omphacite by phengite, albitic plagioclase and epidote as well as partial re-equilibration and/or overgrowths in phengite and sodic amphibole, producing characteristic step-like compositional zoning patterns. The observed textures, together with the map-scale distribution of the samples, suggest open-system, pervasive and reactive fluid flux across large rock volumes above the subducted slab. Thermodynamic modelling indicates a minimum amount of fluid of 0·1–0·5 wt % interacting with the wall-rocks. Phase relations and reaction textures indicate mobility of K, Ca, Fe and Mg, whereas Al is relatively immobile in these medium-temperature–high-pressure fluids. Furthermore, the thermodynamic models show that recycling of previously fractionated material, such as in the cores of garnet porphyroblasts, largely controls the compositional re-equilibration of the exhumed rock body.

Description: Remote sensing has become the most important data source for the digital elevation model (DEM) generation. DEM analyses can be applied in various fields and many of them require appropriate DEM visualization support. Analytical hill-shading is the most frequently used relief visualization technique. Although widely accepted, this method has two major drawbacks: identifying details in deep shades and inability to properly represent linear features lying parallel to the light beam. Several authors have tried to overcome these limitations by changing the position of the light source or by filtering. This paper proposes a new relief visualization technique based on diffuse, rather than direct, illumination. It utilizes the sky-view factor—a parameter corresponding to the portion of visible sky limited by relief. Sky-view factor can be used as a general relief visualization technique to show relief characteristics. In particular, we show that this visualization is a very useful tool in archaeology as it improves the recognition of small scale features from high resolution DEMs.

Description: The need for integrated methodological framework for sustainability assessment has been widely discussed and is urgent due to increasingly complex environmental system problems. These problems have impacts on ecosystems and human well-being which represent a threat to economic performance of countries and corporations. Integrated assessment crosses issues; spans spatial and temporal scales; looks forward and backward; and incorporates multi-stakeholder inputs. This study aims to develop an integrated methodology by capitalizing the complementary strengths of different methods used by industrial ecologists and biophysical economists. The computational methodology proposed here is systems perspective, integrative, and holistic approach for sustainability assessment which attempts to link basic science and technology to policy formulation. The framework adopts life cycle thinking methods—LCA, LCC, and SLCA; stakeholders analysis supported by multi-criteria decision analysis (MCDA); and dynamic system modelling. Following Pareto principle, the critical sustainability criteria, indicators and metrics (i.e., hotspots) can be identified and further modelled using system dynamics or agent based modelling and improved by data envelopment analysis (DEA) and sustainability network theory (SNT). The framework is being applied to development of biofuel supply chain networks. The framework can provide new ways of integrating knowledge across the divides between social and natural sciences as well as between critical and problem-solving research.

Description: SUMMARY We present a new 3-D transversely isotropic shear wave velocity model of the European and Mediterranean upper mantle obtained by analysis of surface waves. Data used are fundamental-mode Rayleigh and Love group velocity measurements in the period range 35–170 s, taken on seismograms recorded by European stations for regional earthquakes. The tomographic inversion to map the 3-D earth structure is split into two steps. First, we regionalize the group velocity dispersion measurements, obtaining distinct geographical group velocity maps at different periods; then, each local dispersion curve is inverted separately to find the shear wave velocity structure at depth. The inversion benefits from using a priori information from a 3-D global mantle model (S20RTS) and a new detailed European crustal model (EPcrust) to constrain the shallower layers. The inversion scheme follows a non-linear iterative algorithm by which Rayleigh and Love group slowness are inverted simultaneously for the best-fitting isotropic Voigt shear wave speed and radial anisotropy parameter ( v SH − v SV ) . Final merging of the v S profiles results in a new higher resolution 3-D model of European upper mantle. We find that Western Europe and Mediterranean Sea are mainly characterized by relatively low velocities, strongly contrasting with the fast roots of the Eastern European Craton. Many regional scale structures are also evident in the model, thus providing insights into the complex geodynamic framework of the European continent. Most prominent are the low-velocity West Mediterranean spreading basins and European Cenozoic rift system, and seismically fast features connected to subduction of Adria microplate, Hellenic Arc and Calabrian Arc. Radial anisotropy does not vary very significantly with respect to the PREM profile, as available data only resolve lateral variations to a limited degree due to trade-off with velocity. EPmantle has the potential to provide a reliable seismological reference for the upper-mantle structure in the broad European region.

Description: SUMMARY The main objective of this paper is to present a methodology for predicting spectrum compatible seismograms at sites where main shock is not recorded by other sites from which the main event has been observed using recorded aftershocks. The genetic algorithm and multitaper method are employed to minimize differences between the calculated synthetic seismogram and observed data. The multitaper spectrum approach being used permits the reduction of spectral leakage providing the possibility of matching the Fourier spectral amplitudes corresponding to synthetic and observed waveforms. The proposed technique is applied to the 2006 Silakhor earthquake ( M w = 6.1, Iran) as a case study. The six component waveforms (L, T and V) at two stations were simulated by minimizing errors between the synthetic response spectra, the Fourier spectral amplitudes and those of the observed. The results are shown to be in a good agreement with those of the observed data. The three components of waveforms at the other three stations were predicted incorporating the calculated parameters. The validity of the technique is proven through demonstrating good agreement of the response spectra and Fourier spectral amplitudes corresponding to the predicted and observed seismograms at other three stations. It is concluded that, the methodology is applicable for predicting realistic acceleration time histories to be used in seismic performance evaluation of existing structures in the region under study.

Description: SUMMARY Crustal density models derived from seismic velocity models by means of velocity–density conversions typically reproduce the main features of the observed gravity anomaly over the area but often show significant misfits. Given the uncertainty in the relationship between velocity and density, seismically derived density models should be regarded as an initial estimate of the true subsurface density structure. In this paper, we present a method for estimating the adjustments necessary to a seismically derived density model to improve the fit to gravity data. The method combines the Genetic Algorithm paradigm with linear inversion as a way to approach the non-linear and linear aspects of the problem. The models are divided into three layers representing the sedimentary column, the crystalline crust and the lithospheric mantle; the depths of these layers are determined from the seismic velocity model. Each of the layers is divided into a number of provinces and a density adjustment (Δρ) value is found for each province so that the residual gravity (difference between the observed gravity anomaly and the anomaly calculated for the seismically derived model) is minimized while keeping Δρ between predefined bounds. The preferred position of the province boundaries is found through the artificial evolution of a population of solutions. Given the stochastic nature of the algorithm and the non-uniqueness of the problem, different realizations can yield different solutions. By performing multiple realizations we can analyse a set of solutions by taking their mean and standard deviation, providing not only an estimate of the Δρ distribution in the subsurface but also an estimate of the associated uncertainty. Synthetic tests prove the ability of the algorithm to accurately recover the location of province boundaries and the Δρ values for a known model when using noise-free synthetic data. When noise is added to the data, the algorithm broadly recovers the features that define the known model despite greater standard deviations of the solutions and the occurrence of artefacts in the mean solutions. The algorithm was applied to four profiles across the Caribbean–South America Plate boundary. Some general patterns in the distribution of Δρ were observed consistently in the profiles and are correlated with the interpretations of the velocity models. Positive Δρ values in the sedimentary layer, negative Δρ values for island arc and extended island arc crust with an abrupt change to positive values in South American crust, and positive values in the mantle under the continent and island arc with a transition to negative values under the Caribbean oceanic crust.

Description: SUMMARY Active high-temperature (〉150 °C) geothermal areas like the Krafla caldera, NE-Iceland, often show distinct magnetic lows in aeromagnetic anomaly maps suggesting a destruction of magnetic minerals by hydrothermal activity. The main alteration processes in such an environment are low-temperature oxidation (〈350 °C, maghemitization) and fluid–rock interactions. We investigated the rock magnetic properties [natural remanent magnetization (NRM) magnetic susceptibility and their temperature and field variation] and the mineralogy, using X-ray diffraction, microscopic methods and electron microprobe analyses, of two drill cores (KH1 and KH3) from the rim of the Krafla caldera. The drill cores have distinctly lower NRM values (average 〈3 A m –1 ) compared to younger surface basalts (average 20 A m –1 ) along with a large variation in magnetic susceptibility (1.3 × 10 −7 – 4.9 × 10 −5 m 3 kg –1 ). The secondary mineral assemblage (sulphides, sphene, rutile and chlorite) indicates an alteration within the chlorite–smectite zone for both cores without depth zoning. Optical miscroscopy in combination with the Bitter technique and backscatter electron microscopy along with the thermomagnetic analyses allow distinguishing two different magnetomineralogical groups of titanomaghemite: (1) titanomaghemite with intermediate titanium concentration and probably high vacancy concentration, and (2) titanomaghemite with low titanium concentration and low vacancy concentration. The mineral assemblages, textures and magnetic properties deduced from the mentioned magnetic measurements indicate two-stage transformation mechanism: (1) Dissolution of titanium at low pH under oxidizing conditions. The ulvöspinel component of titanomagnetite and ilmenite forms rutile or sphene, and Fe 2 + migrates out of the spinel lattice forming titanomaghemite. (2) Formation of pyrite and dissolution of remaining titanomaghemite under reducing and acidic conditions. The latter mechanism produces ghost textures (all titanomaghemite is transformed and only their former grain shapes are preserved), with only paramagnetic minerals left and ferrimagnetic minerals nearly dissolved. This mechanism could explain the significant magnetization loss, which is seen in many local magnetic anomaly lows within the oceanic crust and volcanic islands like Iceland or Hawaii. The production of nanoporous textures in titanomaghemites is suggested as a mechanism for the enhancement of the magnetic susceptibility values related to the hydrothermal alteration of Krafla.

Description: SUMMARY We determine centroid-moment tensor (CMT) solutions by minimizing waveform differences between observed and simulated seismograms based on an adjoint method. Synthetic seismograms and Fréchet derivatives are calculated based on a spectral-element method. The non-linear adjoint CMT inversion algorithm requires three simulations for each iteration: one ‘forward’ simulation to obtain synthetics for the current source parameters, one ‘adjoint’ simulation which involves injecting time-reversed differences between observed and simulated seismograms as simultaneous virtual sources at each of the receivers, and an extra forward simulation to compute the step length in the conjugate-gradient direction. Whereas the vertical component of the adjoint wavefield reflects the radiation pattern near the centroid location, the components of the adjoint strain tensor capture the elements of the moment tensor. We use the method to determine adjoint CMT solutions for two representative southern California earthquakes using recent 3-D crustal model CVM-6.2. The adjoint CMT solutions are in good agreement with classical Hessian-based CMT solutions involving 3-D Green's functions. In general, adjoint CMT inversions require fewer numerical simulations than traditional Hessian-based inversions. This faster convergence holds promise for multiple moment-tensor and kinematic rupture inversions in 3-D earth models.

Description: SUMMARY We present a parallel goal-oriented adaptive finite element method that can be used to rapidly compute highly accurate solutions for 2.5-D controlled-source electromagnetic (CSEM) and 2-D magnetotelluric (MT) modelling problems. We employ unstructured triangular grids that permit efficient discretization of complex modelling domains such as those containing topography, dipping layers and multiple scale structures. Iterative mesh refinement is guided by a goal-oriented error estimator that considers the relative error in the strike aligned fields and their spatial gradients, resulting in a more efficient mesh refinement than possible with a previous approach based on the absolute errors. Reliable error estimation is accomplished by a dual weighted residual method that is carried out via hierarchical basis computations. Our algorithm is parallelized over frequencies, wavenumbers, transmitters and receivers, where adaptive refinement is performed in parallel on subsets of these parameters. Mesh sharing allows an adapted mesh generated for a particular frequency and wavenumber to be shared with nearby frequencies and wavenumbers, thereby efficiently reducing the parallel load of the adaptive refinement calculations. We demonstrate the performance of our algorithm on a large cluster computer through scaling tests for a complex model that includes strong seafloor topography variations and multiple thin stacked hydrocarbon reservoirs. In tests using up to 800 processors and a realistic suite of CSEM data parameters, our algorithm obtained run-times as short as a few seconds to tens of seconds.

Description: SUMMARY Geomagnetic induction responses such as geomagnetic depth sounding (GDS), magnetotelluric (MT), and horizontal geomagnetic transfer function (HTF) at long periods are used to estimate the electrical conductivity in the deep mantle. The responses in the period range that are shorter than 10 5 s (about 1 day) are in many cases considered to be local or regional induction problems in which the source field is approximated by plane waves and therefore the sphericity of Earth is not taken into account. In the period range between 10 4 and 10 5 s, the most dominant signature of the magnetic field variation is the solar quiet daily ( Sq ) variation and its higher harmonics. Therefore, when we obtain the responses due to the quasi-white background spectrum composed of plane waves, we regard the Sq field variations as noises in estimating the responses, and line spectra of the variations are removed from the observed time-series before the responses are calculated. However, with this approach, the calculated responses tend to possess a discontinuity at a period of about 10 4 s, and the response functions show common features at longer periods irrespective of the location of the observation site. It is particularly well known that the imaginary part of the induction vector tends to have a significant westward component for periods ranging between 10 4 and 10 5 s. Such features cannot be easily explained by the effect of the electrical conductivity structure alone. Examination of the phase of the HTF implies that the responses in the same period range are affected by the signature of sources of finite wavelength moving westward. Thus, it is suggested that the response functions in this period range were under the effect of the Sq field variations, even though the line spectra of them were removed and the responses were estimated at periods separate from the harmonic periods of Sq field variation. We examined how the influences of the Sq field appear on the responses numerically by a forward modelling. Results show most of the characteristic features in observed response functions can be ascribed to Sq source effects.

Description: The physiological principles predicting growth (3-PG) model is generally used to estimate gross primary productivity (GPP) in forest plantations. All existing parameter values in the 3-PG model for GPP estimation have been set as the standard values for eucalyptus and pine plantations. We propose that the 3-PG model can be applied to deciduous broadleaf forests dominated by Betula platyphylla via appropriate parameterization of their structure and functions. The allometric relationships between stem biomass and stem diameter, and between foliage biomass and stem biomass, were determined for the biomass partitioning ratio. Additionally, a temperature modifier was considered appropriate because it affected canopy quantum efficiency. After parameterization, the model showed a good correlation between the estimated results and the data from experimental plots in central and northern Japan. At both sites, GPP peaked around August and was 0 during the winter, when the canopy is bare of leaves. Furthermore, a sensitivity analysis was conducted to determine the most influential parameter relative to the output. GPP was sensitive to changes in canopy quantum efficiency and optimum temperature. Among the meteorological data used, solar radiation and temperature had great impacts on GPP, therefore, these parameters should be carefully considered to produce accurate results.

Description: Based on predicted changes in the magnitude and distribution of global precipitation, temperature and river flow under the A1B and A2 scenarios of the Intergovernmental Panel on Climate Change Special Report on Emissions Scenarios (IPCC SRES), this study assesses the potential impacts of climate change and CO2 fertilization on global agriculture, and its interactions with trade liberalization, as proposed for the Doha Development Round. The analysis uses the new version of the GTAP-W model, which distinguishes between rainfed and irrigated agriculture and implements water as an explicit factor of production for irrigated agriculture. Significant reductions in agricultural tariffs lead to modest changes in regional water use. Patterns are non-linear. On the regional level, water use may go up for partial liberalization, and down for more complete liberalization. This is because different crops respond differently to tariff reductions, and because trade and competition matter too. Moreover, trade liberalization tends to reduce water use in water scarce regions, and increase water use in water abundant regions, even though water markets do not exist in most countries. Considering impacts of climate change, the results show that global food production, welfare and GDP fall over time while food prices increase. Larger changes are observed under the SRES A2 scenario for the medium term (2020) and under the SRES A1B scenario for the long term (2050). Combining scenarios of future climate change with trade liberalization, countries are affected differently. However, the overall effect on welfare does not change much.

Description: The purpose of the paper is two-fold. First, to identify the main sustainability issues that Scottish food supply chain actors are concerned with and any differences that exist between primary producers, processors and distributors and consumers; and second, to explore the implications of respondents’ views for the direction of food and drink policy in Scotland. The analysis was based on a dataset assembled from the written responses to the National Food Policy discussion in Scotland, which contains opinions on the different dimensions of sustainability (economic, environmental and social) from a broad range of individuals and organizations representing different segments of the Scottish population. The empirical analyses involved comparing the responses according to two criteria: by food supply chain stakeholder and by geographical region. The results indicated that whilst there were differences among the studied groups, the importance of social and economic sustainability were strongly evident in the foregoing analysis, highlighting issues such as diet and nutrition, the importance of local food, building sustainability on sound economic performance, the market power of supermarkets, and regulation and support in building human and technical capabilities.

Description: As a member of the Annex 1 countries to the Kyoto Protocol of the United Nations Framework Convention on Climate Change, Japan is committed to reducing 6% of the greenhouse gas emissions. In order to achieve this commitment, Japan has undertaken several major mitigation measures, one of which is the domestic measure that includes ecologically friendly lifestyle programs, utilizing natural energy, participating in local environmental activities, and amending environmental laws. Mitigation policies could be achieved if public responses were strong. As the internet has increasingly become an online platform for sharing environmental information, public responses to the need for reducing greenhouse gas emissions may be assessed using available online tools. We used Google Insights for Search, Google AdWords Keyword Tool, and Google Timeline View to assess public responses in Japan based on the interest shown for five search terms that define global climate change and its mitigation policies. Data on online search interests from January 04, 2004 to July 18, 2010 were analyzed according to locations and categories. Our study suggests that the search interests for the five chosen search terms dramatically increased, especially when new mitigation policies were introduced or when climate change related events were organized. Such a rapid increase indicates that the Japanese public strongly responds to climate change mitigation policies.