Evolution of the late Pleistocene Mojanda–Fuya Fuya volcanic complex (Ecuador), by progressive adakitic involvement in mantle magma sources

The Mojanda–Fuya Fuya Volcanic Complex consists of two nearby volcanoes, Mojanda and Fuya Fuya. The older one, Mojanda volcano (0.6 to 0.2 Ma), was first constructed by andesites and high-silica andesites forming a large stratovolcano (Lower Mojanda). This edifice was capped by a basaltic andesite and andesitic cone (Upper Mojanda), which collapsed later to form a 3-km-wide summit caldera, after large phreatomagmatic eruptions. The Lower Fuya Fuya edifice was constructed by the extrusion of viscous Si-rich andesitic lavas and dacitic domes, and the emission of a thick sequence of pyroclastic-flow and fallout deposits which include two voluminous rhyolitic layers. An intermediate construction phase at Fuya Fuya is represented by a mainly effusive cone, andesitic in composition (San Bartolo edifice), the construction of which was interrupted by a major sector collapse in the Late Pleistocene. Finally, a complex of thick siliceous lavas and domes was emplaced within the avalanche amphitheatre, forming the Upper Fuya Fuya volcanic centre. This paper shows that the general evolution from an effusive to an explosive eruptive style is related to a progressive adakitic contribution to the magma source. Although all the rocks of the complex are included in the medium-K field of continental arcs, the Fuya Fuya suite (61–75 wt.% SiO₂) shows depletion in Y and HREE and high Sr/Y and La/Yb values, compared to the less silicic Mojanda suite (55–66.5 wt.% SiO₂). The Mojanda calc-alkaline suite was generated by partial melting of an adakite-metasomatised mantle source that left a residue with 2% garnet, followed by fractional crystallization of dominant plagioclase + pyroxene + olivine at shallow, intra-crustal depths. For Fuya Fuya, geochemical and mineralogical data suggest either (1) partial melting of a similar metasomatised mantle with more garnet in the residue (4%), followed by fractional crystallization involving plagioclase, amphibole and pyroxene, or (2) mixing of mafic mantle-derived magma from the Mojanda suite and slab melts, followed by the same fractional crystallization process.     

An inter-comparison of tidal solutions computed with a range of unstructured grid models of the Irish and Celtic Sea Regions

Three finite element codes, namely TELEMAC, ADCIRC and QUODDY, are used to compute the spatial distributions of the M₂, M₄ and M₆ components of the tide in the sea region off the west coast of Britain. This region is chosen because there is an accurate topographic dataset in the area and detailed open boundary M₂ tidal forcing for driving the model. In addition, accurate solutions (based upon comparisons with extensive observations) using uniform grid finite difference models forced with these open boundary data exist for comparison purposes. By using boundary forcing, bottom topography and bottom drag coefficients identical to those used in an earlier finite difference model, there is no danger of comparing finite element solutions for “untuned unoptimised solutions” with those from a “tuned optimised solution”. In addition, by placing the open boundary in all finite element calculations at the same location as that used in a previous finite difference model and using the same M₂ tidal boundary forcing and water depths, a like with like comparison of solutions derived with the various finite element models was possible. In addition, this open boundary was well removed from the shallow water region, namely the eastern Irish Sea where the higher harmonics were generated. Since these are not included in the open boundary, forcing their generation was determined by physical processes within the models. Consequently, an inter-comparison of these higher harmonics generated by the various finite element codes gives some indication of the degree of variability in the solution particularly in coastal regions from one finite element model to another. Initial calculations using high-resolution near-shore topography in the eastern Irish Sea and including “wetting and drying” showed that M₂ tidal amplitudes and phases in the region computed with TELEMAC were in good agreement with observations. The ADCIRC code gave amplitudes about 30 cm lower and phases about 8° higher. For the M₄ tide, in the eastern Irish Sea amplitudes computed with TELEMAC were about 4 cm higher than ADCIRC on average, with phase differences of order 5°. For the M₆ component, amplitudes and phases showed significant small-scale variability in the eastern Irish Sea, and no clear bias between the models could be found. Although setting a minimum water depth of 5 m in the near-shore region, hence removing wetting and drying, reduced the small-scale variability in the models, the differences in M₂ and M₄ tide between models remained. For M₆, a significant reduction in variability occurred in the eastern Irish Sea when a minimum 5-m water depth was specified. In this case, TELEMAC gave amplitudes that were 1 cm higher and phases 30° lower than ADCIRC on average. For QUODDY in the eastern Irish Sea, average M₂ tidal amplitudes were about 10 cm higher and phase 8° higher than those computed with TELEMAC. For M₄, amplitudes were approximately 2 cm higher with phases of order 15° higher in the northern part of the region and 15° lower in the southern part. For M₆ in the north of the region, amplitudes were 2 cm higher and about 2 cm lower in the south. Very rapid M₆ tidal-phase changes occurred in the near-shore regions. The lessons learned from this model inter-comparison study are summarised in the final section of the paper. In addition, the problems of performing a detailed model–model inter-comparison are discussed, as are the enormous difficulties of conducting a true model skill assessment that would require detailed measurements of tidal boundary forcing, near-shore topography and precise knowledge of bed types and bed forms. Such data are at present not available.     

The line-of-sight velocity distributions of simulated merger remnants

We use Gauss–Hermite functions to study the line-of-sight velocity distributions in simulated merger remnants. Our sample contains 16 remnants; eight produced by mergers between disc galaxies of equal mass, and eight produced by mergers between disc galaxies with mass ratios of 3:1. The equal-mass mergers display a wide range of kinematic features, including counterrotation at large radii, orthogonally rotating cores and misaligned rotational axes. Most of the unequal-mass remnants exhibit fairly regular disc-like kinematics, although two have kinematics more typical of the equal-mass remnants. Our results may be compared with observations of early-type objects, including ellipticals with misaligned kinematic axes, counterrotating systems and S0 galaxies.     

Analytical perspective on trace element species of interest in exploration

Analysis of soil and sediment samples, using selective extraction methods to distinguish different phases, is of particular interest in exploration geochemistry to locate deeply buried mineral deposits. There are various mechanisms of binding labile elements in the secondary environment, including physical and chemical sorption, precipitation, chelation and complexation. Phases present in soils and sediments which are likely to scavenge ‘free' elements include amorphous Mn and Fe oxides, the humic and fulvic components of humus, and clays. This paper reviews these forms of trace elements and the methods in current use to quantify them. Examples of precision data, both for control and survey samples, are given with respect to trace elements dissolved from the ‘soluble organic' component of humus, Mn oxides and amorphous Fe oxides. The high sensitivity of inductively coupled plasma mass spectrometry (ICP–MS) is required to measure accurately and precisely a large suite of trace elements, especially where only small fractions of elements are dissolved by such leaches as the commercially available Enzyme and MMI (Mobile Metal Ion) extractions. The relative standard deviations (RSD) obtained for 33 elements (e.g. Ag, Cd, In, I) in the standard reference sample (SRM), TILL-2, are in the range 0.5–8% for the hydroxylamine hydrochloride (NH₂OH·HCl) leach designed to extract hydrous Fe and Mn oxides. The corresponding RSDs for elements in the reactive Mn oxide phase extracted by the Enzyme leach are in the range 3–19% except for some trace elements at levels close to detection limit (e.g. Cd, Bi). The RSDs obtained for field duplicates are inferior to those for analytical replicates (i.e. sample splits), probably a reflection of different concentrations of the host phase. In one soil survey, the Fe extracted by a 0.25 M NH₂OH·HCl leach ranged conservatively from 0.2 to 1.7% whereas the Mn extracted by the Enzyme leach varied extensively, from 0.3 to >999 ppm. In contrast, precision, at 1–7% RSD, for field duplicates was found to be comparable with that for both analytical duplicates and the SRM, LKSD-4, for elements associated with the humic and fulvic component of humus samples sieved to <177 μm.     

The Vulnerability of the Australian Beef Industry to Impacts of the Cattle Tick (Boophilus microplus) under Climate Change

An integrated assessment is presented of the potential impacts of the cattle tick (Boophilus microplus Canestrini) on the Australian beefindustry under climate change. The project was carried out as a case study to test an impact assessment approach that was designed to integrate biological, production and socio-economic impacts on managed and natural systems. A climate-driven, tick population model was run for European, zebu and crossbred cattle breeds having different levels of resistance to cattle ticks. A geographical information system (GIS) was used to organise spatial data on climate scenarios and industry statistics and to undertake regional analyses.A comparison was made of the two available approaches to conducting impact assessments, namely a bottom-up approach using sensitivity analysis and a top-down approach using climate change scenarios from a global circulation model (GCM) (CSIRO, 1996). The output, in terms of the abundance of tick populations and reductions in cattle productivity for each breed showed significant expansions in potential geographical impacts. In the absence of any adaptation measures, the results indicated changes in the losses in live weight gain of cattle tick ranging from 7780 tonnes per year by 2030 to 21637 tonnes per year by 2100, in comparison with estimates for current losses of 6594 tonnes per year.The principal adaptation options available to the beef industry are to switch to breeds that are more resistant to cattle ticks, or to increase the frequency of treatments with various tick control products. In this paper we focus on switching breeds as an adaptive measure when appropriate damage thresholds are triggered under the climate change scenarios. When adaptation measures were put in place, the losses ranged from 4962 tonnes in 2030 to 5619 tonnes in 2100 compared with 2636 tonnes at present if all producers adopted the optimal breed structure. Optimal breed structure was defined as one that would prevent tick numbers per animal exceeding 100 ticks per animal for European and 700 ticks per animal for crossbred breeds of cattle in any week of the year under a tick control strategy that was suitable for present climatic conditions. The lower threshold for European breeds reflects their vulnerability to explosive increases in numbers because of their low resistance to ticks. The results of the analyses using the GCM scenarios were used in an economic model to calculate costs of lost live-weight gain for 2030, 2070 and 2100. The greatest increases in costs were incurred in the southern parts of the current distribution in Queensland and potentially in northern New South Wales if the present quarantine barrier failed.Given the great uncertainty of the nature of possible regional changes in climate, analyses of the sensitivity of losses in live weight gain to changes in climatic variables were also undertaken. The analyses included a measure of likely impacts of cattle tick on the beef cattle industry, in the absence of adaptation measures, as a baseline measure of sensitivity. The likely impacts on crossbred cattle were insensitive to the climatic variables.When adaptive breed changes were allowed, the economic impacts on the industry were insensitive to the GCM scenarios. This suggests that, at least in this instance, reducing the uncertainties in climate change scenarios is not a priority if the adaptation strategies can be implemented in a cost-effective manner. Finally we made a qualitative assessment of the sustainability and robustness of alternative approaches to adaptation and assessed regional vulnerability to cattle tick under climate change. The conclusions were so strongly dependent on assumptions about the future of other global changes, in particular the ability to maintain quarantine barriers and to retain effective acaricides at comparable costs to the present, that we strongly recommend that risk assessments of climate change extend to all relevant variables in involved in global change where possible.     

Bulk density effects on soil hydrologic and thermal characteristics: A numerical investigation

Soil bulk density (ρ_b) is commonly treated as static in studies of land surface dynamics. Magnitudes of errors associated with this assumption are largely unknown. Our objectives were to (a) quantify ρ_b effects on soil hydrologic and thermal properties and (b) evaluate effects of ρ_b on surface energy balance and heat and water transfer. We evaluated 6 soil properties, volumetric heat capacity, thermal conductivity, soil thermal diffusivity, water retention characteristics, hydraulic conductivity, and vapour diffusivity, over a range of ρ_b, using a combination of 6 models. Thermal conductivity, water retention, hydraulic conductivity, and vapour diffusivity were most sensitive to ρ_b, each changing by fractions greater than the associated fractional changes in ρ_b. A 10% change in ρ_b led to 10–11% change in thermal conductivity, 6–11% change in saturated and residual water content, 49–54% change in saturated hydraulic conductivity, and 80% change in vapour diffusivity. Subsequently, 3 field seasons were simulated with a numerical model (HYDRUS‐1D) for a range of ρ_b values. When ρ_b increased 25% (from 1.2 to 1.5 Mg m^?3), soil temperature variation decreased by 2.1 °C in shallow layers and increased by 1 °C in subsurface layers. Surface water content differed by 0.02 m³ m^?3 for various ρ_b values during drying events but differences mostly disappeared in the subsurface. Matric potential varied by >100 m of water. Surface energy balance showed clear trends with ρ_b. Latent heat flux decreased 6%, sensible heat flux increased 9%, and magnitude of ground heat flux varied by 18% (with a 25% ρ_b increase). Transient ρ_b impacted surface conditions and fluxes, and clearly, it warrants consideration in field and modelling investigations.     

A new approach for measuring surface hydrological connectivity

The development of surface hydrological connectivity is a key determinant of flood magnitude in drylands. Thresholds in runoff response may be reached when isolated runoff-generating areas connect with each other to form continuous links to river channels, enabling these areas to contribute to flood hydrographs. Such threshold behaviour explains observed nonlinearities and scale dependencies of dryland rainfall–runoff relationships and complicates attempts at flood prediction. However, field methods for measuring the propensity of a surface to transmit water downslope are lacking, and conventional techniques of infiltration measurement are often inappropriate for use on non-agricultural drylands. Here, we argue for a reconceptualization of the dryland surface runoff process, suggesting that the downslope transfer of water should be considered alongside surface infiltration; that is, there is a need for the “aggregated” measurement of infiltration and overland flow hydraulics. Surface application of a set volume of water at a standardized rate generates runoff that travels downslope; the distance it travels downslope is determined by infiltration along the flow, integration of flow paths, and flow resistance. We demonstrate the potential of such a combined measurement system coupled with structure-from-motion photogrammetry to identify surface controls on runoff generation and transfer on dryland hillslopes, with vegetation, slope, surface stone cover, and surface roughness all having a significant effect. The measurement system has been used on slopes up to 37° compared with the flat surface typically required for infiltration methods. On average, the field workflow takes ~10–15 min, considerably quicker than rainfall simulation. A wider variety of surfaces can be sampled with relative ease, as the method is not restricted to stone and vegetation-free land. We argue that this aggregated measurement represents surface connectivity and dryland runoff response better than standard hydrological approaches and can be applied on a much greater variety of dryland surfaces.     

Radiative Shock Experiments At Luli

We present the set-up and the results of a supercritical radiative shock experiment performed with the LULI nanosecond laser facility. Using specific designed targets filled with xenon gaz at low pressure, the propagation of a strong shock with a radiative precursor is evidenced. The main measured quantities related to the shock (electronic density, propagation velocities, temperature, radial dimension) are presented and compared with various numerical simulations.     

Impact of India–Asia collision on SE Asia: The record in Borneo

Borneo occupies a central position in the Sundaland promontory of SE Asia. It has a complex Cenozoic geological history of sedimentation and deformation which began at about the same time that India is commonly suggested to have started to collide with Asia. Some tectonic reconstructions of east and SE Asia interpret a large SE Asian block with Borneo at its centre which has been rotated clockwise and displaced southwards along major strike–slip faults during the Cenozoic due to the indentation of Asia by India. However, the geological history of Borneo is not consistent with the island simply forming part of a large block extruded from Asia. The large clockwise rotations and displacements predicted by the indentor model for Borneo are incompatible with palaeomagnetic evidence and there is no evidence that the major strike–slip faults of the Asian mainland reach Borneo. Seismic tomography shows there is a deep high velocity anomaly in the lower mantle beneath SE Asia interpreted as subducted lithosphere but it can be explained just as well by alternative tectonic models as by the indentor model. Very great thicknesses of Cenozoic sediments are present in Borneo and circum-Borneo basins, and large amounts of sediment were transported to the Crocker turbidite fan of north Borneo from the Eocene to the Early Miocene, but all evidence indicates that these sediments were derived from local sources and not from distant sources in Asia elevated by India–Asia collision. The Cenozoic geological history of Borneo records subduction of the proto-South China Sea and Miocene collision after this ocean lithosphere was eliminated, and a variety of effects resulting from long-term subduction beneath SE Asia. There is little to indicate that India–Asia collision has influenced the Cenozoic geological record in Borneo.