Quick Clay Research Articles

Spatial quick-clay predictions using multi-criteria evaluation in SW Sweden

The transformation of marine and glaciomarine clay deposits into high sensitive and quick clays is largely dependent on the influence of local and regional geologic history and the resulting stratigraphy. The general conditions that facilitate quick-clay development are well known from numerous laboratory investigations during the last century, but their local and regional in-field variation is less understood. In this study, the geographic distribution of quick clay in SW Sweden is predicted using a multicriteria evaluation model that incorporates both qualitative information (established theory and expert judgment concerning the influences on both quick-clay development and the stratigraphic and geomorphologic distribution of sediment types) and observational data (maps of surficial deposits, geotechnical records and digital elevation data). This information duality cannot be avoided if knowledge from different disciplines is utilized. Considering this, model transparency is important for improvements and for characterizing its reliability for risk analysis. The model was constructed stepwise by an initial parameterization with subsequent hierarchical structuring, weighting and standardization of criteria, before running the full analysis. Comparisons between regional model results and geotechnically documented localities have yielded promising results concerning the model’s ability to predict general trends. However, the large natural and site-specific variability of clay sensitivities is not always captured by the model. These deviations are examined and suggestions are given for minimizing their effect. Applications of model methodology and results are briefly discussed.

High-resolution 3D reflection seismic investigation over a quick-clay landslide scar in southwest Sweden

Quick-clay landslides often occur in the northern hemisphere in areas that were covered by Pleistocene glaciations. They are particularly common along the shorelines of the Göta River in southwestern Sweden. Characterization of potential landslide areas and identification of features that indicate high risk are necessary to better understand the triggering mechanisms of these events. Therefore, an intensive characterization project was initiated at the Fråstad landslide in Sweden. Part of the characterization program included the acquisition of 3D reflection seismic data to image structures in the normally consolidated sediments, as well as the bedrock topography below the landslide scar. Two seismic horizons within the glacial and postglacial sediments were observed. The shallowest seismic horizon (here, referred to as S1) corresponds to a coarse-grained layer that was previously detected by eight geotechnical boreholes located within the 3D survey area. Discontinuities in S1, mapped by the 3D reflection seismic data, occur across a zone that correlates with the landslide scar boundary, suggesting that this zone may have played a role in triggering and/or in limiting the extension of the landslide. If S1 is truncated by or mixed with clays in this zone, then the outflow of water from the permeable S1 into the clays above may have increased the amount of quick clays above this zone. The increased outflow of water may also have caused a higher pore-water pressure south of the zone, which in turn could have acted as a trigger for the landslide. We evaluated the potential of using the 3D reflection seismic method as a complement to drilling and other geophysical methods when performing landslide site investigations. We also demonstrated the importance of further investigating the relationship between 3D subsurface geometries and landslide development.

Towards geophysical and geotechnical integration for quick‐clay mapping in Norway

ABSTRACTQuick clay is a known hazard in formerly‐glaciated coastal areas in e.g., Norway, Sweden and Canada. In this paper, we review the physical properties of quick clays in order to find a suitable, integrated and multi‐disciplinary approach to improve our possibilities to accurately identify the occurrence of quick clay and map its extent both vertically and laterally. As no single geophysical method yields optimal information, one should combine a variety of geophysical methods with geotechnical data (in situ measurements using Cone Penetration Testing (CPTU), Seismic CPTU (SCPTU) and Resistivity CPTU (RCPTU); laboratory tests) for an in‐depth quick‐clay assessment at a given site. In this respect, geophysical data are used to fill the gaps between geotechnical boreholes providing ground‐truth. Such an integrated and multi‐disciplinary approach brings us closer to 2D or pseudo‐3D site characterization for quick clays and as such, an improved assessment of the potential hazard they pose. The integrated approach is applied in practice on two Norwegian quick‐clay sites. The first site, Hvittingfoss, was remediated against potential landslides in 2008 whereas the second one, Rissa, was the scene of a major quick‐clay landslide in 1978, quick clays being still present over a large area. The collected data and preliminary site characterizations illustrate the high diversity as well as the complexity and clearly emphasize the need for higher resolution, careful imaging and calibration of the data in order to accomplish the assessment of a quick‐clay hazard.

Constitutive Model for Cemented Quick Clay in Deep-Mix Stabilization

Based on experimental results of a laboratory study on lime and cement stabilized quick clay, carried out at NTNU, Norway, a two yield surface constitutive model has been formulated. The model reproduces the mechanical behavior of the stabilized soil in the p’-q stress space. In the stress space the two yield surfaces appear as a “cap” and a “wedge” respectively. The formulated model “QUICKSTAB” has satisfactorily been able to simulate the overall geomechanical behaviour of stabilized quick clay and it agrees very well with the results established by the laboratory study. Typical model parameter values are also presented herein for various soil-binder mixtures. An extensive parametric study has been recommended to thoroughly explain the behavior of the model, including tests following different stress paths and extension of the study to include anisotropic cap-surface.

Relationship between electrical resistivity and basic geotechnical parameters for marine clays

Recently, considerable efforts have been made in the attempt to map quick clay areas using electrical resistivity measurements. However there is a lack of understanding regarding which soil parameters control the measured resistivity values. To address this issue, inverted resistivity values from 15 marine clay sites in Norway have been compared with basic geotechnical index properties. It was found that the resistivity value is strongly controlled by the salt content of the pore fluid. Resistivity decreases rapidly with increasing salt content. There is also a relatively clear trend of decreasing resistivity with increasing clay content and plasticity index. Resistivity values become very low (≈5 Ω·m) for high clay content (>50%), medium- to high-plasticity (Ip ≈ 20%) materials with salt content values greater than about 8 g/L (or corresponding remoulded shear strength values greater than 4 kPa). For the range of values studied, there is poor correlation between resistivity and bulk density and between resistivity and water content. The data studied suggest that the range of resistivity values corresponding to quick clay is 10 to 100 Ω·m, which is consistent with other published limits. A comparison is made between two-dimensional electrical resistivity tomography (ERT) and resistivity cone penetration test (RCPTU) data for two of the sites and the two sets of data show similar trends and values irrespective of scale effect.

A hypothesis of the Senoumi submarine megaslide in Suruga Bay in Japan—based on the undrained dynamic-loading ring shear tests and computer simulation

The distinctive bathymetric feature exists in the Suruga Bay, Japan. It has been called as Senoumi (Stone flower sea) from old times. Senoumi is a 30 km wide and 20 km long concave feature. Its origin has not been explained yet; however, the feature might be a combined consequence of intensive tectonic activity in the plate border, landslides, and a submarine flow coming from the Oi River. If the Senoumi was caused by a landslide, the latter would be larger than any on-land landslide in Japan. The downshelf “exit” from this feature is much narrower than its central part. This is not usual shape of landslides, but it is similar to the liquefied landslides such as those in quick clays which mobilize great strength reduction after failure. To study Senoumi as a landslide, the shear behaviors of the following three soil samples were investigated by the cyclic and seismic undrained stress control ring shear tests. One sample is volcanic ash taken from the base of landslide deposits (mass transport deposits), from 130 to 190 m deep layer below the submarine floor which was drilled and cored by the Integrated Ocean Drilling Program Expedition 333. Another two samples are the Neogene silty–sand and silt taken from the Omaezaki hill adjacent to the Senoumi, because the shear zone might have been formed in Neogene layers extending from on-land to the continental shelf. The largest strength reduction from peak to steady-state shear resistance in the undrained cyclic loading test was found in volcanic ash. The strength reduction in Neogene silty–sand was smaller than volcanic ash, while the Neogene silt mobilized the least post-failure strength reduction. An integrated model simulating the initiation and motion of earthquake-induced rapid landslides (landslide simulation (LS)-RAPID, Sassa et al. Landslides 7–3:219–236, 2010) was applied to this study. The steady-state shear resistance and other geotechnical parameters measured by the undrained ring shear tests and the greatest strong motion record in the 2011 off-the-Pacific Coast of Tohoku earthquake (M w 9.0), also known as “2011 Tohoku Earthquake” at the observation point MYG004 (2,933 gal) were input to this model. As the result, it was found that landslides would be triggered by 0.30–1.0 times of MYG004 in volcanic ash, 0.4–1.0 times of MYG004 in Neogene silty–sand and Neogene silt, though the depth and area of triggered landslides were different in soils and intensity of shaking. Feature, created by LS-RAPID using the parameters of volcanic ash, was most similar to the Senoumi in depth and extent. The result obtained from this study includes a hypothesis to be proved, but presents the strong need to investigate the risk of the large-scale submarine landslides which could enhance tsunami wave and possibly enlarge the submarine landslide retrogressively into the adjacent coastal plain by the upcoming mega earthquake in the Nankai Trough.

Mapping of quick clay using geoelectrical imaging and CPTU‐resistivity

ABSTRACTQuick clay has a major impact on landslide risk and it is therefore of considerable interest to map its presence and extent. In Sweden, quick clay has been involved in most landslides in soft clay with serious consequences. The predominant method for detection of quick clay in Sweden has been to take undisturbed samples and to perform fall‐cone tests on the clay in its undisturbed and remoulded state. Originally deposited in saltwater in a marine environment, the salt maintains the stability of the clay. When the salt is leached out, the clay can become quick. When salt is leached from clay of marine origin the resistivity increases. In this study the intention was to calibrate electrical imaging with cone penetration tests with resistivity measurements (CPTU‐R) and measurement of the total penetration resistance, i.e. the total rod friction, together with both geotechnical and chemical analyses on specimens in the laboratory.The results show that electrical imaging can be used for separation of leached soil volumes in marine clays that may form quick clay, from those where the salt content remains too high for this. In the dry crust and thin weathered zone at the top, the resistivity is high but the clay is non‐quick. Also soils with less clay content will have higher resistivity without being quick. The technique may thus be used as a screening tool in order to delimit areas where further investigations are needed from areas that do not require more attention. This has a potential of saving significant resources if used in a relatively early stage of the survey process. It can also increase the overall quality and reliability of the survey results.The induced polarization (IP) results are consistent and seem to be geologically realistic, and appear to contain additional information to the resistivity that is related to material or electrochemical properties, although it is not clear how due to lack of sufficiently detailed reference data.The electrical imaging gives a general picture of the variation in resistivity along soil sections. The CPTU‐R gives variations at depth with very high vertical resolution in one specific location. There is generally good agreement between the models based on electrical imaging and the CPTU‐R. The CPTU‐R results may be used to calibrate the electrical imaging results with quick clay estimations based on rod friction.

Shear Behavior of Sensitive Marine Clay-Concrete Interfaces

AbstractThis paper presents the results of an experimental study on Leda clay-concrete interface shear behavior. Leda clay specimens with various dry densities and salt contents are used to study the influence of these parameters on the interface shear behavior. The results indicate the following: (1) interface shear resistance increases as the concrete surface roughness increases; (2) increasing the dry density of remolded Leda clay enhances its interface shear resistance; and (3) a higher salinity of Leda clay pore water results in higher interface shear strength.

Modeling parameters of structured clays as a multivariate normal distribution

This study explores the possibility of modeling liquidity index, undrained shear strength, remolded undrained shear strength, preconsolidation stress, and vertical effective stress of structured clays (sensitive or quick clays) as a multivariate normal distribution. The literature is replete with correlation equations between two soil parameters. Consistent synthesis of more than two soil parameters through construction of a multivariate probability distribution function is rare, despite obvious practical usefulness of such an approach. This study compiles a large database of structured clays to construct the multivariate probability distribution among the aforementioned five soil parameters. This multivariate distribution is then used to simulate the correlations between soil parameters of interest and to derive useful equations for Bayesian inference. This constructed multivariate distribution and equations are further validated by another independent database of structured clays as well as by empirical equations proposed in the literature.

Flow of clays

Recent experimental and theoretical research into physical phenomena in clays is reviewed. Clays’ present place in the context of modern materials science is briefly discussed, and illustrated through the rich behavior recently displayed in this physical model system. We will show that in order to understand macroscopic flow behaviors in these systems, it is crucial to know the underlying nanostructures in detail. With the clay nanostructural basis at hand, we will review recent advances in clay systems from the geological example of quick clay flows and avalanches, to materials science and the stability, strength and flow of smart electrorheological clay structures. In the case of natural quick clay, there is now hope of establishing a protocol for avalanche preditction based on rheological sample data. In materials science, the use of electric fields together with flow in order to improve the processing of clay composite materials may open new unexplored avenues. We will finally discuss that due to the interplay of van der Waals and electrostatic forces screened by ions at the nanoscale, clays may either form a glass, or a gel and thus give fundamental insights into the elusive questions related to materials universal aging flow dynamics.

Multi‐method geophysical mapping of quick clay

ABSTRACTMarine clay deposits in coastal, post‐submarine areas of Scandinavia and North America may be subjected to quick clay landslides and hence significant efforts are being taken to map their occurrence and extent. The purpose of this paper is to assess the use of a number of geophysical techniques for identifying quick clay. The investigated area, Smørgrav, located in southern Norway has a history of quick clay sliding, the most recent event occurring in 1984. Geophysical techniques that are used include electromagnetic conductivity mapping, electrical resistivity tomography, seismic refraction and multichannel analysis of surface waves. These results are compared to geotechnical data from bore samples, rotary pressure soundings and cone penetration testing. A number of these approaches have proved promising for identifying quick clay, in particular electrical resistivity tomography and electromagnetics, which delineated a zone of quick clay that had previously been confirmed by rotary pressure soundings and sampling. Seismic refraction was useful for determining the sediment distribution as well as for indicating the presence of shallow bedrock whereas the multichannel analysis of surface‐waves approach suggested differences between the intact stiffness of quick and unleached clay. It is observed that quick clay investigations using discrete rotary pressure soundings can be significantly enhanced by using, in particular, electrical resistivity tomography profiles to link together the information between test locations, perhaps significantly reducing the need for large numbers of soundings.

‘We are not afraid of flooding’ …but what about landslide? The effects of assumed and perceived hazards on the value of residential locations

The issue of interest is as to whether it is possible to use the added value generated by the positive effect of having housing developments close to water to abate the negative effects that may arise in some of the same locations. This paper is based on a critical literature review together with expert interviews. It comments on the methodology of spatial benefit-cost analysis (benefit for price premium ‘generated’ by the amenity; costs for price discounts ‘generated’ by the risk) in relation to design of housing and hazard management mechanisms within the context of urban land use in contemporary Trondheim, Norway. The particular issue at stake concerns potential quick clay landslide hazard areas. Prior research suggests that the role of situation by a coast, lake or river should not be overlooked when assessing possibilities for financing general water related hazard abatement schemes. This is potentially a win-win-situation: safety enhancing public works could be financed based on the added value of new developments at the given location. The purported kind of innovative financial mechanism however requires close cooperation between private developers and local authorities – in other words, governance instead of government – as well as flexible legal codes. Unfortunately the study area to a great extent lacks such institutional possibilities due to an outmoded conceptualization of the role of government intervention.

Monofrequency in situ damping measurements in Ottawa area soft soils

In Ottawa, Canada, unusually high amplification ratios have recently been measured in clayey silts (called ‘Leda Clays’) at low levels of earthquake-induced ground shaking. However, the contribution of seismic Q, or material damping ( ξ=1/2 Q), to the overall ground motion at soft soil sites across the city is not well understood. This research investigates attenuation measurements in soft soils ( V s <250 m/s) for ongoing seismic hazard evaluation in the Ottawa area. The work focuses on in situ measurements of damping in two deep boreholes drilled into Leda Clay. To investigate the possibility of frequency-dependent dynamic properties of these materials at low strains, a new approach to the spectral ratio technique has been developed for the measurement of Q s in the field using a mono-frequency vibratory source (generating signals between 10 and 100 Hz), and two identical downhole 3-component geophones. Monofrequency signals also allowed for the measurement of dispersion (variation of velocity with frequency). Analysis of the data show that dynamic properties are, for the most part, independent of frequency in the homogenous silty soils, yielding negligible variation in shear wave velocity (<2 m/s) across the frequency test band, and small strain Q s 's ranging from 170 to 200 (damping of 0.25–0.30%) over soil thickness intervals ranging from 10 to 60 m. At intervals within 20 m of the ground surface, laminated silt and clay beds of elevated porosity are found to have slight influence on the frequency dependence of damping for frequencies greater than 70 Hz (damping increase to 0.6%).

Combined geophysical and geotechnical approach to ground investigations and hazard zonation of a quick clay area, mid Norway

Mapping of quick clay is important for hazard zonation, planning and protection purposes. The present study focuses on an area prone to quick clay landslides in mid Norway, which is investigated through a combination of geophysical and geotechnical methods. The following classes are suggested for a first-order interpretation of resistivity profiles in areas with few or no previous investigations: Unleached clay deposits: 1–10 Ωm; Leached clay deposits, possibly quick: 10–100 Ωm; Dry crust clay deposits and coarse sediments: >100 Ωm. In the study area, 14–80 Ωm was found as the main resistivity interval for quick clay. The resistivity values from the present study are compared to previously published values. Classification of material from resistivity values is influenced by local conditions, and there is an overlap between the classes. Resistivity profiles can give valuable information for hazard zonation and may assist in maximising subsequent intrusive investigations.

Expanding collapse in partially submerged granular soil slopes

The traditional approach to stability analysis of granular slopes leads to the safety factor that is associated with a planar failure surface approaching the slope face, whether the slope is “dry” or submerged. However, for partially submerged slopes, a more critical, nonplanar failure surface can be formed. A family of geometrically similar surfaces can be found that is characterized by the same safety factor. If the safety factor drops down to unity and the slope becomes unstable, then a mechanism of any size can form. Alternatively, the failure may start at some small region and then the volume of the mechanism of failure can expand, giving rise to a progressive failure of a different kind that is typically associated with slopes. This progression has the character of a “disturbance” or a shock-like kinematic discontinuity propagating into the soil at rest. A quantitative analysis is presented and it is demonstrated that the soil dilates while the mechanism expands, leaving the slope weakened and susceptible to a deep failure. This is a plausible mode of failure of partially submerged slopes, the type that is most likely responsible for large subaqueous landslides, and is similar to the well-documented instability propagation in “quick clay.”

Quick Clay and Landslides of Clayey Soils

We study the rheology of quick clay, an unstable soil responsible for many landslides. We show that above a critical stress the material starts flowing abruptly with a very large viscosity decrease caused by the flow. This leads to avalanche behavior that accounts for the instability of quick clay soils. Reproducing landslides on a small scale in the laboratory shows that an additional factor that determines the violence of the slides is the inhomogeneity of the flow. We propose a simple yield stress model capable of reproducing the laboratory landslide data, allowing us to relate landslides to the measured rheology.

Integrated database for rapid mass movements in Norway

Abstract. Rapid gravitational slope mass movements include all kinds of short term relocation of geological material, snow or ice. Traditionally, information about such events is collected separately in different databases covering selected geographical regions and types of movement. In Norway the terrain is susceptible to all types of rapid gravitational slope mass movements ranging from single rocks hitting roads and houses to large snow avalanches and rock slides where entire mountainsides collapse into fjords creating flood waves and endangering large areas. In addition, quick clay slides occur in desalinated marine sediments in South Eastern and Mid Norway. For the authorities and inhabitants of endangered areas, the type of threat is of minor importance and mitigation measures have to consider several types of rapid mass movements simultaneously. An integrated national database for all types of rapid mass movements built around individual events has been established. Only three data entries are mandatory: time, location and type of movement. The remaining optional parameters enable recording of detailed information about the terrain, materials involved and damages caused. Pictures, movies and other documentation can be uploaded into the database. A web-based graphical user interface has been developed allowing new events to be entered, as well as editing and querying for all events. An integration of the database into a GIS system is currently under development. Datasets from various national sources like the road authorities and the Geological Survey of Norway were imported into the database. Today, the database contains 33 000 rapid mass movement events from the last five hundred years covering the entire country. A first analysis of the data shows that the most frequent type of recorded rapid mass movement is rock slides and snow avalanches followed by debris slides in third place. Most events are recorded in the steep fjord terrain of the Norwegian west coast, but major events are recorded all over the country. Snow avalanches account for most fatalities, while large rock slides causing flood waves and huge quick clay slides are the most damaging individual events in terms of damage to infrastructure and property and for causing multiple fatalities. The quality of the data is strongly influenced by the personal engagement of local observers and varying observation routines. This database is a unique source for statistical analysis including, risk analysis and the relation between rapid mass movements and climate. The database of rapid mass movement events will also facilitate validation of national hazard and risk maps.

Development of the submarine channel in front of the Nidelva River, Trondheimsfjorden, Norway

The Nidelva submarine channel is a 4 km long and up to 100 m deep erosional feature seaward of a former outlet of the river Nidelva, in Trondheimsfjorden, Norway. Mechanisms of channel evolution and sedimentary processes associated to its development are revealed by analyzing the sea bed morphology and the stratigraphy of the fjord sediments, based on high resolution seismic data, swath bathymetry imaging and gravity cores. Seismic data suggest that the submarine channel is eroded into Late Pleistocene clays filling a bedrock depression. A thick Holocene deltaic succession overlay the glacio-marine clays. Seismic facies suggest that turbidity currents likely initiated downcutting of the channel in the Early-Holocene. These flows were guided by sea bed morphology and bedrock exposures, favouring confluence and acceleration. The area was subjected to glacio–isostatic rebound and rapid fall of relative sea-level in the Holocene. This resulted in a fast delta progradation likely with an enhanced activity of delta derived sediment gravity flows promoting incision. In Mid- to Late-Holocene, such flows, including hyperpycnal flows, were responsible for deep incision of the channel. Laminated clay-rich beds are interpreted as deposited by the overspill of turbiditic flows along the channel and likely originate from large terrestrial quick clay slides. In the last 2000 years, erosive flows initiated by mass wasting along the coast controlled the evolution of the channel after the mouth of the Nidelva River moved to a more easterly position in the embayment. Present slope failure along the channel flanks seems to be controlled by the distribution of the laminated clay-rich layers in the stratigraphy.

Mapping of quick clay formations using geotechnical and geophysical methods

Quick clay has been involved in most serious, large clay slides in Sweden, Norway, and Canada. This paper describes geotechnical and geophysical methods that can be used to locate and map the extent of quick clay formations. Surface resistivity measurements and four different sounding methods have been tested. The results have been compared with sensitivities determined using fall-cone tests. The investigation shows that there is a correlation between sensitivity and electrical resistivity, which can be used to discriminate between marine clays that have been leached sufficiently to possibly form quick clay and those for which the salt content remains sufficiently high to prevent this. Although the most reliable evaluation of the variation in sensitivity was obtained by the CPT with additional measurement of total penetration force, this investigation suggests that any sounding method that uses a constant rate of advance into the ground and in which the penetration force applied on the top of the rods is measured may be used for quick clay mapping.

Probability risk assessment of landslides: A case study at Finneidfjord

Probabilistic risk assessments are increasingly being considered the most appropriate framework for engineers to systematically base decisions on hazard mitigation issues. This paper aims to show the advantages of a quantitative risk assessment by application to a historical case study. The generalized integrated risk assessment framework has been applied retrospectively to a submarine landslide that occurred in 1996 near the village of Finneidfjord in northern Norway. Over 1 million cubic metres of predominantly quick clay was displaced. Even though it was triggered underwater on the embankment of the Sørfjord, the retrogressive nature of the slide resulted in it encroaching 100–150 m inland. The triggering mechanism is believed to have been the placement of fill, from a nearby tunnelling project, on the foreshore of the embankment. This paper is a retrospective quantitative evaluation of the risk to the neighbouring houses, the persons in those houses, and the persons in open spaces caused by the placement of increasing levels of embankment fill. A probabilistic approach, making use of second-moment modelling and first-order second-moment approximation is adopted. It aims to demonstrate the advantages of this type of risk assessment in understanding complex and integrated hazards, particularly those in populated environments.