ACTA CARSOLOGICA 34/2 2 317-348 LJUBLJANA 2005 COBISS: 1.01 CONDENSATION CORROSION: A THEORETICAL APPROACH KONDENZACIJSKA KOROZIJA: TEORETIČNI PRISTOP WOLFGANG DREYBRODT12, FRANCI GABROVŠEK2 & MATIJA PERNE3 1 Karst Processes Research Group, University of Bremen, Germany 2 Karst Research Institute, ZRC SAZU, Titov trg 2, SI-6230 Postojna, Slovenia 3 Student, Faculty of Mathematics and Physics, University of Ljubljana, Slovenia Prejeto / Received: 20. 9. 2005 Abstract UDC: 551.435.8:551.3.053 Wolfgang Dreybrodt, Franci Gabrovšek & Matija Perne: Condensation corrosion: A theoretical approach Condensation of water from warm, humid air to cold rock walls in caves is regarded to play a significant role in speleogenesis. The water condensing to the cave walls quickly attains equilibrium with the carbon dioxide in the surrounding air, and consequently dissolves limestone or gypsum forming various types of macro- ,meso-, and micromorphologies. In this paper we present the basic physical principles of condensation and give equations, which allow a satisfactory estimation of condensation rates. Water condensing to a cooler wall releases heat of condensation, which raises the temperature of the wall thus reducing the temperature difference AT between the warm air and the cave wall. Furthermore one has to take into account the heat flux from the air to the cave wall. This defines the boundary conditions for the equation of heat conduction. For a constant temperature of the air initial condensation rates are high but then drop down rapidly by orders of magnitude during the first few days. Finally constant condensation rates are attained, when the heat flux into the rock is fully transmitted to the surface of the karst plateau. For spherical and cylindrical conduits these can be obtained as a function of the depth Z below the surface. When diurnal or seasonal variations of the air temperature are active as is the case close to cave entrances, condensation rates can become quite significant, up to about 10-6 m/year. The theoretical results are applied also to corrosion of speleothems and the formation of »röhrenkarren« as described by Simms (2003). To convert condensation rates into retreat of bedrock the saturation state of the solution must be known. In the appendix we present experiments, which prove that in any case the solution flowing off the rock is saturated with respect to limestone or gypsum, respectively.. Key words: condensation, corrosion, speleogenesis, heat transfer. Izvleček UDK: 551.435.8:551.3.053 Wolfgang Dreybrodt, Franci Gabrovšek & Matija Perne: Kondenzacijska korozija: teoretični pristop O speleogenetskem pomenu vode, ki iz toplega vlažnega zraka kondenzira na hladne jamske stene je napisanega veliko. Kondenzirana voda se hitro uravnoteži z ogljikovem dioksidom v jamski atmosferi, zato raztaplja apnenec in pri tem tvori različne skalne oblike. V članku predstavimo fizikalne osnove in podamo enačbe, ki omogočajo približno oceno hitrosti kondenzacije v različnih pogojih. Zaradi kondenzacijske toplote in prenosa toplote iz zraka na steno, se temperatura stene viša, pri čemer se zmanjšuje razlika temperature med zrakom in steno. To predstavlja robni pogoj za prevajanje toplote iz jame. Pri konstantni temperaturi zraka, hitrost kondenzacije v nekaj dneh pade za več velikostnih redov, dokler ne doseže končne vrednosti, pri kateri se ves toplotni tok prenese na površje. Slednjo za primer krogelnih in valjastih prostorov opišemo kot funkcijo globine rova Z. Pri dnevnih in sezonskih spemembah temperature, značilnih v bližini jamskih vhodov, je hitrost kondenzacije in posledično korozije do 1^m na leto. Teoretične rezultate uporabimo tudi za izračun korozije kapnikov in tipa škrapelj (röhrenkarren), ki jih opisuje Simms (2003). Za pretvorbo hitrosti kondenzacije v hitrost korozije, potrebujemo podatek o stopnji nasičenja vode. V dodatku predstavimo poskus, ki dokazuje, da se voda, ki kondenzira na jamske stene, hitro nasiti tako v primeru apnenca kot v primeru sadre. Ključne besede: kondenzacija, korozija, speleogeneza, prevajanje toplote. INTRODUCTION Water vapor from a cave atmosphere condensing to the walls of a cave creates a water film in equilibrium with the partial pressure pCO2 of the cave atmosphere. This solution is therefore aggressive to limestone and the dissolution process based on it has been termed condensation corrosion. The most recent and comprehensive review on condensation in karst and its role on hydrology and speleogenesis has been published by Dublyansky and Dublyansky (2000) and Klimchouk et al. (1996). They report that in summer condensation supplies a significant amount of water (up to 20% of the total dry reason run-off) to karst systems. They also discuss the role, condensation corrosion plays in sculpturing cupolas in limestone and in gypsum caves, when evaporation from open water surfaces at elevated temperature, above that of the cave wall, produces warm cave air saturated with water vapor. Condensation of water to the cave walls supplies an aggressive solution, which runs down the cave wall and is replaced by fresh, condensing water. Cigna and Forti (1986) and Calafora at al. (1993) have reported on field measurements of condensation. Large cupolas in limestone and gypsum caves are explained by this mechanism which schematically is shown in Fig. 1 (Audra et al. 2002) . Fig. 1: Water evaporates from an open surface of elevated temperature Ta to the cave air in thermal equilibrium with the water. Vapor pressure at temperature Ta is given by p. The partial pressure Pco^ in the air is in equilibrium with Pco^ in the water. The vapor condenses at the cave wall with temperature T^ < Ta and flows back as a thin film, designated by arrows. These water films become saturated with respect to the mineral composing the rock (limestone or gypsum). The free surface of water at elevated temperature Ta could be a lake of hydrothermal water. Water condenses at the cave walls if their temperature Tf is below the dew point of the air. The condensed water flows back in a thin film covering the cave walls. As long as the temperatures Ta and Tf are constant in time this is a continuous, everlasting process. To estimate the dissolution rate one has to know the wall temperature Tf. If initially the cave wall and the surrounding rock is at temperature Tfi, water from the air with temperature Ta>Tfi starts to condense to the wall. This causes a continuous flow of heat into the rock. A small part of it is due to heat transfer from the warmer air into the cold rock. The major part stems from release of heat of condensation from the condensing water. This heat transfer increases the temperature of the cave wall, causing reduction of the rate of condensing water. Eventually a stationary state is reached and the temperature of the cave wall, Tfstat becomes independent of time. It is also possible that the temperature of the open surface changes diurnally or seasonally, e.g. a river flowing into the cave at elevated temperature during day time (summer) and at low temperature during night (winter). In this paper, we address the question, what are the rates of condensation and what average annual retreat of bedrock follows as consequence. The first part of this paper presents a theory of condensation rates for the stationary state of condensation for time independent temperature of the cave air, fully saturated with water vapor. Then we will turn to situations where the air temperature in the cave changes diurnally or seasonally. These results will be also applied to corrosion of speleothems, which is an important topic in conservation of tourist caves (Avramidis et al. 2001). A recently observed form of small scale cupolas of a few cm diameter and a length of about 10 cm, growing upwards from bedding planes of limestone at a lake shore (Simms 2003) is also explained by our theoretical findings. To translate rates of condensation into bedrock retreat in m/year one has to know the saturation state of the water, when it flows off the rock. In the appendix we present laboratory experiments on limestone and gypsum, which prove that those waters are saturated with respect to limestone and gypsum, respectively, in all relevant situations. This work shall provide a theoretical basis to the interpretation of field data. Therefore we assume two scenarios of boundary conditions: a) stable, time independent temperature of the cave air or b) periodic variations in temperature, diurnally or seasonally. We do not ask the question, how these conditions arise in detail. This question deals with cave climate and is a complex subject (Badino 1995; Wigley & Brown 1976). Our question to be answered is: How effective under given boundary conditions is water condensation to the cave walls or to speleothems, and what are the rates of condensation corrosion? BASIC THEORY The amount of water condensing to an exposed rock surface per unit of time and surface area determines the retreat of rock by dissolution. Therefore, it is of utmost importance to give a reliable estimation. We assume that the exposed surface of the rock is covered by a thin film of condensed water, about 10-2 cm deep. Such a thin film will quickly come to thermal equilibrium with the temperature Tf of the cave wall. It keeps its constant depth since there is flow from the rock surface down to the cave floor, which compensates for condensing water. Water condenses to the film, if the partial pressure Pa of vapor in the cave air exceeds the vapor pressure Pf at the temperature Tf of the water film. In the following, we assume that the cave air is saturated, i.e. relative humidity is 100%. For condensation, water molecules must be transported to the cave wall, where they attach to the water film. We assume that the cave air is well mixed, but close to the rock a diffusion boundary of thickness sD [m] exists. Transport of water molecules is effected through this layer by molecular diffusion. The flux of vapor Fv to the surface of the film is given by Fick's law as [gmV] 1 where D is the constant of molecular diffusion of water molecules in air (D = 2.5 ■ 10-5 m2/s), c m ^ m ' a and cf are the concentrations in g/m3 of water vapor at temperatures Ta and Tf respectively. By use of the equation for ideal gases this can be written as [gmV] 2 Pa and Pf are the vapor pressures measured in Pa, Ta is the temperature of the cave air in K. In all following calculations we use Ta=300 K. M = 18 g/mol is the molecular weight of water and R is the universal gas constant R = 8.314 J mol-1K-1. Pa-Pf can be approximated by AP [Pa] 3 where AP/ AT = 100 Pa/K for the temperatures of interest. The vapor condensing to the water film releases heat of condensation qc = 2450 J/g. This is a very high amount of energy. A condensation rate of 1 mm/day = 1 kg/(m2day) causes a heat flux of 28 W/m2, equivalent to the heat produced by lightening the cave by electrical bulbs of 28 W, with one bulb on each square meter of the cave walls. The flux of heat Fq released by condensation is given by D„ M AP 4.42 10' _ An additional flux from the warm cave air to the cooler water film is given by Fc (T^.T^) [Wm-'] ka is the thermal conductivity of air ( 2.6^ 10-2 Wm-1K-1) and sT is the thickness of the thermal boundary layer, which is related to sD by sD = sT(Sc/Pr)"3. Sc is the Schmidt number for diffusion and Pr the Prandtl number for heat convection(Beek & Muttzall 1975). For air, Sc=Pr=1. Therefore sD= sT=s. The total heat flux to the water film is the sum of Eqs. 4 and 5 Pu,, =-{T,-Tf), k=0.0702 Wm'K' 6 S This heat flux causes an increase of the temperature of the water film, until the heat flux transported into the rock equals the heat flux from the cave air into the water film. If one knows the temperature Tf of the water film, equal to that of the rock surface the amount of condensing water can be obtained by inserting Eq. 3 into Eq. 2, as F = AT L-Tf RZ S In a first approach, we reduce the problem to a simple one-dimensional setting. Fig. 2 shows this. An extended large cave is located at depth Z below the surface. We regard only the vertical heat flux towards the surface, thus reducing the problem to one dimension in z-direction. This gives the following situation. At z = 0 the temperature T0 is kept constant for t > 0. The cave roof at depth Z experiences a heat flow Ftot given by Eq. 6. A solution to this problem is given in Carslaw & Jaeger (1959, p.125). We will not give the complete solution here. We extract what we need for further discussion: The solution has a stationary part and the surface temperature of the rock at the cave wall becomes To ♦ ♦ ♦ ♦♦ ♦ 1 F,„.=|(VT,(t)) T=T(t=0)=T„ Z ♦ ♦ # ♦♦ ♦ 1 t T. ka — , Eq. 10 can be approximated by K 1- 1 ka 1 k 4t 11 The rates of condensing water for the the fast initial rise of temperature can now be obtained from Eqns. 7, 10 and 11 Py- AP V / M RT„ exp V V aV vKt erfc ka 4Kt 12 For times t> /7 \ fM^ k V — , a slow decline of the temperature difference Ta-Tf(t) towards stationary state yields rates M r 1 n AT > RT, ■sItoc k -^t ^ [Ta-T,] gm'^s ' 13 It is important to note that this is independent of s. The minimum rate of condensation is reached at stationary state after time 5ts. 1.0 0.8 p= 0.6 fei ^ 0.4 e 0.2 0.0 a) t,=1.54 h - - Z=10 m - - Z=30m -Z=irTfinity I CO 8=10^m 2 4 6 8 10 12 14 16 18 20 22 24 Time [hours] 1.000 0.998 3 0.996 0.994 0.992 0.990 1 ' I'll 1 ' 7 -►T^=15.3y A/(1+A)=0.9991 |-»-T„=1.7y A/(1+A)=0.9974 — Z=10m — Z=^m " — Z= infinite £=10^ m I.I.I. 1 10 20 30 Time [years] 40 50 Fig. 3: Time dependence of the temperature difference T-T(t) normalized to initial temperature difference T-Tfi. a) Results for the first 24 hours. b) Results for long times. Until time Tf normalised temperature difference rise to 0.84. After this time it rises according to Eq.11 until finally it approaches a stationary state after time 5ts. The black curve is calculated with s=1&3m and Z=10 m and the blue one for Z=30 m. The red curve is calculated for Z=> krs/k the condensation rate becomes independent of s. This is true for Z >10 m and s<10-3 m. ' The decay time towards this stationary state is given by T. = ' ß'K [S] 15 At a depth of 100 m this is about 127 years, a short time with respect to the time scales of cave evolution. Fig. 4a shows the rates of condensation for Z = 100 m and various values of s. Fig. 4b depicts the amount of condensed water after time t, which is obtained by integration of the curves in Fig. 4a. In summary, the following picture emerges. For times smaller than 0.5ts the temperature increase of the cave wall can be well approximated by Eqns.10 and 11. Later, for times t > 5ts the stationary state is reached by an exponential approach and the temperature becomes constant. Fig. 3 shows this for s = 0.001 m. The black line shows (Tf - Tfi)/(Ta - Tfi) for a depth of Z = 10 m, the blue line gives the temperature rise for Z = 30 m. The red curve shows the temperature for the semi-infinite plane with Z=ro. CONDENSATION IN STATIONARY BOUNDARY CONDITIONS Our one-dimensional model is idealistic because it requires a cave chamber with horizontal dimension, one order of magnitude larger than its height. Other geometries are more suitable, e.g. a spherical room with diameter Ds or a cylindrical conduit of diameter Dc and length L, both buried at depth Z. These are shown in Fig. 5a-c. For such situations analytical results are not available. But the general behavior is similar to the idealistic one-dimensional case, and it is possible to obtain the temperature at stationary state, reached after t > 5 Ts, by the theory of conduction shape factors (Incropera & DeWitt 2002). At stationary state, the wall of the rock is at temperature Tfs'a' and the surface temperature at z = 0 is To. See Fig. 2. The total amount of heat flowing from the conduit to the surface is then given by [W] 16 where S [m] is the shape factor. The expressions for some shape factors are given in Table 1. System Restriction Shape factor Sphere in a semi infinite z > D/2 2nD medium (Fig. 5a) \-DIAZ Cylinder in a semi infinite medium (Fig. 5b) L >> D z > 3D/2 2kL ln(4Z//)) Slab (Fig. 5c) o/Z Table 1: Shape factors for some typical geometries as shown in Figs. 5a-c. 0.05 i 0.04 ro.o3 =S0.02 10.01 0.00 a) Z=100m stationary 50 100 150 200 Time [years] 250 10 - 0.035 8 J _ 0.030 ® 1 0.025 4 1 2 0.020 0 - 0.015 - 0.010 - 0.005 0.000 300 100 150 200 Time [years] Fig. 4: a) Rates of condensation and retreat of bedrock as function of time divided by the initial temperature difference Ta-Tfi. Values of s are written at the curves. Note that for large times all curves coincide, i.e the rates become independent of s. b) Amount of condensed water per m2 and K of temperature difference and total retreat of bedrock for the cases of Fig. 4a. Note that although the variation of s covers two orders of magnitude, the amount of condensed water only weakly depends on s. The total heat transferred to the wall of the cave is 17 where a is the surface area of the cave. Conservation of energy requires Qw = Qr. Therefore one finds 1+ J_ ka k-sS 18 From this the rate of condensed water can be calculated by use of Eq. 7. Fig. 5d shows these rates per degree K for circular cave rooms and cylindrical conduits, buried at depth Z. Note that for s << 0.038a/S the rates are independent of s, which is true for all practical applications. a) b) z p * ■ fTi")/ ^-MTd T, CT C) 5=: 1 > 0.1 0.01 1E-3 50 100 Z, Depth [m] .d) ' ' ' ' ^__ '8=0.001 m - " \ ^ ^ ~ -Sphere -_____ ■-D(0)=1m,- - D=5m,- ■ ............ Condurt ■-D(0)=1m,- - D=5m,- ■ ■ D=10m Plane , i , i 1 , " 1 „ 0.1 'S 0.011 I 1E-3 150 200 Fig. 5: a-c) Geometrical configurations for the shape factors in Table 1. d) Condensation rates and retreat of bedrock in the stationary state in dependence on depth Z for spheres and cylindrical conduits with various diameters. D(0) is the initial diameter. The blue curve represents one dimensional scenario (slab with thickness Z) of Fig. 5c. See also Fig. 2. Spheres show the highest condensation rates, almost independent of their depth Z below ground, but dependent on their diameter DS. As we will show in the appendix, films of water running down rock-walls of gypsum or limestone can be regarded as saturated with respect to gypsum or limestone. Therefore at a pCO2 = 0.00035 atm, for limestone with saturation concentration of 60 mg/L, 1 g/m-2s-1 of condensed water corresponds to retreat of rock by 7.510-4 m/year using p = 2500 m^ as density of a slightly porous limestone. For gypsum with a saturation concentration of 2.5 g/L (independent on pCO2) and density 2300 kg/m3, the retreat of rock can be approximately obtained by multiplying the values for limestone (pCO2 = 0.00035 atm) by a factor of 50. Now we assume the following situation. A spherical cave of initial diameter Ds(0) is invaded by geothermal waters with temperature Ta, higher than the initial rock temperature Tfi. Water evaporates from a lake, stable for long times. We assume the cave 100 m below ground. Then after about 100 years the stationary state is attained. The change of diameter Ds is then approximated quite accurately by ^«2xl.8.10-=x7.5-10^x-^[t^-T^) [m/year] 19 S H * K 2 For limestone one gets dD. dt D. Note that growth rates decrease with 1/DS. MO^—-(TL-T^) [m/year] 20 Integration of Eq. 20 yields t-iT^-Tj:) 21 where t is in years. A cave with initial diameter of 1 m evolves into a cupola of 10 m diameter independent of its depth Z for D/(4Z)>10 in 5 ■107 years for Ta - Tfi = 1 K. A sphere with initial diameter of 0.1 m needs, however, only 530000 years to reach a diameter of 1 m for T -T= 1 K. fi For a cylindrical conduit buried at depth Z, [m/year] 22 Dc which can be integrated -Z)^(0|ln^-i| + Z)(0)|ln®-i| = 2.10-x/x(7;-r^) [m/year] 23 For Z = 100 m growth from initially 1 m to 10 m takes 2108 years for Ta - Tfi = 1 K, and correspondingly 2107 years if the temperature difference is 10 K. Fig. 6 shows growth times for spheres and cylindrical conduits at various depths nd with various diameters. A final statement must be given. In this section we have assumed that temperature difference between the cave and the surface and 100% relative humidity of the cave air are independent of time. This may be the case when hydrothermal waters form open surfaces, from which evaporating waters condense at cave walls, delivering a constant flow of water in equilibrium with the pCO2 of the cave atmosphere. These waters attain equilibrium with respect to the soluble rock, either limestone or gypsum (see appendix). The saturation concentration for limestone is c = 60mg/L ■ ^l^^jJJ^ mg/L, where pCO2 is the partial pressure of CO2 in the cave atmosphere, measured in atm. Our findings of Eq. 21 explain the existence of spherical niches and cupola in caves of Hungary (Muller, 1974) in Italy, (Cigna & Forti 1986), and in hypogenetic caves in France (Audra et al. 2002). Sarbu & Lascu (1997) report on the measurements of active condensation corrosion in Movile Cave, Romania. Mostly in nature external boundary conditions are not constant in time. Annual fluctuations in the temperature of the cave air might occur when a river flows into a cave, with warm water in summer, which evaporates and condenses at the cave walls. In the winter time, however, the water is colder than the temperature at the cave walls and condensation stops. In view, that to attain stationary state under time independent boundary conditions takes several years, such cases cannot be described by the considerations above. 0.0 £=10"'m CONDUIT (Z=25 m) INITIAL DIAMETER 0.1 m 0.5 m 1.0 m _L 5.0x10® 1.0x10' Time[y] 1.5x10' 2.0x10' Fig. 6: Diameters of spheres and cylindrical conduits with initial diameter of 0.1, 0.5 and 1 m at depth of 25 m as a function of time. The curves are calculated from equations 21 and 23, respectively. CONDENSATION AT PERIODIC BOUNDARY CONDITIONS As already shown Fig. 4a presents the time dependence of the rate of condensing water. A first rapid decline until time t=Tf (see Eq.12) / I \2 1 I k£) 1 ^ 24 is followed by a slow decline (see Eq. 13) exhibiting a t-"2- time dependence. After the time t = 5Tg stationary state is attained (see Eq. 15). Note that the time to reach stationary state is controlled by the depth Z of the cave below ground and the thermal diffusivity K, whereas the decline toward stationary state depends solely on the thermal properties k , k, and kr and is independent of the depth Z. (See Eq.13). Variations of the cave temperature can be caused by diurnal fluctuations when warm humid air flows through the cave during day time, but during the night cold dry air enters. If the initial temperature of the rock is Tfi and the temperature of the cave air with a humidity of 100%, is Ta within the time t = 9(krs/k)2/K the temperature Tf of the rock has increased such that Ta-Tf = 0.18(Ta-Tfi). From then on Ta-Tf is given by Eq. 13. Note that for s = 0.001 m, t = 9(krs/k)2/K = 5.5-103 s whereas the equilibration time 5ts= 9105 s for Z = 1 m. Therefore at depths Z > 1m daily variations in Ta cause temperature changes of the rock, which are governed by the entity (kr/k). When cold air, at temperature Tn colder than the actual temperature Tf at the cave walls enters into the cave during the night, condensation stops. Two possibilities can be envisaged. The water, which has condensed during day time evaporates during the night. In this case, in a first approximation, the temperature of the cave wall drops towards the temperature of the cold air in about the same time as is needed to approach to the temperature of the warm air during day time. This is shown by Fig. 7a. As a result calcium carbonate dissolved during the day, will be precipitated during the night, and this process disintegrates the texture of the rock, leaving a weathered rind of corroded material (Auler & Smart 2004). For s = 0.001 m, the amount of water condensing within 10 hours is about 100 g/(m2K) as can be obtained by integrating Eqs. 12 and 13. See Fig. 4b, which shows the amount of condensed water per square meter as a function of time for various values of s and a temperature difference Ta-Tfi of 1 K. In the other extreme, when the condensed water flows sufficiently fast from the rock, evaporation is excluded. Heat transfer from the rock through the boundary layer is affected only by thermal conduction. The heat transfer coefficient ka/sT (see Eq.5) is therefore lower by a factor of 2.7 compared to that, when condensation is present (see Eq. 6). The time dependence of cooling or heating is given by the dimensionless variable k(Kt)"2 /(krs). With k lower by a factor of 2.7 cooling takes a time, longer by a factor of 7 than heating. Therefore during the cooling period the wall of the rock cools down slower. This way, after several cycles a stationary amplitude is established with constant temperature differences of the rock wall, and condensation rates are reduced. Fig. 7b shows this schematic concept. These two scenarios are extremes. Depending on the relative humidity of the air during the day and during the night, rates will be between these two extreme limits. We will discuss this latter in detail. Seasonal variations can be caused when warm water from the surface enters into the cave in summer time and condensation takes place, whereas in winter time, when cold water flows into the 16 -i—1—I—1—I—1—I—1—I—1—I—1—I—1—I—1— 14 ........"7 12 10 ■ikj ......r ..... .....r .....n ......T £=0.001 m ...ikj.....iy....iy...jy....!k a) With evaporation 4 5 6 Time [days] 10 ----T £=0.001 m No evaporation b) 1_I_1_I_1_I_1_I_1_lL 8 10 12 14 16 18 20 Time [days] Fig. 7: Wall temperature for diurnal temperature variations of Ta in the cave. For simplicity an abrupt change between day temperature Tay and lower night temperature T^^^^^ is assumed. The initial temperature of the cave wall Tfi 4104 s, where tc is the period of condensation in seconds. For tc = 8.4106 s ^100 days this corresponds to about 1500 g/(m2K), equivalent to a retreat of cave wall by 3.510-8 m/(yK). During winter time condensation stops and the cave walls have sufficient time to cool to low temperatures. Then in summer condensation starts again. Note that so far all calculations of retreat of rock are based on the assumption that pCO2 in the cave is at atmospheric level with 3.510-4 atm, and that the solution flowing off the wall has attained saturation with respect to calcite. For elevated values of pCO2 in the cave all numbers given so far must be multiplied by (pCO2/0.00035)"3 to account for elevated pressure pCO2 in the cave. CONDENSATION CORROSION ON SPELEOTHEMS Many researchers have observed features of surface corrosion on speleothems, which they interpret as a results of condensation corrosion (Auler & Smart 2004; Dublyansky & Dublyansky 2000; Tarhule-Lips & Ford 1998). In this section the physical background is discussed. Fig. 8 depicts the thermal boundary conditions for a stalagmite. At the outer surface heat flux is given by the action of condensation. At the base heat is transferred into the colder base rock. As can be calculated from Fig. 8 this heat flow is negligible. The decay time which is needed in such a case to approach thermal equilibrium (Luikov 1968), page 217) is TStal = D/(4^j2k), where is '2tSk the first root of the equatio^g n =-- ■ /j, . For D = 1 m, is n2/4 and TStal = 2 days. When thermal equilibrium is reached after 5ts condensation stops. Therefore to keep condensation active, diurnal variations in cave temperature are necessary. When the equilibration time 5TStal is longer than a day, and when during the night evaporation takes place a stationary state with high average rates is reached after some days. Lower rates are effective when evaporation is excluded during the night. We restrict to a more simple approximation. We assume stalagmites with diameters less than 0.3 m corresponding to TStal = 0.2 day. In this case the stalagmite practically reaches temperature Ta during day time. Even, if no evaporation is present during the night, the time to cool down is longer only a factor of 2. This results from the dependence of the root ^ on the heat transfer coefficient k (Luikov 1968, page 217). The maximal amount of condensing water during one day can be estimated from conservation of energy. The total heat transferred to the stalagmite must be equal to the increase of internal energy due to heating from the initial temperature T at t = 0 to the stationary temperature Ta. According to Eqs. 4 and 5 about 70% of the total heat transfer results from condensation. Therefore to a sufficiently good approximation one has where Mc is the mass of condensed water during 1 day, c is the specific heat of limestone (0.88 kJ/kg K) and AT is the change of temperature until the end of heating. It is close to AT = (Tday-Tn.ght). Ms is the mass of the stalagmite with diameter D and length L. From Eq. 26 one finds the amount of water which condenses during one period of heating. Dividing this by the surface area of the stalagmite gives the amount of water condensed per m2 during the heating period as F„ = — • ^ pDAT « 0.11 DAT 4 qc kg/m' this corresponds to retreat of surface by =3.8-10"'DAT [m/day] or 5^=1.4-10^Z)Ar [m/year] for D < 0.3 m. Note that for compact stalagmites we use the density of 2700 kg/m3. 27 28 n— D=2r oz Fig. 8: Heat rates transferred via corresponding surface of a stalagmite with length L and radius r. k=0.07 Wm-'K-', kr=1.3 Wm-'K-'. Even at a temperature gradient of 10K/m at the base, F3<0.1K. T is the surface temperature of the stalagmite. For larger diameter D > 1 m the equilibration times for heating and cooling are longer than 5 days. If such a stalagmite experiences condensational heating with equilibration time Tj= TStal during the time tj and cooling with equilibration time t2, during the time t2 after several cycles the maximum temperature and the minimum temperature become independent on time. The temperature difference Tmax - Tmin can be estimated in the following way. During heating the temperature Th(t) of stalagmite to an acceptable approximation is given by 1-exp 29 where Tmin is the minimum temperature at the end of prior cooling, and Th is the temperature of the cave air during heating. For the cooling period t2 with equilibration time t2 the temperature Tc of the stalagmite is given by / t 30 where Tc is the air temperature in cave during cooling and Tmax the temperature of the stalagmite at the end of prior heating. From 29 and 30 we get Solving this equation one finds (l-expH,/ir,))(l-expH/T,)) 31 32 l-exp(-?i/Ti)exp(-?2/Tj) To estimate the amount of condensation corrosion, this value of AT must be used in Eq. 28. AT in Eq. 32 becomes 0.25^^7) for i = -^ = 0.5 and 0.12^^?;) for _i = .^ = o.25 . Ti Tj For ^^»tj expansion of the exponents and regarding t^ = TStal yields Te*^/ U Therefore one finds from Eq. 28 tKJl D =3.3 10 .7 (TWD D [m/year] [m/year], with tj = 12 h 33 34 It should be noted that this number is valid for D > l m and presents an upper limit, because we have assumed tacitly that the temperature inside the stalagmite is homogenous and equal to the surface temperature. Eqs. 28 and 34 show that for stalagmites with diameters D between 0.1 m to several meters and temperature differences during night and day of 10°C, corrosion rates are on the order of 10-4 to 10-5 cm/year. Auler and Smart (2004) have estimated rates of corrosion on stalagmites by measuring the depth of the weathered rind and determining the age of the unaltered speleothem calcite below. The highest values they observed was 410-5 cm/year. In most cases the rates were lower by one to two orders of magnitude. FIELD EXPERIMENTS TO MEASURE CONDENSATION CORROSION Tarhule-Lips and Ford (1998) suspended gypsum plates of about 1 cm in thickness on nylon strings for about one year in flank-margin caves of the Caribbean. From the measured weight loss they report corrosion rates of 2.4 ■ 10-2 cm/year. The thermal behavior of such isolated plates is very similar to that of stalagmites with diurnal variation of temperature and Eq. 28 remains valid if one replaces the diameter D by the thickness of the plate. From this one obtains ^^ ^ • 10"' cm/year, assuming AT=10 K. This is three orders of magnitude lower than the experimental findings. It should be noted that these gypsum plates reach thermal equilibrium after a time of only a few minutes. After this time condensation stops and renewed cooling and subsequent heating are necessary to revive it. Sarbu and Lascu (1997) report condensation rates of water in Movile-Cave, Romania, where they collected water from a 10 cm by 10 cm glass plate suspended at a distance of 10 cm from the cave wall. Movile cave is heated by a hydrothermal lake with water temperature of 21°C. In its upper level, where the glass plate was suspended air temperature is 21°C, but the temperature of the cave walls is between 19.4° and 15.7°C. In this case the glass plate achieves thermal equilibrium after a few minutes. The observed high run-off rates of about 20 g/month cannot be explained by condensation. From Eq. 27 one obtains only 2.5 g/month for a plate of 1 cm thickness. Summarizing, we state that condensation to small scale objects (0.1-1 cm) is subject to diurnal variations of temperature. In thermal stable cave environments, or where variations of temperature change seasonally, it can be excluded. RÖHRENKARREN, A SMALL SCALE EXAMPLE OF CONDENSATION CORROSION Recently Simms (2003) reported on vertical, upward tapering tubes in limestone exposed in the epiphreatic zone at the shores of several lakes in Ireland. The dimensions of these solutional features comprising almost perfect circular tubes are between 1 to 5 cm in diameter and up to 30 cm in length. Simms suggests that these "Röhrenkarren" originate from condensation corrosion within air pockets trapped by seasonal high stands of the lake. During March 2000 water temperature TW stayed nearly constant, whereas the temperature of the overlying rock showed diurnal variations, fluctuating around the water temperature with amplitude of about 5°C. In winter time, when surface temperatures become low extended periods of rock temperatures TR below water temperature are likely. Therefore conditions for condensation corrosion to operate are valid. The lake water is saturated with respect to calcite such that dissolution of the rock in contact with water is excluded. Since the region of the lakes was glaciated about 15 ky ago, rates of dissolution were estimated about 2.10-3 cm/year, if one assumes continuous dissolution to present. With this information it is possible to test our theoretical predictions. The distance between the surface of the rock to the apex of the Röhrenkarren is on the order of meters. Therefore thermal equilibrium is attained in a few days. During winter time an extended period of rock temperatures TR several degrees below that of the lake water supports condensation in a steady state. To estimate the rates we use Eq. 19. In the initial state of the evolution of the Röhrenkarren air is entrapped on irregularities of the rock. The geometry of this bubble could be approximated by a half-sphere with radius R as depicted in Fig. 9. The air entrapped is stagnant, due to the closed conditions. Therefore, diffusion of vapor from the water surface to the rock wall is through this stagnant air. s in Eq. 19 must therefore be replaced by the length L of the tube. The retreat of bedrock is focused to the apex, because heat flow is highest there, whereas the walls exhibit lower heat flow. This way, a circular tube can propagate upwards into the bedrock. Using Eq. 19 one finds — = 1.3510-^-dt [m/year] 35 Initially L = R, therefore to a good approximation dt = 1.3510 [m/year] 36 Initially with L = 0.01 m the growth rate is 1.310-4 cm/year K. Since condensation is active only during cold weather this must be reduced by a factor of 0.3 or so, to about 510-5 cm/(year K). When the length increases the growth rate drops with 1/L. During the growth of the tube the entrapped air is in contact with the lake water and consequently pCO2 in this air is in equilibrium with the pCO2 in the water, which can enhance dissolution rates by a factor of two for pCO2=210-3 atm, likely in lake water. Therefore with a temperature difference of 10 K maximal rates are 10-3 cm/year, dropping to 10-4 cm/year when a depth L of 10 cm is reached. Rock temperature T, Lake water with temperature Tw>Tp Fig. 9: Röhrenkarren: A slab of rock with a depth of a few meters is exposed to the surface. Its bottom is flooded by the water from the lake. This way air is entrapped. An initial irregularity of radius R (dotted line) grows into a cylindrical shape of length L. The temperature of the air inside is close to the water temperature. In winter time the surface temperature is lower than that of the lake water and condensation is active. DISCUSSION To illustrate our theoretical results we give a numerical solution of a representative example. Fig. 10 represents a rectangular cave conduit parallel to the surface of a limestone plateau. It is located at a depth 25 m below ground and its cross section is 10 m x 10 m. Inside this cave the air temperature is 10°C and relative humidity is 100%. The temperature at the surface is 0°C. Note that only the temperature difference is significant. The other boundaries are assumed to be adiabatic, and cannot transmit heat. This is an approximation, which is valid for t < y2/K, where y is the distance of the cave to the outer limit. t = y2/K is the time when the thermal front reaches the adiabatic boundaries. We solve the differential heat conduction equation dt = K d/ ^ dz' 37 with the boundary conditions, as discussed above by a finite difference program. Fig. 11 illustrates the results. It shows isotherms as they evolve in time. At the beginning (1106 s, Fig. 11a) an almost circular temperature field has developed, symmetrical around the conduit. The temperature has changed only in the vicinity of the cave. After 107 s (Fig. 11b) the distance the temperature front has propagated is about, (107k)"2 = 2.4 m, close to what we see in Fig. 11a. After 2108 s (Fig. 11c). After 1109 s (Fig. 11e) a thermal gradient develops, directed from ceiling of the cave toward the surface. Most of the heat from the ceiling flows vertically to the surface. This can be seen from the flow lines depicted in Fig. 11. In this region the temperature distribution becomes stable in time, as can be visualized in Figs. 11d, e, f. Keeping in mind that t = ^^ (confer Eq. 15) this is a reasonable result. Constant temperature (T=TJ E o To=0 °C T=10°C I CO