G-Probe-1 - An International Proficiency Test for Microprobe Laboratories - Report on Round 1 : February 2002 (TB-1 Basaltic Glass)

Results are presented for round one of a new international proficiency testing programme designed for microprobe laboratories involved in the routine analysis of silicate minerals. The sample used for this round was TB-1, a basaltic glass fused and prepared by the USGS. Thirty nine laboratories contributed data to this round, the majority of major element results being undertaken by EPMA and the majority of trace elements by LA-ICP-MS. Assigned values were derived from the median of results produced by nine selected laboratories that analysed powdered material by conventional ICP-MS, INAA and XRF techniques using bulk powders of the sample. Submitted microprobe results were evaluated using a target precision calculated using the Horwitz function, adopting the same criteria as those used for "applied" geochemistry laboratories in the companion GeoPT proficiency testing programme for laboratories involved in the routine bulk analysis of silicate rocks. An evaluation of results from participating microprobe laboratories indicated that overall, data were compatible with this precision function. A comparison between the performance of bulk and microprobe techniques used in the analysis of the basaltic glass showed remarkably good agreement, with significant bias only observed for the major oxide MgO.     

Field Survey and Numerical Simulations: A Review of the 1998 Papua New Guinea Tsunami

— The Papua New Guinea (PNG) tsunami of 1998 is re-examined through a detailed review of the field survey as well as numerous numerical computations. The discussion of the field survey explores a number of possible misinterpretations of the recorded data. The survey data are then employed by a numerical model as a validation tool. A Boussinesq model and a nonlinear shallow water wave (NLSW) model are compared in order to quantify the effect of frequency dispersion on the landslide-generated tsunami. The numerical comparisons indicate that the NLSW model is a poor estimator of offshore wave heights. However, due to what appears to be depth-limited breaking seaward of Sissano spit, both numerical models are in agreement in the prediction of maximum water elevations at the overtopped spit. By comparing three different hot-start initial profiles of the tsunami wave, it is shown that the initial shape and orientation of the tsunami wave is secondary to the initial displaced water mass in regard to prediction of water elevations on the spit. These numerical results indicate that agreement between numerical prediction of runup values with field recorded values at PNG cannot be used to validate either a NLSW tsunami propagation model or a specific landslide tsunami hot-start initial condition. Finally, with the use of traditional tsunami codes, a new interpretation of the PNG runup measurements is presented.     

Climatic and volcanic forcing revealed in a 50,000-year diatom record from Lake Massoko, Tanzania

The interactions between climatic and volcanic forcing on diatom communities contained in a 50,000-year sedimentary sequence from Lake Massoko, Tanzania, were examined. At the century scale, 19 discrete tephra inputs to the lake isolated the sedimentary nutrient supply and shifted the diatom communities to those tolerant of low phosphorus levels, whereas at the millennial scale, diatom-inferred shifts in precipitation–evaporation based on conductivity optima and diatom life-form ratios were broadly similar to lake-level reconstructions from Lake Rukwa, Lake Malawi, and others in the region. Some fluctuations of Lake Massoko are consistent with the precession-driven changes in insolation, but the major climate shifts do not relate directly to orbital forcing of summer insolation south of the equator and show more consistency with records from the equatorial and northern tropics that receive rainfall from the passing of the intertropical convergence zone. Sea surface temperatures are strongly correlated to multimillennial-scale climate patterns over this region of Africa.     

Calculating the two-dimensional magnetotelluric Jacobian in finite elements using reciprocity

To speed up the calculation of the field Jacobian for 2-D magnetoteliuric inversion using finite elements, the principle of electromagnetic reciprocity is applied. The governing relationship for the Jacobian of the field along strike is obtained by differentiating the Helmholtz equation with respect to the resistivity of each region in the finite-element mesh. The result is a similar Helmholtz equation for the Jacobian, with new sources distributed over all nodes within the parameter medium. However, according to the principle of electromagnetic reciprocity, the roles of sources and receivers are interchangeable. Utilizing reciprocity, the field values obtained from the original forward problem and for new unit sources imposed at the receivers are then utilized in the calculation of the Jacobian by simple multiplication and summation with finite-element terms at each rectangle in the mesh. For the auxiliary (across-strike) fields, the Jacobian terms are obtained by solving source vectors loaded with parabola coefficients used in the approximation to Maxwell's equations. Jacobian terms for the apparent resistivity ( p _a), the impedance phase (φ) and the vertical magnetic field ( K _zy) are then calculated utilizing the parallel- and auxiliary-field Jacobians. Comparison of Jacobian values obtained from reciprocity calculations and by differencing two forward solutions show that the reciprocity method is accurate and can be used to decrease the number of calculations required to obtain sensitivities by one to two orders of magnitude.