Subsurface Bedrock Research Articles

Novel method to determine bedrock from a single sensor based seismic reflection signal using Haar wavelet transform

The determination of bedrock depth is crucial across various earth sciences and related fields. Geophysical techniques, notably the continuous wavelet transform, are increasingly employed to map subsurface bedrock structures. This method allows the transformation of 1D time domain seismic reflection signals into a 2D image, providing time–frequency representations for feature extraction. Leveraging the advantages of wavelet theory, particularly the Haar wavelet known for its simplicity and compact support, this study utilizes continuous Haar wavelet local maxima lines to enhance bedrock analysis. Seismic reflection data is acquired using an extremely sensitive geophone as a receiver along with a high dynamic range data logger. The Common Midpoint (CMP) approach, coupled with optimized offset distances, is utilized to ensure robust data quality and signal strength. Simple sources, such as a sledgehammer, generate sound waves for data collection prior to borehole drilling, correlating acquired data with lithogeological information. This research emphasizes the significance of adjusting parameters like sampling frequency, choice of wavelet, and optimal offset distances. The methodology employed, involving a single geophone, a one-channel data recorder, and a basic sledgehammer as a sound source, offers an inexpensive, straightforward, and non-invasive approach. The study demonstrates the effectiveness of Haar wavelet local maxima lines in accurately determining bedrock depth, producing results closely aligned with both predicted and actual depths.

Imaging an Enigmatic Paleovalley with Passive Seismic Methods (Unaweep Canyon, Colorado, United States)

Abstract High-resolution passive seismic imaging of shallow subsurface structures is often challenged by the scarcity of coherent body-wave energy in ambient noise recorded at surface stations. We show that the autocorrelation (AC) of teleseismic P-wave coda extracted from just one month of continuous recording at 5 Hz geophones can overcome this limitation. We apply this method to investigate the longitudinal subsurface bedrock structure of Unaweep Canyon—a paleovalley in western Colorado (United States) with complex evolution. Both fluvial and glacial processes have been proposed to explain the canyon’s genesis and morphology. The teleseismic P-wave coda AC retrieves zero-offset reflections from the shallow (200–500 m depth) basement interface at 120 stations along a 5 km long profile. In addition, we invert interferometrically retrieved surface-wave dispersion for the shear-wave structure of the sedimentary fill. Combined interpretation of these results and other geophysical and well data suggests an overdeepened basement geometry most consistent with glacial processes.

Built-Up Area Extraction from GF-3 SAR Data Based on a Dual-Attention Transformer Model

Built-up area (BA) extraction using synthetic aperture radar (SAR) data has emerged as a potential method in urban research. Currently, typical deep-learning-based BA extractors show high false-alarm rates in the layover areas and subsurface bedrock, which ignore the surrounding information and cannot be directly applied to large-scale BA mapping. To solve the above problems, a novel transformer-based BA extraction framework for SAR images is proposed. Inspired by SegFormer, we designed a BA extractor with multi-level dual-attention transformer encoders. First, the hybrid dilated convolution (HDC) patch-embedding module keeps the surrounding information of the input patches. Second, the channel self-attention module is designed for dual-attention transformer encoders and global modeling. The multi-level structure is employed to produce the coarse-to-fine semantic feature map of BAs. About 1100 scenes of Gaofen-3 (GF-3) data and 200 scenes of Sentinel-1 data were used in the experiment. Compared to UNet, PSPNet, and SegFormer, our model achieved an 85.35% mean intersection over union (mIoU) and 94.75% mean average precision (mAP) on the test set. The proposed framework achieved the best results in both mountainous and plain terrains. The experiments using Sentinel-1 shows that the proposed method has a good generalization ability with different SAR data sources. Finally, the BA map of China for 2020 was obtained with an overall accuracy of about 86%, which shows high consistency with the global urban footprint. The above experiments proved the effectiveness and robustness of the proposed framework in large-scale BA mapping.

Spatial distribution modeling of subsurface bedrock using a developed automated intelligence deep learning procedure: A case study in Sweden

Due to associated uncertainties, modelling the spatial distribution of depth to bedrock (DTB) is an important and challenging concern in many geo-engineering applications. The association between DTB, the safety and economy of design structures implies that generating more precise predictive models can be of vital interest. In the present study, the challenge of applying an optimally predictive three-dimensional (3D) spatial DTB model for an area in Stockholm, Sweden was addressed using an automated intelligent computing design procedure. The process was developed and programmed in both C++ and Python to track their performance in specified tasks and also to cover a wide variety of different internal characteristics and libraries. In comparison to the ordinary Kriging (OK) geostatistical tool, the superiority of the developed automated intelligence system was demonstrated through the analysis of confusion matrices and the ranked accuracies of different statistical errors. The results showed that in the absence of measured data, the intelligence models as a flexible and efficient alternative approach can account for associated uncertainties, thus creating more accurate spatial 3D models and providing an appropriate prediction at any point in the subsurface of the study area. • Two automated deep learning intelligence models through C++ and Python were developed. • 3D pattern of subsurface geo-spatial DTB distributions for Stockholm, Sweden was presented. • The visualised predicted model using C++ showed more accurate results than Python and OK. • Flaws of OK in modelling the limited and heterogeneous distributed data were demonstrated. • The capacity of intelligence models for subsurface DTB characterizing was improved.

Subsurface Topographic Modeling Using Geospatial and Data Driven Algorithm

Infrastructures play an important role in urbanization and economic activities but are vulnerable. Due to unavailability of accurate subsurface infrastructure maps, ensuring the sustainability and resilience often are poorly recognized. In the current paper a 3D topographical predictive model using distributed geospatial data incorporated with evolutionary gene expression programming (GEP) was developed and applied on a concrete-face rockfill dam (CFRD) in Guilan province- northern to generate spatial variation of the subsurface bedrock topography. The compared proficiency of the GEP model with geostatistical ordinary kriging (OK) using different analytical indexes showed 82.53% accuracy performance and 9.61% improvement in precisely labeled data. The achievements imply that the retrieved GEP model efficiently can provide accurate enough prediction and consequently meliorate the visualization insights linking the natural and engineering concerns. Accordingly, the generated subsurface bedrock model dedicates great information on stability of structures and hydrogeological properties, thus adopting appropriate foundations.

Geo-Electric Assessment of Kali Failed Dam Project Aliero, North Western Nigeria

An integrated geophysical investigation using very low frequency electro-magnetic (VLF-EM) and electrical resistivity methods using (VES) technique has been carried out for investigation of possible causes of dam failure in Kali, Aliero Northwestern Nigeria. A total of Three 3 (VLF-EM) traverses and Ten 10 vertical electrical soundings (VES) techniques have been carried out to identify the fracture zones, seepage and possible threat that lead to the collapse of the dam. The interpretation of VLF-EM data is carried out qualitatively using Fraser and Karous-Hjelt filters and the interpretation of sounding data is done using IPI2 win software. VLF-EM data signifies the exact fracture zones founded on higher values of relative current density that coincide to the subsurface fracture. The survey was conducted on an interval of 10 m interval along traverses at South-North direction ranging from 620 - 640 m in length. The plot of filtered real and filtered imaginary identified noticeable fracture zone close to the surface which was later further delineated by (VES) technique. VES soundings were carried out at locations of prominent VLF anomalies presumably typical of basement fracture using AB/2 of 80 m. Results interpreted show layers thickness and their corresponding resistivity, top soil (0.5 m to 3.45 m, 51.9 Ωm to 318 Ωm), underlain by clay to confined clay (4.56 m to 54 m, 7.88 Ωm to 60.2 Ωm) which is underlain by fractured/fresh basement (5.76 m to 73.9 m, 40.3 Ωm to 197 Ωm). The overburden materials underlain are far away to subsurface (bed rock) in most of the area to support the dam. The clay materials which could serve as conduit, exist within the shallow overburden which was not excavated during the construction. Therefore it was concluded that the dam cannot retain its strength to oppose the dam against collapse.

Major flood events recorded in the Holocene sedimentary sequence of the uplifted Ladiko and Makrisia basins near ancient Olympia (western Peloponnese, Greece)

Detailed palaeoenvironmental studies were conducted in the Ladiko and Makrisia basins near the Alpheios River and ancient Olympia (western Peloponnese, Greece) to assess major landscape changes during the Holocene. Previous studies and literature data document that the area experienced crust uplift of minimum 13 m to 30 m since the mid-Holocene. Geological archives were sampled along a vibracore transect connecting the Ladiko and Makrisia basins. Sediment cores were analyzed using sedimento-logical, geochemical and micropalaeontological methods. Geochronological reconstruction of major landscape changes is based on a set of 24 radiocarbon dates. Geophysical studies were carried out using electrical resistivity tomography (ERT) and Direct Push-Electrical Conductivity (DP-EC) measurements to detect stratigraphic changes and subsurface bedrock structures. The stratigraphic record of the uplifted lake basins of Ladiko and Makrisia revealed two major lithostratigraphic units. Unit I, predominantly composed of clay, silt and silty fine sand, reflects prevailing low-energy sedimentary conditions typical of quiescent (fluvio-)limnic waterbodies. Unit II is made out of fine to coarse sand and documents repeated interferences of unit I associated with abrupt and temporary high-energy flood type (= heft) events. We found signals of four different heft events (H1 to H4) showing strong stratigraphic and geochronological consistencies along the vibracore transect. The following age ranges were determined: H1 – between 4360 – 4330 cal BC and 4320 – 4080 cal BC; H2 – be- tween 2830 – 2500 cal BC and 2270 – 2140 cal BC; H3 – between 1220 –1280 cal AD and 1290 –1390 cal AD; H4 – between 1640 –1800 cal AD and 1650 –1800 cal AD. Different hypotheses concerning the characteristics, potential trigger mechanisms and causes of the flood events were tested against the background of strong Holocene crust uplift and using a variety of different methodological approaches: Geomorphological and granulometric aspects, micropalaeontological contexts, geochronological data sets, numerical simulation of flooding events, local tectonic uplift, and the palaeoclimate background were taken into account. We hypothesize that, during the mid-Holocene, the study area was affected by tsunami events, namely between 4360 – 4330 cal BC and 4320 – 4080 cal BC (H1) and between 2830 – 2500 cal BC and 2270 – 2140 cal BC (H2). These ages are very well consistent with the supra-regional and regional tsunami event signal retrieved from many coastal archives in large parts of western Greece. The timing of flood events H1 and H2 is highly consistent with ages of (supra-)regional tectonic events known from literature and is not consistent with increased flood indices of palaeoclimate data available for western Greece. Tsunami inundation scenarios based on numerical simulation are highly consistent with vibracoring and geophysical (ERT, DP-EC) data. In contrast, heft events H3 and H4 are possibly related to phases of increased precipi- tation and flooding activity in the Mediterranean or to land-based geomorphological processes triggered by regional tectonic events (RTE). Neolithic, Chalcolithic as well as Early and Middle Helladic human activities documented at ancient Olympia were most probably affected by tsunami heft events H1 and H2. Sandy deposits of tsunami event H2, covering the prehistorical tumulus, seem to have been used as a higher and dry base to construct the apsidal houses in the center of the later sanctuary at Olympia. The site, already abandoned, must have again been subject to major flood events during the 13/14th cent. AD and the 17–19th cent. AD associated with heft events H3 and H4.

Unravelling the shape and stratigraphy of a glacially-overdeepened valley with reflection seismic: the Lienz Basin (Austria)

We reveal the subsurface bedrock topography and sedimentary succession of one of the deepest glacially-formed basins in the Eastern Alps: the Lienz Basin in the Upper Drau Valley (Tyrol), by means of seismic reflection. A dense source-receiver spacing, supplied by autonomous receivers, and a prestack depth-migration processing scheme were essential to distinguish the various deposits in fine detail, such as slumping, fan delta deposits, and a modified monocline on the basin flank. These details support our interpretation of the seismic stratigraphy that consists of, e.g., subglacial till of last glacial maximum (LGM) age and possibly older, laminated basin fines, and gravel/coarse sand. The maximum depth of the basin is 622 m, at the junction of two major basement faults that are not clearly visible in the seismic reflections. We regard the overdeepening in this longitudinal valley as the result of glacier confluence during the LGM. Subglacial meltwaters utilized the higher erodibility of faulted rocks, as indicated by channel structures. The adverse slope (2.6%) along the valley axis exceeds the gradient ice-surface slope (0.4–0.5%) during the LGM by more than fivefold. We thus suggest this feature is a product of a pre-LGM phase, since adverse slopes greater than ~ 1.2 times the ice surface slope promote the freezing of water in subglacial channels and prevent efficient water flushing of sediments. Integrating other studies allows us to estimate the local overdeepening of the Lienz Basin and that of the whole Upper Drau Valley to be 146 m and 530 m, respectively. At the beginning of lacustrine sedimentation, we estimate the paleo-water depth to be at least 216 m.

Stratigraphy and geothermal assessment of Mesozoic sandstone reservoirs in the Øresund Basin – exemplified by well data and seismic profiles

The Øresund Basin in the transnational area between Sweden and Denmark forms a marginal part of the Danish Basin. The structural outline and stratigraphy of the Mesozoic succession is described, and a novel interpretation and description of the subsurface geology and geothermal potential in the North Sjælland Half-graben is presented. The subsurface bedrock in the basin includes several Mesozoic intervals with potential geothermal sandstone reservoirs. Parts of the succession fulfill specific geological requirements with regard to distribution, composition and quality of the sandstones. A characterisation of these is presently of great interest in the attempt to identify geothermal reservoirs suitable for district heating purposes. The results presented in this paper include for the first time a comprehensive description of the stratigraphic intervals as well as the characteristics of the potential Mesozoic geothermal reservoirs in the Øresund region, including their distribution, composition and physical properties. This is illustrated by seismic cross-sections and well sections. In addition, results from analyses and evaluations of porosity, permeability, formation fluids and temperature are presented. Six potential geothermal reservoirs in the Mesozoic succession are described and assessed. Primary focus is placed on the characteristics of the reservoirs in the Lower Triassic and Rhaetian–Lower Jurassic succession. The study shows that the Mesozoic reservoir sandstones vary considerably with respect to porosity and permeability. Values range between 5–25% for the pre-Rhaetian Triassic sandstones and are commonly >25% for the Rhaetian–Lower Jurassic and Lower Cretaceous sandstones. The corresponding permeability rarely reaches 500 mD for the pre-Rhaetian Triassic reservoirs, while it is commonly above one Darcy for the Rhaetian–Lower Jurassic and the Lower Cretaceous sandstones. The interpreted formation temperatures are 45–50°C at 1500 m, 60–70°C at 2000 m and 70–90°C at 2500 m depth. The combined results provide a geological framework for making site-specific predictions regarding appraisal of viable geothermal projects for district heating purposes in the region as well as reducing the risk of unsuccessful wells.

Use of electromagnetic–terrain conductivity and DC–resistivity profiling techniques for bedrock characterization at the 15th-of-May City extension, Cairo, Egypt

A joint multi-spacing electromagnetic–terrain conductivity meter and DC–resistivity horizontal profiling survey was conducted at the anticipated eastern extensional area of the 15th-of-May City, southeastern Cairo, Egypt. The main objective of the survey was to highlight the applicability, efficiency, and reliability of utilizing such non-invasive surface techniques in a field like geologic mapping, and hence to image both the vertical and lateral electrical resistivity structures of the subsurface bedrock. Consequently, a total of reliable 6 multi-spacing electromagnetic–terrain conductivity meter and 7 DC–resistivity horizontal profiles were carried out between August 2016 and February 2017. All data sets were transformed–inverted extensively and consistently in terms of two-dimensional (2D) electrical resistivity smoothed-earth models. They could be used effectively and inexpensively to interpret the area's bedrock geologic sequence using the encountered consecutive electrically resistive and conductive anomalies. Notably, the encountered subsurface electrical resistivity structures, below all surveying profiles, are correlated well with the mapped geological faults in the field. They even could provide a useful understanding of their faulting fashion. Absolute resistivity values were not necessarily diagnostic, but their vertical and lateral variations could provide more diagnostic information about the layer lateral extensions and thicknesses, and hence suggested reliable geo-electric earth models. The study demonstrated that a detailed multi-spacing electromagnetic–terrain conductivity meter and DC–resistivity horizontal profiling survey can help design an optimal geotechnical investigative program, not only for the whole eastern extensional area of the 15th-of-May City, but also for the other new urban communities within the Egyptian desert.

Investigations of Changjing earth fissures, Jiangyin, Jiangsu, China

Fissuring of the earth began at Changjing, Jiangyin City, Jiangsu Province, in 1995. The earth fissure zone consists of a main fissure with a number of secondary fissures on both sides, which range from 20 to 60 m wide. The main fissure is about 500 m long, trends 15° to 30°, and has a scarp 100–200 mm high. Extensive field investigations, including shallow seismic surveys and geological drilling, have been conducted over 20 years. Deposits of a Quaternary paleoriver at 150–160 m depth curve around a bedrock spur which is situated beneath, and lies parallel to, the fissure at about 90 m. A regional trough has formed due to land subsidence over the paleoriver channel, with Changjing at its eastern edge. The fissures resulted from differential subsidence of the Quaternary deposits as a result of extended over-exploitation of groundwater. Groundwater extraction facilitates compaction of the Quaternary aquifers and aquitards from the increase of the effective in situ stresses, which leads to subsidence. Positive relief of the subsurface bedrock defines the location and formation of the earth fissures. On one hand, differential subsidence causes the Quaternary sediments to rotate over the ridge, leading to tensional cracking. On the other hand, temporal differential in subsiding aggravates the vertical shearing or faulting of the fissure, forming a scarp.

Geological mapping for the urban area of Tarragona

A detailed and systematic geological characterization of the urban area of Tarragona has been conducted under the project ‘1:5000 scale Urban geological map of Catalonia’ of the Institut Cartografic i Geologic de Catalunya. The results of this characterization are organized into: (1) a geological information system that includes all the information acquired; (2) a geological model focused on the identification, characterization and correlation of the geological materials and structures that are present in the area, and (3) a detailed geological map that represents the synthesis of all the collected information. The mapping project integrates in a GIS environment pre-existing cartographic documentation (geological and topographical), core data from compiled boreholes, descriptions of geological outcrops within the urban network and neighbouring areas, physico-chemical characterization of representative samples of geological materials, detailed geological mapping of Quaternary sediments, subsurface bedrock and artificial deposits, and 3D modelling of the main geological surfaces. The geological model is structured in a system of geological units that form a chronostratigraphic point of view are structured in Preneogene (from Mesozoic to Paleogene), Neogene, Quaternary and Anthropocene. The description of the geological units is guided by a systematic procedure. It includes the main lithological and structural features of the units that constitute the geological substratum and represents the conceptual base of the published map sheets of the urban geological map of the Tarragona metropolitan area. These map sheets are composed of a principal map, geological cross sections, and several complementary maps, charts and tables. Regardless of the geological map units, the principal map also represents the main artificial deposits, features related to the geohistorical processes, contours of outcrop areas, information obtained in stations, borehole data, and contour lines of the top of the pre-quaternary basement surface. The most representative complementary maps are the Quaternary map, the subsurface bedrock map, and the isopach map of thickness of superficial and anthropogenic deposits. The map also includes charts and tables of relevant physico-chemical parameters of the geological materials, harmonized downhole lithological columns from selected boreholes, stratigraphic columns, and photographs and figures illustrating the geology of the mapped area and how urbanization has changed the natural environment. The development of systematic urban geological mapping projects, such as the example of Tarragona’s case, which provides valuable resources to address targeted studies related to urban planning, geoengineering works, soil pollution and other important environmental issues that society should deal with in the future.

Floodplain Change in the Glenmalure Valley, South-east Leinster

The Glenmalure Valley is one of several U-shaped valleys which drain the eastern flanks of the Wicklow Mountains, 45km south west of Dublin. Five channel and floodplain reaches are recognised along a 25km stretch of the Glenmalure River on the basis of spatial changes in stream power and sediment load characteristics. The downvalley sequence of floodplain types within the river is atypical of that commonly associated with decreasing slope and sediment calibre downstream. Radiocarbon dating together with the use of anomalous concentrations of heavy metals, which are used as stratigraphic indicators of alluvial deposition post 1800, also reveal temporal changes in the rate of Holocene floodplain development within the valley. Downstream changes in floodplain processes within the valley are explained, therefore, in terms of spatial and temporal changes in the river's ability to adjust its slope, form and bed texture within the constraints of geological and Pleistocene geomorphic inheritance. The presence of a bedrock basin in the valley floor profile, infilled by Pleistocene lacustrine deposits, is believed to have caused a local change in base level leading to a significant reduction in valley slope and local stream power in the middle reaches of the valley. The subsurface bedrock form and the Pleistocene palaeolake also had a significant impact on coarse sediment transport, storage and transfer through the river system and hence on the ability of the Holocene river to regrade its slope and laterally rework Pleistocene valley fills for floodplain construction.

Radiomagnetotelluric mapping of marlstone and limestone in the Silurian bedrock of Gotland

The Silurian of Gotland is composed of a series of stacked carbonate platforms. Airborne electromagnetic measurements show areas with high electrical resistivity, dominated by dense limestone, and areas with lower resistivity, dominated by marlstone. On Gotland, three main regional, up to 80 m thick, limestone-dominated bedrock packages overlie marlstone-dominated strata. The thickness of the limestone and the subsurface relationship between limestone and marlstone was investigated in the area of Ale–Ardre on central Gotland by use of the radiomagnetotelluric (RMT) method. A significant difference in resistivity between the limestone (1000–10 000 Ωm) and the marlstone ( < 100 Ωm) renders favorable RMT conditions. Four profiles, up to 4.3 km long, were measured across the upper part of the Klinteberg Formation and lower part of the Hemse Group (late Wenlock–early Ludlow). The results show a distinct change in resistivity at 50–80 m depth, correlating with the transition between biostromal–biohermal limestone and shelf marlstone–wackestone lithofacies, as observed in boreholes. The measurements were used to create resistivity models that make it possible to map the subsurface lateral relationship between limestone- and marlstone-dominated sequences. The results also provide important information on distribution and bedding, which give supporting information regarding the sequence stratigraphic model for the Klinteberg Fm–Hemse Group in the Ala–Ardre–Sigsarve area. The RMT method utilized in this study proved to be a very useful method for mapping subsurface bedrock in carbonate-dominated geological settings such as those on Gotland.

Erosion of a granite inselberg, Gross Spitzkoppe, Namib Desert

Namibia, with its passive margin setting, long-term tectonic stability, and long-lasting arid climate is a typical granitic landscape characterized by numerous dominating inselbergs. The Namib has been the focus of quantitative studies on the overall rates of bedrock exhumation, escarpment retreat, and sediment generation, transport, and deposition. Results from these studies indicate steady bedrock erosion rates ranging between 1 and 5mmkyr−1 over the last 105–106yr. This rate was determined from samples collected mostly from small inselbergs with low relief. How fast the large and “classic” inselbergs, which dominate the Namib Desert landscape, erode, generate sediment, and contribute to the overall sediment supply has yet to be systematically examined.We document erosion rates of the Gross Spitzkoppe (GS), one of the largest inselbergs in the Namib Desert, Namibia, based on measured concentrations of 10Be in samples collected from bedrock, boulders, and surface grus. The results suggest a slow lowering rate (1–2mmkyr−1), while the overall slope processes deliver sediment to the base of the cliff with 5–6×105atomsg−1 quartz. This concentration is equivalent to an average erosion rate of ~8mmkyr−1, 2–3 times faster than the bedrock lowering rate. Thus, the inselberg has been eroding by cliff retreat. The concentration of cosmogenic isotopes in slope sediment is controlled by exposure time on the slope, weathering of exposed bedrock, and weathering of subsurface bedrock.Sediment continues to accumulate cosmogenic isotopes while being transported on the flat pediments surrounding the inselberg. Within 1–2km from the base of the GS, the concentration doubles to ~10×105atomsg−1 quartz. The increase is achieved by dosing due to exposure and mixing with highly dosed sediment eroded from near subsurface bedrock. From this area, grus is transported through a network of low gradient channels to the Atlantic Ocean, ~100km away. Cosmogenic isotope concentration in sediment increases to ~30×105 atoms of 10Beg−1 quartz, which is 3–4 times higher than that measured in the grus close to the GS and is controlled by the erosion rate of stable bedrock throughout the entire drainage basin. Dosing due to exposure during transport obviously contributes to the overall measured concentration in the sediment as it approaches the ocean.

Locating and characterising groundwater storage areas within an alpine watershed using time‐lapse gravity, GPR and seismic refraction methods

AbstractUnconsolidated sediments in alpine watersheds can store glacier melt and snowmelt as groundwater, which helps sustain flow in mountain rivers during dry periods. However, the amount and distribution of groundwater storage in rugged alpine terrain is not well understood, hindering our ability to predict the rate and timing of groundwater discharge into alpine streams. We show how non‐invasive time‐lapse microgravity surveys can be used to gauge the spatial distribution of groundwater storage changes within a large (ca 1500 × 1000 m) moraine–talus field of the Lake O'Hara alpine watershed of the Canadian Rockies. Additional ground‐penetrating radar (GPR) and seismic refraction surveys provide complementary information on subsurface bedrock topography and reveal the location of a major northwest–southeast trending depression that likely controls groundwater flow to an alpine lake contiguous with the moraine–talus field. Repeat relative gravity measurements made on a network of 80 gravity stations over and around the moraine–talus field during the summers of 2009 and 2010 reveal gravity changes of up to 25 µgal. Although the small gravity changes associated with groundwater flowing out of storage areas are noisy, significant changes are evident on the eastern side of the moraine–talus field. © Her Majesty the Queen in Right of Canada 2012.

Is climate change responsible for changing landslide activity in high mountains?

ABSTRACTClimate change, manifested by an increase in mean, minimum, and maximum temperatures and by more intense rainstorms, is becoming more evident in many regions. An important consequence of these changes may be an increase in landslides in high mountains. More research, however, is necessary to detect changes in landslide magnitude and frequency related to contemporary climate, particularly in alpine regions hosting glaciers, permafrost, and snow. These regions not only are sensitive to changes in both temperature and precipitation, but are also areas in which landslides are ubiquitous even under a stable climate. We analyze a series of catastrophic slope failures that occurred in the mountains of Europe, the Americas, and the Caucasus since the end of the 1990s. We distinguish between rock and ice avalanches, debris flows from de‐glaciated areas, and landslides that involve dynamic interactions with glacial and river processes. Analysis of these events indicates several important controls on slope stability in high mountains, including: the non‐linear response of firn and ice to warming; three‐dimensional warming of subsurface bedrock and its relation to site geology; de‐glaciation accompanied by exposure of new sediment; and combined short‐term effects of precipitation and temperature. Based on several case studies, we propose that the following mechanisms can significantly alter landslide magnitude and frequency, and thus hazard, under warming conditions: (1) positive feedbacks acting on mass movement processes that after an initial climatic stimulus may evolve independently of climate change; (2) threshold behavior and tipping points in geomorphic systems; (3) storage of sediment and ice involving important lag‐time effects. Copyright © 2011 John Wiley &amp; Sons, Ltd.

A deep rock laboratory in the Dellen impact crater

The Deep Rock Lab is a platform to establish a comprehensive subsurface bedrock characterization approach, by integrating site characterization techniques applied from different disciplines of geo-mechanics, geochemistry, hydrogeology, structural geology, lithology and geophysics, with consideration of the effects of coupled geological processes of importance for the understanding of groundwater renewal, continental shield deformations, engineering issues related to geological disposal of nuclear waste and CO2, and geothermal energy retrieval in crystalline rocks. The approach will focus on the physics and chemistry of crystalline rocks and groundwater with down-the-hole measurements of relevant variables, using and developing more efficient geo-scientific site investigation techniques for deep boreholes at a chosen site, and develop more advanced down-the-hole measurements and numerical modelling methods with more advanced inversion algorithms to help integrate data interpretations and object representations. The goal is to develop this platform into a long-term research facility that can be readily used by the scientific community for both subsurface fundamental and engineering-oriented research. Such a platform will be especially important for the education of PhD students for generations to come. The integrated drilling and research facility is suggested to be located at the Dellen site. This site has an impact crater with a large range of expected physical property changes with depth, complex and multiple thermal processes that have affected the bedrock, a favorable infrastructure and local supporting activities, and a large body of existing geo-scientific data.

Erosion of an Ancient Mountain Range, The Great Smoky Mountains, North Carolina and Tennessee

Analysis of 10Be and 26Al in bedrock (n=10), colluvium (n=5 including grain size splits), and alluvial sediments (n=59 including grain size splits), coupled with field observations and GIS analysis, suggest that erosion rates in the Great Smoky Mountains are controlled by subsurface bedrock erosion and diffusive slope processes. The results indicate rapid alluvial transport, minimal alluvial storage, and suggest that most of the cosmogenic nuclide inventory in sediments is accumulated while they are eroding from bedrock and traveling down hill slopes. Spatially homogeneous erosion rates of 25 - 30 mm Ky−1 are calculated throughout the Great Smoky Mountains using measured concentrations of cosmogenic 10Be and 26Al in quartz separated from alluvial sediment. 10Be and 26Al concentrations in sediments collected from headwater tributaries that have no upstream samples (n=18) are consistent with an average erosion rate of 28 ± 8 mm Ky−1, similar to that of the outlet rivers (n=16, 24 ± 6 mm Ky−1), which carry most of the sediment out of the mountain range. Grain-size-specific analysis of 6 alluvial sediment samples shows higher nuclide concentrations in smaller grain sizes than in larger ones. The difference in concentrations arises from the large elevation distribution of the source of the smaller grains compared with the narrow and relatively low source elevation of the large grains. Large sandstone clasts disaggregate into sand-size grains rapidly during weathering and downslope transport; thus, only clasts from the lower parts of slopes reach the streams. 26Al/10Be ratios do not suggest significant burial periods for our samples. However, alluvial samples have lower 26Al/10Be ratios than bedrock and colluvial samples, a trend consistent with a longer integrated cosmic ray exposure history that includes periods of burial during down-slope transport. The results confirm some of the basic ideas embedded in Davis’ geographic cycle model, such as the reduction of relief through slope processes, and of Hack’s dynamic equilibrium model such as the similarity of erosion rates across different lithologies. Comparing cosmogenic nuclide data with other measured and calculated erosion rates for the Appalachians, we conclude that rates of erosion, integrated over varying time periods from decades to a hundred million years are similar, the result of equilibrium between erosion and isostatic uplift in the southern Appalachian Mountains.

Space weathering on Eros: Constraints from albedo and spectral measurements of Psyche crater

Abstract— We present combined multi‐spectral imager (MSI) (0.95 μm) and near‐infrared spectrometer (NIS) (0.8–2.4 μm) observations of Psyche crater on S‐type asteroid 433 Eros obtained by the Near‐Earth Asteroid Rendezvous (NEAR)—Shoemaker spacecraft. At 5.3 km in diameter, Psyche is one of the largest craters on Eros which exhibit distinctive brightness patterns consistent with downslope motion of dark regolith material overlying a substrate of brighter material. At spatial scales of 620 m/ spectrum, Psyche crater wall materials exhibit albedo contrasts of 32–40% at 0.946 μm. Associated spectral variations occur at a much lower level of 4–8% (±2–4%). We report results of scattering model and lunar analogy investigations into several possible causes for these albedo and spectral trends: grain size differences, olivine, pyroxene, and troilite variations, and optical surface maturation. We find that the albedo contrasts in Psyche crater are not consistent with a cause due solely to variations in grain size, olivine, pyroxene or lunar‐like optical maturation. A grain size change sufficient to explain the observed albedo contrasts would result in strong color variations that are not observed. Olivine and pyroxene variations would produce strong band‐correlated variations that are not observed. A simple lunar‐like optical maturation effect would produce strong reddening that is not observed. The contrasts and associated spectral variation trends are most consistent with a combination of enhanced troilite (a dark spectrally neutral component simulating optical effects of shock) and lunar‐like optical maturation. These results suggest that space weathering processes may affect the spectral properties of Eros materials, causing surface exposures to differ optically from subsurface bedrock. However, there are significant spectral differences between Eros' proposed analog meteorites (ordinary chondrites and/or primitive achondrites), and Eros' freshest exposures of subsurface bright materials. After accounting for all differences in the measurement units of our reflectance comparisons, we have found that the bright materials on Eros have reflectance values at 0.946 μm consistent with meteorites, but spectral continua that are much redder than meteorites from 1.5 to 2.4 μm. Most importantly, we calculate that average Eros surface materials are 30–40% darker than meteorites.