Review: Saltwater intrusion in fractured crystalline bedrock

During the past few years, the number of regional and national assessments of groundwater quality in regard to saltwater intrusion in coastal aquifers has increased steadily. However, most of the international literature on saltwater intrusion is focused on coastal plains with aquifers in unconsolidated material. Case studies, modelling approaches and parameter studies dealing with saltwater intrusion in those systems are abundant. While the hydrogeology of fractured rock has been intensively studied with both modelling approaches and parameter studies—mainly in relation to deep-laying fractured crystalline bedrock as potential waste repositories—case studies on saltwater intrusion in shallow fractured rocks are still an exception. This review summarizes the actual knowledge on saltwater intrusion in fractured crystalline rock. In combination with short overviews of the processes of saltwater intrusion, flow in fractured systems and the genesis of these systems, the review highlights the importance of the fracture systems and its specific characteristics. Fracture properties are a direct consequence of the geological history as well as the current situation of the coastal area. A holistic assessment of water quality in coastal areas hosting fractured crystalline bedrock therefore requires the combination of different approaches in order to investigate the impact of saltwater intrusion through the fractured system.

Development characteristic and formation mechanism analysis of collapse sinkholes in Wugaishan town, Chenzhou city, Hunan province, China

The cover collapse sinkholes occurred and concentrated in Wugaishan town, Chen zhou city since 1996. The results are combined with results of site investigation, geophysical prospecting and in situ groundwater monitoring data, allowing the development characteristics and formation mechanism of surficial collapse incidents to be summarized. Collapse sinkholes are significantly active in recent years and mostly develop in the rainy season ranging from April to June and generally show a zonal distribution along the topography of study area from SW to NE. 92.31 % of total collapse events occurred in the thickness of overburden material ranged from 0 to 15 m, which indicated that overlying material less than 15 m was easier to collapse. The results show that collapse sinkholes have strong relationship with characteristic of overburden material, which sharply decrease in internal physical and mechanical property of bottom layer. Furthermore, substantial cavities formed within bedrock are the best transport channels and storage spaces for the unconsolidated material, especially under the condition of dynamic undulation of groundwater level. The formation mechanism of collapse sinkhole is divided into three types: infiltration erosion, coupling air implosion with vacuum cavitation and saturation erosion. Each formation mechanism is related to changes of groundwater level. When groundwater level rose above the soil-bedrock interface, saturated subsoil were easier to disintegrate into small particles and migrate downward as the vertical seepage of groundwater. The hydraulic gradient increased and became the predominant factor for the development of soil cavity as groundwater level dropped below the soil-bedrock interface. Moreover, when groundwater level sharply surged up at the relative sealed environment, the upward erosion roof of cavity would be more likely to collapse by the entrapped air blasting.

The Equilibrium Concept, or…(Mis)concept in Beaches

Beaches, as deposits of unconsolidated material at the land/water interface, are open systems where input and output items constitute the sediment budget. Beach evolution depends on the difference between the input/output to the system; if positive the beach advances, if negative the beach retreats. Is it possible that this difference is zero and the beach is stable? The various processes responsible for sediment input and output in any beach system are here considered by taking examples from the literature. Results show that this can involve movement of a volume of sediments ranging from few, to over a million cubic meters per year, with figures continuously changing so that the statistical possibility for the budget being equal can be considered zero. This can be attributed to the fact that very few processes are feedback-regulated, which is the only possibility for a natural system to be in equilibrium. Usage of the term “beach equilibrium” must be reconsidered and used with great caution.

GEOPHYSICAL INVESTIGATIONS OF A POTENTIAL LANDSLIDE AREA IN MAYOON, HUNZA DISTRICT, GILGIT-BALTISTAN, PAKISTAN

The Mayoon landslide in the Hunza District is a slowly developed, non-catastrophic landslide that has gained its importance in the last few years after its rapid activation and fast slip rate. The area is characterized by high earthquake hazards (zone 3 with a peak ground acceleration value of 2.4–3.2 m/s2) by the Building Code of Pakistan due to frequent earth quakes. The past high earthquake activity in the area has displaced the foliated rocks towards the south and is responsible for opening the bedrock joints. The head and body of the landslide are covered by unconsolidated material and have fractures of varying lengths and widths. The non-invasive geophysical techniques, including Ground Penetrating Radar (GPR) and Electrical Resistivity Soundings (ERS), are deployed to evaluate the Mayoon landslide subsurface. The subsurface is interpreted into a two-layer model. Bright reflectors and highly variable resistivity characterize the top layer (Layer-1). This layer is associated with a loose, highly heterogeneous, fragmented material deposited under glacial settings over the existing bedrock. Hyperbolic reflections and intermediate resistivity characterize the bottom layer (Layer-2). This layer is associated with foliated metamorphic bedrock. The hyperbolic reflections show faults/fractures within the bedrock. The extension of these fractures/faults with depth is uncertain due to decay in the GPR signal with depth. The intermediate resistivity shows the bedrock is weathered and foliated. Reflections within Layer-1 have disrupted directly above the fractures/faults suggesting a possible movement. A bright reflection between the two layers highlights the presence of the debonded surface. Loose material within Layer-1 coupled with debonding possesses a significant hazard to generate a landslide under unfavourable conditions, such as an intense rainstorm or earthquake activity.

Lithologic, geomorphic, and permafrost controls on recent landsliding in the Alaska Range

Because landslide regimes are likely to change in response to climate change in upcoming decades, the need for mechanistic understanding of landslide initiation and up-to-date landslide inventory data is greater than ever. We conducted surficial geologic mapping and compiled a comprehensive landslide inventory of the Denali National Park road corridor to identify geologic and geomorphic controls on landslide initiation in the Alaska Range. The supplemental geologic map refines and improves the resolution of mapping in the study area and adds emphasis on surficial units, distinguishing multiple glacial deposits, hillslope deposits, landslides, and alluvial units that were previously grouped. Results indicate that slope angle, lithology, and thawing ice-rich permafrost exert first-order controls on landslide occurrence. The majority (84%) of inventoried landslides are <0.01 km2 in area and occur most frequently on slopes with a bimodal distribution of slope angles with peaks at 18° and 28°. Of the 85 mapped landslides, a disproportionate number occurred in unconsolidated sediments and in felsic volcanic rocks. Weathering of feldspar within volcanic rocks and subsequent interactions with groundwater produced clay minerals that promote landslide initiation by impeding subsurface conductivity and reducing shear strength. Landslides also preferentially initiated within permafrost, where modeled mean decadal ground temperature is −0.2 ± 0.04 °C on average, and active layer thickness is ∼1 m. Landslides that initiated within permafrost occurred on slope angles ∼7° lower than landslides on seasonally thawed hillslopes. The bimodal distribution of slope angles indicates that there are two primary drivers of landslide failure within discontinuous permafrost zones: (1) atmospheric events (snowmelt or rainfall) that saturate the subsurface, as is commonly observed in temperate settings, and (2) shallow-angle landslides (<20° slopes) in permafrost demonstrate that permafrost and ice thaw are also important triggering mechanisms in the study region. Melting permafrost reduces substrate shear strength by lowering cohesion and friction along ice boundaries. Increased permafrost degradation associated with climate change brings heightened focus to low-angle slopes regionally as well as in high-latitude areas worldwide. Areas normally considered of low landslide potential will be more susceptible to shallow-angle landslides in the future. Our landslide inventory and analyses also suggest that landslides throughout the Alaska Range and similar climatic zones are most likely to occur where low-cohesion unconsolidated material is available or where alteration of volcanic rocks produces sufficient clay content to reduce rock and/or sediment strength. Permafrost thaw is likely to exacerbate slope instability in these materials and expand areas impacted by landslides.

Assessment of the Dnieper Alluvial Riverbed Stability Affected by Intervention Discharge Downstream of Kaniv Dam

Along the middle reaches of the Dnieper River in central Ukraine, braided riverbeds with many islands have developed in alluvial valleys. In the 1970s, six dams were commissioned, and respective monitoring infrastructure was installed. Riverbanks and valley floors composed of unconsolidated material have much lower bank strengths and are susceptible to fluvial erosion and bank collapse, particularly during the release of high flow volumes from hydropower dams. The regulation of the Dnieper River along a cascade of storage reservoirs caused significant changes in its active river channel and hydrological regime. In order to estimate channel stability downstream of the Kaniv reservoir, we conducted an analysis of the hydraulic conditions in terms of changes in flow velocity and propagation of waves caused by intervention water discharges from the Kaniv Hydroelectric Power Plant (HPP). In this paper, we assess the hydromorphological parameters of the studied river reach as well as the characteristics of the related erosion and deposition zones. Therefore, a monitoring framework for channel processes (MCP) downstream of the Kaniv HPP was installed. The analysis of the intervention discharge parameters was conducted based on measurements from July 2015. Channel stability was expressed by the following factors: Lohtin’s number (L), Makkaveev’s (Kc) factor of stability, and a complex index of stability (Mx) by Grishanin. This study shows that the velocity of artificial wave propagation may reach a speed of up to 74.4 km·h−1. The wave propagates for a distance of approx. 45 km within 65 min at a mean velocity of 37.4 km·h−1. The L, Kc, and Mx indicators used in this work showed that when water discharge increased (e.g., during typical peak-capacity operation), the channel becomes unstable and sediments are subject to erosion processes. The riverbed stability indicators clearly illustrate that an increase in parameter values is not dependent on the distance to the dam. The results are valuable for sustainable sediment management at catchment scale and hence, directly applicable in water management.

Geomorphological properties of the island of Hvar beaches (Croatia, Eastern Adriatic Coast)

<p>Beaches are sedimentary forms, situated at the land-sea contact formed of unconsolidated material which can range in size from sand up to cobbles and boulders. They are one of the most dynamic coastal forms and particularly sensitive to changes (natural and/or anthropogenic).</p><p>Geomorphological properties of the island of Hvar beaches were analyzed by means of field mapping, ortho-photo and GIS analysis. All the beaches have been mapped and measured at four different points in time by means of Ortho-photos from State Geodetic Administration (2011, 2014 and 2017) and HERE maps (2019). Larger beaches have been also measured in the field with a GPS receiver (from July 2018 until July 2019). GIS and statistical calculations and visualization were done in ArcGIS 10.4 software.</p><p>Hvar is the longest Croatian island with a length of 67,8 km, and the fourth in size with a surface of 297,4 km<sup>2</sup>. Along its 254 km long coastline 247 beaches have been mapped which make up 3,8 % of total coastal length. The beaches are rather small relating generally to pocket beaches. Only 14,7 % of beaches are larger than 500 m<sup>2</sup> and 59,95 % are smaller than 200 m<sup>2</sup>. According to the sediment size gravely beaches predominate with 95,5 %, while only 4,5 % relates to sand beaches.</p><p>This study revealed that on the island of Hvar four major morphological types of beaches can be distinguished: beaches formed in fan material at the gully mouth (82,6 %), beaches under the cliff (9,3 %), beaches formed in Aeolian deposits (4,45 %) and artificial or anthropogenic beaches (2,4 %). 1,2 % are undefined. The majority of beaches, 75%, are today under anthropogenic impact while only 25% is completely natural.</p><p>Along the eastern Adriatic coast most of the beaches are formed in torrential material derived from the land accumulated at the gully mouth. Here we revealed that this is also the case on the Island of Hvar (82,6 %). Those beaches are parts of a larger geomorphological system which links the backward drainage basin with the beach. Consequently, here we test if the surface of the beaches correlates with the surface of the drainage basins. Taking into account all the beaches of that morphological type (204 beaches) the correlation revealed to be rather low (r<sup>2</sup>=0,37). However, taking into account only the beaches without any anthropogenic impacts the correlation becomes more significant (r<sup>2</sup>=0,64). This probably points to the disturbing effects of the anthropogenic activity on beaches sediment budget of the island of Hvar.</p><p>A part of this research was made with the support of the Croatian Science foundation (HRZZ-IP-2019-04-9445).</p>

Sediment Supply and Hydrogeological Hazard in the Quebrada De Humahuaca (Province of Jujuy, Northwestern Argentina)—Rio Huasamayo and Tilcara Area

This paper describes the hydrogeological hazard in a reach of Quebrada de Humahuaca, (Upper Valley of Rio Grande de Jujuy, in the Argentine Andes), elected a World Heritage Site by UNESCO in 2003. Along the Quebrada, the Rio Huasamayo, flowing into the Rio Grande, formed a large alluvial fan where the village of Tilcara was built. The final reach of Rio Huasamayo is an artificial channel, embanked by unconsolidated material removed from the riverbed. The village is located in an area at a lower elevation with respect to the riverbed, still in aggradation; so it is affected by an evident hydrogeological hazard. The main cause of the riverbed aggradation is the enormous sediment supply from the slopes of the Rio Huasamayo basin. The aim of the paper is to estimate the soil loss on the slopes and the consequent sediment supply to the main stream, identifying the areas of the basin mostly affected by erosion processes that cause the aggradation of the Rio Huasamayo riverbed. In this case, due to the lack of hydrometeorological stations (monitoring rainfalls, temperature, flow rates, etc.), soil loss and sediment supply to the main stream cannot be estimated through the application of commonly used models in the literature (e.g., USLE, RUSLE, USPED). Here the Gavrilovic method (EPM) was applied in combination with the data of the CORINE Project, allowing the estimation of the volume of material exiting from the catchment. So the main supplying areas (sub-basins of the Rio Huasamayo) can be identified where focused interventions for the control of solid transport could be realized, to mitigate the process of riverbed aggradation.

Preliminary Assessment on Riverine Soil Properties of Sg. Lebir, Kelantan, Peninsular Malaysia

Soil is the unconsolidated material in which plants are rooted. Soil texture is determined by the amount of sand, silt and clay in the soil. Understanding soil properties can help to better management and minimize soil disturbance. These soils were collected for soil properties and bulk density along the 55.2km of Sg. Lebir. Soil on Earth commonly contains of various particles such as sand, silt and clay. These particles commonly classes by refer to soil textural triangle. 40g of air-dried sample diluted into 1000ml of deionised water were left until the percentage of sand, silt and clay observed then calculated as percentage. As calculated, the mean percentage of sand contained is 66%, clay 17% and silt 17%. Besides that, it was concluded that most of the soils texture along 55.2km of Sg. Lebir had a sandy loam soil with pH slightly acidic (6.08) with densities the mean is 0.24 gcm-3.

Thermal properties variations in unconsolidated material for very shallow geothermal application (ITER project)

In engineering, agricultural and meteorological project design, sediment thermal properties are highly important parameters, and thermal conductivity plays a fundamental role when dimensioning ground heat exchangers, especially in very shallow geothermal systems. Herein, the first 2 m of depth from surface is of critical importance. However, the heat transfer determination in unconsolidated material is difficult to estimate, as it depends on several factors, including particle size, bulk density, water content, mineralogy composition and ground temperature. The performance of a very shallow geothermal system, as a horizontal collector or heat basket, is strongly correlated to the type of sediment at disposal and rapidly decreases in the case of dry-unsaturated conditions. The available experimental data are often scattered, incomplete and do not fully support thermo-active ground structure modeling. The ITER project, funded by the European Union, contributes to a better knowledge of the relationship between thermal conductivity and water content, required for understanding the very shallow geothermal systems behaviour in saturated and unsaturated conditions. So as to enhance the performance of horizontal geothermal heat exchangers, thermally enhanced backfilling material were tested in the laboratory, and an overview of physical-thermal properties variations under several moisture and load conditions for different mixtures of natural material was here presented.