A study of erosion rates on salt diapir surfaces in the Zagros Mountains, SE Iran

Salt exposures and weathering residuum on several salt diapirs in different geographic/climatic settings were studied. Anhydrite, gypsum, hematite, calcite, dolomite, quartz, and clay minerals are the main constituents of the weathering residuum covering the salt diapirs in various thicknesses. Erosion rates of residuum as well as of rock salt exposures were measured at selected sites for a period of 5 years by plastic pegs as benchmarks. Recorded data were standardized to a horizontal surface and to long-term mean precipitation. For the rock salt exposures the following long-term denudation rates were determined of 30–40 mm a⁻¹ for coastal diapirs and up to 120 mm a⁻¹ for mountain salt diapirs. Long-term mean superficial denudation rate measured on weathering residuum of low thickness reached 3.5 mm a⁻¹ on coastal diapirs. The total denudation rate estimated for the thin residuum is close to 4–7 mm a⁻¹ based on apparent correlation with the uplift rate on Hormoz and Namakdan diapirs. Denudation of rock salt exposures is much faster compared to parts of diapirs covered by weathering residuum. The extent of salt exposures is an important factor in the morphological evolution of salt diapirs as it can inhibit further expansion of the diapir. Salt exposures produce huge amounts of dissolved and clastic load, thus affecting the surrounding of the diapir.     

Three-dimensional forward and backward modelling of diapirism: numerical approach and its applicability to the evolution of salt structures in the Pricaspian basin

Numerical studies of ductile deformations induced by salt movements have, until now, been restricted to two-dimensional (2D) modelling of diapirism. This paper suggests a numerical approach to model the evolution of three-dimensional (3D) salt structures toward increasing maturity. This approach is also used here to restore the evolution of salt structures through successive earlier stages. The numerical methodology is applied to study several model examples. We analyse a model of salt diapirs that develop from an initial random perturbation of the interface between salt and its overburden and restore the evolved salt diapirs to their initial stages. We show that the average restoration errors are less than 1%. An evolutionary model of a 2D salt wall loaded by a 2D pile of sediments predicts a decomposition of the salt wall into 3D diapiric structures when the overburden of salt is supplied by 3D synkinematic wedge of sediments. We model salt extrusion feeding a gravity current over the depositional surface and estimate an average rate of extrusion and horizontal velocity of salt spreading. Faulting of the overburden to salt overhangs initiates new secondary diapirs, and we analyse the growth of these secondary diapirs. We also study how lateral flow effects the evolution of salt diapirs. The shape of a salt diapir can be very different if the rate of horizontal flow is much greater than the initial rate of diapiric growth solely due to gravity. We discuss the applicability of the results of the models to the evolution of Late Permian salt structures in the Pricaspian basin (Russia and Kazakhstan). These structures are distinguishable into a variety of styles representing different stages of growth: salt pillows, diapirs, giant salt massifs, 2D diapiric walls and 3D stocks complicated by large overhangs. The different sizes, shapes and maturities of salt structures in different parts of the Pricaspian basin reflect areal differences in salt thickness and loading history. Our results suggest that the numerical methodology can be employed to analyse the evolution of all salt structures that have upbuilt through younger ductile overburdens.     

Geochemical assessement of actinide isolation in a German salt repository environment

A significant criterion in evaluating disposal strategies for high-level nuclear waste is the assessment of the isolation capacity for the most radiotoxic radionuclides, the actinides. Important processes pertinent to potential mobilization from the waste forms, retention in secondary phases and migration of actinides in the geochemical environment of the near field of disposal locations are summarized. Criteria are formulated for assessing engineered barrier performance as a geochemical barrier for actinide long-term retention.     

Occurrence of salt water above fresh water in dynamic equilibrium in a coastal groundwater flow system near De Panne, Belgium

A salt water lens is found above fresh water under the shore between Dunkerque (France) and Nieuwpoort (Belgium). This inverse density distribution is in a dynamic equilibrium. It develops due to the infiltration of salt water on the back shore during high tide. Under this salt water lens, water infiltrated in the adjacent dune area flows towards the sea and discharges at the seabed. This water quality distribution differs from the classic salt water wedge under fresh water described in the literature. Here, the evolution to this water quality distribution is simulated with a density dependent numerical model. A large tidal range, shore morphology and a permeable groundwater reservoir are the main conditions for the observed water quality distribution. By altering these conditions, intermediate water quality distributions between the classic salt water wedge and the one discussed here develop. Based on these simulations, it is expected that similar kinds of inverse density distribution could be present in a number of coastal areas, which have tides, a gently sloping shore and a permeable substratum.     

库车坳陷中段基底断裂特征及其对盐构造变形的影响

基底断裂对盐构造的形成演化过程有重要控制作用。针对库车坳陷盐构造最为发育的克拉苏构造带和秋里塔格构造带,运用最新的二维和三维地震资料解释方案,分析了基底断裂的几何学特征、形成过程,及其对后期盐构造变形的影响。结果表明：库车坳陷不同地区基底断裂和盐构造变形样式有较大不同;基底断裂几何样式的差异,是导致盐构造复杂多变的重要原因。基底断裂的活动一定程度上影响了原始盐层的沉积厚度,还控制了盐构造的发育位置和变形样式。物理模拟实验表明：基底断裂在后期变形中相对于无断裂区可优先活动,使盐下层构造变形的传播距离变大,而盐上层变形传播距离减小。此外基底断裂的活动还可促使断裂顶部盐上层的褶皱作用增强,并发育新的断层。有限元数值模拟实验表明：在构造挤压初期,即使是远离山前的秋里塔格构造带,基底断裂顶部也可出现明显的剪应力和应变增大区域,这是该地区盐构造规模宏大的重要原因。     

The salt decay of medieval bricks at a vault in Brarup Church,Denmark

This paper presents an investigation of the decay of bricks in the chancel vault of Brarup Church located close to the Baltic Sea at the island of Falster in Denmark. The aim of the work was to study a peculiar phenomenon in order to prescribe appropriate treatment. Although protected by a lime plaster, some bricks were pulverized up to 50 mm deep from the topside. The decay occurred in a random pattern over the structure, with undamaged bricks positioned next to deteriorated bricks. The brick structure was investigated by Mercury Intrusion Porosiometry. All bricks had a majority of pores in the range 1–10 μm, but the decayed bricks also had a fraction of pores with a diameter less than 100 nm. The difference may be due to a lower firing temperature, or it may be a consequence of the salt decay. Salt analysis proved that the deteriorated bricks were contaminated with sodium chloride, which could have degraded the brick structure. However, the powdering of the bricks may be caused by the precipitation of gypsum, which was found in the brick powder but not in the deteriorated brick itself. According to calculations with the computer program RUNSALT, the precipitation of gypsum is much influenced by the presence of sodium chloride. The gypsum precipitates due to changes in temperature between 0 and 30°C and variations in relative humidity between 75 and 95%. Seasonal climatic changes in the attic were measured to cover this interval. Perhaps sodium chloride acted as a catalyst that facilitated the accumulation and precipitation of gypsum at the topside of the bricks. A sacrificial plaster, which is the traditional treatment for salt contaminated structures, would not prevent such decay.     

Influence of a low-viscosity layer between rigid inclusion and viscous matrix on inclusion rotation and matrix flow: A numerical study

We have used 2-D finite element modelling to investigate the influence of a permanent low-viscosity layer between matrix and inclusion on matrix flow and inclusion rotation under viscous simple shear flow. Rigid inclusions of different shape (circle, square, ellipse, lozenge, rectangle and skewed rectangles) and aspect ratio (R) were used. The calculated matrix flow pattern is neither bow tie nor eye-shaped. It is a new flow pattern that we call cat eyes-shaped, which is characterized by: (i) straight streamlines that slightly bend inwards at the inclusion's crests; (ii) elongate eye-shaped streamlines on each side of the inclusion; (iii) stagnation points in the centre of the eyes; (iv) absence of closed streamlines surrounding the inclusion; (v) changes in flow configuration with inclusion orientation; the lines of flow reversal bend and tilt, closed streamline circuits may disappear, and streamlines may bend outwards at the inclusion's crests.Concerning inclusion rotation, the numerical results show that: (i) a low-viscosity layer (LVL) makes inclusions with R = 1 rotate synthetically, but the rotation rate depends upon shape (circle or square) and orientation. Therefore, shape matters in the slipping mode. (ii) All studied shapes with R > 1 rotate antithetically when starting with the greatest principal axis (e₁) parallel to the shear direction ( = 0°); (iii) rotation is limited because there is a stable equilibrium orientation (_se) for all studied shapes with R > 1. (iii) There is also an unstable equilibrium orientation (_ue), and both _se and _ue depend upon inclusion's R and shape.The present numerical results closely agree with previous results of analogue experiments with a permanent low viscosity interface. Only minor deviations related with small shape differences were detected.     

琼东南盆地深水区泥-流体底辟发育特征与天然气成藏

琼东南盆地深水区发育泥-流体底辟构造,为油气垂向运移提供通道。文章基于琼东南区域构造与沉积演化认识,利用高精度三维地震资料,对深水区底辟构造进行了系统研究,识别出了底辟的分布范围,并分析了其发育机制以及对天然气成藏的控制作用。琼东南盆地深水区底辟主要发育在乐东—陵水凹陷结合部的陵南斜坡带及凸起之上,向凹陷中心规模逐渐减小。新近纪时期盆地快速沉降和生烃作用共同导致了深部异常高压的存在,其为底辟发育提供了原始的驱动力;陵南斜坡带古近纪发育的断裂为区域内的构造薄弱带,且在新近系发生活化为异常高压释放提供突破口,进而导致了底辟的形成,即底辟形成时间为新构造运动时期。底辟构造可沟通崖城组烃源岩及黄流组、莺歌海组储集体,是深水西区油气垂向运移的主要通道,控制了天然气藏的分布。根据底辟构造与烃源、储层的配置关系,提出了陵南斜坡带和松南低凸起区是深水区下步勘探的重点区带。     

Installation of a vertical slurry wall around an Italian quarry lake: complications arising and simulation of the effects on groundwater flow

Slurry walls are non-structural barriers that are constructed underground to impede groundwater flow or manage groundwater control problems. The study area is in the Piemonte plain (Italy), close to the River Po. Quarrying works carried out below the piezometric surface created two big quarry lakes. The local groundwater system is characterized by a lower semi-confined aquifer, which is overlain by a semi-permeable bed of clayey peat (aquitard) and an upper unconfined aquifer. Locally, the peat fades away and the granulometry of this horizon becomes silty sandy. A planned enlargement of the quarry will increase the size and depth of the quarry lakes. So the aquitard bed between the two aquifers will be damaged, creating a mixing rate of groundwater. Such a procedure would not be compatible with the presence of two municipal wells upstream from the quarries. Consequently, the installation of a vertical diaphragm (slurry wall) is recommended to separate the aquifers and to act as a filter for the groundwater flowing from the unconfined to the semi-confined aquifer. To predict the consequences caused by the installation of the vertical diaphragm separating the unconfined aquifer and the semi-confined one, a specifically adjusted finite-difference model was used. The model showed a maximum rising of the water table equal to 12 cm, just upstream of the diaphragm and for a distance of about 100 m, and a maximum lowering of 2 cm just downstream of the diaphragm. However, the slurry wall would not cause any change in the piezometric head in the area where there are municipal wells and, hence, will not have any negative effect on the functionality of the municipal wells. Moreover, the migration of water from the unconfined aquifer through the vertical diaphragm will stimulate a series of attenuation and auto-depuration processes of eventual contaminants. These processes are due to the higher crossing time that the groundwater flow takes to go through the vertical barrier (t _a = 96.5 days, whereas for the horizontal semi-permeable layer t _a = 9.6 days). So, the vertical diaphragm can be a resolutive element, representing a mediation and separation factor between the unconfined and the semi-confined aquifers along the border of the quarrying areas, and a protective barrier for the water quality of the quarry lake and the semi-confined aquifer.     

基于供水风险分析的汛限水位控制范围研究

以大连市碧流河水库为例,根据大连市供用水系统的风险及其变化过程,分析超蓄洪水资源的利用情况及其效益,并在对抬高水库汛限水位造成的淹没损失分析的基础上,利用边际成本分析方法确定控制的范围。计算结果表明,随着大连市用水量的增加,碧流河水库汛限水位可适当抬高。     

The role of saline solution properties on porous limestone salt weathering by magnesium and sodium sulfates

Saline solution properties, viscosity in particular, are shown to be critical in salt weathering associated with sodium and magnesium sulfate crystallization in porous limestone. The crystallization of sodium and magnesium sulfate within a porous limestone has been studied at a macro- and microscale using different techniques, including mercury intrusion porosimetry, environmental scanning microscopy and X-ray computed tomography. Such analysis enabled the visualization of the crystallization process in situ, and at high magnification, yielding critical information as to where and how salts crystallize. Sodium sulfate decahydrate (mirabilite) tends to crystallize in large pores as euhedral micron-sized crystals formed at low supersaturation near to the surface of the stone. In contrast, magnesium sulfate heptahydrate (epsomite) tends to precipitate as anhedral wax-like aggregates formed at high supersaturation and distributed homogeneously throughout the stone pore system filling large and small pores. While the former crystallization behavior resulted in scale formation, the latter led to crack development throughout the bulk stone. Ultimately, the contrasting weathering behavior of the two sulfates is explained by considering differences in flow dynamics of solutions within porous materials that are mainly connected with the higher viscosity of magnesium sulfate saturated solution (7.27 cP) when compared with sodium sulfate saturated solution (1.83 cP). On the basis of such results, new ways to tackle salt weathering, particularly in the field of cultural heritage conservation, are discussed.     

The influence of petrophysical properties on the salt weathering of porous building rocks

The influence of pore structure, water transport properties and rock strength on salt weathering is evaluated by means of a thorough rock characterisation and a statistical analysis. The pore structure was described in terms of its porosity, pore size distribution (quantified by mean pore radius) and specific surface area, density and water transport was characterised by means of water permeability (saturated flow) and capillary imbibition (unsaturated flow); whilst the rock strength test was carried out using uniaxial compressive strength, compressional and shear wave velocities, dynamic elastic constants and waveform energy and attenuation were obtained from the digital analysis of the transmitted signal. A principal component analysis and a stepwise multiple regression model was carried out in order to examine the direct relationships between salt weathering and petrophysical properties. From the principal component analysis, two main components were obtained and assigned a petrophysical meaning. The first component is mostly linked to mechanical properties, porosity and density whereas the second component is associated with the water transport and pore structure. Salt weathering, quantified by the percentage of weight loss after salt crystallisation, was included in both principal components, showing its dependence on their petrophysical properties. The stepwise multiple regression analysis found that rock strength has a predominant statistical weight in the prediction of salt weathering, with a minor contribution of water transport and pore structure parameters.     

Numerical modeling of the development of a preferentially leached layer on feldspar surfaces

Compositional depth profiles in the leached layer of feldspar surfaces are usually interpreted by using analytical solutions which introduce oversimplifying assumptions. Here we present a general multicomponent interdiffusion numerical model for simulating cation release from a preferentially leached layer on feldspar surfaces in acid solutions. The numerical model takes into account interdiffusion, dissolution of the solid phase (represented by a moving boundary problem), and adsorption in the leached layer. Effective diffusion coefficients of ions vary with concentration along the leached layer. Governing equations of ions diffusion in the leached layer are solved numerically with a finite element method implemented in a multicomponent reactive transport code, CORE^3D, previously verified against analytical solutions of compositional depth profiles. The numerical model is tested with published X-ray photoelectron spectroscopy (XPS) data on early development of compositional profiles of labradorite leached in pH 2 HCl solutions. Model parameters are estimated by fitting depth profiles of Ca and Al measured at 12, 26, 48, 72, and 143 h. The best fit is achieved with tracer diffusion coefficients of 4 × 10⁻¹⁸, 8 × 10⁻¹⁷, 3.4 × 10⁻¹⁷, and 7 × 10⁻¹⁸ cm²/s for H, Na, Ca, and Al, respectively, which fall within the range of values reported in the literature. Our estimate of the retreat velocity corresponding to the dissolution rate is 3 × 10⁻¹³ cm/s. Results of sensitivity runs show that computed compositional profiles are sensitive to most model parameters.     

Numerical prediction of the creep behaviour of an unstabilised and a chemically stabilised soft soil

This paper examines the ability of volumetric and deviatoric creep laws associated with constitutive models to simulate the creep behaviour of a soft soil in its natural state or chemically stabilised state. Initially, the models/laws are validated by oedometer and triaxial creep tests, for the stabilised and unstabilised soils. Finally, the long-term behaviour of an embankment built on a soft soil reinforced with deep mixing columns is predicted based on the properties for a curing time of 28 days. The results show that the creep phenomenon should be considered in a long-term analysis of deep mixing columns.     

The Influence of Floodplain Compartmentalization on Flood Risk within the Rhine–Meuse Delta

The present compartmentalization layout within the river polders in the Dutch Rhine–Meuse delta is the result of abandonment and partial removal of secondary dikes and the construction of modern infrastructure embankments. These structures will guide the flow of water in case the polder would inundate. Through the application of a 2-D flood propagation model in the polder Land van Maas en Waal this study explores whether restoration or removal of old dike remnants would contribute to a reduction of the risk and damage during an inundation. A systematic set of 28 flood scenarios was simulated and for each scenario an additional damage and risk assessment was carried out. It is concluded that a simple removal or total restoration will not reduce flood damage, but that this must be achieved by a strategic compartment plan. With such a plan old dike remnants and present embankments can be used to keep water away from vulnerable and valuable areas for as long as possible and to guide the floodwater to areas that are considered less vulnerable.     

Control of rheological stratification on rifting geometry: a symmetric model resolving the upper plate paradox

Numerical experiments reproduce the fundamental architecture of magma-poor rifted margins such as the Iberian or Alpine margins if the lithosphere has a weak mid-crustal channel on top of strong lower crust and a horizontal thermal weakness in the rift center. During model extension, the upper crust undergoes distributed collapse into the rift center where the thermally weakened portion of the model tears. Among the features reproduced by the modeling, we observe: (1) an array of tilted upper-crustal blocks resting directly on exhumed mantle at the distal margin, (2) consistently oceanward-dipping normal faults, (3) a mid-crustal high strain zone at the base of the crustal blocks (S-reflector), (4) new ocean floor up against a low angle normal fault at the tip of the continent, (5) shear zones consistent with continentward-dipping reflectors in the mantle lithosphere, (6) the mismatch frequently observed between stretching values inferred from surface extension and bulk crustal thinning at distal margins (upper plate paradox). Rifting in the experiment is symmetric at a lithospheric scale and the above features develop on both sides of the rift center. We discuss three controversial points in more detail: (1) weak versus strong lower crust, (2) the deformation pattern in the mantle, and (3) the significance of detachment faults during continental breakup. We argue that the transition from wide rifting towards narrow rifting with a pronounced polarity towards the rift center is associated with the advective growth of a thermal perturbation in the mantle lithosphere.     

Modelling the settlement behaviour of a strip footing on sloping sandy fill under cyclic loading conditions

This paper investigates the settlement behaviour of a strip footing seating on the crest of an embankment and subjected to cyclic loading. The embankment fill is a dense sand and the issue is the gradual accumulation of settlement over a large number of load repetitions. Cyclic triaxial tests were first conducted to develop a consistent but simple material model for numerical implementation. Particular emphasis was placed on linking the stress-strain behaviour of an unload-reload cycle to the accumulation of permanent strain, with only five input parameters required to model the cyclic behaviour. The material model was implemented in a numerical analysis to compute the settlement behaviour obtained from model tests conducted by another researcher. It is pertinent to highlight that the same soil, compacted to same density at same moisture content, was used for both the cyclic triaxial tests and model tests. Reasonable to good agreement between the experimental and numerical results was achieved.     

The behaviour of an impermeable flexible raft on a deep layer of consolidating soil

Analytic solutions to the problem of the time-settlement behaviour of raft foundations have been limited in the past to flexible or rigid loadings, and have treated the foundation as being completely permeable. In this paper, solutions are presented for smooth circular rafts of any flexibility causing consolidation of a deep homogeneous clay layer, where the raft may be considered permeable or impermeable. Results for the time-dependent behaviour of contact stresses, pore pressures, raft displacements and moment in the raft are presented.     

The mean and turbulence structure simulation of the monsoon trough boundary layer using a one-dimensional model withe-l ande-ε closures

An attempt has been made here to study the sensitivity of the mean and the turbulence structure of the monsoon trough boundary layer to the choice of the constants in the dissipation equation for two stations Delhi and Calcutta, using one-dimensional atmospheric boundary layer model withe-ε turbulence closure. An analytical discussion of the problems associated with the constants of the dissipation equation is presented. It is shown here that the choice of the constants in the dissipation equation is quite crucial and the turbulence structure is very sensitive to these constants. The modification of the dissipation equation adopted by earlier studies, that is, approximating the Tke generation (due to shear and buoyancy production) in theε-equation by max (shear production, shear + buoyancy production), can be avoided by a suitable choice of the constants suggested here. The observed turbulence structure is better simulated with these constants. The turbulence structure simulation with the constants recommended by Aupoixet al (1989) (which are interactive in time) for the monsoon region is shown to be qualitatively similar to the simulation obtained with the constants suggested here, thus implying that no universal constants exist to regulate dissipation rate. Simulations of the mean structure show little sensitivity to the type of the closure parameterization betweene-l ande-ε closures. However the turbulence structure simulation withe-ε. closure is far better compared to thee-l model simulations. The model simulations of temperature profiles compare quite well with the observations whenever the boundary layer is well mixed (neutral) or unstable. However the models are not able to simulate the nocturnal boundary layer (stable) temperature profiles. Moisture profiles are simulated reasonably better. With one-dimensional models, capturing observed wind variations is not up to the mark.