Bathymetric inversion of South China Sea from satellite altimetry data

1　ＩｎｔｒｏｄｕｃｔｉｏｎＴｏｄｅｖｅｌｏｐｔｈｅｏｃｅａｎｗｉｄｅｌｙａｎｄｄｅｅｐｌｙ ,ｗｅｎｅｅｄａｂｕｎｄａｎｔｏｃｅａｎｉｎｆｏｒｍａｔｉｏｎ .Ａｓａｎｅｓｓｅｎｔｉａｌｐａｒｔｏｆｓｕｃｈｉｎｆｏｒｍａｔｉｏｎ ,ｓｅａｆｌｏｏｒｔｏｐｏｇｒａｐｈｙｐｌａｙｓａｖｅｒｙｉｍｐｏｒｔａｎｔｒｏｌｅｉｎａｖａｒｉｅｔｙｏｆｍａｒｉｎｅａｃｔｉｖｉｔｉｅｓ .Ｈｏｗｅｖｅｒ,ｔｈｅｈｉｇｈｃｏｓｔｆｏｒｏｃｅａｎｂａｔｈｙｍｅｔｒｉｃｓｕｒｖｅｙｉｎｇｌｉｍｉｔｓｔｈｅａｐｐｌｉｃａｔｉｏｎｏ…     

SAINF: A toolbox for analyzing the effect of point-like,line-like and polygon-like infrastructural features on the distribution of point-like non-infrastructural features

This paper introduces a GIS-based toolbox, called SAINF, for analyzing the effects of infrastructural features on the distribution of non-infrastructural features. One of the distinct characteristics of SAINF is that it can deal not only with point-like infrastructural features but also line-like (e.g. railways) and polygon-like (e.g. big parks) features. SAINF analyzes these effects with the conditional nearest neighbor distance method and the cross K-function method. First, the paper briefly refers to these methods. Second, the paper illustrates how to use SAINF with an actual example. SAINF can be downloaded via the Internet without charge for non-profit uses.We express our thanks to Exceed Co. Ltd. for helping us program SAINF, and to Miki Arimoto for making the data of high-class apartment buildings and parks in Kiba. This development was partly supported by Grant-in-aid for Scientific Research No. 10202201 (Spatial information science for human and social sciences) and No. 14350327 (Regional supporting network planning for the elderly) by the Ministry of Education, Culture, Sports, Science and Technology of Japan.     

Statistical Downscaling Based on Dynamically DownscaledPredictors： Application to Monthly Precipitation in Sweden

A prerequisite of a successful statistical downscaling is that large-scale predictors simulated by the General Circulation Model (GCM) must be realistic. It is assumed here that features smaller than the GCM resolution are important in determining the realism of the large-scale predictors. It is tested whether a three-step method can improve conventional one-step statistical downscaling. The method uses predictors that are upscaled from a dynamical downscaling instead of predictors taken directly from a GCM simulation. The method is applied to downscaling of monthly precipitation in Sweden. The statistical model used is a multiple regression model that uses indices of large-scale atmospheric circulation and 850-hPa specific humidity as predictors. Data from two GCMs (HadCM2 and ECHAM4) and two RCM experiments of the Rossby Centre model (RCA1) driven by the GCMs are used. It is found that upscaled RCA1 predictors capture the seasonal cycle better than those from the GCMs, and hence increase the reliability of the downscaled precipitation. However, there are only slight improvements in the simulation of the seasonal cycle of downscaled precipitation. Due to the cost of the method and the limited improvements in the downscaling results, the three-step method is not justified to replace the one-step method for downscaling of Swedish precipitation.     

Dip and anisotropy effects on flow using a vertically skewed model grid

Darcy flow equations relating vertical and bedding-parallel flow to vertical and bedding-parallel gradient components are derived for a skewed Cartesian grid in a vertical plane, correcting for structural dip given the principal hydraulic conductivities in bedding-parallel and bedding-orthogonal directions. Incorrect-minus-correct flow error results are presented for ranges of structural dip (0 < or = theta < or = 90) and gradient directions (0 < or = phi < or = 360). The equations can be coded into ground water models (e.g., MODFLOW) that can use a skewed Cartesian coordinate system to simulate flow in structural terrain with deformed bedding planes. Models modified with these equations will require input arrays of strike and dip, and a solver that can handle off-diagonal hydraulic conductivity terms.     

Surveillance for previously unmonitored organic contaminants in the San Francisco Estuary

The San Francisco Estuary Regional Monitoring Program initiated surveillance monitoring to identify previously unmonitored synthetic organic contaminants in the San Francisco Estuary. Organic extracts of water samples were analyzed using gas chromatography-mass spectrometry in full scan mode. The major contaminant classes identified in the samples were fire retardants, pesticides, personal care product ingredients, and plasticizers. Evidence from the literature suggests that some of these contaminants can persist in the environment, induce toxicity, and accumulate in marine biota and in higher food chain consumers. The major sources of these contaminants into the marine environment are the discharge of municipal and industrial wastewater effluents, urban stormwater, and agricultural runoff. As a proactive effort, it is suggested that surveillance studies be used routinely in monitoring programs to identify and prevent potential problem contaminants from harming the marine environment.     

Volumetric characteristics of lava flows from interferometric radar and multispectral satellite data: the 1995 Fernandina and 1998 Cerro Azul eruptions in the western Galápagos

We have used a suite of remotely sensed data, numerical lava flow modeling, and field observations to determine quantitative characteristics of the 1995 Fernandina and 1998 Cerro Azul eruptions in the western Galápagos Islands. Flank lava flow areas, volumes, instantaneous effusion rates, and average effusion rates were all determined for these two eruptions, for which only limited syn-eruptive field observations are available. Using data from SPOT, TOPSAR, ERS-1, and ERS-2, we determined that the 1995 Fernandina flow covers a subaerial area of 6.5×10⁶ m² and has a subaerial dense rock equivalent (DRE) volume of 42×10⁶ m³. Field observations, ATSR satellite data, and the FLOWGO numerical model allow us to determine that the effusion rate declined exponentially from a high of ~60–200 m³ s^-1 during the first few hours to <5 m³ s^-1 prior to ceasing after 73 days, with a mean effusion rate of 4–16 m³ s^-1. Integrating the ATSR-derived, exponentially declining effusion rate over the eruption duration produces a total (subaerial + submarine) DRE volume of between 27 and 100×10⁶ m³, the range in values being due to differing assumptions about heat loss characteristics; only values in the higher part of this range are consistent with the independently derived subaerial volume. Using SPOT, TOPSAR, ERS-1, and ERS-2 data, we determine that the 1998 Cerro Azul flow is 16 km long, covers 16 km², and has a DRE volume of 54×10⁶ m³. FLOWGO produces at-vent velocity and effusion rate values of 11 m s^-1 and ~600 m³ s^-1, respectively. The velocity value agrees well with the 12 m s^-1 estimated in the field. The mean effusion rate (total DRE volume/duration) was 7–47 m³ s^-1. Dike dimensions, fissure lengths, and pressure gradients along the conduit based on magma chamber depth estimates of 3–5 km produce mean effusion rates for the two eruptions that range over nearly four orders of magnitude, the range being due to uncertainty in the magma viscosity, dike dimensions, and pressure gradient between magma chamber and vent. Although somewhat consistent with mean effusion rates from other techniques, their wide range makes them less useful. The exponentially declining effusion rates during both eruptions are consistent with release of elastic strain being the driving mechanism of the eruptions. Our results provide independent input parameters for previously published theoretical relationships between magma chamber pressurization and eruption rates that constrain chamber volumes and increases in volume prior to eruption, as well as time constants of exponential decay during the eruption. The results and theoretical relationships combine to indicate that at both volcanoes probably 25–30% of the volumetric increase in the magma chamber erupted as lava onto the surface. In both eruptions the lava flow volumes are less than 1% of the magma chamber volume.     

Was a dacite dome implicated in the 9,500 <Emphasis Type="SmallCaps">b.p.</Emphasis> collapse of Mt Ruapehu? A palaeomagnetic investigation

Constraining the process by which volcanoes become unstable is difficult. Several models have been proposed to explain the driving forces which cause volcanic edifices to catastrophically collapse. These include models for destabilisation of volcanic flanks by wedging due to dyke intrusion and the weakening of mechanical properties by pressurisation of pore fluids. It is not known which, if any, of the models are relevant to particular sector collapse events. Recent developments in the palaeomagnetic estimation of emplacement temperatures of volcaniclastic rocks have shown that even relatively low emplacement temperatures can be recorded by volcaniclastics with high fidelity. We have carried out a palaeomagnetic study of emplacement temperatures to investigate the role of igneous activity in the initiation of the 9,500 b.p. Murimotu sector collapse of Mt Ruapehu, New Zealand. This debris avalanche deposit has three fades which are stratigraphically superimposed, and the lowermost fades contains three lithological assemblages representing different segments of the edifice which were transported with little internal mixing within the flow. We have determined that some of the dacite-bearing assemblage 1, fades 1 was hot (∼350 °C) during transport and emplacement, whereas none of the other lithological assemblages of fades contained hot material. Our interpretation is that a dacite dome was active on the ancient Ruapehu edifice immediately prior to the Murimotu sector collapse. The partially cooled carapace of the dome and material shed from this part was incorporated into the avalanche deposit, along with cold lavas and volcaniclastics. We have not found evidence for incorporation of material at or close to magmatic temperatures, at least in the sampled locations. Our palaeomagnetic work allows us to develop a comprehensive, new palaeomagnetic classification of volcaniclastics. Published online: 25 January 2003 Editorial responsibility: D. Dingwell     

Seismicity and crustal deformation preceding the January 1996 eruptions at Karymsky Volcanic Center,Kamchatka

Two explosive eruptions occurred on 2 January 1996 at Karymsky Volcanic Center (KVC) in Kamchatka, Russia: the first, dacitic, from the central vent of Karymsky volcano, and the second, several hours later, from Karymskoye lake in the caldera of Akademia Nauk volcano. The main significance of the 1996 volcanic events in KVC was the phreatomagmatic eruption in Karymskoye lake, which was the first eruption in this lake in historical time, and was a basaltic eruption at the acidic volcanic center. The volcanic events were associated with the 1 January Ms 6.7 (Mw 7.1) earthquake that occurred at a distance of about 9–17 km southeast from the volcanoes just before the eruptions. We study the long-term (1972–1995) and short-term (1–2 January 1996) characteristics of crustal deformations and seismicity before the double eruptive event in KVC. The 1972–1995 crustal deformation was homogeneous and characterized by a gradual extension with a steady velocity. The seismic activity in 1972–1995 developed at the depth interval from 0 to 20 km below the Akademia Nauk volcano and spread to the southeast along a regional fault. The seismic activity in January 1996 began with a short sequence of very shallow microearthquakes (M ~0) beneath Karymsky volcano. Then seismic events sharply increased in magnitude (up to mb 4.9) and moved along the regional fault to the southeast, culminating in the Ms 6.7 earthquake. Its aftershocks were located to the southeast and northwest from the main shock, filling the space between the two active volcanoes and the ancient basaltic volcano of Zhupanovsky Vostryaki. The eruption in Karymskoye lake began during the aftershock sequence. We consider that the Ms 6.7 earthquake opened the passageway for basic magma located below Zhupanovsky Vostryaki volcano that fed the eruption in Karymskoye lake.