Sediment Thickness Research Articles

Seismic vulnerability distribution of the earthquake prone area in Central Aceh

Abstract Central Aceh is known as an earthquake-prone area due to the presence of the Sumatra Fault and secondary fault crossing the area resulting the M 6.1 earthquake in 2013. The objective of the research is to investigate seismic vulnerability of Central Aceh. Microtremor observations were conducted at 25 locations in Central Aceh analysed using the Horizontal to Vertical Spectral Ratio (HVSR). Seismic amplification is represented by the maximum values of H/V. The dominant frequency f0 representing the sediment thickness was obtained when the H/V curve is maximum. The frequency values ranged from 14.9 to 25.9 Hz, and the seismic vulnerability values varied from 0.02 to 3.1. The top soil consists of thin layers of sedimentary material dominated by tertiary rocks. The rocks are classified as Type IV and Class I based on the dominant frequency values. This result is consistent with the geology, consisting of sedimentary and volcanic rocks. Areas with a seismic vulnerability value > 1.0 is considered to have a high risk of earthquake damage because they are likely to experience earthquake amplification. Areas with seismic vulnerability value > 1.0 are found in West Aceh and the Geureudong volcano area.

Correlation between flood couplet-based sediment yield and rainfall patterns in a small watershed on the Chinese Loess Plateau

Understanding how rainfall patterns impact sediment yields during rainfall events is important to explore soil erosion changes and to prevent watershed flood disasters. However, few studies identified sediment yields variation affected by rainfall patterns at continuous deposited flood couplet scale. This study analyzed sediment yields under different rainfall patterns and correlated particle depositional features on the condition that the erosive environment was obviously changed by anthropogenic activities in a check dam-controlled small watershed on the Chinese Loess Plateau. Results showed that erosive rainfall events, having formed 40 flood couplets behind the check dam between 1978 and 2010, can be grouped into three rainfall patterns (i.e., rainfall patterns I, II, and III) according to changes in rainfall amount, erosivity, and duration. Rainfall pattern III, characterized by high rainfall amount and erosivity within a short duration (i.e., rainstorm) accounted for only 15.0 % of erosive rainfall events while yielding 44.7 % of the total sediment. Conversely, rainfall pattern I (the most common, 55.0 %), characterized by its long duration and relatively low rainfall amount and erosivity, only contributed 25.3 % of the total sediment. Both sediment yields and thickness proportions of coarse particle layer in flood couplets exhibited a decreasing trend under all three rainfall patterns. Meanwhile, the main factor that impacted sediment yields changed from rainfall in the 1980s and 1990s to the “Grain for Green” program dominated anthropogenic activities since 1999. For a small watershed with a few terraces and only one check dam at the outlet, artificial vegetation restoration significantly reduced coarse soil particle loss and sediment yield. Findings showed that sediment yield and deposition feature variations in flood couplets under different rainfall patterns are helpful in improving watershed soil and water conservation measures against flood disasters, particularly under rainstorms.

The spatio-temporal variation of exhumation and its impacts on the resource distribution in the Northern Carnarvon Basin, Australia

Deformation and rock fracture during exhumation can enhance rock porosity and permeability, thereby improving fluid migration through reservoirs and destroying seal integrity in caprocks. Therefore, analysing the spatio-temporal distribution of exhumation in connection to fluid migration models can be crucial for energy resource exploration and storage. Geological and geophysical evidence suggests that parts of the Northern Carnarvon Basin (NCB) were exhumed. However, the spatio-temporal distribution of this exhumation and its impacts on the region's resource distribution is yet to be established. We analysed porosity, sonic logs and paleotemperature data from 210 boreholes using compaction and vitrinite reflectance analysis techniques to constrain the spatio-temporal distribution of exhumation and its correlation to the region's resource distribution. We discovered that multiple exhumation phases occurred since the Late Triassic. The highest exhumation occurred in the Mid Jurassic, where approximately 1000 m of sediment thickness was eroded from the western and southwestern parts of the Exmouth Plateau. Exhumation at this time spatially correlates with magmatic intrusions in the NCB, and paleotemperature profiles indicate that sediments in the exhumed regions were partly heated by increased paleo-basal heat flow, suggesting that magmatic bodies through underplating and crustal intrusions were the primary causes of this exhumation. Another significant exhumation phase occurred in the Early Cretaceous, removing ∼800 m of sediments from the Exmouth, Barrow and Dampier Sub-basins. Depth-dependent extension events and magmatic intrusions possibly caused exhumation at this time. Significant exhumation ceased in the Early Cretaceous, and our study has established that the exhumation episodes in the NCB were part of large-scale regional events that affected the Australian Northwest Shelf. Combining results from different estimation techniques increased the spatio-temporal coverage of the analysis and improved the precision and accuracy of exhumation results.

Impacts of riverbed aggradation on groundwater regime in a lowland area

Abstract In this study, the influence of riverbed silting on the groundwater regime in a lowland area was investigated. The study area is situated at the Rye Island (Žitný Ostrov) in Slovakia, along the Gabčíkovo – Topoľníky canal, which is part of the drainage-irrigation canal system constructed in this locality. The Rye Island is an area with very low slope (0.25 10–4) and good climatic conditions for aquatic vegetation, therefore the canals are influenced by intensive silting processes. The spatial and temporal patterns of surface water – groundwater exchange are significantly influenced by the thickness of riverbed sediments and their permeability. The aim of this study was to evaluate the thickness and hydraulic conductivity of bed sediments in the Gabčíkovo – Topoľníky canal and to examine their influence on the groundwater – surface water interaction in the area. The hydraulic conductivity of the sediments was assessed from undisturbed samples by the falling head method. The obtained data were used for numerical simulations of groundwater heads by the TRIWACO model for different drainage and infiltration resistance conditions in the area of interest. The results of this study can support the planning of canal maintenance.

Cenozoic structure of the Zhu 3 Depression in the Pearl River Mouth Basin and its response to the dynamic background of the South China Sea

The Zhu 3 Depression is located at the transition between the South China Sea and the South China Block. The Cenozoic structure reflects the dynamic background of the South China Sea. We used 3D seismic and log data to study the structural evolution of the Zhu 3 Depression. We defined the characteristics of Cenozoic rift deformation and the migration of the depocentre based on fault activity and the distribution of sediment thickness. The spatial and temporal differences in the Cenozoic structure of the Zhu 3 Depression are due to the influence of pre-existing faults and regional stresses. Our results show that there are three populations of faults in the Zhu 3 Depression, striking NE, east–west and NW. The NE-striking faults are mainly large-scale boundary faults. The east–west-striking faults are small and were activated in the later stages of deformation. The NW-striking faults are continuous at depth, but are en echelon and parallel to each other at shallow depths. Based on these results, we suggest that the NE- and NW-striking pre-existing faults divide the Zhu 3 Depression into different structural zones and controlled the differences in deformation in the basin. The regional stress direction, which is controlled by regional plate interactions and the Red River Fault Zone, changed clockwise from NW-trending during the Paleocene to NE-trending during the Mid-Miocene. These findings will contribute to a better understanding of the evolution of the entire South China Sea. Thematic collection: This article is part of the Emerging knowledge on the tectonics of the South China Sea collection available at: https://www.lyellcollection.org/topic/collections/south-china-sea

Cumulative sedimentation hazard map of urban areas subject to hyperconcentrated flash flood: A case study of Suide County in the Wuding River basin, China

AbstractFlash floods can carry substantial sediment, posing significant sedimentation hazards in hilly cities. The sedimentation hazard map can reproduce the sediment thickness and extent of an extreme events scenario, playing an important role in sediment risk management. However, current research primarily focuses on modeling the inundation area and depth of floods, while studying sedimentation hazard caused by flash floods in urban areas remains insufficient. This paper aims to address this gap by utilizing a numerical model that simulates hyperconcentrated flow in hilly urban areas using the two‐dimensional hydro‐sediment‐morphological model to compile the cumulative sedimentation hazard map. The model, built upon the open‐source TELEMAC‐MASCARET framework, incorporates Zhang Hongwu's formula to simulate sediment‐carrying capacity, particularly suitable for hyper‐sediment concentration near the riverbed. This paper uses the data of extreme flash flood events in the Wuding River basin in 2017 to simulate and compile the cumulative sedimentation hazard map. The hazard map delineates the sedimentation hazard extent and level attributable to overbank floodplain sedimentation. Notably, the sediment thickness is highest in areas near the levees on both sides of the Dali River. Moreover, the map illustrates the extent of channel erosion resulting from hyperconcentrated floods, which could jeopardize bank stability.

Groundwater Exploration Within Shallow Depths Around Distinct Litho-Petrological Contact Zones

The world’s increasing demand for freshwater has put a strain on existing reserves in recent years, necessitating further exploration of new groundwater reserves. However, geologic provinces located within the contact boundary of Sedimentary Basin-Basement Complex region are somewhat faced with the problem of shallow sediment thickness. As a result, exploring for groundwater poses a challenge as the sediment thickness might be insufficient to host productive aquifers. Hence, exploring for faults and fractures zones that are embedded on the hard rock underneath the sediment is very essential and complimentary. In this study, high resolution aeromagnetic data over a litho-petrological contact zone of the Middle Benue Trough (MBT) and the Adamawa Massif were analyzed to delineate groundwater exploration areas using a simplified qualitative and quantitative approach. Visual inspection of isolated and aggregate anomalies is performed by analyzing their shapes, dimensions, lateral extents, and discontinuities. Lineaments were then extracted from the residual magnetic map to determine faults, fractures, and joints. A two-source-depth model is indicated by the Spectral Analysis technique used to determine the thickness of the sediments in the area. The deeper source (Z1) has a thickness of 0.9 - 3.6 km, while the shallow source (Z2) has a thickness of 0.1 - 0.8 km. While in some areas the overburden thickness may be sufficient for groundwater exploration, in others it was not, as indicated by the shallow magnetic depth, Z2. Nonetheless, areas of shallow thickness but with aggregates of faults, fractures, and joints were identified as possible locations for groundwater accumulation through magnetic lineament exploration.

Inference of source signatures of merchant ships in shallow ocean environmentsa).

An ocean acoustics experiment in 2017 near a shipping lane on the New England continental shelf in about 75 m of water provided an opportunity to evaluate a methodology to extract source signatures of merchant ships in a bottom-limited environment. The data of interest are the received acoustic levels during approximately 20 min time intervals centered at the closest position of approach (CPA) time for each channel on two 16-element vertical line arrays. At the CPA ranges, the received levels exhibit a frequency-dependent peak and null structure, which possesses information about the geophysical properties of the seabed, such as the porosity and sediment thickness, and the characterization of the source, such as an effective source depth. The modeled seabed is represented by two sediment layers, parameterized with the viscous grain shearing (VGS) model, which satisfies causality, over a fixed deep layered structure. Inferred estimates of the implicit source levels require averaging an error function over the full 20 min time intervals. Within the 200-700 Hz band, the Wales-Heitmeyer model captures the inferred frequency dependence of the source levels.

Calibration of horizontal seismic design spectra of offshore earthquake ground motions by applying improved differential evolution method

Design Response Spectrum (DRS) serves as a crucial cornerstone for seismic design, encapsulating the characteristic parameters of the design spectrum. To investigate the offshore design spectrum required for seismic design of offshore structures, improved differential evolution (DE) methods are first used to determine the characteristic parameters of the horizontal design spectra of offshore ground motions. In assessing the stability and precision of the chosen improved DE methods, we scrutinize the applicability of various renowned DE methodologies in the calibration process. Our findings indicate that differential evolution with global and local neighborhoods (GLDE) exhibits superior stability and accuracy compared to other DE methods, rendering it suitable for calibrating the design spectrum to delineate optimal characteristic parameters of the DRS. Subsequently, the properties of the offshore horizontal design spectra are investigated by applying the recommended GLDE to 6756 seismic data collected from three offshore regions. The site conditions of water depth and sedimentary thickness of the ocean-bottom stations were investigated by the finely divided seismic data with classified earthquake types. The ocean-bottom stations are classified as A-D by the water depth and sedimentary thickness, and the offshore DRS for four site classes with different earthquake types are presented in this study. According to the observed offshore response spectra with different water depth groups, it is found that the water depth should be considered in the offshore design spectra mainly due to the differences between buried and unburied seafloor stations, and the spectral values with deep water layer are considerably larger than those with shallow water. Offshore DRSs are compared with onshore DRSs of seismic codes and different peak ground acceleration (PGA) groups are presented in this study.

Correction to: Mapping of geological structures and sediment thickness from analysis of aeromagnetic data over the Obudu Basement Complex of Nigeria

Correction to: Mapping of geological structures and sediment thickness from analysis of aeromagnetic data over the Obudu Basement Complex of Nigeria

Sedimentary facies and mineral provenance of Upper Triassic sandstones offshore Kvitøya, Svalbard: implications for palaeogeographic interpretations in the northern Barents Shelf area

Upper Triassic (Carnian) sandstones of the De Geerdalen Formation cored south of the island of Kvitøya (80°N), north-easternmost Svalbard, are described in terms of sedimentary facies and petrography and compared regionally in the northern Barents Shelf. The succession off Kvitøya is characterized by its great thickness and is dominated by deltaic deposits with high sand content of lithic–feldspathic compositions. Comparison of sediment facies and sandstone compositions with adjacent areas suggest that the succession off Kvitøya is part of a larger delta system with its main sediment source from the east. The delta sedimentation was terminated by marine transgression in the earliest Norian. The sandstone compositions off Kvitøya differ from nearby locations by the higher content of cherty rock fragments and reworked volcanic debris in the Kvitøya sandstone, which is most distinct in the lower part of the succession. Provenance signatures are investigated by mineral–chemical analysis of detrital feldspars, rock fragments, garnet and Cr-spinel, characterizing a wide variety of igneous, metamorphic and sedimentary terranes, including palaeo-Urals and areas farther to east. Additional, more proximal sediment source areas may also have existed that could explain the increased sediment thickness and the mineralogical immature sandstone compositions of the Carnian sediments off Kvitøya.

Time and mechanism of formation of mud volcanoes in the South Caspian basin according to seismic data

There are 277 mud volcanoes located in the South Caspian basin and the Apsheron-Gobustan periclinal trough and the Lower Kura depression that open to it. On land there are 177 mud volcanic structures, from which steam, water, mud, gas and oil are constantly erupting at approximately 500 points. Mud volcanoes located in the Caspian Sea (about 100) also play an important role in the degassing (mainly methane) of the sedimentary layer of the earth’s crust. Based on the results of geological and geophysical research, it has been established that almost all mud volcanoes discovered on the territory of Azerbaijan are associated with hydrocarbon deposits. In the Oligocene – Lower Miocene, mainly in the Maikop time, a thick clayey stratum (Maikop) - source rocks – accumulated in the South Caspian depression. At the end of the Miocene, there was a sharp drop in erosion levels in this basin. As a result of the activity of several river systems (Pleo-Volga, Paleo-Kura, Paleo-Pirsagat, Paleo-Amu Darya) over the Maikop clay deposits, a system of alternating layers of low and high density was formed with a total thickness of sediments of eight kilometers or more. The Maikop clay formations buried under this multilayer system, possessing the properties of a non-Newtonian fluid, created intrusions in the overlying environment, which at subsequent stages of development of the sedimentation basin became eruptive channels of mud volcanoes. The results of research using 3D seismic data made it possible to identify the mechanism of formation of mud volcanism, which, according to the authors, is associated with the migration of hydrocarbons from oil and gas source rocks to reservoir rocks, with the channels of mud volcanoes being considered as the main migration routes. Geological interpretation of seismic data shows that the formation of mud volcanoes in the South Caspian basin supposedly began from the beginning of the Miocene and continued in parallel during subsequent periods of development of the sedimentation basin. As the thickness of the sedimentary cover increased, mud volcanoes periodically became active.

Spatial and temporal variation in long-term sediment accumulation in a back-barrier salt marsh

In situ persistence of salt marshes in the face of sea-level rise relies on their ability to maintain substrate elevation through sufficient vertical accretion of sediment. However, sedimentation rates in salt marshes vary spatially and temporally, which complicates the assessment of their ability to keep up with sea-level rise. Here, we explore the spatial and temporal variation in sediment accumulation in a single back-barrier salt marsh site. Using one-time in situ measurements at the landscape scale, we obtained synoptic information of elevation and sediment thickness over the entire salt marsh in a chronosequence over centuries. Repeated measurements along short elevation transects (0.3–0.9 m +MHT) revealed decadal changes, complementing the broader marsh data with detailed information on elevation, thickness of the marsh deposits and accumulation rates. Thickness of the deposits was largely related to the elevation gradient: the sediment layer was thinner at the higher marsh (near the dunes and far away from the intertidal flats), and thicker at the lower marsh (near the intertidal flats). Moreover, the thickness of the layer increased with salt marsh age along the chronosequence, and age accounted for 72 % of variability in sediment accumulation. The rate of sediment accumulation was higher than the local rate of sea-level rise in the younger marsh, whereas it was equal to the rate of sea-level rise in the older marsh. In the older salt marsh, sediment accumulation was lower, possibly due to autocompaction in the thicker, older layers. Both at the landscape scale and along short elevation transects within individual drainage basins, sediment accumulation decreased with distance to sediment supply routes. However, their relative importance depended on the scale of observation. Distance to creeks accounted for 17 % of the variability in sediment accumulation at the landscape scale, compared to 4 % at the smaller scale. Similarly, the influence of distance to intertidal flats varied from 1 % at the landscape scale to 13 % at the smaller scale. Our main findings indicate that lower-elevation older marshes and higher-elevation younger marshes far away from sediment sources are at risk of not keeping pace with the local rate of sea-level rise and are potentially vulnerable to increased flooding.

Sounding out horse mussel sediment thickness: an integrated data approach

Horse mussel beds are listed as a Priority Marine Feature (PMF) in Scotland for their influence in the creation of diverse benthic communities and provision of ecosystem services. In Scotland, horse mussel beds are also recognised for their importance in marine carbon sequestration. Unfortunately, there is a paucity of data on horse mussel bed carbon stocks and sediment thickness. There are also significant financial and logistical barriers with which to address these data gaps. This study investigates the robustness of Sub-Bottom Profiling (SBP) as a cost-effective method to quantify horse mussel sediment thickness across a landscape. Integrating SBP and Drop-Down Video (DDV) data, this study also details an integrated approach to investigate the links between horse mussel habitat condition and sediment thickness. With the addition of abiotic data, this study uses Structural Equation Modelling (SEM) to elucidate key relationships between abiotic factors, biotic variables, and sediment thickness. There is a significant positive correlation of horse mussel habitat condition and sediment thickness. Average horse mussel total sediment thickness, across all measures of habitat condition, was 1.37m. These findings substantially increase previous estimates of horse mussel sediment thickness, and potential value to climate change mitigation through blue carbon frameworks. This study highlights the importance of both abiotic and biotic factors on marine carbon sediment quantification.

THE STRUCTURE IF THE QUATERNARY SHEET IN THE PETROZAVODSK BAY OF LAKE ONEGA ACCORDING TO SEISMOACOUSTICS

We present the results of seismoacoustic works carried out in the Petrozavodsk Bay of Lake Onega, the goal was to study a near-surface geological section. Seismoacoustic investigations were presented by the reflection method in a modification ultra-high resolution seismic survey. During the office work, the data were processed and interpreted, in which five seismic complexes were identified and correlated with the information about the lake structure. The result of the work was the maps construction of all five identified reflecting horizons covering all Petrozavodsk Bay. It was possible to estimate the sediments abundance and thickness corresponding to each complex.

Exploring inferred geomorphological sediment thickness as a new site proxy to predict ground-shaking amplification at regional scale: application to Europe and eastern Türkiye

Abstract. To test whether a globally inferred sediment thickness value from geomorphological studies can be used as a proxy to predict earthquake site amplification, we derive site-amplification models from the relation between empirical amplification for sites in Europe and Türkiye and the geomorphological sediment thickness. The new site-amplification predictions are then compared to predictions from site-amplification models derived using the traditional site proxies, VS30 inferred from slope, slope itself, and geological era and slope combined. The ability of each proxy to capture the site amplification is evaluated based on the reduction in site-to-site variability caused by each proxy. The results show that the highest reduction is caused by geological era and slope combined, while the geomorphological sediment thickness shows a slightly larger or equal reduction in site-to-site variability as inferred VS30 and slope. We therefore argue that including geology and geomorphology in site-amplification modelling on regional scale can give an important added value and that globally or regionally inferred models for soil and sediment thickness from fields beyond engineering seismology can have a great potential in regional seismic hazard and risk assessments. Furthermore, the differences between the site-amplification maps derived from different proxies capture the epistemic uncertainty of site-amplification modelling. While the different proxies predict similar features on a large scale, local differences can be large. This shows that using only one proxy when predicting site amplification does not capture the full epistemic uncertainty, which is demonstrated by looking into detail on the site-amplification maps predicted for eastern Türkiye and Syria, where the devastating Kahramanmaraş earthquake sequence occurred in February 2023.

Deciphering the sedimentary imprint of tsunamis and storms in the Lesser Antilles (Saint Martin): A 3500-year record in a coastal lagoon

Our study aims at identifying and characterizing tsunami and storm deposits by combining sedimentological, geochemical, and radiocarbon dating analyses. X-ray computed microtomography (micro-CT) is used to characterize the sedimentary fabric of the deposits. The study was conducted on a transect of oriented short sediment cores from a small lagoon at Saint Martin (Lesser Antilles), where both event types could be observed. Using micro-CT data of different event deposits, we were able to derive paleo-flow orientations and directions and compare the dynamics of the events to their environmental setting. The results showed that seven extreme-wave event (EWE) deposits, originating from hurricanes or tsunamis, had the highest Ca/Fe ratios, distinctive for marine input to the lagoon, and were characterized by coarser grains. Some also presented above the sandy part of the event deposit a siliciclastic geochemical signal rich in fine siliciclastic sediment with organic matter, which might correspond to backwash deposits. The thickness and frequency of these sandy EWE deposits in our cores decreased from the proximal to the distal zone in the center of the lagoon, with high lateral and vertical variability in sediment thickness. The youngest EWE deposit relates to the unprecedentedly powerful Category 5 Hurricane Irma, which barely left a fine layer of sand in the lagoon, demonstrating the lagoon's resilience to tropical storms. Five of the six other EWEs, characterized by micro-CT-derived sedimentary fabric, likely originated from tsunamis and occurred over the last 3500 years BP with a recurrence interval of 300–400 years, based on radiocarbon dating. A well-studied tsunami deposit from Pre-Columbian times (∼400 years BP) showed oriented sedimentary fabric, which could be used for paleo-flow reconstruction and EWE-type distinction in the Lesser Antilles, suggesting the use of an additional method for identifying EWE-type deposits in the geological record rather than other methods alone. To develop effective strategies for mitigating the natural hazards faced by communities in the coastal areas of the Lesser Antilles, it will be crucial to examine EWEs in the Caribbean region beyond historic times successfully.

Horizontal ground-motion model for subduction slab earthquakes using offshore ground motions in the Japan Trench area

The S-net, a large-scale network of permanent ocean-bottom seismographs in the Japan Trench, consisting of 150 stations, has recorded a large number of offshore ground motions that can be used to establish the empirical attenuation relationship of offshore ground motions in the subduction zone. Due to the different attenuation characteristics among different earthquake types in subduction zones, the earthquakes are classified into four tectonic types of subduction slab, interface, shallow crustal, and upper-mantle earthquakes according to the classification scheme of Zhao et al. (2015). Predictive models and uncertainties in offshore ground motions were investigated for different earthquake types. This study establishes a horizontal offshore subduction slab earthquake ground-motion model (GMM) and compares the offshore and onshore slab earthquake GMMs. As the site conditions at ocean-bottom stations are different from those at land stations, the effects of water depth and sediment thickness are taken into account in the offshore slab earthquake GMM. Due to differences in the burial methods of ocean-bottom stations, stations were divided into buried and unburied to investigate the effects of the burial methods. Therefore, regression analysis was used to propose an offshore slab earthquake GMM considering the magnitude, focal depth, distance, water depth, sediment thickness, and burial method. By separating the within-event residuals, the single-station standard deviation is presented. Compared to the onshore GMM, the predicted spectra of the offshore GMM are significantly larger at long periods. For the attenuation rate, the offshore attenuation rate is lower than that of the onshore attenuation rate for short periods, but is basically consistent with the onshore attenuation rate for long periods. The proposed GMM can be used to predict the offshore ground motions for slab earthquakes with rupture distances less than 300 km, focal depths less than 110 km, and moment magnitude between 4 and 7.4.

Influence of gravity on methane hydrate dissociation characteristics by depressurization in marine hydrate reservoirs

Since the thickness of hydrate-bearing sediments is usually tens to hundreds of meters, the influence of gravity on hydrate dissociation characteristics can no longer be ignored. In this study, experimental simulations were designed to investigate the influences of gravity on hydrate dissociation characteristics by depressurization. The experimental results showed that for sealed water-rich marine hydrate reservoir, gravity could cause more water to migrate downward. The free pore space in sediments was thus rapidly vacated, resulting in faster hydrate dissociation and gas production. The average gas production rate increased from 0.018 mol/min to 0.038 mol/min. In addition, the free pore space at the top of the sediments was vacated first. Hydrate dissociation there was thus faster. For unsealed water-rich marine hydrate reservoir, gravity could not only cause more and faster migration of water, but also more easily cause the formation of dominant channels for water migration, leading to an effective reduction in reservoir pressure. The cumulative water yield increased from 1380g to 8950g. Hydrate thus dissociated. Hydrate dissociated faster at positions close to the mining well where the pressure was lower. These findings can provide deep insights into the spatial evolution of multiphase flows during the exploitation of natural gas hydrates.

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.