Pore Water Research Articles

Experimental study on the influence of cyclic loading frequency on liquefaction characteristics of saturated sand

A series of undrained cyclic torsional shear tests was conducted to investigate the effect of cyclic loading frequency on the liquefaction characteristics of saturated sand using a hollow cylinder apparatus. The test results show that the dilative and contractive tendencies of various saturated sands are not only related to the physical properties of the sand, but are also affected by loading frequency. Under low-frequency loading, the saturated sand has a dilative behaviour, excess pore water pressure fluctuates after initial liquefaction and the soil maintains the ability to resist liquefaction to some extent after the initial liquefaction. The liquefaction mode in terms of the stress–strain relationship generally performs as the ‘cyclic mobility’. However, under high-frequency loading, the saturated sand has a contractive behaviour, and excess pore water pressure generally remains stable after the initial liquefaction. The liquefaction mode in terms of stress–strain relationship generally exhibits as ‘cyclic instability’. The deformation caused by low-frequency loading is significantly larger compared with that caused by high-frequency loading. At higher loading frequencies, the phase transformation stress ratio increases with the increase of loading frequency, and gradually approaches the failure stress ratio.

Investigations of the Effect of Artificial Rainfall on the Pore Water Pressure and Slope Surface Displacement of Loess Slopes

Investigations of the Effect of Artificial Rainfall on the Pore Water Pressure and Slope Surface Displacement of Loess Slopes

Energy-Based Assessment of Liquefaction Behavior of a Non-Plastic Silt Based on Cyclic Triaxial Tests

Earthquakes cause cyclic shear deformations in soil and build-up of excessive pore water pressure as a result of undrained loading, accompanied with rearrangement of soil particles and degradation in stiffness of the soil due to decrease in effective stresses. During loading, the onset of soil liquefaction is defined as a stress state in which the excess pore water pressure is equalized to the total stress. From this point of view, assessment of the pore water pressure development pattern under cyclic loading has been one of the most salient research topics in geotechnical and earthquake engineering. In this study, results of a series of cyclic triaxial tests on non-plastic silt specimens consolidated under 100 kPa effective isotropic consolidation pressure were used to question the modelling ability of pore pressure development models previously proposed for sands. Tests were performed on specimens of 6 different initial relative densities (Dr) ranging between 30-80% and 10 different cyclic stress ratios (CSR). The key parameters of pore water pressure development and shear deformation in the energy-based model used are relative density, cyclic stress ratio and number of cycles. The results revealed that, these energy-based models have a strong potential in evaluation of pore water pressure development pattern of non-plastic silts. Test results also show that the increase in relative density and decrease in CSR causes a ladderlike behavior among pore water pressure and cyclic shear strain, which is relevantly rendered by energy-based models.

Numerical Simulation of Rainfall‐induced Xianchi Reservoir Landslide in Yunyang, Chongqing, China

AbstractA calamitous landslide happened at 22:00 on September 1, 2014 in the Yunyang area of Chongqing City, southwest China, enforcing the evacuation of 508 people and damaging 23 buildings. The landslide volume comprised 1.44 million m3 of material in the source area and 0.4 million m3 of shoveled material. The debris flow runout extended 400 m vertically and 1600 m horizontally. The Xianchi reservoir landslide event has been investigated as follows: (1) samples collected from the main body of landslide were carried out using GCTS ring shear apparatus; (2) the parameters of shear and pore water pressure have been measured; and (3) the post‐failure characteristics of landslide have been analyzed using the numerical simulation method. The excess pore‐water pressure and erosion in the motion path are considered to be the key reasons for the long‐runout motion and the scale‐up of landslides, such as that at Xianchi, were caused by the heavy rainfall. The aim of this paper is to acquired numerical parameters and the basic resistance model, which is beneficial to improve simulation accuracy for hazard assessment for similar to potentially dangerous hillslopes in China and elsewhere.

Vallisneria spiralis L. adaptive capacity improves pore water chemistry and increases potential nitrification in organic polluted sediments

BackgroundMacrophytes may modify benthic biodiversity and biogeochemistry via radial oxygen loss from roots. This condition contrasts sediments anoxia, allows roots respiration, and facilitates aerobic microbial communities and processes in the rhizosphere. Simultaneously, the rhizosphere can stimulate anaerobic microorganisms and processes via exudates or by favoring the build-up of electron acceptors as nitrate. As eutrophication often results in organic enrichment in sediments and large internal nutrients recycling, an interesting research question is to investigate whether plants maintain the capacity to stimulate aerobic or anaerobic microbial communities and processes also under elevated organic pollution.MethodsA manipulative experiment was carried out under laboratory-controlled conditions. Microcosms containing bare sediments and sediments transplanted with the macrophyte Vallisneria spiralis L. were created. The effect of the plant was investigated on sediments with moderate (8%) and elevated (21%) organic matter content, after an acclimatization period of 30 days. Chemical and physical parameters, microbial community composition and the potential rates of nitrification, denitrification and nitrate ammonification were measured at two different depths (0–1 and 1–5 cm) after the acclimatization period to evaluate the role of roots.ResultsVallisneria spiralis grew and assimilated pore water nutrients at the two organic matter levels and vegetated sediments had always nutrient-depleted porewaters as compared to bare sediments. Nitrifying microbes had a lower relative abundance and diversity compared to denitrifying bacteria. However, regardless of the organic content, in vegetated sediments nitrifiers were detected in deeper horizons as compared to bare sediments, where nitrification was confined near the surface. In contrast, potential denitrification rates were not affected by the presence of roots, but probably regulated by the presence of nitrate and by root-dependent nitrification. Potential nitrate ammonification rates were always much lower (< 3%) than potential denitrification rates.ConclusionsVallisneria spiralis affects N-related microbial diversity and biogeochemistry at moderate and elevated organic matter content, smoothing bottom water–pore water chemical gradients and stimulating nitrification and nitrogen loss via denitrification. These results suggest the possibility to deploy V. spiralis as a nature-based solution to counteract eutrophication in freshwater systems impacted by high loads of organic matter, for example, downstream of wastewater treatment plants.

Evaluating the influence of Cynodon dactylon on the wave force and wave erosion in the water-level fluctuation zone of the Three Gorges Reservoir Area

Wave erosion is the main erosion type in the water-level fluctuation zone (WLFZ) of the Three Gorges Reservoir Area (TGRA). Despite vegetation can effectively mitigate wave erosion in the WLFZ, its influence on the wave force and wave erosion remains unclear. Therefore, the wave experiments were conducted under 3 Cynodon dactylon coverage rates (0, 30% and 60%) and 9 wave conditions (3 wave heights of 4, 6 and 8 cm combined with 3 wave periods of 1, 2 and 3 s) to analyse the wave force (expressed as the wave pressure on the slope surface and the pore water pressure in the slope) and wave erosion rate, and the factors influencing wave erosion were identified. The results indicated that the wave pressure, pore water pressure and wave erosion rate increased by 19.14%–104.75%, 16.84%–65.04% and 23.33%–91.64%, respectively, as wave height increases. The wave pressure decreased by 1.50%–31.23% followed by an increase by 22.05% to 87.10% with the increase of wave period, whereas the pore water pressure and wave erosion rate decreased by 28.33%–53.59% and 20.46%–63.59%, respectively. However, these quantities decreased by 2.10%–50.84%, 17.06%–40.23% and 17.28%–82.18%, respectively, with the increase of Cynodon dactylon coverage rate. It was also discovered that the pore water pressure and Cynodon dactylon coverage rate attained the highest positive and negative correlation coefficients with the wave erosion rate, respectively. In addition, pore water pressure accumulation is the most critical influence factor on wave erosion, and Cynodon dactylon could effectively reduce the pore water pressure via its roots, thus improving the slope wave erosion resistance. This study could be useful to understand the mechanism of plants on controlling wave erosion and could provide a scientific reference for wave erosion control and the ecological construction in the WLFZ.

Analytical model for geotextile-enhanced horizontal drain vacuum consolidation of slurries

Horizontal drains are gradually introduced to the vacuum preloading method to improve dredged slurries by adding geotextiles to alleviate the blockage in the consolidation process. This study considers the consolidation of slurries enhanced by the vacuum preloading method with geotextile combined horizontal drains based on a double-layered consolidation model. The model approximates geotextile as a special soil layer possessing an equivalent consolidation factor. An analytical solution of the layered consolidation model is obtained using the Laplace transform and the finite Fourier transform method. The effectiveness of the solution is verified by comparing it with the one-dimensional double-layered consolidation solution and the one-dimensional consolidation with a partially permeable boundary. Through comparison with laboratory experiments, the model shows good fitness with the test results in the literature. The influences of related parameters, including the drain arrangement densities, soil parameters, and geotextile parameters, are discussed on average consolidation degree and pore water pressure. The influence mechanism is explained regarding drainage path and vacuum pressure transfer. Findings demonstrate that geotextile facilitates vacuum transfer and promotes soil consolidation, especially when the smaller density of drains’ deposition and lower soil permeability are applied.

Spatial distribution, sources, and human health risk assessment of elevated nitrate levels in groundwater of an agriculture-dominant coastal area in Hainan Island, China

In order to comprehend the groundwater nitrate contamination in Hainan Island, which is recognized as China's sole “tropical agricultural production base,” the distribution characteristics of nitrate in different aquifers and land-use types were elucidated through hydrochemical and geostatistical analysis methods. Principal component analysis (PCA), hierarchical cluster analysis (HCA), and absolute principal component score-multiple linear regression (APCS-MLR) were employed for qualitative and quantitative analysis of the nitrate source. The results showed that NO3−-high (>88.57 mg/L) groundwater accounted for about 10 % of the total distribution area of the western coastal area in Hainan Island. Moreover, the proportions of NO3−-high groundwater in pore water and fissure water were 10 % and 7 %, respectively. The order of groundwater nitrate concentration in different land-use types was as follows: bare land, cultivated land, construction land, grassland, and woodland. The main sources of elevated nitrate concentrations in pore water were animal feces and nitrogen/potassium fertilizers, accounting for 45.4 % of the overall contribution. These substances infiltrated through the vadose zone, leaching out and mobilizing selenium within this zone, ultimately leading to the characteristic enrichment patterns observed in pore water, including nitrate, arsenic, potassium, and selenium. Furthermore, the occurrence of high nitrate in fissure water was mainly attributed to the leaching of cultured manure, septic tank leakage, and infiltration of domestic sewage, accounting for a contribution rate of 32.9 %. Simultaneously, under oxidizing conditions, selenium became activated and entered groundwater through water flow, thereby exhibiting common enrichment characteristics with nitrate and potassium in fissure water. In addition, the assessment of human health risks revealed that 16 % and 40 % of groundwater samples posed unacceptable non-carcinogenic risks to adults and children respectively, with a higher risk for local children. The study recommends the reasonable application of chemical fertilizers, as well as the strengthening of sewage treatment capacity and standardized disposal of manure and feces, to reduce nitrate contamination of groundwater at source and ensure the safety of drinking water.

Experimental investigation of thawing behavior of saline soils using resistivity method

Abstract Electrical resistivity method has been widely used to study permafrost and to monitor the process of freezing-thawing. However, a thorough understanding of the mechanism of electrical response during thawing is missing. In this study, we investigated the thawing behavior of saline soils in the temperature range from roughly −10 to 15°C considering the effects of soil type and salinity. A total of nine experiments were performed with three soil types (silica sand, sandy soil, and silt) and three salinities (0.01, 0.1, and 1 S m−1). The results show that resistivity variations with temperature can be divided into three stages. In Stage I, tortuosity and unfrozen water content play major roles in the decrease of resistivity. In Stage Ⅱ, which is an isothermal or near isothermal process, resistivity still decreases slightly due to the thawing of residual ice and pore water movement. In Stage III, ionic mobility plays an important impact on decreasing resistivity. In addition, the isothermal process is found to only occur in silica sand that can be explained by latent heat effect. Exponential and linear models linking temperature with resistivity are used to fit the experimental data in Stages I and III. The fitting parameter in different models shows great correlation with soil type and salinity. Furthermore, unfrozen water content below 0°C is also estimated and uncertainty of estimation is analyzed.

From Chingford to Carsington: impact of soil mechanics on dam construction, 1937–1987

In 1937 the investigation of the failure during construction of Chingford number 2 embankment dam brought the concepts of modern soil mechanics to bear on the construction of British embankment dams. Effective stress theory highlighted the crucial significance of pore-water pressure in understanding embankment stability and it became standard practice to monitor pore pressures during embankment construction using twin-tube hydraulic piezometers. By the 1980s, embankment stability might have seemed assured, with embankment design utilising limit equilibrium stability analyses that incorporated strength parameters of fill and foundation materials determined by laboratory testing. However, in 1984 the upstream slope of Carsington dam suffered a major slip in the final stage of construction. Recently developed advanced finite-element analyses were able to demonstrate that progressive failure was a major factor in the slope instability, but the failure was also a reminder of the continuing need for experience-based judgement in design.

River–Spring Connectivity and Hydrogeochemical Processes in a Karst Water System of Northern China: A Case Study of Jinan Spring Catchment

Frequent surface water–groundwater interactions and prevalent anthropogenic inputs make karst water systems vulnerable to human disturbance. As a typical karst region in North China, the Jinan Spring Catchment has become increasingly threatened due to rapid population growth and urban expansion. In this study, the local river–spring interaction and its interference with the hydrogeochemical evolution of groundwater are evaluated based on water stable isotopes and hydrochemistry. Twenty-two karst groundwater, eleven Quaternary pore water, sixteen spring water, and thirty-two surface water samples were collected during low- and high-flow conditions over the course of a year. The isotopic signatures of four different water types display significant differences, reflecting the recharge–discharge relationship of the karst water system. Mountainous springs feature lighter isotopes, whereas urban springs have significantly heavier isotopes. The result of end-member mixing analysis shows that the surface–groundwater interaction varies spatially and temporally within the spring catchment. Urban springs receive considerable replenishment from the surface water, especially after rainy episodes (up to 50%), while mountainous springs show little hydraulic dependence on surface water leakage (4~6%). Local mineral dissolution (including calcite, dolomite, gypsum, and halite), CO2 dissolution/exsolution, and cation exchange are the main hydrogeochemical processes constraining water chemistry in the spring catchment. The deterioration of water quality can be attributed to anthropogenic influences involving the discharge of domestic effluents, agricultural activities, and irrigation return flow. The findings of this work can improve our understanding of the complex karst water system and serve as a reference for sustainable groundwater management in other karst areas of northern China.

Observed soil arching-induced ground deformation and stress redistribution behind braced excavation

This paper presents observed arching-induced ground deformation and stress redistribution behind braced excavation using the top-down construction method. The soil properties around the excavation were determined by laboratory and field tests. The ground deformation, soil displacement vector, strain path, principal strain, maximum shear strain, lateral earth pressure, pore water pressure, and effective stress path are presented based on the measured data. The majority of soil behind the wall is under volumetric expansion, indicating consolidation, creep behavior, or a combination of both. Besides, two periods of increases in pore pressure are observed, due to stress transfer from the lower to the upper parts (i.e., soil arching effect). The deep inward movement of the wall and the nearby soil accounts for the distribution of lateral earth pressure acting on the wall. The soil located behind the area of maximum wall deformation and adjacent to the wall, as well as the soil below the excavation base intersected by the shear plane, is in an active stress state. The lateral earth pressure at 5 m from the left excavation wall showed minimal changes, due to the combined effects of soil arching from lateral excavation and shield tunneling.

The Possibility of Estimating the Permafrost’s Porosity In Situ in the Hydrocarbon Industry and Environment

Global warming firstly influences the permafrost regions where numerous and rich world hydrocarbon deposits are located. Permafrost thawing has caused severe problems in exploring known hydrocarbon deposits and searching for new targets. This process is also dangerous for any industrial and living regions in cold regions. Knowledge of permafrost’s ice and unfrozen water content is critical for predicting permafrost behavior during the water–ice transition. This is especially relevant when ice and permafrost are melting in many regions under the influence of global warming. It is well known that only part of the formation’s pore water turns into ice at 0 °C. After further lowering the temperature, the water phase transition continues, but at gradually decreasing rates. Thus, the porous space is filled with ice and unfrozen water. Laboratory data show that frozen formations’ mechanical, thermal, and rheological properties strongly depend on the moisture content. Hence, porosity and temperature are essential parameters of permafrost. In this paper, it is shown that by combining research in three fields, (1) geophysical exploration, (2) numerical modeling, and (3) temperature logging, it is possible to estimate the porosity of permafrost in situ. Five examples of numerical modeling (where all input parameters are specified) are given to demonstrate the procedure. This investigation is the first attempt to quantitatively analyze permafrost’s porosity in situ.

Early Diagenetic Processes in the Sediments of the Krka River Estuary

To study the processes that govern the post-depositional mobility of metals in the estuarine sediment, five sediment cores were sampled in the Krka River estuary (Croatia). The obtained concentration ranges in the pore water were 0.057–49.7 μM for Fe, 0.310–100 μM for Mn, 0.068–26.8 nM for Co, 0.126–153 nM for Cu, 11.5–2793 nM for Zn, 0.222–31.3 nM for Pb, 4.09–59.4 nM for U, 38.8–2228 nM for Mo, and 0.065–2.29 nM for As. The vertical distribution of metals in the dissolved and solid fraction of the sediment, coupled with other diagenetic tracers (e.g., dissolved sulphide), demonstrate the importance of early diagenetic reactions, in particular Fe and Mn oxyhydroxide and sulphate reduction, for the cycling of metals in the sediment. The redox zonation in the sediment was compressed, and the suboxic zone occurs immediately below the sediment–water interface. The estimated benthic fluxes in the estuary were 5220 kg y−1 for Fe, 27,100 kg y−1 for Mn, 6.00 kg y−1 for Co, 20.5 kg y−1 for Cu, 5.16 kg y−1 for Pb, 111 kg y−1 for Mo, and 87.3 kg y−1 for As. The riverine input was more important than the benthic flux, except in the case of Mn and Fe.

Study on the influence of pore water pressure on shear mechanical properties and fracture surface morphology of sandstone

To further investigate the weakening effect of pore water pressure on intact rock mechanics properties and characteristics of fracture surface after failure, direct shear tests of sandstone were conducted under different pore pressure. A 3D scanner was employed to digitize the morphology of the post-shear fracture surface. The variogram function was applied to quantify the anisotropic characteristics of post-shear fracture surface. The relationship between deformation during shear failure of intact rock and quantitative parameters of fracture surface after shear failure was initially established. It can be found that amplitudes of the sinusoidal surface determine the maximum value of variogram, and period affect lag distance that reach the maximum value of variogram. Test results revealed that the increase of pore pressure has obvious weakening effect on shear strength and deformation of rock. Moreover, the increase of pore pressure makes the shear fracture surface flatter. It can be obtained that both Sillmax and Rangemax are positively related to shear strain, but negatively related to normal strain.

Weathering Intensity, Paleoclimatic, and Progressive Expansion of Bottom-Water Anoxia in the Middle Jurassic Khatatba Formation, Southern Tethys: Geochemical Perspectives

The Jurassic Period was a significant phase of variable organic matter accumulation in paleo-shelf areas of the southern Tethys (Egypt). Reconstructing the paleoredox conditions, paleoclimate, and weathering intensity, along with the role of terrigenous sediment flux and mineralogical maturity, is important for understanding basin infill history and prevalent paleoenvironmental conditions. Here, inorganic geochemical data are presented from the Middle Jurassic Khatatba Formation and two samples from the underlying Ras Qattara and the overlying Masajid formations in the Jana-1x well, Shushan Basin, Western Desert. Twenty-four (24) whole-rock samples were analyzed for their major and trace element composition and carbonate content. The Khatatba Formation represents one of the major hydrocarbon source rocks in the North Western Desert, Egypt. Redox conditions were assessed based on enrichment factors of redox-sensitive elements Mo, V, U, and Co. Results revealed that the Khatatba Formation was deposited under predominant anoxic bottom and pore water conditions, in contrast to the oxic settings that were prevalent during the deposition of the Ras Qattara and Masajid formations. Continental weathering intensity and paleoclimate were reconstructed based on several proxies, such as the chemical index of alteration (CIA), K2O/Rb, Rb/Sr, Ln(Al2O3/Na2O), and Al/K ratios, indicating that the studied succession was deposited during alternating phases between weak and moderate weathering intensity under arid and warm-humid climates, respectively. Periods of enhanced continental weathering were associated with high values of clastic ratios such as Si/Al, Ti/Al, and Zr/Al, suggesting increased terrigenous sediment supply during intensified hydrological cycling. These ratios further provided inferences about the changes in sediment grain size, such as a change from shale to coarse silt- and sand-size fractions.

Numerical Study of Pore Water Pressure in Frozen Soils during Moisture Migration

Frost heaving in soils is a primary cause of engineering failures in cold regions. Although extensive experimental and numerical research has focused on the deformation caused by frost heaving, there is a notable lack of numerical investigations into the critical underlying factor: pore water pressure. This study aimed to experimentally determine changes in soil water content over time at various depths during unidirectional freezing and to model this process using a coupled hydrothermal approach. The agreement between experimental water content outcomes and numerical predictions validates the numerical method’s applicability. Furthermore, by applying the Gibbs free energy equation, we derived a novel equation for calculating the pore water pressure in saturated frozen soil. Utilizing this equation, we developed a numerical model to simulate pore water pressure and water movement in frozen soil, accounting for scenarios with and without ice lens formation and quantifying unfrozen water migration from unfrozen to frozen zones over time. Our findings reveal that pore water pressure decreases as freezing depth increases, reaching near zero at the freezing front. Notably, the presence of an ice lens significantly amplifies pore water pressure—approximately tenfold—compared to scenarios without an ice lens, aligning with existing experimental data. The model also indicates that the cold-end temperature sets the maximum pore water pressure value in freezing soil, with superior performance to Konrad’s model at lower temperatures in the absence of ice lenses. Additionally, as freezing progresses, the rate of water flow from the unfrozen region to the freezing fringe exhibits a fluctuating decline. This study successfully establishes a numerical model for pore water pressure and water flow in frozen soil, confirms its validity through experimental comparison, and introduces an improved formula for pore water pressure calculation, offering a more accurate reflection of the real-world phenomena than previous formulations.

Analytical solution for one‐dimensional thaw consolidation model with double moving boundaries

AbstractA one‐dimensional thaw consolidation model considering the density change from pore ice to pore water is established, and the model describes a special type of moving boundary problem with double moving boundaries. An analytical solution for the model under a time‐varying external load is developed using certain form of superposition principle and the similarity type of general solution. Some known solutions in literature can be recovered as special cases of the analytical solution once the density change from pore ice to pore water is neglected. If the thawing front was to cease at certain time, consolidation of thawed soil in a fixed region is then encountered, and an analytical solution for the post‐thawing consolidation problem is developed using the Green's function method. Computational examples of the analytical solutions are presented. First, comparison between our model and classical model of Morgenstern and Nixon (MN model) is conducted, showing the error caused by neglecting the density change from pore ice to pore water. The MN model overestimates the excessive pore water pressure at locations near the soil surface, while underestimates it at locations near the thawing front; the error caused by neglecting the density change becomes more pronounced with increasing ice content of frozen soil. Second, comparison between the thaw consolidation process under an instant load and that under a corresponding exponential load is made. The difference in thaw consolidation behavior between the two situations mainly displays during the early stage, when there is obvious discrepancy between the two external loads; the degree of consolidation and thaw consolidation settlement are more sensitive to the discrepancy in external load than the excessive pore water pressure is.

Effects of pockmark activity on iron cycling and mineral composition in continental shelf sediments (southern Baltic Sea)

AbstractPockmarks are formed as a result of gas (methane) or/and groundwater outflow from the sea bottom. Methane, the second most important (after CO2) greenhouse gas, has a significant impact on biogeochemical processes in the bottom sediments by affecting the cycling of some elements, e.g. C, Fe, and S. Active pockmarks may also lead to changes in water column conditions by causing nutrients release from sediments. In the present study, we have focused on the impact of biogeochemical processes in pockmarks (methanogenesis, anaerobic methane oxidation, and groundwater seepage) on the transformation of iron (Fe) and the mineral composition of the sediment. In pore water, concentrations of hydrogen sulfide, phosphate, ammonia, sulfate, chloride, dissolved inorganic carbon, iron, and methane were analyzed. In the sediment, Fe speciation was performed using sequential extraction. The mineral composition was determined using powder X-Ray diffraction and scanning electron microscopy. The results from two pockmarks (with active gas seepage and groundwater infiltration) and two reference stations in the southern Baltic Sea show that geochemical conditions in pockmark sediments are significantly different from those in the typical muddy sea bottom. Pore water in pockmarks is characterized by lower sulfate and higher dissolved carbon concentrations as compared to areas of the seafloor where such structures are absent. This is due to the outflow of groundwater, which was confirmed by lower chloride concentration. In addition, sulfate is used to oxidize methane diffusing from deeper layers. Sediments in pockmarks are enriched in Fe(II) carbonates and depleted in Fe(III) (oxy)hydroxides, resulting from the anaerobic oxidation of methane with Fe(III) (Fe-AOM). Ferrous iron produced in large quantities during Fe-AOM is precipitated with carbonates.

Seawater contamination by coring and pore water sampling of marine sediments

Coring of sediments with minimal physical disturbance, and sampling of pore water with minimal contamination from seawater, are critical in marine geochemistry and microbiology. Yet, sediment coring generally causes smearing of sediment down along the inside of the core liner and leaves seawater and fluid mud trapped around the core. This seawater is a source of diffusive or advective exchange of solutes with the core, which affects the chemistry of retrieved pore water. The exchange depends on the time elapsed between coring and the end of pore water extraction. The degree of seawater contamination is difficult to determine and has generally been neglected in the geochemical literature. Yet, it is possible to quantify seawater contamination due to the high millimolar concentration of sulfate in seawater and low micromolar concentration in the subsurface methanic zone, where sulfate is largely depleted. Analysis of 540 sediment cores from a global database showed that sulfate concentration data beneath the sulfate-methane transition were often truncated and stated as “0” or “0,0” mM, thereby preventing an evaluation of seawater contamination. Among the 181 sediment cores for which the mM sulfate concentrations were reported with two or three decimal places, concentrations measured in methanic sediments were up to an order of magnitude higher than the estimated in situ concentrations in gravity cores and two orders of magnitude higher in cores recovered by scientific ocean drilling. A comparison between pore water extraction by Rhizon suction samplers and by high-pressure squeezers in cores from Baltic Sea Expedition 347 of the IODP (International Ocean Discovery Program) showed no systematic difference in the degree of contamination. A comparison with a perfluorocarbon (PFC) tracer injected by IODP coring into the drilling fluid showed a seawater contamination similar to that based on sulfate in the methanic zone. In IODP cores, where the true in situ sulfate concentration is generally unknown, we calculated the theoretical sulfate concentration that corresponded to equilibrium between SO42−, Ba2+ and barite (BaSO4) to show that the equilibrium sulfate concentrations in the methanic zone are much lower than the sulfate concentrations measured in pore water samples, again indicating seawater contamination. We used a time-resolved diffusion model to calculate how contaminating seawater ions penetrate from the external liner fluid into sediment cores before pore water is extracted. We conclude that nearly all published pore water samples retrieved from subsurface sediments are contaminated by seawater to some degree. We also conclude that sensitive sulfate analyses in methanic sediments constitute a unique opportunity to quantify this contamination and thereby provide a better possibility to control it.