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.     

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.     

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.     

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.     

The effect of surface moisture on the entrainment of dune sand by wind: an evaluation of selected models

Abstract Reliable predictions of wind erosion depend on the accuracy of determining whether erosion occurs or not. Among the several factors that govern the initiation of soil movement by wind, surface moisture is one of the most significant. Some widely used models that predict the threshold shear velocity for particle detachment of wet soils by wind were critically reviewed and evaluated. Wind‐tunnel experiments were conducted on pre‐wetted dune sand with moisture contents ranging from 0·00 to 0·04 kg kg^?1. Sand samples were exposed to different wind speeds for 2 min. Moisture content was determined gravimetrically before and after each experiment, and the saltation of sand particles was recorded electronically with a saltiphone. Shear velocities were deduced from the wind speed profiles. For each moisture content, the experiments were repeated at different shear velocities, with the threshold shear velocity being determined by least‐squares analysis of the relationships between particle number rates and shear velocity. Within the 2‐min test runs, temporal changes in particle number rates and moisture contents were detected. A steep increase in the threshold shear velocity with moisture content was observed. When comparing the models, large differences between the predicted results became apparent. At a moisture content of 0·007 kg kg^?1, which is half the moisture content retained to the soil matrix at a water tension (or matric potential) of ?1·5 MPa, the increase in ‘wet’ threshold shear velocity predicted with the different models relative to the dry threshold shear velocity ranged from 117% to 171%. The highest care should therefore be taken when using current models to predict the threshold shear velocity of wet sediment. Nevertheless, the models of Chepil (1956; Proc. Soil Sci. Soc. Am., 20, 288–292) and Saleh & Fryrear (1995; Soil Sci., 160, 304–309) are the best alternatives available.     

The effect of hydrodynamic conditions on true flotation and entrainment in flotation of a complex sulphide ore

The separation efficiency and selectivity of flotation are directly proportional to recoveries of the mineral species in the feed due to true flotation and entrainment. In this study, effects of the hydrodynamic conditions on true flotation and entrainment were investigated by using a fractional factorial experimental design. A method previously described in the literature was applied to determine the contributions of true flotation and entrainment in flotation of a complex sulphide ore. In order to apply the method, the kinetic flotation tests were conducted under various hydrodynamic conditions defined by some physical variables. Some of these tests were conducted in the presence and absence of a collector to evaluate the self-induced floatability. The selectivity index of the mineral species for entrainment was seen to be suitable evaluation of the non-selectivity and efficiency of the entrainment. Furthermore, the results of the size-by-size analysis of the froth products indicated that the presence of the self-induced hydrophobic particles in the feed is as important as the presence of very fine particles for accurate estimation of true flotation and entrainment in flotation of a complex sulphide ore. In addition, the estimated results for entrainment in flotation of the complex sulphide ore can be misleading. Therefore, a new approach would be necessary to determine the contributions of true flotation and entrainment in flotation of a complex sulphide ore.     

A three‐dimensional numerical model of sediment transport,erosion and deposition within a network of channel belts,floodplain and hill slope: extrinsic and intrinsic controls on floodplain dynamics and alluvial architecture

A three‐dimensional numerical model of sediment transport, erosion and deposition within a network of channel belts and associated floodplain is described. Sediment and water supply are defined at the upstream entry point, and base level is defined at the downstream edge of the model. Sediment and water are transported through a network of channels according to the diffusion equation, and each channel has a channel belt with a width that increases in time. The network of channels evolves as a result of channel bifurcation and abandonment (avulsion). The timing and location of channel bifurcation is controlled stochastically as a function of the cross‐valley slope of the floodplain adjacent to the channel belt relative to the down‐valley slope, and of annual flood discharge. A bifurcation develops into an avulsion when the discharge of one of the distributaries falls below a threshold value. The floodplain aggradation rate decreases with distance from the nearest active channel belt. Channel‐belt degradation results in floodplain incision. Extrinsic (extrabasinal, allogenic) and intrinsic (intrabasinal, autogenic) controls on floodplain dynamics and alluvial architecture were modelled, and sequence stratigraphy models were assessed. Input parameters were chosen based on data from the Rhine–Meuse delta. To examine how the model responds to extrinsic controls, the model was run under conditions of changing base level and increasing sediment supply. Rises and falls in base level and increases in sediment supply occurred over 10 000 years. Rising base level caused a wave of aggradation to move up‐valley, until aggradation occurred over the entire valley. Frequency of bifurcations and avulsions increased with rate of base‐level rise and aggradation rate. Channel‐belt width varied with water discharge and the lifespan of the channel belt. Wide, connected channel belts (and high channel‐deposit proportion) occurred around the upstream inflow point because of their high discharge and longevity. Less connected, smaller channel belts occurred further down‐valley. Such alluvial behaviour and architecture is also found in the Rhine–Meuse delta. During base‐level fall, valley erosion occurred, and the incised valley contained a single wide channel belt. During subsequent base‐level rise, a wave of aggradation moved up‐valley, filling the incised valley. Bifurcation and avulsion sites progressively moved upstream. Relatively thin, narrow channel belts bordered and cut into the valley fill. These results differ substantially from existing sequence stratigraphy models. The increase in sediment supply from upstream resulted in an alluvial fan. Most bifurcations and avulsions occurred at the fan apex (nodal avulsion), and channel belts were the widest and the thickest here (giving high channel‐deposit proportion) due to their high discharge and longevity. The width and thickness of channel belts decreased down‐valley due to decreased discharge, longevity and aggradation rate. This behaviour occurs in modern alluvial fans. Intrinsic controls also affect floodplain dynamics and alluvial architecture. Variation of aggradation rate, bifurcation frequency and number of coexisting channel belts occurred over periods of 500 to 2000 years, compared with 10 000 years for extrinsic controls. This variation is partly related to local aggradation and degradation of channel belts around bifurcation points. Channel belts were preferentially clustered near floodplain margins, because of low floodplain aggradation rate and topography there.     

Adaptive hierarchical upscaling of flow in heterogeneous reservoirs based on an a posteriori error estimate

This paper treats the upscaling of the absolute permeability in a heterogeneous reservoir. By replacing the fine scale permeability tensor with an upscaled, or effective permeability tensor, a modelling error is introduced. An a posteriori error estimate on this modelling error is formulated and tested. An implementation of the theory, based on domain decomposition coupled with a hierarchical representation of the absolute permeability field, is given. As hierarchical basis functions we have chosen the Haar system, which leads to a wavelet representation of the permeability. The wavelet representation offers a natural upscaling technique which resembles the highcut filters commonly used in signal analysis. This procedure represents an adaptive upscaling method. The numerical results show that this method conserves both the dissipation and the mean velocity in the problem fairly well. The a posteriori error estimate on the modelling error coupled with domain decomposition methods constitutes a powerful modelling tool.