Submarine pyroclastic deposits formed during the 20th May 2006 dome collapse of the Soufrière Hills Volcano, Montserrat

The 20th May 2006 lava dome collapse of the Soufrière Hills Volcano, Montserrat, had a total non-dense rock equivalent (non-DRE) collapse volume of approximately 115 × 10 6 m 3. The majority of this volume was deposited into the ocean. The collapse was rapid, 85% of the mobilized volume being removed in just 35 min, giving peak pyroclastic flow flux of 66 × 10 3 m 3 s -1. Channel and levee facies on the submarine flanks of the volcano and formation of a thick, steep-sided ridge, suggest that the largest and most dense blocks were transported proximally as a high concentration granular flow. Of the submerged volume, 30% was deposited from the base of this granular flow, forming a linear, high-relief ridge that extends 7 km from shore. The remaining 70% of the submerged volume comprises the finer grain sizes, which were transported at least 40 km by turbidity currents on gradients of <2°. At several localities, the May 2006 distal turbidity currents ran up 200 m of topography and eroded up to 20 cm of underlying substrate. Multiple turbidites are preserved, representing current reflection from the graben margins and deflection around topography. The high energy of the May 2006 collapse resulted in longer submarine run out than the larger (210 × 10 6 m 3) Soufrière Hills dome collapse in July 2003.

A molecular phylogeny for the Tribe Dacini (Diptera: Tephritidae) : systematic and biogeographic implications

With well over 700 species, the Tribe Dacini is one of the most species-rich clades within the dipteran family Tephritidae, the true fruit flies. Nearly all Dacini belong to one of two very large genera, Dacus Fabricius and Bactrocera Macquart. The distribution of the genera overlap in or around the Indian subcontinent, but the greatest diversity of Dacus is in Africa and the greatest diversity of Bactrocera is in south-east Asia and the Pacific. The monophyly of these two genera has not been rigorously established, with previous phylogenies only including a small number of species and always heavily biased to one genus over the other. Moreover, the subgeneric taxonomy within both genera is complex and the monophyly of many subgenera has not been explicitly tested. Previous hypotheses about the biogeography of the Dacini based on morphological reviews and current distributions of taxa have invoked an out-of-India hypothesis; however this has not been tested in a phylogenetic framework. We attempted to resolve these issues with a dated, molecular phylogeny of 125 Dacini species generated using 16S, COI, COII and white eye genes. The phylogeny shows that Bactrocera is not monophyletic, but rather consists of two major clades: Bactrocera s.s. and the ‘Zeugodacus group of subgenera’ (a recognised, but informal taxonomic grouping of 15 Bactrocera subgenera). This ‘Zeugodacus’ clade is the sister group to Dacus, not Bactrocera and, based on current distributions, split from Dacus before that genus moved into Africa. We recommend that taxonomic consideration be given to raising Zeugodacus to genus level. Supportive of predictions following from the out-of-India hypothesis, the first common ancestor of the Dacini arose in the mid-Cretaceous approximately 80 mya. Major divergence events occurred during the Indian rafting period and diversification of Bactrocera apparently did not begin until after India docked with Eurasia (50–35 mya). In contrast, diversification in Dacus, at approximately 65 mya, apparently began much earlier than predicted by the out-of-India hypothesis, suggesting that, if the Dacini arose on the Indian plate, then ancestral Dacus may have left the plate in the mid to late Cretaceous via the well documented India–Madagascar–Africa migration route. We conclude that the phylogeny does not disprove the predictions of an out-of-India hypothesis for the Dacini, although modification of the original hypothesis is required.

Evidence for carbonate platform failure during rapid sea-level rise; ca 14 000 year old bioclastic flow deposits in the Lesser Antilles

Bioclastic flow deposits offshore from the Soufrie`re Hills volcano on Montserrat in the Lesser Antilles were deposited by the largest volume sediment flows near this active volcano in the last 26 kyr. The volume of these deposits exceeds that of the largest historic volcanic dome collapse in the world, which occurred on Montserrat in 2003. These flows were most probably generated by a large submarine slope failure of the carbonate shelf comprising the south west flank of Antigua or the east flank of Redonda; adjacent islands that are not volcanically active. The bioclastic flow deposits are relatively coarse-grained and either ungraded or poorly graded, and were deposited by non cohesive debris flow and high density turbidity currents. The bioclastic deposit often comprises multiple sub-units that cannot be correlated between core sites; some located just 2 km apart. Multiple sub-units in the bioclastic deposit result from either flow reflection, stacking of multiple debris flow lobes, and/or multi-stage collapse of the initial landslide. This study provides unusually precise constraints on the age of this mass flow event that occurred at ca 14 ka. Few large submarine landslides have been well dated, but the slope failures that have been dated are commonly associated with periods of rapid sea-level change.

Anatomy of a submarine pyroclastic flow and associated turbidity current : July 2003 dome collapse, Soufrière Hills volcano, Montserrat, West Indies

The 12 to 13 July 2003 andesite lava dome collapse at the Soufrière Hills volcano, Montserrat, provides the first opportunity to document comprehensively both the sub-aerial and submarine sequence of events for an eruption. Numerous pyroclastic flows entered the ocean during the collapse, depositing approximately 90% of the total material into the submarine environment. During peak collapse conditions, as the main flow penetrated the air–ocean interface, phreatic explosions were observed and a surge cloud decoupled from the main flow body to travel 2 to 3 km over the ocean surface before settling. The bulk of the flow was submerged and rapidly mixed with sea water forming a water-saturated mass flow. Efficient sorting and physical differentiation occurred within the flow before initial deposition at 500 m water depth. The coarsest components (∼60% of the total volume) were deposited proximally from a dense granular flow, while the finer components (∼40%) were efficiently elutriated into the overlying part of the flow, which evolved into a far-reaching turbidity current.

Reconstructing the event stratigraphy from the complex structural - stratigraphic architecture of an Archaean volcanic – intrusive – sedimentary succession : the Boorara Domain, Eastern Goldfields Superterrane, Western Australia

Two main deformational phases are recognised in the Archaean Boorara Domain of the Kalgoorlie Terrane, Eastern Goldfields Superterrane, Yilgarn Craton, Western Australia, primarily involving southover- north thrust faulting that repeated and thickened the stratigraphy, followed by east northeast – west-southwest shortening that resulted in macroscale folding of the greenstone lithologies. The domain preserves mid-greenschist facies metamorphic grade, with an increase to lower amphibolite metamorphic grade towards the north of the region. As a result of the deformation and metamorphism, individual stratigraphic horizons are difficult to trace continuously throughout the entire domain. Volcanological and sedimentological textures and structures, primary lithological contacts, petrography and geochemistry have been used to correlate lithofacies between faultbounded structural blocks. The correlated stratigraphic sequence for the Boorara Domain comprises quartzo-feldspathic turbidite packages, overlain by high-Mg tholeiitic basalt (lower basalt), coherent and clastic dacite facies, intrusive and extrusive komatiite units, an overlying komatiitic basalt unit (upper basalt), and at the stratigraphic top of the sequence, volcaniclastic quartz-rich turbidites. Reconstruction of the stratigraphy and consideration of emplacement dynamics has allowed reconstruction of the emplacement history and setting of the preserved sequence. This involves a felsic, mafic and ultramafic magmatic system emplaced as high-level intrusions, with localised emergent volcanic centres, into a submarine basin in which active sedimentation was occurring.

Tephrochronology of marine sediments around the Island of Montserrat, Lesser Antilles volcanic arc

The recent history of the Soufrière Hills Volcano, Montserrat, Lesser Antilles volcanic arc, is reconstructed using data obtained from recently drilled submarine cores.Tephra layers in these cores preserve a record of the volcanic history of Montserrat back to ~250 ka on the basis of micropaleontology and stable isotope stratigraphy. Stratigraphic relationships identified in the cores collected in 2002 and 2005 document the fate of both pyroclastic flows entering the ocean to the east of Montserrat and carbonate-rich turbidites sourced from the carbonate platformssurrounding the islands of the Lesser Antilles. Using oxygen isotope stratigraphy, micropalaeontological analysis and Carbon-14 dating, it can be shown that three significant volcanic events, including the on-going eruption, have occurred over the last 12 ka. Preceding this was a time of volcanic quiescence, with three carbonate-rich turbidite events being documented in many of the cores. Our data suggest that these events occurred during Marine Isotope Stage 2, following the Last Glacial Maximum (LGM) and onset of post-glacial sea level rise.

Assessment of flood-related recharge to alluvial aquifers of an irrigated catchment following extended drought conditions using 3D visualisation and environmental tracers, Lockyer Valley, southeast Queensland, Australia.

The Lockyer Valley in southeast Queensland, Australia, hosts an economically significant alluvial aquifer system which has been impacted by prolonged drought conditions (~1997 to ~ 2009). Throughout this time, the system was under continued groundwater extraction, resulting in severe aquifer depletion. By 2008, much of the aquifer was at <30% of storage but some relief occurred with rains in early 2009. However, between December 2010 and January 2011, most of southeast Queensland experienced unprecedented flooding, which generated significant aquifer recharge. In order to understand the spatial and temporal controls of groundwater recharge in the alluvium, a detailed 3D lithological property model of gravels, sands and clays was developed using GOCAD software. The spatial distribution of recharge throughout the catchment was assessed using hydrograph data from about 400 groundwater observation wells screened at the base of the alluvium. Water levels from these bores were integrated into a catchment-wide 3D geological model using the 3D geological modelling software GOCAD; the model highlights the complexity of recharge mechanisms. To support this analysis, groundwater tracers (e.g. major and minor ions, stable isotopes, 3H and 14C) were used as independent verification. The use of these complementary methods has allowed the identification of zones where alluvial recharge primarily occurs from stream water during episodic flood events. However, the study also demonstrates that in some sections of the alluvium, rainfall recharge and discharge from the underlying basement into the alluvium are the primary recharge mechanisms of the alluvium. This is indicated by the absence of any response to the flood, as well as the observed old radiocarbon ages and distinct basement water chemistry signatures at these locations. Within the 3D geological model, integration of water chemistry and time-series displays of water level surfaces before and after the flood suggests that the spatial variations of the flood response in the alluvium are primarily controlled by the valley morphology and lithological variations within the alluvium. The integration of time-series of groundwater level surfaces in the 3D geological model also enables the quantification of the volumetric change of groundwater stored in the unconfined sections of this alluvial aquifer during drought and following flood events. The 3D representation and analysis of hydraulic and recharge information has considerable advantages over the traditional 2D approach. For example, while many studies focus on singular aspects of catchment dynamics and groundwater-surface water interactions, the 3D approach is capable of integrating multiple types of information (topography, geological, hydraulic, water chemistry and spatial) into a single representation which provides valuable insights into the major factors controlling aquifer processes.

Orthogonality defect and reduced search-space size for solving integer least-squares problems

In the context of ambiguity resolution (AR) of Global Navigation Satellite Systems (GNSS), decorrelation among entries of an ambiguity vector, integer ambiguity search and ambiguity validations are three standard procedures for solving integer least-squares problems. This paper contributes to AR issues from three aspects. Firstly, the orthogonality defect is introduced as a new measure of the performance of ambiguity decorrelation methods, and compared with the decorrelation number and with the condition number which are currently used as the judging criterion to measure the correlation of ambiguity variance-covariance matrix. Numerically, the orthogonality defect demonstrates slightly better performance as a measure of the correlation between decorrelation impact and computational efficiency than the condition number measure. Secondly, the paper examines the relationship of the decorrelation number, the condition number, the orthogonality defect and the size of the ambiguity search space with the ambiguity search candidates and search nodes. The size of the ambiguity search space can be properly estimated if the ambiguity matrix is decorrelated well, which is shown to be a significant parameter in the ambiguity search progress. Thirdly, a new ambiguity resolution scheme is proposed to improve ambiguity search efficiency through the control of the size of the ambiguity search space. The new AR scheme combines the LAMBDA search and validation procedures together, which results in a much smaller size of the search space and higher computational efficiency while retaining the same AR validation outcomes. In fact, the new scheme can deal with the case there are only one candidate, while the existing search methods require at least two candidates. If there are more than one candidate, the new scheme turns to the usual ratio-test procedure. Experimental results indicate that this combined method can indeed improve ambiguity search efficiency for both the single constellation and dual constellations respectively, showing the potential for processing high dimension integer parameters in multi-GNSS environment.