Hazard assessment of volcanic ballistic impacts at Mt Chihshin,Tatun Volcano Group,northern Taiwan

The assessment of vulnerability provides valuable knowledge in the risk assessment steps of a risk governance process. Given the multiscale, multilevel, and multisectorial aspects of flood risk, the diversified entities that directly and indirectly intervene in risk management require specific outputs from the assessment studies. Urban areas in estuarine margins are particularly exposed and vulnerable to flooding. Such interface conditions are found in the Old City Centre of the Seixal, located in the Tagus estuary, Portugal. Here, two distinct methodologies were applied for the assessment of territorial vulnerability. A regional, lower-scale, methodology explores the application of the statistical procedure based on the SoVI^® at the statistical block level. A second, local and higher-scale, methodology is based in data collected through field matrices at the building and statistical sub-block level. Comparison of results revealed that the lower-scale assessment provides information on the vulnerability drivers at the regional and municipal level. Nevertheless, only at a higher-scale, it is possible to characterize and differentiate the smaller geographical units of analysis that compose the Old City Centre of Seixal. The lower-scale vulnerability assessment allows a strategic response, based on adaptation measures such as spatial planning, institutional capacity building and public awareness. The local level assessment provides more accurate knowledge to support local emergency planning and the allocation of operational and material resources at the urban level. Nevertheless, rather than antagonistic, both models can be considered as complementary, having in mind the requirements of an holistic flood risk governance model.     

Nonlinear Time Series Analysis of Volcanic Tremor Events Recorded at Sangay Volcano,Ecuador

The assumption that volcanic tremor may be generated by deterministic nonlinear source processes is now supported by a number of studies at different volcanoes worldwide that clearly demonstrate the low-dimensional nature of the phenomenon. We applied methods based on the theory of nonlinear dynamics to volcanic tremor events recorded at Sangay volcano, Ecuador in order to obtain more information regarding the physics of their source mechanism. The data were acquired during 21–26 April 1998 and were recorded using a sampling interval of 125 samples s^–1 by two broadband seismometers installed near the active vent of the volcano. In a previous study Johnson and Lees (2000) classified the signals into three groups: (1) short duration (<1 min) impulses generated by degassing explosions at the vent; (2) extended degassing chugging events with a duration 2–5 min containing well-defined integer overtones (1–5 Hz) and variable higher frequency content; (3) extended degassing events that contain significant energy above 5 Hz. We selected 12 events from groups 2 and 3 for our analysis that had a duration of at least 90 s and high signal-to-noise ratios. The phase space, which describes the evolution of the behavior of a nonlinear system, was reconstructed using the delay embedding theorem suggested by Takens. The delay time used for the reconstruction was chosen after examining the first zero crossing of the autocorrelation function and the first minimum of the Average Mutual Information (AMI) of the data. In most cases it was found that both methods yielded a delay time of 14–18 samples (0.112–0.144 s) for group 2 and 5 samples (0.04 s) for group 3 events. The sufficient embedding dimension was estimated using the false nearest neighbors method which had a value of 4 for events in group 2 and was in the range 5–7 for events in group 3. Based on these embedding parameters it was possible to calculate the correlation dimension of the resulting attractor, as well as the average divergence rate of nearby orbits given by the largest Lyapunov exponent. Events in group 2 exhibited lower values of both the correlation dimension (1.8–2.6) and largest Lyapunov exponent (0.013–0.022) in comparison with the events in group 3 where the values of these quantities were in the range 2.4–3.5 and 0.029–0.043, respectively. Theoretically, a nonlinear oscillation described by the equation ++g(x)=fcost can generate deterministic signals with characteristics similar to those observed in groups 2 and 3 as the values of the parameters ,,f, are drifting, causing instability of orbits in the phase space.     

Three-dimensional shallow velocity structure beneath Taal Volcano,Philippines

Based on its numerous historical explosive eruptions and high potential hazards to nearby population of millions, Taal Volcano is one of the most dangerous “Decade Volcanoes” in the world. To provide better investigation on local seismicity and seismic structure beneath Taal Volcano, we deployed a temporary seismic network consisting of eight stations from March 2008 to March 2010. In the preliminary data processing stage, three periods showing linear time-drifting of internal clock were clearly identified from noise-derived empirical Green’s functions. The time-drifting errors were corrected prior to further data analyses. By using VELEST, 2274 local earthquakes were manually picked and located. Two major earthquake groups are noticed, with one lying beneath the western shore of Taal Lake showing a linear feature, and the other spreading around the eastern flank of Taal Volcano Island at shallower depths. We performed seismic tomography to image the 3D structure beneath Taal Volcano using the LOTOS algorithm. Some interesting features are revealed from the tomographic results, including a solidified magma conduit below the northwestern corner of Taal Volcano Island, indicated by high Vp, Vs, and low Vp/Vs ratio, and a large potential hydrothermal reservoir beneath the center of Taal Volcano Island, suggested by low Vs and high Vp/Vs ratio. Furthermore, combining earthquake distributions and tomographic images, we suggest potential existence of a hydrothermal reservoir beneath the southwestern corner of Taal Lake, and a fluid conduit extending to the northwest. These seismic features have never been proposed in previous studies, implying that new hydrothermal activity might be formed in places away from the historical craters on Taal Volcano Island.     

Empirical relationships between aftershock area dimensions and magnitude for earthquakes in the Mediterranean Sea region

Fault dimension estimates derived from the aftershock area extent of 36 shallow depth (≤ 31 km) earthquakes that occurred in the Mediterranean Sea region have been used in order to establish empirical relationships between length, width, area and surface-wave/moment magnitude. This dataset consists of events whose aftershock sequence was recorded by a dense local or regional network and the reported location errors did not exceed on average 3–5 km. Surface-wave magnitudes for these events were obtained from the NEIC database and/or published reports, while moment magnitudes as well as focal mechanisms were available from the Harvard/USGS catalogues. Contrary to the results of some previously published studies we found no evidence in our dataset that faulting type may have an effect on the fault dimension estimates and therefore we derived relationships for the whole of the dataset. Comparisons, by means of statistical F-tests, of our relationships with other previously published regional and global relationships were performed in order to check possible similarities or differences. Most such comparisons showed relatively low significance levels (< 95%), since the differences in source dimension estimates were large mainly for magnitudes lower than 6.5, becoming smaller with increasing magnitude. Some degree of similarity, however, could be observed between our fault length relationship and the one derived from aftershock area lengths of events in Greece, while a difference was found between our regional and global fault length relationships. A calculation of the ratio defined as the fault length, derived from our relationships, to the length estimated from regional empirical relationships involving surface ruptures showed that it can take a maximum value of about 7 for small magnitudes while it approaches unity at Ms 7.2. When calculating the same ratio using instead global empirical relationships we see the maximum value not exceeding 1.8, while unity is reached at Mw 7.8, indicating the existence of a strong regional variation in the fault lengths of earthquakes occurring in the Mediterranean Sea region. Also, a relationship between the logarithms of the rupture area and seismic moment is established and it is inferred that there is some variation of stress drop as a function of seismic moment. In particular, it is observed that for magnitudes lower than 6.6 the stress drop fluctuates around 10 bar, while for larger magnitudes the stress drop reaches a value as high as 60 bar.     

Calcite twinning strain variations across the Proterozoic Grenville orogen and Keweenaw-Kapuskasing inverted foreland,USA and Canada

We report the calcite twinning strain results of a traverse across the Grenville orogen from Parry Sound, Ontario (NW) to Ft. Ann, New York (SE), including the younger, adjacent Ordovician Taconic allochthon. Fifty four carbonates (marbles, calcite veins, Ordovician limestone) were collected resulting in 68 strain analyses on mechanically twinned calcite (n = 2337 grains) across the Central Gneiss Belt (CGB; 3 samples), the Central Metasedimentary Belt (CMB; 27 samples), the Central Granulite Terrane (CGT; Adirondack's; 13 samples) and the Ottawan Orogenic Lid (OOL; 11 samples). Twinning strains in the greenschist-grade OOL marbles preserve N–S shortening and U-Pb titanite ages (～1150 Ma; n = 4) document these marbles formed during the Shawinigan (1190–1140 Ma) part of the Grenville orogen. From northwest to southeast, the Ottawan (1095–1020 Ma) twinning strain is dominantly a layer-parallel shortening fabric oriented N–S (Parry Sound), then becomes parallel to the Grenville thrust direction (NW–SE) across the CMB to the Adirondack Highlands where the sub-horizontal shortening strain becomes margin-parallel (SW–NE). Within the regional sample suite there are two areas studied in detail, the Bancroft shear zone (n = 11) and a roadcut on the southeast side of the Adirondack Mountains (Ft. Ann, NY; n = 8). Marbles from the Bancroft shear zone contain calcite grains with 2 sets of twin lamellae (e₁ and e₂). The better-developed e₁ sets (n = 406) record a horizontal fabric oriented NW–SE whereas the younger e₂ lamellae (n = 146) preserve a margin-parallel (SW–NE) horizontal fabric. Both the e₁ and e₂ strains record an overprint vertical shortening strain (NEV), perhaps related to extensional orogenic collapse. We also report an Ottawan orogen-aged granoblastic mylonite (1093 Ma, U-Pb zircon; 1102 Ma Ar-Ar biotite) in the Keweenaw thrust hanging wall 500 km inboard of the Grenville front and interpret the relations of Grenville-Keweenaw far-field dynamics.     

Use of microfluidic experiments to optimize MICP treatment protocols for effective strength enhancement of MICP-treated sandy soils

Microbially induced calcium carbonate (CaCO₃) precipitation (MICP) has been extensively studied for soil improvement in geotechnical engineering. The quantity and size of calcium carbonate crystals affect the strength of MICP-treated soil. In this study, microfluidic chip experiments and soil column experiments were conducted to optimize MICP treatment protocols for effective strength enhancement of MICP-treated sandy soils. The microscale experiments reveal that, due to Ostwald ripening, longer injection intervals allow crystals to dissolve and reprecipitate into larger crystals regardless of the concentration of cementation solution. Even though a cementation solution input rate of 0.042 mol/l/h is sufficient to maintain a high chemical transformation efficiency, a further reduction in the input rate by about four times resulted in an increase in the size of crystals produced by the end of treatment from about 40 to 60 μm. These findings were applied in soil column experiments. Results showed that significantly larger crystals and higher soil strength were achieved when the normalized rate of cementation solution injection was reduced from 0.042 to 0.021 mol/l/h. Crystal size and soil strength increased slightly more when the normalized input rate was further reduced from 0.021 to 0.010 mol/l/h. This study demonstrates how data from microscale microfluidic experiments that examine the effects of injection intervals and concentration of cementation solution on the properties of calcium carbonate crystals can be used to optimize MICP treatment in macroscale sand soil column experiments for effective strength enhancement.

     

Soil stabilization with microbially induced desaturation and precipitation (MIDP) by denitrification: a field study

Acta Geotechnica - A field trial was performed to test the feasibility of microbially induced desaturation and precipitation (MIDP) by denitrification for the stabilization of a silty soil....     

Seismicity characteristics of a potentially active Quaternary volcano: The Tatun Volcano Group, northern Taiwan

The Tatun Volcano Group (TVG) is located at the northern tip of Taiwan, near the capital Taipei and close to two nuclear power plants. Because of lack of any activity in historical times it has been classified as an extinct volcano, even though more recent studies suggest that TVG might have been active during the last 20 ka. In May 2003 a seismic monitoring project at the TVG area was initiated by deploying eight three-component seismic stations some of them equipped with both short-period and broadband sensors. During the 18 months observation period local seismicity mainly consisted of high frequency earthquakes either occurring as isolated events, or as a continuous sequence in the form of spasmodic bursts. Mixed and low frequency events were also present during the same period, even though they occurred only rarely. Arrival times from events with clear P-/S-wave phases were inverted in order to obtain a minimum 1D velocity model with station corrections. Probabilistic nonlinear earthquake locations were calculated for all these events using the newly derived velocity model. Most high frequency seismicity appeared to be concentrated near the areas of hydrothermal activity, forming tight clusters at depths shallower than 4 km. Relative locations, calculated using the double-difference method and utilising catalogue and cross-correlation differential traveltimes, showed insignificant differences when compared to the nonlinear probabilistic locations. In general, seismicity in the TVG area seems to be primarily driven by circulation of hydrothermal fluids as indicated by the occurrence of spasmodic bursts, mixed/low frequency events and a b-value (1.17 ± 0.1) higher than in any other part of Taiwan. These observations, that are similar to those reported in other dormant Quaternary volcanoes, indicate that a magma chamber may still exist beneath TVG and that a future eruption or period of unrest should not be considered unlikely.