Soil liquefaction during the Arequipa Mw 8.4, June 23, 2001 earthquake, southern coastal Peru

The Arequipa June 23, 2001, earthquake with a moment magnitude of Mw 8.4 struck southern Peru, northern Chile and western Bolivia. This shallow (29 km deep) interplate event, occurring in the coupled zone of the Nazca subduction next to the southeast of the subducting Nazca ridge, triggered very localized but widely outspread soil liquefaction. Although sand blows and lateral spreading of river banks and road bridge abutments were observed 390 km away from the epicenter in the southeast direction (nearing the town of Tacna, close to the Chile border), liquefaction features were only observed in major river valleys and delta and coastal plains in the meizoseismal area. This was strongly controlled by the aridity along the coastal strip of Southern Peru. From the sand blow distribution along the coastal area, a first relationship of isolated sand blow diameter versus epicentral distance for a single event is ever proposed. The most significant outcome from this liquefaction field reconnaissance is that energy propagation during the main June 23, 2001, event is further supported by the distribution and size of the isolated sand blows in the meizoseismal area. The sand blows are larger to the southeast of the epicenter than its northwestern equivalents. This can be stated in other words as well. The area affected by liquefaction to the northwest is less spread out than to the southeast. Implications of these results in future paleoliquefaction investigations for earthquake magnitude and epicentral determinations are extremely important. In cases of highly asymmetrical distribution of liquefaction features such as this one, where rupture propagation tends to be mono-directional, it can be reliably determined an epicentral distance (between earthquake and liquefaction evidence) and an earthquake magnitude only if the largest sand blow is found. Therefore, magnitude estimation using this uneven liquefaction occurrence will surely lead to underrating if only the shortest side of the meizoseismal area is unluckily studied, which can eventually be the only part exhibiting liquefaction evidence, depending on the earthquake location and the distribution of liquefaction-prone environments.     

Analysis of earth dams affected by the 2001 Bhuj Earthquake

An earthquake of magnitude of 7.6 (M_w 7.6) occurred in Bhuj, India on January 26, 2001. This event inflicted damages of varying extents to a large number of small to moderate size multi-zone earth dams in the vicinity of the epicenter. Some of the distress was due to the liquefaction of saturated alluvium in foundation. Liquefaction was relatively localized for the majority of these dams because the earthquake struck in the middle of a prolonged dry season when the reservoirs behind these dams were nearly empty and shallow alluvium soils underneath the downstream portions of the dams were partly dry. Otherwise, liquefaction of foundation soils would have been more extensive and damage to these dams more significant. Six such dams have been examined in this paper. Four of these facilities, Chang, Shivlakha, Suvi, and Tapar were within the 50 km of epicenter region. These dams underwent free-field ground motion with peak ground accelerations between 0.28g to 0.52g. Of these Chang Dam underwent severe slumping, whereas Shivlakha, Suvi, and Tapar Dams were affected severely especially over the upstream sections. Fatehgadh Dam and Kaswati Dam were affected relatively less severely. Foundation conditions underneath these dams were first examined for assessing liquefaction potential. A limited amount of subsurface information available from investigations undertaken prior to the earthquake indicates that, although the foundation soils within the top 2.0 to 2.5 m underneath these dams were susceptible to liquefaction, Bhuj Earthquake did not trigger liquefaction because of lack of saturation of these layers underneath the downstream portions of these dams. These dams were then analyzed using a simple sliding block procedure using appropriate estimates of undrained soil strength parameters. The results of this analysis for these structures were found to be in general agreement with the observed deformation patterns.     

An interpretation and catalogue of paleoseismicity in Sweden

Paleoseismic data provide a long-term record of seismic activity to predict hazards for periods longer than one to a few centuries. In Sweden, the analysis reveals there was a drastic change in dominant seismic mode from a high to super-high deglacial mode to a low to moderately low mode in present and Late Holocene time. Paleoseismic criteria and characteristics include numerous different sources of information; viz. primary faults, bedrock deformation, sedimentary deformation, rock and sediment slides, liquefaction, sorting by shaking, tsunamis, differing geomorphic expressions, disordering and ordering of magnetic particles. By applying multiple criteria, it was possible to identify 44 paleoseismic events, including 23 events of estimated M 6–7, 12 events of M 7–8 and 6 events of M > 8. Varve-dating often allows a precision as to a single year, in one case even to the season of a year. The key for paleoseismic reconstruction and testing is the application of multiple criteria.     

Late Pleistocene and Holocene paleoseismology of an intraplate seismic zone in a large alluvial valley, the New Madrid seismic zone, Central USA

The New Madrid seismic zone (NMSZ) is an intraplate right-lateral strike-slip and thrust fault system contained mostly within the Mississippi Alluvial Valley. The most recent earthquake sequence in the zone occurred in 1811–1812 and had estimated moment magnitudes of 7–8 (e.g., [Johnston, A.C., 1996. Seismic moment assessment of stable continental earthquakes, Part 3: 1811–1812 New Madrid, 1886 Charleston, and 1755 Lisbon. Geophysical Journal International 126, 314–344; Johnston, A.C., Schweig III, E.S, 1996. The enigma of the New Madrid earthquakes of 1811–1812. Annual Reviews of Earth and Planetary Sciences 24, 339–384; Hough, S.E., Armbruster, J.G., Seeber, L., Hough, J.F., 2000. On the modified Mercalli intensities and magnitudes of the New Madrid earthquakes. Journal of Geophysical Research 105 (B10), 23,839–23,864; Tuttle, M.P., 2001. The use of liquefaction features in paleoseismology: Lessons learned in the New Madrid seismic zone, central United States. Journal of Seismology 5, 361–380]). Four earlier prehistoric earthquakes or earthquake sequences have been dated A.D. 1450 ± 150, 900 ± 100, 300 ± 200, and 2350 B.C. ± 200 years using paleoliquefaction features, particularly those associated with native American artifacts, and in some cases surface deformation ([Craven, J. A. 1995. Paleoseismology study in the New Madrid seismic zone using geological and archeological features to constrain ages of liquefaction deposits. M.S thesis, University of Memphis, Memphis, TN, U.S.A.; Tuttle, M.P., Lafferty III, R.H., Guccione, M.J., Schweig III, E.S., Lopinot, N., Cande, R., Dyer-Williams, K., Haynes, M., 1996. Use of archaeology to date liquefaction features and seismic events in the New Madrid seismic zone, central United States. Geoarchaeology 11, 451–480; Guccione, M.J., Mueller, K., Champion, J., Shepherd, S., Odhiambo, B., 2002b. Stream response to repeated co-seismic folding, Tiptonville dome, western Tennessee. Geomorphology 43(2002), 313–349; Tuttle, M.P., Schweig, E.S., Sims, J.D., Lafferty, R.H., Wolf, L.W., Haynes, M.L., 2002. The earthquake potential of the New Madrid seismic zone, Bulletin of the Seismological Society of America, v 92, n. 6, p. 2080–2089; Tuttle, M.P., Schweig III, E.S., Campbell, J., Thomas, P.M., Sims, J.D., Lafferty III, R.H., 2005. Evidence for New Madrid earthquakes in A.D. 300 and 2350 B.C. Seismological Research Letters 76, 489–501]). The two most recent prehistoric and the 2350 B.C. events were probably also earthquake sequences with approximately the same magnitude as the historic sequence.Surface deformation (faulting and folding) in an alluvial setting provides many examples of stream response to gradient changes that can also be used to date past earthquake events. Stream responses include changes in channel morphology, deviations in the channel path from the regional gradient, changes in the direction of flow, anomalous longitudinal profiles, and aggradation or incision of the channel ([Merritts, D., Hesterberg, T, 1994. Stream networks and long-term surface uplift in the New Madrid seismic zone. Science 265, 1081–1084.; Guccione, M.J., Mueller, K., Champion, J., Shepherd, S., Odhiambo, B., 2002b. Stream response to repeated co-seismic folding, Tiptonville dome, western Tennessee. Geomorphology 43 (2002), 313–349]). Uplift or depression of the floodplain affects the frequency of flooding and thus the thickness and style of vertical accretion or drowning of a meander scar to form a lake. Vegetation may experience trauma, mortality, and in some cases growth enhancement due to ground failure during the earthquake and hydrologic changes after the earthquake ([VanArdale, R.B., Stahle, D.W., Cleaveland, M.K., Guccione, M.J., 1998. Earthquake signals in tree-ring data from the New Madrid seismic zone and implications for paleoseismicity. Geology 26, 515–518]). Identification and dating these physical and biologic responses allows source areas to be identified and seismic events to be dated.Seven fault segments are recognized by microseismicity and geomorphology. Surface faulting has been recognized at three of these segments, Reelfoot fault, New Madrid North fault, and Bootheel fault. The Reelfoot fault is a compressive stepover along the strike-slip fault and has up to 11 m of surface relief ([Carlson, S.D., 2000. Formation and geomorphic history of Reelfoot Lake: insight into the New Madrid seismic zone. M.S. Thesis, University of Arkansas, Fayetteville, Arkansas, U.S.A]) deforming abandoned and active Mississippi River channels ([Guccione, M.J., Mueller, K., Champion, J., Shepherd, S., Odhiambo, B., 2002b. Stream response to repeated co-seismic folding, Tiptonville dome, western Tennessee. Geomorphology 43 (2002), 313–349]). The New Madrid North fault apparently has only strike-slip motion and is recognized by modern microseismicity, geomorphic anomalies, and sand cataclasis ([Baldwin, J.N., Barron A.D., Kelson, K.I., Harris, J.B., Cashman, S., 2002. Preliminary paleoseismic and geophysical investigation of the North Farrenburg lineament: primary tectonic deformation associated with the New Madrid North Fault?. Seismological Research Letters 73, 393–413]). The Bootheel fault, which is not identified by the modern microseismicity, is associated with extensive liquefaction and offset channels ([Guccione, M.J., Marple, R., Autin, W.J., 2005, Evidence for Holocene displacements on the Bootheel fault (lineament) in southeastern Missouri: Seismotectonic implications for the New Madrid region. Geological Society of America Bulletin 117, 319–333]). The fault has dominantly strike-slip motion but also has a vertical component of slip. Other recognized surface deformation includes relatively low-relief folding at Big Lake/Manila high ([Guccione, M.J., VanArdale, R.B., Hehr, L.H., 2000. Origin and age of the Manila high and associated Big Lake “Sunklands”, New Madrid seismic zone, northeastern Arkansas. Geological Society of America Bulletin 112, 579–590]) and Lake St. Francis/Marked Tree high ([Guccione, M.J., VanArsdale, R.B., 1995. Origin and age of the St. Francis Sunklands using drainage patterns and sedimentology. Final report submitted to the U. S. Geological Survey, Award Number 1434-93-G-2354, Washington D.C.]), both along the subsurface Blytheville arch. Deformation at each of the fault segments does not occur during each earthquake event, indicating that earthquake sources have varied throughout the Holocene.     

The 17 March 2005 Kuzulu landslide (Sivas, Turkey) and landslide-susceptibility map of its near vicinity

Landslides are common natural hazards in the seismically active North Anatolian Fault Zone of Turkey. Although seismic activity, heavy rainfall, channel incisions, and anthropogenic effects are commonly the main triggers of landslides, on March 17, 2005, a catastrophic large landslide in Sivas, northeastern of Turkey, the Kuzulu landslide, was triggered by snowmelt without any other precursor. The initial failure of the Kuzulu landslide was rotational. Following the rotational failure, the earth material in the zone of accumulation exhibited an extremely rapid flow caused by steep gradient and high water content. The Agnus Creek valley, where Kuzulu village is located, was filled by the earth-flow material and a landslide dam was formed on the upper part of Agnus Creek. The distance from the toe of the rotational failure down to the toe of the earth flow measured more than 1800 m, with about 12.5 million m³ of displaced earth material. The velocity of the Kuzulu landslide was extremely fast, approximately 6 m/s. The main purposes of this study are to describe the mechanism and the factors conditioning the Kuzulu landslide, to present its environmental impacts, and to produce landslide-susceptibility maps of the Kuzulu landslide area and its near vicinity. For this purpose, a detailed landslide inventory map was prepared and geology, slope, aspect, elevation, topographic-wetness index and stream-power index were considered as conditioning factors. During the susceptibility analyses, the conditional probability approach was used and a landslide-susceptibility map was produced. The landslide-susceptibility map will help decision makers in site selection and the site-planning process. The map may also be accepted as a basis for landslide risk-management studies to be applied in the study area.     

午后多重极光弧观测研究

本文利用2001-2003年南极中山站175天全天空摄像机观测,对午后多重极光弧的出现率及其与Kp指数的关系进行了统计分析,结果表明午后多重极光弧出现率呈一单峰分布,最大发生率出现在1445UT(1645MLT),其位置在1500MLT极光热点(1300-1700MLT)近夜侧的部分。与地磁活动指数Kp的相关统计分析表明,Kp值为2-3之间时多重极光弧有较大的出现率,这说明中等地磁活动情形下午后多重极光弧有较高的出现率。事件分析表明多重极光弧的强度变化与地磁Pc5脉动具有较高的相关性,并且有类似的频谱特征,这说明午后多重极光弧可能与同时出现的Pc5地磁脉动有关。     

Development of a soil moisture‐based distributed hydrologic model for determining hydrologically based critical source areas

A simple grid cell‐based distributed hydrologic model was developed to provide spatial information on hydrologic components for determining hydrologically based critical source areas. The model represents the critical process (soil moisture variation) to run‐off generation accounting for both local and global water balance. In this way, it simulates both infiltration excess run‐off and saturation excess run‐off. The model was tested by multisite and multivariable evaluation on the 50‐km² Little River Experimental Watershed I in Georgia, U.S. and 2 smaller nested subwatersheds. Water balance, hydrograph, and soil moisture were simulated and compared to observed data. For streamflow calibration, the daily Nash‐Sutcliffe coefficient was 0.78 at the watershed outlet and 0.56 and 0.75 at the 2 nested subwatersheds. For the validation period, the Nash‐Sutcliffe coefficients were 0.79 at the watershed outlet and 0.85 and 0.83 at the 2 subwatersheds. The per cent bias was less than 15% for all sites. For soil moisture, the model also predicted the rising and declining trends at 4 of the 5 measurement sites. The spatial distribution of surface run‐off simulated by the model was mainly controlled by local characteristics (precipitation, soil properties, and land cover) on dry days and by global watershed characteristics (relative position within the watershed and hydrologic connectivity) on wet days when saturation excess run‐off was simulated. The spatial details of run‐off generation and travel time along flow paths provided by the model are helpful for watershed managers to further identify critical source areas of non‐point source pollution and develop best management practices.     

Liquefaction potential evaluations by energy-based method and stress-based method for various ground motions

An energy-based liquefaction potential evaluation method (EBM) previously developed was applied to a uniform sand model shaken by seismic motions recorded at different sites during different magnitude earthquakes. It was also applied to actual liquefaction case histories in Urayasu city during the 2011 M9.0 Tohoku earthquake and in Tanno-cho during the 2003 M8.0 Tokachi-oki earthquake. In all these evaluations, the results were compared with those by the currently used stress-based method (SBM) under exactly the same seismic and geotechnical conditions. It was found that EBM yields similar results with SBM for several ground motions of recent earthquakes but has easier applicability without considering associated parameters. In Urayasu city, the two methods yielded nearly consistent results by using an appropriate coefficient in SBM for the M9.0 earthquake, though both overestimated the actual liquefaction performance, probably because effects of plasticity and aging on in situ liquefaction strength were not taken into account. In Tanno-cho, EBM could evaluate actual liquefaction performance due to a small-acceleration motion during a far-field large magnitude earthquake while SBM could not.