Determination of critical slip surfaces using mutative scale chaos optimization

Limit equilibrium is a common method used to analyze the stability of a slope, and minimization of the factor of safety or identification of critical slip surfaces is a classical geotechnical problem in the context of limit equilibrium methods for slope stability analyses. A mutative scale chaos optimization algorithm is employed in this study to locate the noncircular critical slip surface with Spencer’s method being employed to compute the factor of safety. Four examples from the literature—one homogeneous slope and three layered slopes—are employed to identify the efficiency and accuracy of this approach. Results indicate that the algorithm is flexible and that although it does not generally provide the minimum FS, it provides results that are close to the minimum, an improvement over other solutions proposed in the literature and with small relative errors with respect to other minimum factor of safety (FS) values reported in the literature.

(Table 1) Experimental shear strength values for IODP Exp315 and Exp316 holes

We investigate the mechanics of slope failures on the Nankai accretionary complex offshore Japan in the vicinity of a major out-of-sequence thrust fault (termed the "megasplay"). Incorporating laboratory-measured shear strength of slope sediments sampled during Integrated Ocean Drilling Project (IODP) Expeditions 315 and 316 with local seafloor slope angles from bathymetric data and constraints on in-situ effective stress conditions from drilling, we find that slopes in the study area are stable and submarine landslides are not expected to occur under static conditions. In order to assess the possibility of slope failure triggered by coseismic rupture of the megasplay fault, we use empirical relations for strong ground motion attenuation from earthquakes with Mw 6-9. We find that the slope sediments should be stable based on computations from one model, developed from a catalog of worldwide subduction zone earthquakes (Youngs et al., 1997, doi:10.1785/gssrl.68.1.58). However, using a different model developed primarily from a catalog of crustal earthquakes in Japan (Kanno et al., 2006, doi:10.1785/0120050138), we find that slopes should be unstable for earthquakes 8 <= Mw <= 9, and possibly unstable for events with 6 <= Mw < 8, depending on the proximity of rupture to the seafloor. Considering limitations of the models and geologic observations of slope failure recurrence, the true slope stability is likely to be in between the predictions of the two models, and we suggest that it may be modulated by long-term pore pressure fluctuations.

I-82 Prosser Landslide Investigation Benton City, WA

The scope of this technical report is to establish the mechanisms by which the eastbound lanes of Interstate 82 at mile post (MP) 91.9 near Benton City continue to deform. Within the Washington State Department of Transportation (WSDOT), the area is known as the Prosser Landslide and has been an ongoing concern since the 1980s. Results from previous technical investigations have been conflicted or inconclusive as to whether landslide movement persists beneath or through the shear key-buttress or that pavement distress is related to swelling of a clay-rich unit that underlies the slope and interstate. For this report, the following steps were taken. First, I conducted a desk review of archived reports, memos, data, and drill logs from the original construction of I-82 and previous geotechnical investigations commissioned by WSDOT. Findings of this desk review are reported in Part III. Second, WSDOT drillers drilled two new boreholes at the Prosser Landslide site above the buttress and instrumentation was installed within the boreholes. Borehole logs produced from the 2013 drilling can be found in Appendix A of this report. Material retrieved from the suspected failure zone during drilling was tested at the WSDOT Materials Lab by WSDOT personnel for its mechanical properties including Atterberg limits, grain-size analysis, and residual shear strength (Appendix B). Samples were also analyzed for mineral content using X -ray powder diffraction (XRD). These data and observations are reported in Part III and Appendix C. Finally, using drill logs produced by WSDOT from the latest drilling and from historic drilling campaigns, I constructed a 2-D geologic model of the landslide site. This model is the basis for slope stability analysis reported in Part IV and Appendix D. This study concludes that the deformation observed in the eastbound lanes of I-82 could be the result of continued landslide movement, despite previous remediation efforts.

The Mineralogy and Strength Characteristics of Selected Glaciolacustrine Clays in the Puget Sound Region

The combination of rainy climate, glaciolacustrine clays, and steep topography of the Puget Lowland creates slope stability issues for the regional population. Several glaciolacustrine deposits of laminated silt and clay of different ages contribute to the likelihood of slope failure. The glaciolacustrine deposits are generally wet, range in thickness from absent to >30m, and consist of laminated silt and clay with sand interbeds at the tops and bottoms, sandy laminae throughout the deposits, occasional dropstones and shear zones. The glaciolacustrine deposits destabilize slopes by 1) impeding groundwater flow percolating through overlying glacial outwash sediments, 2) having sandy laminae that lower strength by increasing pore pressure during wet seasons, and 3) increasing the potential for block-style failure because of secondary groundwater pathways such as laminae and vertical fractures. Eight clay samples from six known landslide deposits were analyzed in this study for their mineralogy, clay fraction and strength characteristics. The mineralogy was determined using X-ray Diffractometry (XRD) which revealed an identical mineralogic suite among all eight samples consisting of chlorite, illite and smectite. Nonclay minerals appearing in the X-ray diffractogram include amphibole and plagioclase after removal of abundant quartz grains. Hydrometer tests yielded clay-size fraction percentages of the samples ranging from 10% to 90%, and ring shear tests showed that the angle of residual shear resistance (phi_r) ranged from 11° to 31°. Atterberg limits of the samples were found to have liquid limits ranging from 33 to 83, with plastic limits ranging from 25 to 35 and plasticity indices ranging from 6 to 48. The results of the hydrometer and residual shear strength tests suggest that phi_r varies inversely with the clay-size fraction, but that this relationship was not consistent among all eight samples. The nature of the XRD analysis only revealed the identity of the clay minerals present in the samples, and provided no quantitative information. Thus, the extent to which the mineralogy influenced the strength variability among the samples cannot be determined given that the mineral assemblages are identical. Additional samples from different locations within each deposit along with quantitative compositional analyses would be necessary to properly account for the observed strength variability.

Shallow Landslide Risk Assessment of the Westerly Bluffs of the 7.5-minute Freeland Quadrangle, Whidbey Island, Washington

Landslides often occur on slopes rendered unstable by underlying geology, geomorphology, hydrology, weather-climate, slope modifications, or deforestation. Unfortunately, humans commonly exacerbate such unstable conditions through careless or imprudent development practices. Due to local geology, geography, and climatic conditions, Puget Sound of western Washington State is especially landslide-prone. Despite this known issue, detailed analyses of landslide risks for specific communities are few. This study aims to classify areas of high landslide risk on the westerly bluffs of the 7.5 minute Freeland quadrangle based on a combined approach: mapping using LiDAR imagery and the Landform Remote Identification Model (LRIM) to identify landslides, and implementation of the Shallow Slope Stability Model (SHALSTAB) to establish a landslide exceedance probability. The objective is to produce a risk assessment from two shallow landslide scenarios: (1) minimum bluff setback and runout and (2) maximum bluff setback and runout. A simple risk equation that takes into account the probability of hazard occurrence with physical and economic vulnerability (van Westen, 2004) was applied to both scenarios. Results indicate an possible total loss as much as $32.6b from shallow landslides, given a setback of 12 m and a runout of 235 m.

A comparison of landslide inventories along the shoreline of Steamboat Island Peninsula, Thurston County, WA: office based protocol vs. field based protocol

The mountain ranges and coastlines of Washington State have steep slopes, and they are susceptible to landslides triggered by intense rainstorms, rapid snow melts, earthquakes, and rivers and waves removing slope stability. Over a 30-year timespan (1984-2014 and includes State Route (SR) 530), a total of 28 deep-seated landslides caused 300 million dollars of damage and 45 deaths (DGER, 2015). During that same timeframe, ten storm events triggered shallow landslides and debris flows across the state, resulting in nine deaths (DGER, 2015). The loss of 43 people, due to the SR 530 complex reactivating and moving at a rate and distance unexpected to residents, highlighted the need for an inventory of the stateís landslides. With only 13% of the state mapped (Lombardo et al., 2015), the intention of this statewide inventory is to communicate hazards to citizens and decision makers. In order to compile an accurate and consistent landslide inventory, Washington needs to adopt a graphic information system (GIS) based mapping protocol. A mapping protocol provides consistency for measuring and recording information about landslides, including such information as the type of landslide, the material involved, and the size of the movement. The state of Oregon shares similar landslide problems as Washington, and it created a GIS-based mapping protocol designed to inform its residents, while also saving money and reducing costly hours in the field (Burns and Madin, 2009). In order to determine if the Oregon Department of Geology and Mineral Industries (DOGAMI) protocol, developed by Burns and Madin (2009), could serve as the basis for establishing Washingtonís protocol, I used the office-based DOGAMI protocol to map landslides along a 40-50 km (25-30 mile) shoreline in Thurston County, Washington. I then compared my results to the field-based landslide inventory created in 2009 by the Washington Division of Geology and Earth Resources (DGER) along this same shoreline. If the landslide area I mapped reasonably equaled the area of the DGER (2009) inventory, I would consider the DOGAMI protocol useful for Washington, too. Utilizing 1m resolution lidar flown for Thurston County in 2011 and a GIS platform, I mapped 36 landslide deposits and scarp flanks, covering a total area of 879,530 m2 (9,467,160 ft2). I also found 48 recent events within these deposits. With an exception of two slides, all of the movements occurred within the last fifty years. Along this same coastline, the DGER (2009) recorded 159 individual landslides and complexes, for a total area of 3,256,570 m2 (35,053,400 ft2). At a first glance it appears the DGER (2009) effort found a larger total number and total area of landslides. However, in addition to their field inventory, they digitized landslides previously mapped by other researchers, and they did not field confirm these landslides, which cover a total area of 2,093,860 m2 (22,538,150 ft2) (DGER, 2009). With this questionable landslide area removed and the toes and underwater landslides accounted for because I did not have a bathymetry dataset, my results are within 6,580 m2 (70,840 ft2) of the DGERís results. This similarity shows that the DOGAMI protocol provides a consistent and accurate approach to creating a landslide inventory. With a few additional modifications, I recommend that Washington State adopts the DOGAMI protocol. Acquiring additional 1m lidar and adopting a modified DOGAMI protocol poises the DGER to map the remaining 87% of the state, with an ultimate goal of informing citizens and decision makers of the locations and frequencies of landslide hazards on a user-friendly GIS platform.