GGR Biennial Review: Atomic Absorption,Inductively Coupled Plasma-Atomic Emission Spectrometry,Neutron Activation Analysis and X-Ray Fluorescence Spectrometry Review for 2008–2009

This review describes developments major and in trace element determination using atomic absorption spectroscopy (AAS), inductively coupled plasma-atomic emission spectrometry (ICP-AES), instrumental neutron activation analysis (INAA), and X-ray fluorescence (XRF) spectrometry that were reported in 2008 and 2009. Publication levels were comparable to those of previous years, except for XRF which showed an increase in the number of articles published in the past 2 years. In terms of number of publications and impact, the most active field was AAS, while INAA was the least active. As expected for well-established and mature methods, novel developments for all four analytical techniques were relatively sparse. For AAS, the most notable publications concerned developments in sample introduction, particularly focussing on solid samples, increasing analyte sensitivity, and higher sample throughput. For ICP-AES, publication of developments in sample introduction, remediation of matrix effects and calibration continues. Compared with past years, there was a clear focus on sample preconcentration methods, but very few publications reporting new hyphenated speciation methods. For INAA, there were several publications exploring the accuracy and robustness of the method, as well as the requirements for INAA to meet criteria for a primary method of measurement. Two other related techniques, delayed neutron activation analysis, and prompt gamma activation analysis, were also described. Hyphenated XRF techniques showed interesting developments in enabling XRF and XRD analysis on the same spot, and further work characterising and calibrating three-dimensional micro-XRF shows promising results for investigating sample heterogeneity.     

The Proterozoic,albitite-hosted,Valhalla uranium deposit,Queensland, Australia: a description of the alteration assemblage associated with uranium mineralisation in diamond drill hole V39

The Valhalla uranium deposit, located 40 km north of Mount Isa, Queensland, Australia, is an albitite-hosted, Mesoproterozoic U deposit similar to albitite-hosted uranium deposits in the Ukraine, Sweden, Brazil and Guyana. Uranium mineralisation is hosted by a thick package of interbedded fine-grained sandstones, arkoses and gritty siltstones that are bound by metabasalts belonging to the ca. 1,780 Ma Eastern Creek Volcanics in the Western Succession of the Mount Isa basin. Alteration associated with U mineralisation can be divided into an early, main and late stage. The early stage is dominated by laminated and intensely altered rock comprising albite, reibeckite, calcite, (titano)magnetite ± brannerite. The main stage of mineralisation is dominated by brecciated and intensely altered rocks that comprise laminated and intensely altered rock cemented by brannerite, apatite, (uranoan)-zircon, uraninite, anatase, albite, reibeckite, calcite and hematite. The late stage of mineralisation comprises uraninite, red hematite, dolomite, calcite, chlorite, quartz and Pb-, Fe-, Cu-sulfides. Brannerite has U–Pb and Pb–Pb ages that indicate formation between 1,555 and 1,510 Ma, with significant Pb loss evident at ca. 1,200 Ma, coincident with the assemblage of Rodinia. The oldest ages of the brannerite overlap with ⁴⁰Ar/³⁹Ar ages of 1,533 ± 9 Ma and 1,551 ± 7 Ma from early and main-stage reibeckite and are interpreted to represent the timing of formation of the deposit. These ages coincide with the timing of peak metamorphism in the Mount Isa area during the Isan Orogeny. Lithogeochemical assessment of whole rock data that includes mineralised and unmineralised samples from the greater Mount Isa district reveals that mineralisation involved the removal of K, Ba and Si and the addition of Na, Ca, U, V, Zr, P, Sr, F and Y. U/Th ratios indicate that the ore-forming fluid was oxidised, whereas the crystal chemistry of apatite and reibeckite within the ore zone suggests that F⁻ and were important ore-transporting complexes. δ¹⁸O values of co-existing calcite and reibeckite indicate that mineralisation occurred between 340 and 380°C and involved a fluid having δ¹⁸O_fluid values between 6.5 and 8.6‰. Reibeckite δD values reveal that the ore fluid had a δD_fluid value between −98 and −54‰. The mineral assemblages associated with early and main stages of alteration, plus δ¹⁸O_fluid and δD_fluid values, and timing of the U mineralisation are all very similar to those associated with Na–Ca alteration in the Eastern Succession of the Mount Isa basin, where a magmatic fluid is favoured for this style of alteration. However, isotopic data from Valhalla is also consistent with that from the nearby Mount Isa Cu deposit where a basinal brine is proposed for the transport of metals to the deposit. Based on the evidence to hand, the source fluids could have been derived from either or both the metasediments that underlie the Eastern Creek Volcanics or magmatism that is manifest in the Mount Isa area as small pegmatite dykes that intruded during the Isan Orogeny.     

Transient deformation following the 30 January 1997 dike intrusion at Kīlauea volcano,Hawai’i

On 30 January 1997 an intrusion on Kīlauea volcano opened a new fissure within the East Rift Zone (ERZ) at Nāpau Crater, 3 km uprift from the ongoing eruptions at Pu’u ’Ō’ō. The fissure eruption lasted 22 h and opened a 5.1 km long, nearly vertical dike 1.9 m, extending from the surface to a depth of 2.4 km (Owen et al. 2000b). During the eruption, the lava pond at Pu’u ’Ō’ō drained, and eruptions ceased there. Pu’u ’Ō’ō eventually refilled in late February and eruptions resumed there on 28 March 1997. Continuous GPS data show a large transient following the 30 January 1997 dike intrusion. After lengthening 40 cm during the initial eruption, the baseline between two stations spanning the ERZ lengthened an additional 10 cm over the following 6 months. A coastal station KAEP also exhibited transient deformation, as it continued to move southward (5 cm) over the same 6-month period. The baseline between two stations spanning Kīlauea’s summit caldera contracted sharply during the eruption, but gradually recovered to slightly longer than its previous length 2 months after the intrusion. We use the extended network inversion filter (McGuire and Segall 2003) to invert continuous GPS data for volume change of a spherical pressure source under Kīlauea’s summit, opening distribution on a nearly vertical dike in the ERZ and potential slip on a decollement 9 km beneath the south flank. Following the 30 January intrusion, rift extension continued below the initial dike intrusion for the duration of the transient. Decollement slip, regardless of its assumed depth, is not required to fit the data. The modeled transient summit reinflation and rift opening patterns under Nāpau crater coincide with changes in observed behavior of Pu’u ’Ō’ō’s lava pond. Rift opening accelerated while Pu’u ’Ō’ō eruptions paused and began to decelerate after the lava pond reappeared nearly a month after the Nāpau eruption. The transient deformation is interpreted as resulting from shallow accommodation of the new dike volume.     

Forecasting Wave Amplitudes after the Arrival of a Tsunami

The destructive Pacific Ocean tsunami generated off the east coast of Honshu, Japan, on 11 March 2011 prompted the West Coast and Alaska Tsunami Warning Center (WCATWC) to issue a tsunami warning and advisory for the coastal regions of Alaska, British Columbia, Washington, Oregon, and California. Estimating the length of time the warning or advisory would remain in effect proved difficult. To address this problem, the WCATWC developed a technique to estimate the amplitude decay of a tsunami recorded at tide stations within the Warning Center’s Area of Responsibly (AOR). At many sites along the West Coast of North America, the tsunami wave amplitudes will decay exponentially following the arrival of the maximum wave (Mofjeld et al., Nat Hazards 22:71–89, 2000). To estimate the time it will take before wave amplitudes drop to safe levels, the real-time tide gauge data are filtered to remove the effects of tidal variations. The analytic envelope is computed and a 2 h sequence of amplitude values following the tsunami peak is used to obtain a least squares fit to an exponential function. This yields a decay curve which is then combined with an average West Coast decay function to provide an initial tsunami amplitude-duration forecast. This information may then be provided to emergency managers to assist with response planning.     

Imaging the ramp–décollement geometry of the Chelungpu fault using coseismic GPS displacements from the 1999 Chi-Chi, Taiwan earthquake

We use coseismic GPS data from the 1999 Chi-Chi, Taiwan earthquake to estimate the subsurface shape of the Chelungpu fault that ruptured during the earthquake. Studies prior to the earthquake suggest a ramp–décollement geometry for the Chelungpu fault, yet many finite source inversions using GPS and seismic data assume slip occurred on the down-dip extension of the Chelungpu ramp, rather than on a sub-horizontal décollement. We test whether slip occurred on the décollement or the down-dip extension of the ramp using well-established methods of inverting GPS data for geometry and slip on faults represented as elastic dislocations. We find that a significant portion of the coseismic slip did indeed occur on a sub-horizontal décollement located at 8 km depth. The slip on the décollement contributes 21% of the total modeled moment release. We estimate the fault geometry assuming several different models for the distribution of elastic properties in the earth: homogeneous, layered, and layered with lateral material contrast across the fault. It is shown, however, that heterogeneity has little influence on our estimated fault geometry. We also investigate several competing interpretations of deformation within the E/W trending rupture zone at the northern end of the 1999 ground ruptures. We demonstrate that the GPS data require a 22- to 35-km-long lateral ramp at the northern end, contradicting other investigations that propose deformation is concentrated within 10 km of the Chelungpu fault. Lastly, we propose a simple tectonic model for the development of the lateral ramp.