the wegener-bergeron-findeisen process – its …...prepub article, 2015 t. storelvmo & i. tan:...

7
proof B Meteorologische Zeitschrift, PrePub DOI 10.1127/metz/2015/0626 MetZet Classic Papers © 2015 The authors The Wegener-Bergeron-Findeisen process – Its discovery and vital importance for weather and climate Trude Storelvmo and Ivy Tan Department of Geology and Geophysics, Yale University, New Haven, USA (Manuscript received July 11, 2014; in revised form October 16, 2014; accepted October 20, 2014) Abstract The Wegener-Bergeron-Findeisen process refers to the rapid growth of ice crystals at the expense of surrounding cloud droplets, which frequently occurs in atmospheric mixed-phase clouds. The process is a result of the difference in saturation vapor pressures with respect to liquid and ice, and may in some circumstances lead to abrupt and complete cloud glaciation at temperatures between 40 °C and 0 °C in the Earth’s atmosphere. The process is named after three eminent scientists who were active in the first half of the 20 th century, among them being German meteorologist Walter Findeisen (1909–1945). In his classical paper published in 1938, Findeisen described the contemporary understanding of the Wegener-Bergeron-Findeisen process and other key cloud microphysical processes. Here, we compare the understanding of aforementioned processes at the time with that of the present, and find that they are remarkably similar. We also discuss how the Wegener-Bergeron-Findeisen process is implemented in state-of-the-art numerical models of the atmosphere, and highlight its importance for both weather and climate. Keywords: Wegener-Bergeron-Findeisen process, mixed-phase clouds, weather and climate 1 Introduction 1 The importance of the Wegener-Bergeron-Findeisen 2 (WBF) process, or simply the Bergeron-Findeisen pro- 3 cess, is well known to meteorologists and climatologists 4 alike. By abruptly transforming non-precipitating liquid 5 clouds to heavily precipitating ice clouds and dramati- 6 cally changing cloud radiative properties, it can have a 7 profound impact on both weather and climate. Its dis- 8 covery, which dates back almost a century ago, should 9 be equally accredited to three eminent scientists of the 10 time: Alfred Wegener, Tor Bergeron and Walter 11 Findeisen. 12 The German scientist Alfred Wegener (1880– 13 1930), well-known for his then-controversial theory on 14 continental drift, first laid the theoretical foundation for 15 the WBF process (Wegener, 1911), by showing that 16 the co-existence of liquid and ice is a thermodynami- 17 cally unstable state. This revelation allegedly came to 18 Wegener while studying the formation of hoarfrost. A 19 decade later, in the winter of 1922, the Swede Tor Berg- 20 eron (1891–1977) found himself pondering the theory 21 put forth in Wegener’s book during a stay at a health 22 resort in Voksenkollen (430m above sea level) outside 23 of Oslo in Norway. Observant as he was, Bergeron 24 had noticed that when the temperature was below freez- 25 ing, nearby forest roads were clear of fog while trees 26 were covered in frost. Fog, however, would typically be 27 present and extend all the way to the ground when tem- 28 peratures were above 0 °C. As an active member of the 29 prestigious Bergen School of Meteorology, Bergeron 30 Corresponding author: Trude Storelvmo, Department of Geology and Geophysics, Yale University, New Haven, CT06511, USA, e-mail: [email protected] became immersed in his duties in Bergen in the years 31 following his discovery in Voksenkollen, so much to the 32 extent that he did not further pursue his ideas on the mat- 33 ter until 1928, when the topic became one of the chapters 34 in his PhD thesis (Bergeron, 1928). 35 It wasn’t until 1938 that Walter Findeisen entered 36 the scene, contributing to the previous work of We- 37 gener and Bergeron by providing additional theoret- 38 ical calculations, as well as cloud chamber experiments 39 to further develop their theories. Findeisen’s PhD the- 40 sis (1931) focused on cloud droplet size distributions, 41 and included cloud chamber experiments, a novel ap- 42 proach at the time. Findeisen’s cloud chamber was ap- 43 proximately 2 m 3 in volume, and was connected to a 44 vacuum pump, allowing the process of adiabatic expan- 45 sion and atmospheric cloud formation to be mimicked 46 in the chamber. After World War II, Findeisen’s cloud 47 chamber was recovered from the ruins of Prague, where 48 Findeisen had his last appointment as director of the 49 Prague branch of the German Meteorological Office. 50 Fig. 1 shows the rebuilt cloud chamber, as it appeared 51 in Podzimek (1957). 52 Findeisen frequently cites the work of Wegener 53 and Bergeron in his seminal paper from 1938 (Find- 54 eisen, 1938, hereafter F38), whose work allowed him 55 to present a coherent and comprehensive overview of 56 the most recent understanding of atmospheric cloud and 57 precipitation formation at the time. As such, the paper 58 goes beyond the WBF process that Findeisen later be- 59 came known for, and can in many ways be considered 60 the first complete description of cloud microphysics as 61 we understand it today. 62 Atmospheric scientists today are privileged to con- 63 duct research in an age as data-rich as the present, with 64 © 2015 The authors DOI 10.1127/metz/2015/0626 Gebrüder Borntraeger Science Publishers, Stuttgart, www.borntraeger-cramer.com

Upload: others

Post on 27-Jun-2020

1 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: The Wegener-Bergeron-Findeisen process – Its …...PrePub Article, 2015 T. Storelvmo & I. Tan: The Wegener-Bergeron-Findeisen process 3 Figure 2: Left: Saturation vapor pressure

proof

BMeteorologische Zeitschrift, PrePub DOI 10.1127/metz/2015/0626 MetZet Classic Papers© 2015 The authors

The Wegener-Bergeron-Findeisen process – Its discovery andvital importance for weather and climateTrude Storelvmo∗ and Ivy Tan

Department of Geology and Geophysics, Yale University, New Haven, USA

(Manuscript received July 11, 2014; in revised form October 16, 2014; accepted October 20, 2014)

AbstractThe Wegener-Bergeron-Findeisen process refers to the rapid growth of ice crystals at the expense of surrounding cloud droplets,which frequently occurs in atmospheric mixed-phase clouds. The process is a result of the difference in saturation vapor pressureswith respect to liquid and ice, and may in some circumstances lead to abrupt and complete cloud glaciation at temperatures between−40 °C and 0 °C in the Earth’s atmosphere. The process is named after three eminent scientists who were active in the first half of the20th century, among them being German meteorologist Walter Findeisen (1909–1945). In his classical paper published in 1938,Findeisen described the contemporary understanding of the Wegener-Bergeron-Findeisen process and other key cloud microphysicalprocesses. Here, we compare the understanding of aforementioned processes at the time with that of the present, and find that theyare remarkably similar. We also discuss how the Wegener-Bergeron-Findeisen process is implemented in state-of-the-art numericalmodels of the atmosphere, and highlight its importance for both weather and climate.

Keywords: Wegener-Bergeron-Findeisen process, mixed-phase clouds, weather and climate

1 Introduction1

The importance of the Wegener-Bergeron-Findeisen2

(WBF) process, or simply the Bergeron-Findeisen pro-3

cess, is well known to meteorologists and climatologists4

alike. By abruptly transforming non-precipitating liquid5

clouds to heavily precipitating ice clouds and dramati-6

cally changing cloud radiative properties, it can have a7

profound impact on both weather and climate. Its dis-8

covery, which dates back almost a century ago, should9

be equally accredited to three eminent scientists of the10

time: Alfred Wegener, Tor Bergeron and Walter11

Findeisen.12

The German scientist Alfred Wegener (1880–13

1930), well-known for his then-controversial theory on14

continental drift, first laid the theoretical foundation for15

the WBF process (Wegener, 1911), by showing that16

the co-existence of liquid and ice is a thermodynami-17

cally unstable state. This revelation allegedly came to18

Wegener while studying the formation of hoarfrost. A19

decade later, in the winter of 1922, the Swede Tor Berg-20

eron (1891–1977) found himself pondering the theory21

put forth in Wegener’s book during a stay at a health22

resort in Voksenkollen (430 m above sea level) outside23

of Oslo in Norway. Observant as he was, Bergeron24

had noticed that when the temperature was below freez-25

ing, nearby forest roads were clear of fog while trees26

were covered in frost. Fog, however, would typically be27

present and extend all the way to the ground when tem-28

peratures were above 0 °C. As an active member of the29

prestigious Bergen School of Meteorology, Bergeron30

∗Corresponding author: Trude Storelvmo, Department of Geology andGeophysics, Yale University, New Haven, CT06511, USA, e-mail:[email protected]

became immersed in his duties in Bergen in the years 31

following his discovery in Voksenkollen, so much to the 32

extent that he did not further pursue his ideas on the mat- 33

ter until 1928, when the topic became one of the chapters 34

in his PhD thesis (Bergeron, 1928). 35

It wasn’t until 1938 that Walter Findeisen entered 36

the scene, contributing to the previous work of We- 37

gener and Bergeron by providing additional theoret- 38

ical calculations, as well as cloud chamber experiments 39

to further develop their theories. Findeisen’s PhD the- 40

sis (1931) focused on cloud droplet size distributions, 41

and included cloud chamber experiments, a novel ap- 42

proach at the time. Findeisen’s cloud chamber was ap- 43

proximately 2 m3 in volume, and was connected to a 44

vacuum pump, allowing the process of adiabatic expan- 45

sion and atmospheric cloud formation to be mimicked 46

in the chamber. After World War II, Findeisen’s cloud 47

chamber was recovered from the ruins of Prague, where 48

Findeisen had his last appointment as director of the 49

Prague branch of the German Meteorological Office. 50

Fig. 1 shows the rebuilt cloud chamber, as it appeared 51

in Podzimek (1957). 52

Findeisen frequently cites the work of Wegener 53

and Bergeron in his seminal paper from 1938 (Find- 54

eisen, 1938, hereafter F38), whose work allowed him 55

to present a coherent and comprehensive overview of 56

the most recent understanding of atmospheric cloud and 57

precipitation formation at the time. As such, the paper 58

goes beyond the WBF process that Findeisen later be- 59

came known for, and can in many ways be considered 60

the first complete description of cloud microphysics as 61

we understand it today. 62

Atmospheric scientists today are privileged to con- 63

duct research in an age as data-rich as the present, with 64

© 2015 The authorsDOI 10.1127/metz/2015/0626 Gebrüder Borntraeger Science Publishers, Stuttgart, www.borntraeger-cramer.com

Page 2: The Wegener-Bergeron-Findeisen process – Its …...PrePub Article, 2015 T. Storelvmo & I. Tan: The Wegener-Bergeron-Findeisen process 3 Figure 2: Left: Saturation vapor pressure

proof

2 T. Storelvmo & I. Tan: The Wegener-Bergeron-Findeisen process Meteorol. Z., PrePub Article, 2015

Figure 1: Walter Findeisen’s cloud chamber, as it appeared inthe publication that reported the first successful measurements withthe chamber after it was recovered and rebuilt after World War II(Podzimek, 1957). Published with permission from ©Springer.

in situ measurements, remote sensing observations, lab-65

oratory work and numerical modeling all contributing to66

this wealth of data. By contrast, the earliest works on67

cloud microphysics by Wegener, Bergeron and Find-68

eisen had very little in situ observations and of course69

no satellite data available to them. It is a testimony to70

their brilliance that they nevertheless came to many of71

the same valid conclusions we find today. In fact, the72

understanding of the microphysical processes involved73

in liquid and ice cloud formation and subsequent cloud74

evolution has changed relatively little relative to that pre-75

sented by F38. To give an example, in F38 Findeisen76

claimed that every precipitation event of at least medium77

intensity, and in particular every event with larger rain-78

drops, is caused by ice crystals. A more recent study79

based on satellite data has confirmed the truth of this80

statement, and found that the majority (∼ 70 %) of trop-81

ical precipitation events indeed originate from the ice82

phase (Lau and Wu, 2003).83

The remainder of this paper is dedicated to compar-84

ing the microphysical processes as they were described85

in F38 with our understanding of them today. Section 286

describes the theoretical foundation for the WBF pro-87

cess, Section 3 discusses the co-existence of supercooled88

liquid (i.e. liquid water existing at temperatures below89

0 °C) and ice crystals required for the WBF process to90

occur, and Section 4 describes how the WBF process 91

is implemented in state-of-the-art numerical models of 92

weather and climate. Finally, Section 5 offers a brief 93

conclusion. 94

2 The Wegener-Bergeron-Findeisen 95

process 96

The WBF process refers to the rapid conversion of liq- 97

uid to ice that may occur when supercooled droplets and 98

ice crystals co-exist (Pruppacher and Klett 2010). 99

The conversion occurs due to the difference in satu- 100

ration vapor pressures over liquid and ice surfaces at 101

temperatures below 273 K (el and ei, respectively, with 102

el > ei), which can be approximated from the Clausius- 103

Clapeyron relation. In other words, an environment that 104

is saturated with respect to liquid water will be highly 105

supersaturated with respect to ice, and the relative dif- 106

ference in supersaturation is exacerbated with decreas- 107

ing temperature (Fig. 2). A common misconception is 108

that the WBF process is automatically activated when 109

liquid and ice co-exist, i.e. ice crystals are guaranteed 110

to grow at the expense of cloud droplets without excep- 111

tion. However, as pointed out for example by Korolev 112

(2007) and Korolev and Mazin (2003), the WBF pro- 113

cess is only one of three possible cases that may occur 114

when a cloud consists of liquid and ice. The WBF pro- 115

cess (i) will occur when the vapor pressure (e) lies be- 116

tween ei and el. The other two possible cases involve 117

either (ii) simultaneous growth of liquid droplets and 118

ice crystals (e > el > ei) or (iii) simultaneous evapo- 119

ration/sublimation of cloud droplets (e < ei < el). Cloud 120

dynamics in the form of small-scale updrafts and down- 121

drafts exert an important control over which case plays 122

out for a given mixture of droplets and crystals. Find- 123

eisen understood this, and described case (ii) in Sec- 124

tion 4 of his 1938 paper: sufficient updrafts and hence 125

adiabatic cooling can result in the counter-intuitive pro- 126

cess of droplet formation and growth in a glaciated 127

cloud. The reason why sufficiently high supersaturations 128

for droplet formation can occur in ice clouds was also 129

addressed in F38: the nuclei on which ice crystals form 130

in the atmosphere are very rare relative to the nuclei that 131

cloud droplets nucleate on. Hence, even though ice crys- 132

tals are present in a cloud, they may not be present in 133

high enough number concentrations for their growth to 134

deplete supersaturation faster than the rate at which high 135

supersaturation is produced via adiabatic cooling. 136

3 Liquid in the supercooled state and 137

the scarcity of ice nuclei 138

The distinction between the two classes of nuclei, cloud 139

condensation nuclei (CCN) and ice nuclei (IN), was of- 140

fered by F38 as an explanation for the frequent obser- 141

vations of supercooled liquid water in the atmosphere 142

that had been reported at the time. While the reported 143

Page 3: The Wegener-Bergeron-Findeisen process – Its …...PrePub Article, 2015 T. Storelvmo & I. Tan: The Wegener-Bergeron-Findeisen process 3 Figure 2: Left: Saturation vapor pressure

proof

Meteorol. Z., PrePub Article, 2015 T. Storelvmo & I. Tan: The Wegener-Bergeron-Findeisen process 3

Figure 2: Left: Saturation vapor pressure over bulk liquid (el) and over bulk ice (ei) as a function of temperature, calculated using theMagnus formula (Magnus, 1844). Right: Absolute (red line) and relative difference (black line) between el and ei, the latter given as(el/ei − 1) · 100 %, all as functions of temperature. The black line corresponds to the supersaturation that would be experienced by an icecrystal forming in a supercooled liquid cloud under the assumption that the water phase is in equilibrium with the liquid phase at the timeof crystal formation.

Figure 3: Supercooled cloud fraction (SCF, %) based on CALIOP retrievals; global mean (green), the dust belt (blue, average over theregion 0–120 ° W and 30–50 ° N) and the Southern Ocean (red, average over the region 0–360 ° W and 60–70 ° S). For further details on howSCFs were calculated based on the CALIOP retrievals, and associated uncertainties, see Tan et al. (2014a).

observations were naturally sporadic and few, they sup-144

ported the existence of liquid at temperatures even be-145

low −20 °C. Now, several decades into the satellite146

era, we are able to take advantage of global datasets147

that provide information on cloud thermodynamic phase148

with relatively high temporal coverage. An example is149

shown in Fig. 3, displaying the observed supercooled150

cloud fraction (SCF, in %) as a function of tempera-151

ture, as retrieved by the Cloud and Aerosol Lidar with152

Orthogonal Polarization (CALIOP) instrument onboard 153

NASA’s Cloud-Aerosol Lidar and Infrared Pathfinder 154

Satellite Observations (CALIPSO) satellite (Hu et al. 155

2009, Tan, Storelvmo and Choi 2014b, Winker et al. 156

2009). These CALIOP retrievals are representative of 157

cloud tops only, and the SCF was calculated by taking 158

the ratio of liquid cloud top pixels to total cloud top pix- 159

els within 2.0 ° longitude by 2.5 ° latitude grid boxes. 160

The cloud top temperatures were determined using the 161

Page 4: The Wegener-Bergeron-Findeisen process – Its …...PrePub Article, 2015 T. Storelvmo & I. Tan: The Wegener-Bergeron-Findeisen process 3 Figure 2: Left: Saturation vapor pressure

proof

4 T. Storelvmo & I. Tan: The Wegener-Bergeron-Findeisen process Meteorol. Z., PrePub Article, 2015

NCEP-DOE Reanalysis 2 dataset (Kanamitsu et al.,162

2002). The satellite observations support a global aver-163

age liquid cloud fraction of ∼ 30 % at −20 °C, in agree-164

ment with the findings from a century ago. Furthermore,165

a comparison of the SCF in the Southern Ocean region166

(0–360 ° W, 60–70 ° S) with those found in the so-called167

Dust belt (0–120 ° W, 30–50 ° N), suggests that SCF is168

spatially heterogeneous and that although certainly im-169

portant, temperature is not the sole factor influencing170

cloud phase.171

Substantial amounts of liquid exist over the Southern172

Ocean, even at temperatures as low as ∼ −30 °C, pre-173

sumably owing to the scarcity of IN that are required174

to initiate freezing at temperatures above approximately175

−40 °C in the atmosphere. Findeisen was aware of this176

and had documented it in Section 1 of F38. He estimated177

that the ratio of CCN to IN number concentrations was178

on the order of 104, a number that matches current obser-179

vations using so-called CCN and IN-counters (Hudson180

and Squires, 1974; Roberts and Nenes, 2005; Rogers181

1988; Stetzer et al., 2008), instrumentation that was182

not available during Findeisen’s time. However, Find-183

eisen did state that it would be possible to design such184

instruments and remarked that they had the potential to185

“clarify the controversial questions” presented in his pa-186

per.187

Given the lack of instrumentation at the time, Find-188

eisen’s description of the properties of CCN versus189

those of IN is remarkably similar to that of any con-190

temporary paper on the subject. In a recent review of191

atmospherically relevant IN, Murray et al. (2012) reaf-192

firmed Findeisen’s description of IN (referred to as sub-193

limation nuclei by Findeisen) in F38 as insoluble par-194

ticles of terrestrial origin, mainly in the form of min-195

eral dust (quartz, according to Findeisen). Beyond this,196

Findeisen also stated that the chemical composition and197

origin of IN are largely unknown, a statement that to198

some degree still holds today. Notwithstanding Find-199

eisen’s pioneering discoveries, atmospheric ice nucle-200

ation is currently a very active field of research, and201

our understanding of what particles are able to act as202

IN under what conditions is rapidly evolving. We now203

know that certain mineral dust types are better at nucle-204

ating ice than others, and that quartz is not a particularly205

good IN (Atkinson et al., 2013). We also know that bio-206

logical particles and potentially anthropogenic particles207

such as soot, ash and metallic particles (Cziczo et al.,208

2009; Hoose and Mohler, 2012) may also be acting as209

IN in the atmosphere.210

In addition to his work on the WBF process, Find-211

eisen had also performed laboratory work on cloud212

droplet formation (F38) at approximately the same time213

that Hilding Köhler was developing the relatively214

straightforward theory of cloud droplet formation by215

the so-called process of “CCN activation” in 1936, now216

sometimes referred to as Köhler Theory. Findeisen’s217

untimely death towards the end of World War II in 1945218

at the tender age of 36 in Prague, meant that Köhler,219

outliving him, was able to influence the field for decades220

thereafter, perhaps resulting in the different legacies of 221

the two scientists. 222

4 Representations of the WBF process 223

in numerical weather and climate 224

models 225

Returning to the WBF process that brought Findeisen 226

fame, its importance for weather and climate has in- 227

creasingly attracted attention in recent years. A real- 228

istic representation of the WBF process in numerical 229

weather prediction (NWP) and global climate models 230

(GCMs) is critical for more accurate simulations of at- 231

mospheric dynamical and radiative processes, and hence 232

the climate system as a whole. The typical horizon- 233

tal resolution of such models is on the order of 10 to 234

100 km, while the WBF process occurs on scales orders 235

of magnitude smaller. Such unresolved processes pose a 236

challenge for numerical models of weather and climate. 237

The impact of these small-scale processes on resolved 238

large-scale processes can be accounted for by includ- 239

ing parameterizations of the small-scale processes that 240

are otherwise unresolved. In recent years, new param- 241

eterizations with various levels of sophistication have 242

been developed. The simplest parameterizations impose 243

a critical threshold of in-cloud ice mixing ratio, above 244

which the WBF process is assumed to become efficient 245

enough to deplete all remaining liquid in the model 246

grid box within a single model time step (typically 247

∼ 30 min) (Storelvmo, Kristjansson and Lohmann 248

2008a, Lohmann and Hoose 2009). However, more re- 249

cent studies have attempted to treat the WBF process in a 250

more rigorous fashion. In a parameterization frequently 251

used in both GCMs and NWP models (Morrison et al. 252

2005), the WBF process is diagnosed based on the rate 253

of depositional growth of ice crystals, A, and the rate of 254

condensation of liquid, Q. If A > Q the WBF process is 255

assumed to deplete liquid water within the model’s time 256

step. While this approach is consistent with the under- 257

standing of how the WBF process operates in the atmo- 258

sphere, it is oversimplified in the sense that it assumes 259

that all cloud properties are uniform within the cloudy 260

portion of each model grid-box, which, in a GCM, typi- 261

cally spans ∼ 100 km in both longitudinal and latitudinal 262

directions. 263

A few studies have sought remedy for the afore- 264

mentioned oversimplified parameterizations of the WBF 265

process by introducing sub-gridscale variability in cloud 266

properties that are key to accurately representing the 267

WBF process (Rotstayn, Ryan and Katzfey 2000, 268

Storelvmo et al. 2008b, Storelvmo et al. 2010, Rot- 269

stayn 1997). In attempt to account for this sub-gridscale 270

variability, Rotstayn (1997) introduced a triangular 271

probability density function (PDF) for the total-water 272

mixing ratio, q, within each model grid box, following 273

Smith (1990). The PDF was centered at the grid box 274

mean total-water mixing ratio. Instead of considering 275

differences in vapor pressure, e, between the two phases 276

Page 5: The Wegener-Bergeron-Findeisen process – Its …...PrePub Article, 2015 T. Storelvmo & I. Tan: The Wegener-Bergeron-Findeisen process 3 Figure 2: Left: Saturation vapor pressure

proof

Meteorol. Z., PrePub Article, 2015 T. Storelvmo & I. Tan: The Wegener-Bergeron-Findeisen process 5

Figure 4: The fraction of cloud in which the WBF process is active, as a function of cloud droplet number concentration (CDNC, 105 m3)and ice crystal number concentration (ICNC, 102m−3), assuming a Gaussian PDF of w centered at 0.1 m/s with a standard deviation of0.2 m/s. Ice crystal and cloud

(as in Section 1), the corresponding difference in satu-277

ration vapor mixing ratio (qs,l and qs,i for liquid and ice,278

respectively) was used to determine the portion of the279

cloud that consists of co-existing ice and liquid. In a grid280

box containing both liquid and ice, coexistence would281

be possible for the portion of the grid box with q > qs,l,282

while the portion with qs,i < q < qs,l would have ice283

clouds only, and q < qs,i would correspond to cloud-free284

conditions. Note that this framework assumes that the285

cloud droplet and ice crystal response to sub-saturation286

is fast, and that complete evaporation occurs within one287

model time step (∼ 30 min) Following up on the work of288

Rotstayn (1997) and Smith (1990), Storelvmo et al.289

(2008b) implemented a normal distribution for the ver-290

tical velocity, w, in place of the triangular PDF for q291

used in Rotstayn (1997). Previously, Korolev (2007),292

Korolev and Mazin (2003) had derived parameteriza-293

tions for the critical updraft above which liquid and ice294

could co-exist (wc,u), and the critical downdraft below295

which both liquid and ice crystals are bound to evaporate296

(wc,d). By combining this with the PDF of w, the evo-297

lution of the thermodynamic phase of clouds can be di-298

vided into three distinct regimes: i) simultaneous growth299

of droplets and ice crystals, ii) growth of ice crystals300

at the expense of cloud droplets (the WBF process),301

and iii) simultaneous evaporation of droplets and ice 302

crystals. wc,u (always positive) and wc,d (always neg- 303

ative) are functions of ice crystal number concentra- 304

tion (ICNC) and cloud droplet number concentration 305

(CDNC), among other variables. Fig. 4 displays the frac- 306

tion of a cloud that will be dominated by the WBF pro- 307

cess (regime ii above) as a function of CDNC and ICNC, 308

calculated according to the formulae in Korolev (2007) 309

(for other assumptions made for the calculation, see the 310

caption of Fig. 4). At high CDNCs, saturation can still 311

be maintained in strong downdrafts by evaporating the 312

many cloud droplets present. The parameterization ac- 313

counts for this by allowing wc,d to become increasingly 314

negative, thereby causing the fraction of the cloud in 315

which ice crystals can grow at the expense of cloud 316

droplets to increase. At high ICNC, ice crystal growth 317

on the many ice crystals present rapidly depletes water 318

vapor and brings the vapor pressure below that of sat- 319

uration with respect to liquid water. In this case, very 320

strong updrafts are required for simultaneous growth of 321

droplets and ice crystals (i.e. wc,u is large). As a result, 322

the fraction of the cloud dominated by the WBF process 323

increases with increasing ICNCs. 324

Independently of how the WBF process is treated 325

in GCMs and/or NWP, the extent to which ice crystals 326

Page 6: The Wegener-Bergeron-Findeisen process – Its …...PrePub Article, 2015 T. Storelvmo & I. Tan: The Wegener-Bergeron-Findeisen process 3 Figure 2: Left: Saturation vapor pressure

proof

6 T. Storelvmo & I. Tan: The Wegener-Bergeron-Findeisen process Meteorol. Z., PrePub Article, 2015

Table 1: Net Cloud Radiative Effect (CRE) evaluated at the top of the atmosphere, Total Water Path (TWP), total precipitation and the ratioof stratiform to convective precipitation for simulations in which i) cloud phase is prescribed according to temperature (i.e. no representationof the WBF process, NO_WBF), ii) a crude critical ice mixing ratio threshold treatment (see above, SIMPLE_WBF) is applied and iii) aWBF treatment which accounts for subgrid scale variability (see above, SUBGRID_WBF) is implemented. Observations are from satelliteretrievals (Loeb et al., 2009; Komurcu et al. 2014).

Simulations NO_WBF SIMPLE_WBF SUBGRID_WBF OBSERVATIONS

Net CRE (Wm−2) −15.7 −24.9 −21.4 −24.5–17.9TWP (gm−2) 148.3 112.2 132.1 47–109Total precipitation (mm/day) 2.84 2.80 2.78 2.74Stratiform/convective precipitation ratio 0.37 0.54 0.51 N/A

and cloud droplets are assumed to be well-mixed with327

each other in the cloud volume, or whether they are as-328

sumed to exist in separate pockets of ice and liquid, is329

of critical importance (Korolev and Isaac, 2006). This330

statement can be illustrated with two sensitivity simu-331

lations1 using the CAM5 GCM (http://www.cesm.ucar.332

edu/models/cesm1.2/cam/) that attempt to mimick pure333

homogeneous (i.e. well-mixed) and pure heterogeneous334

mixing between the two phases (simulations HOM and335

HET, respectively). Between the two simulations, HOM336

yielded three times the amount of vertically integrated337

amount of ice in the atmosphere (Ice Water Path, IWP)338

relative to HET. The contrast between homogeneous and339

heterogeneous mixing was mimicked by reducing the ef-340

ficiency of the WBF process by a factor of 10−6, physi-341

cally corresponding to a situation in which there is only342

contact between the pockets of liquid and ice in rela-343

tively narrow mixing zones. High-frequency in situ mea-344

surements of cloud phase may aid in developing pa-345

rameterizations that realistically represent the degree of346

mixing that occurs within mixed-phase clouds, but such347

measurements are presently scarce.348

To illustrate how sensitive model-simulated cloud,349

radiation and precipitation fields can be to the repre-350

sentation of the WBF process, Table 4 shows global351

mean output from three 5-year simulations with the352

same atmospheric GCM (Storelvmo et al., 2008b). The353

simulations only differ in their treatment of the WBF354

treatment. Drastic changes in the net radiative effect of355

clouds, as well as in the column-integrated amount of356

liquid and ice in the atmosphere, are evident. Since the357

design of atmospheric GCM simulations is such that358

surface evaporation remains relatively constant (because359

of prescribed climatological sea surface temperatures),360

total precipitation does not change much between the361

simulations. However, the partitioning of the precipita-362

tion between the stratiform and convective type can, and363

does, change drastically.364

5 Conclusion365

The WBF process is an extremely powerful microphys-366

ical mechanism that can cause rapid transformation of367

1The simulations were run for one year after a three-month spin-up, at arelatively coarse horizontal resolution of 4 ° × 5 °.

cloud macrophysical and radiative properties. It can play 368

a tremendously important role in both climate forc- 369

ing and feedback mechanisms by amplifying the effect 370

that anthropogenic perturbations in IN have on climate. 371

It may also affect the cloud-climate feedback mecha- 372

nism sometimes referred to as the cloud optical depth 373

feedback (Zelinka, Klein and Hartmann, 2012) by 374

amplifying the effect that warming temperatures has 375

on cloud phase (McCoy, Hartmann and Grosvenor, 376

2014). As such, the significance of the WBF process 377

is gaining attention, as it becomes increasingly clear 378

that realistic representations of aerosol-cloud interac- 379

tions and cloud feedbacks in climate models rely on 380

the accuracy of the representation of the WBF process 381

in these models. In retrospect, Findeisen’s paper from 382

1938 is thus more relevant now than ever before, but 383

ironically for reasons that Findeisen could not have pre- 384

dicted when he wrote his seminal paper. Global warming 385

was not yet detectable at the time, and the early warnings 386

by Arrhenius (1896) had largely been forgotten. 387

References 388

Arrhenius, S., 1896: On the Influence of Carbonic Acid in 389

the Air upon the Temperature of the Ground. – Philosophical 390

Magazine 41, 237–276. 391

Atkinson, J.D., B.J. Murray, M.T. Woodhouse, T.F. Whale, 392

K.J. Baustian, K.S. Carslaw, S. Dobbie, D. O’Sullivan, 393

T.L. Malkin, 2013: The importance of feldspar for ice nucle- 394

ation by mineral dust in mixed-phase clouds. – Nature 498, 395

355–358. 396

Bergeron, T., 1928: Über die dreidimensional verknüpfende 397

Wetteranalyse. – Geophys. Norv. 398

Cziczo, D.J., O. Stetzer, A. Worringen, M. Ebert, S. Wein- 399

bruch, M. Kamphus, S.J. Gallavardin, J. Curtius, 400

S. Borrmann, K.D. Froyd, S. Mertes, O. Mohler, 401

U. Lohmann, 2009: Inadvertent climate modification due to 402

anthropogenic lead. – Nature Geosci. 2, 333–336. 403

Findeisen, W., 1938 [F38]: Kolloid-meteorologische 404

Vorgänge bei Niederschlagsbildung. – Meteorol. Z. 55, 405

121–133. (translated and edited by Volken, E., A.M. Gie- 406

sche, S. Brönnimann. – Meteorol. Z. 24 (2015), DOI: 407

10.1127/metz/2015/0675). 408

Hoose, C., O. Mohler, 2012: Heterogeneous ice nucleation 409

on atmospheric aerosols: a review of results from laboratory 410

experiments. – Atmos. Chem. Phys. 12, 9817–9854. 411

Page 7: The Wegener-Bergeron-Findeisen process – Its …...PrePub Article, 2015 T. Storelvmo & I. Tan: The Wegener-Bergeron-Findeisen process 3 Figure 2: Left: Saturation vapor pressure

proof

Meteorol. Z., PrePub Article, 2015 T. Storelvmo & I. Tan: The Wegener-Bergeron-Findeisen process 7

Hu, Y.X., D. Winker, M. Vaughan, B. Lin, A. Omar,412

C. Trepte, D. Flittner, P. Yang, S.L. Nasiri, B. Baum,413

W.B. Sun, Z.Y. Liu, Z. Wang, S. Young, K. Stamnes,414

J.P. Huang, R. Kuehn, R. Holz, 2009: Calipso/Caliop Cloud415

Phase Discrimination Algorithm. – J. Atmos. Oceanic Tech-416

nol. 26, 2293–2309.417

Hudson, J.G., P. Squires, 1974: Evaluation of a Vertical Contin-418

uous Cloud Chamber. – Bull. Amer. Meteor. Soc 55, 692–692.419

Kanamitsu, M., W. Ebisuzaki, J. Woollen, S.K. Yang,420

J.J. Hnilo, M. Fiorino, G.L. Potter, 2002: Ncep-Doe421

Amip-Ii Reanalysis (R-2). – Bull. Amer. Meteor. Soc. 83,422

1631–1643.423

Köhler, H., 1936: The nucleus in and the growth of hygroscopic424

droplets. –Transactions of the Faraday Society 32, 1152–1161.425

Komurcu, M., T. Storelvmo, I. Tan, U. Lohmann, Y.X. Yun,426

J.E. Penner, Y. Wang, X.H. Liu, T. Takemura, 2014: In-427

tercomparison of the cloud water phase among global climate428

models. – J. Geophys. Res.-Atmos. 119, 3372–3400.429

Korolev, A., 2007: Limitations of the wegener-bergeron-430

findeisen mechanism in the evolution of mixed-phase431

clouds. – J. Atmos. Sci. 64, 3372–3375.432

Korolev, A., G.A. Isaac, 2006: Relative humidity in liquid,433

mixed-phase, and ice clouds. – J. Atmos. Sci. 63, 2865–2880.434

Korolev, A.V., I.P. Mazin, 2003: Supersaturation of water va-435

por in clouds. – J. Atmos. Sci. 60, 2957–2974.436

Lau, K.M., H.T. Wu, 2003: Warm rain processes over tropical437

oceans and climate implications. – Geophys. Res. Lett. 30,438

published online, DOI:10.1029/2003GL018567.439

Loeb, N.G., B.A. Wielicki, D.R. Doelling, G.L. Smith,440

D.F. Keyes, S. Kato, N. Manalo-Smith, T. Wong, 2009:441

Toward Optimal Closure of the Earth’s Top-of-Atmosphere442

Radiation Budget. – J. Climate 22, 748–766.443

Lohmann, U., C. Hoose, 2009: Sensitivity studies of differ-444

ent aerosol indirect effects in mixed-phase clouds. – Atmos.445

Chem. Phys. 9, 8917–8934.446

Magnus, G., 1844: Versuche über die Spannkräfte des Wasser-447

dampfs. – Ann. Phys. Chem. 61, 225–247.448

McCoy, D.T., D.L. Hartmann, D.P. Grosvenor, 2014:449

Observed Southern Ocean Cloud Properties and Short-450

wave Reflection Part 2: Phase changes and low451

cloud feedback. – J. Climate, published online, DOI:452

10.1175/JCLI-D-14–00287.1453

Morrison, H., J.A. Curry, M.D. Shupe, P. Zuidema, 2005:454

A new double-moment microphysics parameterization for ap-455

plication in cloud and climate models. Part II: Single-column456

modeling of arctic clouds. - J. Atmos. Sci. 62, 1678–1693.457

Murray, B.J., D.O’Sullivan, J.D. Atkinson, M.E. Webb,458

2012: Ice nucleation by particles immersed in supercooled459

cloud droplets. – Chem. Soc. Rev. 41, 6519–6554.460

Podzimek, J., 1957: Die Ergebnisse der Ersten Kondensa-461

tionsversuche in der Unterdruckapparatur. – Geofisica pura e462

applicata 36, 76–85.463

Pruppacher, H.R., J.D. Klett, 2010: Microphysics of clouds 464

and precipitation. – Dordrecht; New York: Springer. 465

Roberts, G.C., A. Nenes, 2005: A continuous-flow stream- 466

wise thermal-gradient CCN chamber for atmospheric mea- 467

surements. – Aerosol Sci. Technol. 39, 206–221. 468

Rogers, D.C., 1988: Development of a continuous flow thermal 469

gradient diffusion chamber for ice nucleation studies. – At- 470

mos. Res. 22, 149–181. 471

Rotstayn, L.D., 1997: A physically based scheme for the treat- 472

ment of stratiform clouds and precipitation in large-scale mod- 473

els.1. Description and evaluation of the microphysical pro- 474

cesses. – Quart. J. Roy. Meteor. Soc. 123, 1227–1282. 475

Rotstayn, L.D., B.F. Ryan, J.J. Katzfey, 2000: A scheme 476

for calculation of the liquid fraction in mixed-phase strati- 477

form clouds in large-scale models. – Mon. Wea. Rev. 128, 478

1070–1088. 479

Smith, R.N.B., 1990: A Scheme for Predicting Layer Clouds and 480

Their Water-Content in a General-Circulation Model. – Quart. 481

J. Roy. Meteor. Soc. 116, 435–460. 482

Stetzer, O., B. Baschek, F. Luond, U. Lohmann, 2008: The 483

Zurich Ice Nucleation Chamber (ZINC: – A new instrument to 484

investigate atmospheric ice formation. – Aerosol Sci. Technol. 485

42, 64–74. 486

Storelvmo, T., J.E. Kristjansson, U. Lohmann, 2008a: 487

Aerosol influence on mixed-phase clouds in CAM-Oslo. – 488

J. Atmos. Sci. 65, 3214–3230. 489

Storelvmo, T., J.E. Kristjansson, U. Lohmann, T. Iversen, 490

A. Kirkevag, O. Seland, 2008b: Modeling of the Wegener- 491

Bergeron-Findeisen process-implications for aerosol indirect 492

effects. – Environ. Res. Lett. 3, 045001. 493

Storelvmo, T., J.E. Kristjansson, U. Lohmann, T. Iversen, 494

A. Kirkevag, O. Seland, 2010: Modeling of the Wegener- 495

Bergeron-Findeisen process-implications for aerosol indirect 496

effects. – Environ. Res. Lett. 5, 019801. 497

Tan, I., T. Storelvmo, Y.-S. Choi, 2014a: A comparison of the 498

ice nucleating efficiencies of clean dust, polluted dust, and 499

smoke aerosols in mixed-phase clouds based on spaceborne 500

lidar observations. – J. Geophys. Res. 119, 6653–6665. 501

Tan, I., T. Storelvmo, Y.-S. Choi, 2014b: Spaceborne lidar 502

observations of the ice-nucleating potential of dust, polluted 503

dust, and smoke aerosols in mixed-phase clouds. – J. Geophys. 504

Res. 119, 6653–6665. 505

Wegener, A. 1911. Thermodynamik der Atmosphäre. – Leipzig, 506

Germany: Barth. 507

Winker, D.M., M.A. Vaughan, A. Omar, Y.X. Hu, K.A. Pow- 508

ell, Z.Y. Liu, W.H. Hunt, S.A. Young, 2009: Overview of 509

the CALIPSO Mission and CALIOP Data Processing Algo- 510

rithms. – J. Atmos. Oceanic Technol. 26, 2310–2323. 511

Zelinka, M.D., S.A. Klein, D.L. Hartmann, 2012: Comput- 512

ing and Partitioning Cloud Feedbacks Using Cloud Property 513

Histograms. Part II: Attribution to Changes in Cloud Amount, 514

Altitude, and Optical Depth. – J. Climate 25, 3736–3754. 515