Volcanic and Scientific Activity at Kick ’em Jenny Submarine Volcano 2001–2002: Implications for Volcanic Hazard in the Southern Grenadines,Lesser Antilles

Kick em Jenny submarine volcano, ~8 km north of Grenada, has erupted at least 12 times since it was first discovered in 1939, making it the most frequently active volcano in the Lesser Antilles arc. The volcano lies in shallow water close to significant population centres and directly beneath a major shipping route, and as a consequence an understanding of the eruptive behaviour and potential hazards at the volcano is critical. The most recent eruption at Kick em Jenny occurred on December 4 2001, and differed significantly from past eruptions in that it was preceded by an intensive volcanic earthquake swarm. In March 2002 a multi-beam bathymetric survey of the volcano and its surroundings was carried out by the NOAA ship Ronald H Brown. This survey provided detailed three-dimensional images of the volcano, revealing the detailed morphology of the summit area. The volcano is capped by a summit crater which is breached to the northeast and which varies in diameter from 300 to 370 m. The depth to the summit (highest point on the crater rim) is 185 m and the depth to the lowest point inside the crater is 264 m. No dome is present within the crater. The crater and summit region of Kick em Jenny are located at the top of an asymmetrical cone which is about 1300 m from top to bottom on its western side. It lies within what appear to be the remnants of a much larger arcuate collapse structure. An evaluation of the morphology, bathymetry and eruptive history of the volcano indicates that the threat of eruption-generated tsunamis is considerably lower than previously thought, mainly because the volcano is no longer thought to be growing towards the surface. Of more major and immediate concern are the direct hazards associated with the volcano, such as ballistic ejecta, water disturbances and lowered water density due to degassing.     

δ34S isotope values of dissolved sulfate (SO4 2−) as a tracer for battery acid (H2SO4) contamination in groundwater

The application of sulfur isotope (³⁴S) values of sulfate in groundwater provided the information necessary to evaluate the source, transport and fate of battery acid and associated contaminants at the Gulf Coast Recycling (GCR) facility. The chemical and isotopic composition of groundwater beneath the (GCR) property, a battery recycling facility in east Tampa, Florida, varies more than expected for an area of comparable size. Sulfate (SO₄^2–) values, for example, range from 1.2 to 11,500 mg/L and oxygen and hydrogen isotopes do not attenuate towards the weighted annual mean. Those samples that are high in sulfate generally have a low pH, which immediately indicates battery acid (H₂SO₄) contamination as a potential source for the sulfate. The low pH and high reactivity of the sulfuric acid groundwater cause the formation of hydrogeological microenvironments due to preferential dissolution of carbonate minerals, which in turn causes enhanced recharge and groundwater flow in certain areas; thus, the extreme scatter in the data set. Because of the difficult hydrogeology it is not straightforward to delineate the point-sources of contamination and up to five potential scenarios have to be evaluated: (1) seawater intrusion, (2) upwelling of high-sulfate groundwater, (3) local dissolution of gypsum, (4) an up-gradient contaminant source to the northeast of the GCR property and (5) battery acid contamination.     

An integrated linear/non-linear flow model for the conduit-fissure-pore media in the karst triple void aquifer system

Most karstic aquifer media may be characterized as the triple-void media with highly-varied hydraulic properties, including matrix pore, fissure and conduit, in which liner flow may co-exist with non-linear flow. In this paper, an attempt is made to couple linear flow with non-linear flow in a single unified flow governing equations by introducing the concept of equivalent hydraulic conductivity (EHC) and deriving a general Darcys law for various flow. The expression of EHC in the karst conduit and fissure are also derived. The procedures of numerical implementation are demonstrated via an ideal model and a case study of karst aquifer system in the Beishan Ore Formation area, Guangxi Autonomous Region, China.     

Community preparedness for lava flows from Mauna Loa and Hualālai volcanoes,Kona, Hawai‘i

Lava flows from Mauna Loa and Huallai volcanoes are a major volcanic hazard that could impact the western portion of the island of Hawaii (e.g., Kona). The most recent eruptions of these two volcanoes to affect Kona occurred in a.d. 1950 and ca. 1800, respectively. In contrast, in eastern Hawaii, eruptions of neighboring Klauea volcano have occurred frequently since 1955, and therefore have been the focus for hazard mitigation. Official preparedness and response measures are therefore modeled on typical eruptions of Klauea.The combinations of short-lived precursory activity (e.g., volcanic tremor) at Mauna Loa, the potential for fast-moving lava flows, and the proximity of Kona communities to potential vents represent significant emergency management concerns in Kona. Less is known about past eruptions of Huallai, but similar concerns exist. Future lava flows present an increased threat to personal safety because of the short times that may be available for responding.Mitigation must address not only the specific characteristics of volcanic hazards in Kona, but also the manner in which the hazards relate to the communities likely to be affected. This paper describes the first steps in developing effective mitigation plans: measuring the current state of peoples knowledge of eruption parameters and the implications for their safety. We present results of a questionnaire survey administered to 462 high school students and adults in Kona. The rationale for this study was the long lapsed time since the last Kona eruption, and the high population growth and expansion of infrastructure over this time interval. Anticipated future growth in social and economic infrastructure in this area provides additional justification for this work.The residents of Kona have received little or no specific information about how to react to future volcanic eruptions or warnings, and short-term preparedness levels are low. Respondents appear uncertain about how to respond to threatening lava flows and overestimate the minimum time available to react, suggesting that personal risk levels are unnecessarily high. A successful volcanic warning plan in Kona must be tailored to meet the unique situation there.     

Theory of Basin Scale Dynamics of a Stratified Rotating Fluid

The dynamics of a stratified fluid contained in a rotating rectangular box is described in terms of the evolution of the lowest moments of its density and momentum fields. The first moment of the density field also gives the position of the fluids centre-of-mass. The resulting low-order model allows for fast assessment both of adopted parameterisations, as well as of particular values of parameters. In the ideal fluid limit (neglect of viscous and diffusive effects), in the absence of wind, the equations have a Hamiltonian structure that is integrable (non-integrable) in the absence (presence) of differential heating. In a non-rotating convective regime, dynamically rich behaviour and strong dependence on the single (lumped) parameter are established. For small values of this parameter, in a self-similar regime, further reduction to an explicit map is discussed in an Appendix. Introducing rotation in a nearly geostrophic regime leads through a Hopf bifurcation to a limit cycle, and under the influence of wind and salt to multiple equilibria and chaos, respectively.     

Emplacement mechanisms of the South Kona slide complex,Hawaii Island: sampling and observations by remotely operated vehicle Kaiko

Emplacement of a giant submarine slide complex, offshore of South Kona, Hawaii Island, was investigated in 2001 by visual observation and in-situ sampling on the bench scarp and a megablock, during two dives utilizing the Remotely Operated Vehicle (ROV) Kaiko and its mother ship R/V Kairei. Topography of the bench scarp and megablocks were defined in 3-D perspective, using high-resolution digital bathymetric data acquired during the cruise. Compositions of 34 rock samples provide constraints on the landslide source regions and emplacement mechanisms. The bench scarp consists mainly of highly fractured, vesiculated, and oxidized aa lavas that slumped from the subaerial flank of ancestral Mauna Loa. The megablock contains three units: block facies, matrix facies, and draped sediment. The block facies contains hyaloclastite interbedded with massive lava, which slid from the shallow submarine flank of ancestral Mauna Loa, as indicated by glassy groundmass of the hyaloclastite, low oxidation state, and low sulfur content. The matrix facies, which directly overlies the block facies and is similar to a lahar deposit, is thought to have been deposited from the water column immediately after the South Kona slide event. The draped sediment is a thin high-density turbidite layer that may be a distal facies of the Alika-2 debris-avalanche deposit; its composition overlaps with rocks from subaerial Mauna Loa. The deposits generated by the South Kona slide vary from debris avalanche deposit to turbidite. Spatial distribution of the deposits is consistent with deposits related to large landslides adjacent to other Hawaiian volcanoes and the Canary Islands.     

A technique for modeling the Earth’s gravity field on the basis of satellite accelerations

A technique is proposed for Earths gravity field modeling on the basis of satellite accelerations that are derived from precise orbit data. The functional model rests on Newtons second law. The computational procedure is based on the pre-conditioned conjugate-gradient (PCCG) method. The data are treated as weighted average accelerations rather than as point-wise ones. As a result, a simple three-point numerical differentiation scheme can be used to derive them. Noise in the orbit-derived accelerations is strongly dependent on frequency. Therefore, the key element of the proposed technique is frequency-dependent data weighting. Fast convergence of the PCCG procedure is ensured by a block-diagonal pre-conditioner (approximation of the normal matrix), which is derived under the so-called Colombo assumptions. Both uninterrupted data sets and data with gaps can be handled. The developed technique is compared with other approaches: (1) the energy balance approach (based on the energy conservation law) and (2) the traditional approach (based on the integration of variational equations). Theoretical considerations, supported by a numerical study, show that the proposed technique is more accurate than the energy balance approach and leads to approximately the same results as the traditional one. The former finding is explained by the fact that the energy balance approach is only sensitive to the along-track force component. Information about the cross-track and the radial component of the gravitational potential gradient is lost because the corresponding force components do no work and do not contribute to the energy balance. Furthermore, it is shown that the proposed technique is much (possibly, orders of magnitude) faster than the traditional one because it does not require the computation of the normal matrix. Hints are given on how the proposed technique can be adapted to the explicit assembling of the normal matrix if the latter is needed for the computation of the model covariance matrix.Acknowledgments. Professor R. Klees is thanked for support of the project and for numerous fruitful discussions. The authors are also thankful to Dr. J. Kusche for useful remarks and to Dr. E. Schrama, his solid background in satellite geodesy proved to be very helpful. A large number of valuable comments were made by Dr. S.-C. Han, Dr. P. Schwintzer, and an anonymous reviewer; their contribution is greatly acknowledged. The satellite orbits used in the numerical study were kindly provided by Dr. P. Visser (Aerospace Department, Delft University of Technology). Access to the SGI Origin 3800 computer was provided by Stichting Nationale Computerfaciliteiten (NCF), grant SG-027.     

Occurrence and transport of Irgarol 1051 and its major metabolite in coastal waters from South Florida

Irgarol 1051, a boosting antifouling agent often used to supplement copper based paints was found in surface waters from South Florida at stations collected from the Miami River, Biscayne Bay and selected areas of the Florida Keys. Concentrations of the herbicide ranged from below the method detection limit (1 ng/L) to as high as 182 ng/L in a canal system in Key Largo. The herbicide was present at 93% of the stations and often found in conjunction with its descyclopropyl metabolite (M1) previously reported to be the major degradation product of Irgarol under natural environmental conditions. The 90th percentile concentration calculated for all South Florida samples was 57.6 ng/L. Based on available data on the toxicity of Irgarol to algae and coral, only two stations (approximately 3%) ranked above the LC50 of 136 ng/L reported for the marine algae Naviculla pelliculosa and above the 100 ng/L level reported to reversibly inhibit photosynthesis of intact corals. However, a basic dissipation model for Irgarol using the Key Largo Harbor station as a point source indicated that concentrations of the herbicide decreased rapidly and concentrations below the MDL are observed within 2000 m of the source. No major coral based benthic habitats are documented for all the stations surveyed at distances that Irgarol may pose a substantial risk. However, other types of submerged vegetation like seagrasses are common around the marinas and the effects of Irgarol to such endpoints should be investigated further.     

A simplified method for simulation of strong ground motion using finite rupture model of the earthquake source

We present a simplified method to simulate strong ground motion for a realistic representation of a finite earthquake source burried in a layered earth. This method is based on the stochastic simulation method of Boore (Boore, D. M., 1983, Bull. Seism. Soc. Am. 73, 1865–1894) and the Empirical Greens Function (EFG) method of Irikura (Irikura, K., 1986, Proceedings of the 7th Japan Earthquake symposium, pp. 151–156). The rupture responsible for an earthquake is represented by several subfaults. The geometry of subfaults and their number is decided by the similarity relationships. For simulation of ground motion using the stochastic simulation technique we used the shapping window based on the kinetic source model of the rupture plane. The shaping window deepens on the geometry of the earthquake source and the propagation characteristics of the energy released by various subfaults. The division of large fault into small subfaults and the method for accounting their contribution at the surface is identical to the EGF. The shapping window has been modified to take into account the effect of the transmission of energy released form the finite fault at various boundaries of the layered earth model above the source. In the present method we have applied the correction factor to adjust slip time function of small and large earthquakes. The correction factor is used to simulate strong motion records having basic spectral shape of ² source model in broad frequency range. To test this method we have used the strong motion data of the Geiyo earthquake of 24th March 2001, Japan recorded by KiK network. The source of this earthquake is modelled by a simple rectangular rupture of size 24 × 15 km, burried at a depth of 31 km in a multilayered earth model. This rupture plane is divided into 16 rectangular subfaults of size 6.0 × 3.75 km each. Strong motion records at eight selected near-field stations were simulated and compared with the observed records in terms of the acceleration and velocity records and their response spectrum. The comparison confirms the suitability of proposed rupture model responsible for this earthquake and the efficacy of the approach in predicting the strong motion scenario of earthquakes in the subduction zone. Using the same rupture model of the Geiyo earthquake, we compared the simulated records from our and the EGF techniques at one near-field station. The comparison shows that this technique gives records which matches in a wide frequency range and that too from simple and easily accessible parameters of burried rupture.     

Noril’sk–Talnakh Cu–Ni–PGE deposits: a revised tectonic model

The Norilsk mining district is located at the northwest margin of the Tunguska basin, in the centre of the 3,000×4,000 km Siberian continental flood basalt (CFB) province. This CFB province was formed at the Permo-Triassic boundary from a superplume that ascended into the geometric centre of the Laurasian continent, which was surrounded by subducting slabs of oceanic crust. We suggest that these slabs could have reached the core–mantle boundary, and they may have controlled the geometric focus of the superplume. The resulting voluminous magma intruded and erupted in continental rifts and related extensive flood basalt events over a 2–4 Ma period. Cu–Ni–PGE sulfide mineralization is found in olivine-bearing differentiated mafic intrusions beneath the flood basalts at the northwestern margin of the Siberian craton and also in the Taimyr Peninsula, some 300 km east of a triple junction of continental rifts, now buried beneath the Mesozoic–Cenozoic sedimentary basin of western Siberia. The Norilsk-I and Talnakh-Oktyabrsky deposits occur in the Norilsk–Kharaelakh trough of the Tunguska CFB basin. The Cu–Ni–PGE-bearing mineralized intrusions are 2–3 km-wide and 20 km-long differentiated chonoliths. Previous studies suggested that parts of the magma remained in intermediate-level crustal chambers where sulfide saturation and accumulation took place before emplacement. The 5–7-km-thick Neoproterozoic to Palaeozoic country rocks, containing sedimentary Cu mineralization and evaporites, may have contributed additional metal and sulfur to this magma. Classic tectonomagmatic models for these deposits proposed that subvertical crustal faults, such as the northeast-trending Norilsk–Kharaelakh fault, were major trough-parallel conduits providing access for magmas to the final chambers. However, geological maps of the Norilsk region show that the Norilsk–Kharaelakh fault offsets the mineralization, which was deformed into folds and offset by related reverse faults, indicating compressional deformation after mineralization in the Late Triassic to Early Jurassic. In addition, most of the intrusions are sills, not dykes as should be expected if the vertical faults were major conduits. A revised tectonic model for the Norilsk region takes into account the fold structure and sill morphology of the dominant intrusions, indicating a lateral rather than vertical emplacement direction for the magma into final chambers. Taking into account the fold structure of the country rocks, the present distribution of the differentiated intrusions hosting the Norilsk-I and Talnakh–Oktyabrsky deposits may represent the remnants of a single, >60 km long, deformed and eroded palm-shaped cluster of mineralized intrusions, which are perceived as separate intrusions at the present erosional level. The original direction of sill emplacement may have been controlled by a northeast-trending paleo-rise, which we suggest is present at the southeastern border of the Norilsk–Kharaelakh trough based on analysis of the unconformity at the base of the CFB. The mineralized intrusions extend along this rise, which we interpret as a structure that formed above the extensionally tilted block in the metamorphic basement. Geophysical data indicate the presence of an intermediate magma chamber that could be linked with the Talnakh intrusion. In turn, this T-shaped flat chamber may link with the Yenisei–Khatanga rift along the northwest-trending Pyasina transform fault, which may have served as the principal magma conduit to the intermediate chamber. It then produced the differentiated mineralized intrusions that melted through the evaporites with in situ precipitation of massive, disseminated, and copper sulfide ore. The Norilsk–Kharaelakh crustal fault may not relate to mineralization and possibly formed in response to late Mesozoic spreading in the Arctic Ocean.Editorial handling: P. Lightfoot