Bedrock Elevation Research Articles

Subsurface Geological Profile Interpolation Using a Fractional Kriging Method Enhanced by Random Forest Regression

Analyzing geological profiles is of great importance for various applications such as natural resource management, environmental assessment, and mining engineering projects. This study presents a novel geostatistical approach for subsurface geological profile interpolation using a fractional kriging method enhanced by random forest regression. Using bedrock elevation data from 49 boreholes in a study area in southeast China, we first use random forest regression to predict and optimize variogram parameters. We then use the fractional kriging method to interpolate the data and analyze the variability. We also compare the proposed model with traditional methods, including linear regression, K-nearest neighbors, ordinary kriging, and fractional kriging, using cross-validation metrics. The results indicate that the proposed model reduces prediction errors and enhances spatial prediction reliability compared to other models. The MSE of the proposed model is 25% lower than that of ordinary kriging and 10% lower than that of fractional kriging. In addition, the execution time of the proposed model is slightly higher than other models. The findings suggest that the proposed model effectively captures complex subsurface spatial relationships, offering a reliable and precise solution for performing spatial interpolation tasks.

Evaluating spatially enabled machine learning approaches to depth to bedrock mapping, Alberta, Canada.

Maps showing the thickness of sediments above the bedrock (depth to bedrock, or DTB) are important for many geoscience studies and are necessary for many hydrogeological, engineering, mining, and forestry applications. However, it can be difficult to accurately estimate DTB in areas with varied topography, like lowland and mountainous terrain, because traditional methods of predicting bedrock elevation often underestimate or overestimate the elevation in rugged or incised terrain. Here, we describe a machine learning spatial prediction approach that uses information from traditional digital elevation model derived estimates of terrain morphometry and satellite imagery, augmented with spatial feature engineering techniques to predict DTB across Alberta, Canada. First, compiled measurements of DTB from borehole lithologs were used to train a natural language model to predict bedrock depth across all available lithologs, significantly increasing the dataset size. The combined data were then used for DTB modelling employing several algorithms (XGBoost, Random forests, and Cubist) and spatial feature engineering techniques, using a combination of geographic coordinates, proximity measures, neighbouring points, and spatially lagged DTB estimates. Finally, the results were contrasted with DTB predictions based on modelled relationships with the auxiliary variables, as well as conventional spatial interpolations using inverse-distance weighting and ordinary kriging methods. The results show that the use of spatially lagged variables to incorporate information from the spatial structure of the training data significantly improves predictive performance compared to using auxiliary predictors and/or geographic coordinates alone. Furthermore, unlike some of the other tested methods such as using neighbouring point locations directly as features, spatially lagged variables did not generate spurious spatial artifacts in the predicted raster maps. The proposed method is demonstrated to produce reliable results in several distinct physiographic sub-regions with contrasting terrain types, as well as at the provincial scale, indicating its broad suitability for DTB mapping in general.

Hysteresis and orbital pacing of the early Cenozoic Antarctic ice sheet

Abstract. The hysteresis behaviour of ice sheets arises because of the different thresholds for growth and decline of a continental-scale ice sheet depending on the initial conditions. In this study, the hysteresis effect of the early Cenozoic Antarctic ice sheet to different bedrock elevations is investigated with an improved ice sheet–climate coupling method that accurately captures the ice–albedo feedback. It is shown that the hysteresis effect of the early Cenozoic Antarctic ice sheet is ∼180 ppmv or between 3.5 and 5 ∘C, depending only weakly on the bedrock elevation dataset. Excluding isostatic adjustment decreases the hysteresis effect significantly towards ∼40 ppmv because the transition to a glacial state can occur at a warmer level. The rapid transition from a glacial to a deglacial state and oppositely from deglacial to glacial conditions is strongly enhanced by the ice–albedo feedback, in combination with the elevation–surface mass balance feedback. Variations in the orbital parameters show that extreme values of the orbital parameters are able to exceed the threshold in summer insolation to induce a (de)glaciation. It appears that the long-term eccentricity cycle has a large influence on the ice sheet growth and decline and is able to pace the ice sheet evolution for constant CO2 concentration close to the glaciation threshold.

Simulation of a fully coupled 3D glacial isostatic adjustment – ice sheet model for the Antarctic ice sheet over a glacial cycle

Abstract. Glacial isostatic adjustment (GIA) has a stabilizing effect on the evolution of the Antarctic ice sheet by reducing the grounding line migration following ice melt. The timescale and strength of this feedback depends on the spatially varying viscosity of the Earth's mantle. Most studies assume a relatively long and laterally homogenous response time of the bedrock. However, the mantle viscosity is spatially variable, with a high mantle viscosity beneath East Antarctica and a low mantle viscosity beneath West Antarctica. For this study, we have developed a new method to couple a 3D GIA model and an ice sheet model to study the interaction between the solid Earth and the Antarctic ice sheet during the last glacial cycle. With this method, the ice sheet model and GIA model exchange ice thickness and bedrock elevation during a fully coupled transient experiment. The feedback effect is taken into account with a high temporal resolution, where the coupling time steps between the ice sheet and GIA model are 5000 years over the glaciation phase and vary between 500 and 1000 years over the deglaciation phase of the last glacial cycle. During each coupling time step, the bedrock elevation is adjusted at every ice sheet model time step, and the deformation is computed for a linearly changing ice load. We applied the method using the ice sheet model ANICE and a 3D GIA finite element model. We used results from a regional seismic model for Antarctica embedded in the global seismic model SMEAN2 to determine the patterns in the mantle viscosity. The results of simulations over the last glacial cycle show that differences in mantle viscosity of an order of magnitude can lead to differences in the grounding line position up to 700 km and to differences in ice thickness of the order of 2 km for the present day near the Ross Embayment. These results underline and quantify the importance of including local GIA feedback effects in ice sheet models when simulating the Antarctic ice sheet evolution over the last glacial cycle.

Inlet migration during the turning point between transgressive and regressive stages at the Guaratuba Holocene barrier, Paraná - Southern Brazil

This paper analyzes the sedimentary record of a paleo-inlet that evolved in association with a coastal barrier, using data that provides a high-resolution stratigraphy. During the evolution of the Guaratuba Holocene barrier, autogenic factors associated with the paleo-estuary dynamics compete with allogenic factors of RSL (Relative Sea Level) variation and the influence of climate on sediment input. The integration of geological mapping, vibrocore, GPR (Ground Penetrating Radar), and radiocarbon dating techniques allowed for recognizing of a paleo-inlet migrating east-southeast during the turning point between the transgressive and regressive stages, cyclically affected by episodes of aggradation and reversal in the direction of lateral migration. This dynamic, conditioned by shadow zones of bedrock elevations and the littoral drift of sediments, modifies the stratigraphy and the barrier's depositional architecture, defining diachronous surfaces and lateral variations in the stratigraphic stacking. These features exemplify how barrier-inlet systems are affected by the combination of high frequency, climate-defined events and local coastal processes during their evolution.

Prediction of subglacial lake melt source regions from site characteristics

AbstractSubglacial melt has important implications for ice-sheet dynamics. Locating and identifying subglacial lakes are expensive and time-consuming, requiring radar surveys or satellite methods. We explore three methods to identify source regions for lakes using seven continent-wide environmental characteristics that are sensitive to or influenced by ice-sheet temperature. A simple comparison of environmental properties at lake locations with their continent-wide distributions suggests a statistical relationship (high Kolmogorov-Smirnov statistic) between stable lake locations and ice thickness and surface temperatures, indicating melting under passive conditions. Active lakes, in contrast, show little correlation with direct thermally influenced parameters, instead exhibiting large statistical differences with horizontal velocity and bedrock elevation. More sophisticated techniques, including principal component analysis (PCA) and machine learning (ML) classification, provide better spatial identification of lake types. Positive PCA scores derived from the environmental characteristics correlate with stable lakes, whereas negative values correspond to active lakes. ML methods can also identify regions where subglacial lake melt sources are probable. While ML provides the most accurate classification maps, the combination of approaches adds deeper knowledge of the primary controls on lake formation and the environmental settings in which they are likely to be found.

Basal Melt Patterns around the Deep Ice Core Drilling Site in the Dome A Region from Ice-Penetrating Radar Measurements

Basal melt in the Dome A region will influence the deep-ice-core drilling at Kunlun Station. The melting point (wet bedrock) has a higher reflectivity than the surrounding area, which can be assessed using radar echoes from the bedrock. This paper uses a linear absorption model to determine wet and dry ice–bedrock interfaces around the Kunlun drilling site. In the determination process, an artificial intelligence model was applied to extract the ice–bedrock interface for inferring the ice thickness. Additionally, the various depth-averaged attenuation rates were used to identify the maximal range of basal melting. We mapped the patterns of the wet points on the bottom of the ice sheet and the modeled basal temperature to verify the results of the wet bed conditions. According to these maps of wet bed conditions, the areas with basal melting around the drilling site primarily appear in valley walls with low basal temperatures and are linked with hydraulic potential and bedrock elevation. Identifying the ice–bedrock interface is challenging in the valley bottom area, where the melting points are less than at the valley walls. Additionally, the melting proportions of 11.8% and 3.62% were calculated from two ice-penetrating radar data in this research. The mapped melting points around the site of Kunlun ice core drilling suggest complex ice flow effects and can be used to better interpret archives of old ice for paleoclimate research.

Antecedent geological control on transgressive delta and shoreline preservation: Examples from the SE African shelf

Using a suite of ultra-high resolution geophysical tools, remote operated vehicle dives, legacy cores, and isolated grab samples, we demonstrate at a regional shelf scale, the influence of antecedent geology (basement topography, shelf gradient and submerged shoreline features) on the evolving transgressive shelf stratigraphy of a subaqueous delta. The unconsolidated uppermost delta occurs as isolated remnants scattered across the low-gradient shelf. Seismic data reveal across-shelf heterogeneity in bedrock elevation, with prominent bedrock highs and depressions and several well-preserved aeolianite palaeo-shoreline complexes at water depths of −100, −60 and − 40 m. Analysis of bathymetric and seismic data demonstrates that these pre-existing shoreline complexes exert an overarching control on the distribution patterns of i) deltaic sediments, where they abut the landward flank of the shoreline form, acting as a barrier to seaward dispersal, and ii) shoreface sediments which remain sequestered on the mid-shelf on the seaward flank of the shoreline complex, hampering the landward translation of the shoreface in step with rising sea levels. Moreover, sediment distribution/accumulation is further constrained where elevated portions of the basement topography provide little to no accommodation, and act as zones of sediment bypass, and adjacent basement depressions accommodate sediment to seaward. The gentle antecedent slope, coupled with the gentle shoreline trajectory mediated transgressive erosion directly to landward. This is reinforced where local inflections in the wave ravinement profile form due to the presence of palaeo-shoreline complexes, aiding in the ultimate preservation of these submerged delta facets in the palaeo-shoreline lee. This study shows that low antecedent gradients, and palaeo-shoreline features, lead to development of transgressive coastal profiles that are predisposed to delta overstepping. The pre-existing basement topography and palaeo-shorelines constrain the positioning and morphology of the delta. We suggest that antecedent conditioning has partitioned accommodation for delta accumulation and moderated wave ravinement associated with transgression since the Last Glacial Maximum (LGM). The geological framework has acted as a recurring primary control on the geomorphic evolution of the submarine delta and shelf.

Gravity‐Derived Antarctic Crustal Thickness Based on the Gauss‐FFT Method

AbstractThe information on crustal thickness in Antarctica can provide significant constraints on its crustal deformation and tectonic evolution. To generate reliable images of crustal features, we investigate the model of Moho depth and crustal thickness beneath Antarctica by applying the Bott‐Parker's formulas based on the Gauss‐fast Fourier transform method through a comprehensive analysis of gravity data, ice and sediment thicknesses and bedrock elevation, combined with seismic constraints. Tests with synthetic data indicate that the iterative inversion algorithm can yield a highly accurate Moho topography. Ultimately, inverted crustal thickness reveals a more detailed crustal image by clearly identifying more tectonic elements at different scales than previous results and correlates well with major tectonic provinces. Airy isostasy and flexural isostasy models are used to assess the crustal isostatic compensation. The distinctive negative isostatic anomalies are observed in the Transantarctic Mountains and the East Antarctic areas of the great escarpment in the Dronning Maud Land and Aurora and Wilkes Subglacial Basin, indicating that the low density of the uppermost mantle and lithospheric strength may play important roles in compensating for their elevations. Variations in crustal thickness in interior East Antarctica are analyzed; the thickened crust from Dronning Maud Land to the Gamburtsev Subglacial Mountains may be associated with the collision of continental blocks and is interpreted as the fossil sutures. We compare the relationship between the Moho and Curie interfaces and find that the uppermost mantle is magnetized in some areas of East Antarctica, which may indicate preserved Precambrian cratonic roots.

Evidence of Former Sea Levels from a Passive Seismic Survey at a Sandy Beach; Perranporth, SW England, UK

Since the end of the last glaciation, the United Kingdom’s land surface has been altered by isostatic rebound, rising in the north and sinking in the south. Numerous studies have been published documenting the impact of isostatic rebound on relative sea levels. However, due to the difficulties in acquiring evidence to prove former sea levels, locally, these data can be sparse or absent. In this work, we explored the suitability of the passive seismic survey (PSS) method to estimate the contemporaneous beach thickness in coastal environments where there is a high impedance contrast between the beach deposits and the underlying wave-cut platform. We conducted a three-day survey at Perran Beach, Cornwall, collected 149 measurements using PSS, and interpreted the observations supported by auxiliary topographical, geological, and independent geophysical observation in the study area. The study site is a contemporaneous beach mostly composed of sand underlain by a wave-cut platform composed of igneous and sedimentary rock, therefore high impedance contrast with the sandy beach is anticipated. The elevation of the bedrock relative to the topographical elevation suggests that the bedrock elevation is −15 m ± 5 m below the present day mean sea level, which is coherent with the observation of relative sea level rise along the region of the south-west. The present study contributes to our current limited understanding of land and sea level movements by providing further subsurface information to the coastal geological archive of south-west England, a region currently in need of more data to reconstruct land- and sea-level movements.

Modelling evidence for late Eocene Antarctic glaciations

It is generally believed that a large scale Antarctic ice sheet formed at the Eocene-Oligocene transition (34.44-33.65 Ma). However, oxygen isotope excursions during the late Eocene (38-34 Ma) and geomorphic evidence of glacial erosion suggest that there were ephemeral continental scale glaciations before the Oi-1 event. Here, we investigate the Antarctic ice sheet evolution over a multi-million year timescale during the late Eocene up to the early Oligocene with the most recent estimates of carbon dioxide evolution over this time period and different bedrock elevation reconstructions. A novel ice sheet-climate modelling approach is applied where the Antarctic ice sheet model VUB-AISMPALEO is coupled to the emulated climate from HadSM3 using the coupler CLISEMv1.0. Our modelling results show that short-lived continental scale Antarctic glaciation might have occurred during the late Eocene when austral summer insolation reached a minimum in a narrow range of carbon dioxide concentrations. The Antarctic ice sheet first reached the coast in Prydz Bay and later in the Weddell Sea region, supporting the glaciomarine sediments dated prior to the EOT.

The influence of rock uplift rate on the formation and preservation of individual marine terraces during multiple sea-level stands

AbstractMarine terraces are a cornerstone for the study of paleo sea level and crustal deformation. Commonly, individual erosive marine terraces are attributed to unique sea-level high stands based on the reasoning that marine platforms could only be significantly widened at the beginning of an interglacial. However, this logic implies that wave erosion is insignificant at other times. We postulate that the erosion potential at a given bedrock elevation datum is proportional to the total duration of sea-level occupation at that datum. The total duration of sea-level occupation depends strongly on rock uplift rate. Certain rock uplift rates may promote the generation and preservation of particular terraces while others prevent them. For example, at rock uplift of ~1.2 mm/yr, the Marine Isotope Stage (MIS) 5e (ca. 120 ka) high stand reoccupies the elevation of the MIS 6d–e mid-stand, favoring creation of a wider terrace than at higher or lower rock uplift rates. Thus, misidentification of terraces can occur if each terrace in a sequence is assumed to form uniquely at successive interglacial high stands and to reflect their relative elevations. Developing a graphical proxy for the entire erosion potential of sea-level history allows us to address creation and preservation biases at different rock uplift rates.

Bedrock reconstruction from free surface data for unidirectional glacier flow with basal slip

Glacier ice flow is shaped and defined by several properties, including the bedrock elevation profile and the basal slip distribution. The effect of these two basal properties can present in similar ways in the surface. For bedrock recovery this makes distinguishing between them an interesting and complex problem. The results of this paper show that in some synthetic test cases it is indeed possible to distinguish and recover both bedrock elevation and basal slip given free surface elevation and free surface velocity. The unidirectional shallow ice approximation is used to compute steady state surface data for a number of synthetic cases with different bedrock profiles and basal slip distributions. A simple inversion method based on Newton's method is applied to the known surface data to return the bedrock profile and basal slip distribution. In each synthetic test case, the inversion was successful in recovering both the bedrock elevation profile and the basal slip distribution variables. These results imply that there is a unique bedrock profile and basal slip which give rise to a unique combination of free surface velocity and free surface elevation.

IMPROVEMENT OF STRUCTURAL ANALYSIS BY MODIFICATION OF SITE RESPONSE ANALYSIS AND EARTHQUAKE FORCE DIRECTION

The seismic resistance design of buildings in Indonesia is usually performed based on the most recent Indonesian Seismic Code (ISC). Improvement of the ISC will sometimes produce different (increasing or decreasing) earthquake forces for building design. This paper describes the study of building evaluation by conducting a modified Site Response Analysis (SRA) and a modified incoming earthquake force direction. The modified SRA was performed as a result of in-complete soil boring investigation through the bedrock elevation, using 30 m (soil boring investigation) and 50 m (microtremor investigation) soil profile models. The earthquake force towards the building is computed using a pair (East-West/EW and North-South/NS) of ground motion (acceleration time histories) models. Due to the unpredictability of incoming earthquake forces, these models are improved by adjusting the ground motion direction to the building by rotating a pair of ground motions from 0o to 180o at 10o intervals. The time histories weredeveloped according to the Deterministic (DSHA) and Probabilistic (PSHA) Seismic Hazard Analysis. The analysis was performed by evaluating the moment capacity ratio of the building’s existing column reinforcement in resisting the earthquake force. The surface time histories and spectrum accelerations computed using the 50 m model is slightly higher (10%) than the same output calculated using the 30 m model. In terms of the moment capacity ratio, the existing column reinforcement of the building is strong enough to resist seismic earthquake scenarios, and the minimum value is observed when the ground motion is rotated by 40o.

Inherited periglacial geomorphology of a basalt hill in the Sudetes, Central Europe: Insights from LiDAR‐aided landform mapping

AbstractThis paper presents the results of mapping inherited cold‐climate landforms and deposits at the Muchów Hills site in the Sudetes, Central Europe (51°N). Bedrock supporting the upper slopes is basalt, with well‐developed columnar jointing, and local relief is of the order of 100 m, whereas the middle and lower slopes are underlain by glacigenic deposits. Characteristic cold‐climate landforms include rock cliffs, talus aprons, mid‐slope benches, solifluction sheets and lobes, and are integrated into sediment transport systems linking bedrock elevations with the adjacent plain. Ground identification of these features may be ambiguous, but the interpretation of high‐resolution LiDAR Digital Terrain Modeling data helped to provide a comprehensive spatial picture. Reasons for partial remodeling of inherited periglacial landforms in the Muchów Hills include biogeomorphological activity throughout the Holocene, mainly tree uprooting and resultant regolith mixing.

Time‐transgressive cryoplanation terrace development through nivation‐driven scarp retreat

AbstractCryoplanation terraces are elevated bedrock features resembling giant staircases, with alternating steep scarps and shallow sloping treads. These landscape‐scale features have long been associated with periglacial environments, but the processes involved in their formation remain vaguely specified and contentious. The two leading hypotheses for the formation of cryoplanation terraces are centered on: (1) geologic structure; and (2) nivation‐driven scarp retreat. The nivation‐formation hypothesis invokes scarp retreat under erosion processes associated with late‐lying snowbanks. To test whether cryoplanation terraces develop through scarp retreat, six relative weathering indices (fracture counts, Cailleux roundness, Cailleux flatness, Krumbein sphericity, rebound, and weathering rind thickness) were measured across well‐developed terraces in unglaciated eastern Beringia (central and western Alaska) at Mt Fairplay, Eagle Summit, and Skookum Pass. Statistically significant differences in relative weathering indices detected through chi‐square and multiple‐comparison procedures indicate that material is less weathered closer to scarps, i.e. that these areas were more recently exposed than those distant from the scarp. Based on these findings, a refined model of time‐transgressive cryoplanation terrace development through nivation‐driven scarp retreat is proposed. This new qualitative model addresses the removal of weathered material from terrace treads down side slopes through piping and gravity‐driven mass‐wasting processes. © 2019 John Wiley & Sons, Ltd.

Development of surface ground motion and spectral acceleration based on modified shear wave propagation analysis

The shear wave propagation analysis (SPA) is a method for developing surface ground motion and spectral acceleration by conducting seismic wave propagation from bedrock to surface elevation. The bedrock elevation and soil profile from bedrock through surface elevation are important information needed for SPA. The shear wave velocity, damping ratio, soil density and shear modulus are the four parameters that are needed for developing SPA. This paper describes a modified SPA for developing surface ground motion and spectral acceleration by conducting two different soil deposit model, actual and modified bedrock elevation models. The difference in the two models depends on the bedrock elevation used for SPA. The actual bedrock elevation model is conducted based on the real bedrock elevation observed using microtremor test. However, the second model is performed based on the final depth of soil boring investigation. The study was conducted at the Engineering Faculty of Diponegoro University. Based on the microtremor observation, the bedrock elevation was predicted at 50 metres below the surface level. However, two soil boring investigations conducted in this study area were stopped at 30 metres below the surface level. The surface ground motion and spectral acceleration developed from both models are almost equal.

Indoor and Outdoor Exposure to Radon, Thoron and Thoron Decay Products in a NORM Area with Highly Elevated Bedrock Thorium and Legacy Mines.

Radon (222Rn) and thoron (220Rn), and especially their short-lived decay products, are major contributors to dose received by the public from naturally occurring radioactive material (NORM), particularly in areas with elevated levels of naturally occurring radionuclides. Mining in such areas can involve ventilation of high amounts of these gases, which may influence outdoor levels. In this work, we assessed indoor and outdoor levels of 222Rn, 220Rn and 220Rn decay products (TnDP) in close proximity to an area with elevated bedrock levels of thorium (232Th) and a NORM legacy mining site with high natural ventilation. We assess municipal buildings at distances from a few hundred meters to 2 km from the NORM legacy mines. In some buildings, high indoor levels of 222Rn were observed in winter, as expected for temperate areas. In summer, high indoor levels of 222Rn and 220Rn were observed in some buildings, and very low associated levels of TnDP in actively ventilated buildings may suggest entry by ventilation and an outdoor source. Outdoor levels of TnDP increased with decreased distance from the legacy mines, suggesting dispersal from these during both summer and winter.

Allochthonous sources of iodine and organic carbon in an eastern Ontario aquifer

Regional geochemical characterization of groundwaters in a bedrock aquifer in the Ottawa – St. Lawrence Lowlands of eastern Ontario has identified an iodine (I) anomaly, with values regularly exceeding 150 μg/L and a maximum observed concentration of 10 812 μg/L. The spatial distribution, enrichment mechanisms, and sources of I and organic matter were investigated using geochemical and isotopic data. High-I groundwaters (>150 μg/L) are prevalent in Na–Cl-type groundwaters at low bedrock elevations in areas overlain by thick layers of glacial sediments. I is thought to be linked to massive muds in the glacial sediments overlying the aquifer, deposited during the postglacial incursion of the Champlain Sea 12–10 ka BP. Principal component analysis of I and 18 other chemical parameters revealed correlations among I, salinity, and indicators of microbial oxidation of organic matter, suggesting that the intrusion of saline pore waters affected by decomposition of organic matter such as marine phytoplankton in the massive muds is the dominant process controlling I enrichment in groundwater. 129I/127I ratios in the pre-modern waters vary between near-marine values of 460 × 10−14 and 5 × 10−14, demonstrating that older allochthonous I derived from the surrounding Paleozoic sedimentary terrain also contributed to the I pool in the Champlain Sea basin. 14C ages and δ13C signatures for dissolved organic carbon in groundwater and disseminated organic carbon within the glaciomarine muds highlight an allochthonous source of terrestrial organic carbon predating the Champlain Sea incursion, likely transported via glacial meltwaters in tandem with I to the Champlain Sea basin.

Three-dimensional analysis of unanticipated behavior of a deep excavation

A three-dimensional (3D) effective stress finite element analysis, modified to account for hydrofracturing and gassy soil behavior, is used to examine the potential for the venting of water and gas from a bedrock aquifer and through 13–14 m of low permeability clayey silt between the base of the excavation and the bedrock following excavation to about 24 m in an approximately 40 m thick clayey silt deposit. The clayey deposit contained sand lenses with dissolved gas. The analysis predicts that the exsolution of this dissolved gas, caused by a reduction in total stress due to the excavation, results in liquefaction of the sand in the lenses and consequent lateral deformations of the side slopes. The analysis predicts hydrofracturing through the remaining clayey silt when the excavation reaches its final depth and this explains the venting of water and gas from the underlying aquifer that was observed above a local bedrock high. The presence of gassy sand lenses created weak zones within the clayey deposit that influenced the path of the hydrofracturing. However, the analyses suggest that, for the depth of excavation and bedrock elevation examined, hydrofracturing and subsequent venting would have occurred even if there had been no sand lenses.