Unforced Errors Unforced Errors My mother taught me that in polite - - PowerPoint PPT Presentation

unforced errors unforced errors
SMART_READER_LITE
LIVE PREVIEW

Unforced Errors Unforced Errors My mother taught me that in polite - - PowerPoint PPT Presentation

Unforced Errors Unforced Errors My mother taught me that in polite society, we do not talk about: Unforced Errors My mother taught me that in polite society, we do not talk about: politics, Unforced Errors My mother taught me that in


slide-1
SLIDE 1

Unforced Errors

slide-2
SLIDE 2

Unforced Errors

My mother taught me that in polite society, we do not talk about:

slide-3
SLIDE 3

Unforced Errors

  • politics,

My mother taught me that in polite society, we do not talk about:

slide-4
SLIDE 4

Unforced Errors

  • politics,
  • religion,

My mother taught me that in polite society, we do not talk about:

slide-5
SLIDE 5

Unforced Errors

  • politics,
  • religion,
  • operating systems, or

My mother taught me that in polite society, we do not talk about:

slide-6
SLIDE 6

Unforced Errors

  • politics,
  • religion,
  • operating systems, or
  • cumulus parameterizations.

My mother taught me that in polite society, we do not talk about:

slide-7
SLIDE 7

Unforced Errors

Going against her advice, today I’m going to talk about quasi-equilibrium,
 which is always a good way to start a lively discussion.

  • politics,
  • religion,
  • operating systems, or
  • cumulus parameterizations.

My mother taught me that in polite society, we do not talk about:

slide-8
SLIDE 8

Acknowledgments

Mark Branson Don Dazlich

slide-9
SLIDE 9
slide-10
SLIDE 10

Sources and sinks of buoyancy

∂ ∂t A λ

( )=

K λ, ′ λ

( )M B

′ λ

( )

λmax

d ′ λ + F λ

( )

slide-11
SLIDE 11

Sources and sinks of buoyancy

∂ ∂t A λ

( )=

K λ, ′ λ

( )M B

′ λ

( )

λmax

d ′ λ + F λ

( )

Forcing

slide-12
SLIDE 12

Sources and sinks of buoyancy

∂ ∂t A λ

( )=

K λ, ′ λ

( )M B

′ λ

( )

λmax

d ′ λ + F λ

( )

Response Forcing

slide-13
SLIDE 13

Forcing Response

slide-14
SLIDE 14

ASQE

K λ, ′ λ

( )M B

′ λ

( )

λmax

d ′ λ + F λ

( )≅ 0

slide-15
SLIDE 15

ASQE

K λ, ′ λ

( )M B

′ λ

( )

λmax

d ′ λ + F λ

( )≅ 0

R + F ≅ 0

slide-16
SLIDE 16

The “forcing and response” paradigm

“These prognostic equations involve terms of two types: ‘Cloud terms,’ which depend on the mass flux distribution function…; and ‘large-scale terms,’ such as large-scale advection, surface eddy fluxes, and radiational heating terms, which do not depend on the mass flux distribution function….We call the large-scale terms the large-scale forcing.” —AS74

R + F ≅ 0

slide-17
SLIDE 17

The “forcing and response” paradigm

“These prognostic equations involve terms of two types: ‘Cloud terms,’ which depend on the mass flux distribution function…; and ‘large-scale terms,’ such as large-scale advection, surface eddy fluxes, and radiational heating terms, which do not depend on the mass flux distribution function….We call the large-scale terms the large-scale forcing.” —AS74 “The large-scale forcing can be divided into two parts: … the ‘cloud layer forcing’ and the ‘mixed layer forcing.’” —AS74

R + F ≅ 0

slide-18
SLIDE 18

The mixed-layer forcing

The mixed-layer forcing exerts a powerful influence on the CAPE,
 because what happens in the mixed layer affects an updraft’s buoyancy at all levels.

slide-19
SLIDE 19

The mixed-layer forcing

Surface evaporation Surface sensible heat flux The mixed-layer forcing exerts a powerful influence on the CAPE,
 because what happens in the mixed layer affects an updraft’s buoyancy at all levels.

slide-20
SLIDE 20

The mixed-layer forcing

Surface evaporation Surface sensible heat flux Entrainment into the mixed layer The mixed-layer forcing exerts a powerful influence on the CAPE,
 because what happens in the mixed layer affects an updraft’s buoyancy at all levels.

slide-21
SLIDE 21

The mixed-layer forcing

Surface evaporation Surface sensible heat flux Entrainment into the mixed layer A key missing ingredient:
 Downdrafts The mixed-layer forcing exerts a powerful influence on the CAPE,
 because what happens in the mixed layer affects an updraft’s buoyancy at all levels.

slide-22
SLIDE 22

Next, a few words about Wayne’s 1973 dissertation

slide-23
SLIDE 23

Late June, 1972

Next, a few words about Wayne’s 1973 dissertation

slide-24
SLIDE 24
  • Fig. 38 of Wayne’s dissertation
slide-25
SLIDE 25
  • Fig. 38 of Wayne’s dissertation

“So, Wayne,” I said…

slide-26
SLIDE 26
  • Fig. 38 of Wayne’s dissertation

“So, Wayne,” I said… “Why is the mixed-layer forcing so small?”

slide-27
SLIDE 27

Marshall Islands Data

“Well,” said Wayne…

slide-28
SLIDE 28

AS74 did not include downdrafts, so the moistening by surface evaporation has to be balanced by some combination of horizontal advection and the entrainment of dry air across the top of the mixed layer.

“Well,” said Wayne…

slide-29
SLIDE 29

AS74 did not include downdrafts, so the moistening by surface evaporation has to be balanced by some combination of horizontal advection and the entrainment of dry air across the top of the mixed layer. The vertical resolution of the Marshall Islands data is completely inadequate to reveal the (presumably small) water vapor mixing ratio of the entrained air.

“Well,” said Wayne…

slide-30
SLIDE 30

AS74 did not include downdrafts, so the moistening by surface evaporation has to be balanced by some combination of horizontal advection and the entrainment of dry air across the top of the mixed layer. The vertical resolution of the Marshall Islands data is completely inadequate to reveal the (presumably small) water vapor mixing ratio of the entrained air. An assumption had to be made. Wayne’s thesis doesn’t say what was assumed, and Wayne doesn’t remember.

“Well,” said Wayne…

slide-31
SLIDE 31

AS74 did not include downdrafts, so the moistening by surface evaporation has to be balanced by some combination of horizontal advection and the entrainment of dry air across the top of the mixed layer. The vertical resolution of the Marshall Islands data is completely inadequate to reveal the (presumably small) water vapor mixing ratio of the entrained air. An assumption had to be made. Wayne’s thesis doesn’t say what was assumed, and Wayne doesn’t remember.

“Well,” said Wayne…

He may have assumed that the entrained air was dry enough to balance the surface evaporation.

slide-32
SLIDE 32

“The large-scale forcing can be divided into two parts: … the ‘cloud layer forcing’ and the ‘mixed layer forcing.’” —AS74

slide-33
SLIDE 33
slide-34
SLIDE 34

Our story so far…

ASQE is based on a forcing-and-response paradigm in which the convection responds to “large-scale forcing.” AS74 distinguished between the mixed-layer forcing and the cloud-layer forcing. Simple physical reasoning suggests that the mixed-layer forcing should be strong. Wayne’s thesis includes a figure showing that the mixed-layer forcing is weak. 
 This may have been based on an assumption that the surface evaporation is mostly cancelled by entrainment drying. The tests of QE reported by AS74 are based on the cloud-layer forcing alone. The semi-prognostic tests of Lord (1982) are also based on the cloud-layer forcing alone.

slide-35
SLIDE 35

According to one school of thought, the mixed-layer forcing is dominant.

slide-36
SLIDE 36

Mixed-layer forcing is key. Cumulus downdrafts balance it.

slide-37
SLIDE 37

BLQE

“…convection is regulated by a balance between the respective tendencies of surface fluxes and convective downdrafts to increase and decrease boundary-layer equivalent potential temperature.”

slide-38
SLIDE 38

BLQE

The BLQE hypothesis asserts that the mixed-layer forcing is the primary driver for deep convection.
 The physical argument is that the powerful mixed-layer forcing leads to cumulus downdrafts that cancel it out. 
 This is a hypothetical but explicit and plausible negative feedback of deep convection that regulates the mixed-layer’s properties.

“…convection is regulated by a balance between the respective tendencies of surface fluxes and convective downdrafts to increase and decrease boundary-layer equivalent potential temperature.”

slide-39
SLIDE 39

Inferences from Simple Models of Slow, Convectively Coupled Processes

KERRY EMANUEL

Lorenz Center, Massachusetts Institute of Technology, Cambridge, Massachusetts (Manuscript received 20 March 2018, in final form 28 October 2018) JANUARY 2019 E M A N U E L

195

BLQE is being used in simple models, but as far as I know it’s not being used in any GCM.

Boundary layer quasi equilibrium may be thought

  • f as the limit of (2) as the depth d of the boundary

layer becomes vanishingly small. In that case, (2) may be approximated, after substituting (1) for the sum Md 1 we, as Mu 5 w 1 Fh hb 2 hm . (3) This is our simple way of dealing with deep moist con-

  • vection. While relatively crude, it has been used with

some success in a forecast model of tropical cyclones

slide-40
SLIDE 40

If downdrafts are included in the cumulus parameterization, BLQE can be viewed as a limiting case of ASQE.

slide-41
SLIDE 41

GeophysicalResearchLetters

RESEARCH LETTER

10.1002/2014GL062649

Key Points:

  • The tropical boundary layer is

dried more by entrainment than by downdrafts

  • Downdrafts sometimes inject

high-energy air into the boundary layer

  • Models need better parameteriza-

tions of entrainment at the boundary layer top Correspondence to:

  • K. Thayer-Calder,

katec@ucar.edu

A numerical investigation of boundary layer quasi-equilibrium

  • K. Thayer-Calder1,2 and David Randall3

1Department of Applied Mathematics, University of Wisconsin-Milwaukee, Milwaukee, Wisconsin, USA, 2National Center

for Atmospheric Research, Boulder, Colorado, USA, 3Department of Atmospheric Sciences, Colorado State University, Fort Collins, Colorado, USA

Abstract Despite the large energy input from surface evaporation, the moist static energy (MSE) of the

tropical boundary layer remains relatively constant on large spatial and temporal scales due to lifting of vapor by cloudy updrafts and the addition of dry air from the layers above. Arakawa and Schubert (1974) suggested that drying is due mainly to clear-air turbulent entrainment between cloudy updrafts, while Raymond (1995) described drying due mainly to convective downdrafts. We used cloud-resolving numerical simulations to investigate the transport of MSE into the boundary layer and found turbulent entrainment between clouds to be the dominant process.

Additional support for this conclusion comes from Torri & Kuang (2016) and deSzoeke et al. (2017).

A)

Drying of the mixed layer by: updrafts downdrafts entrainment

slide-42
SLIDE 42

Relevance to BLQE

The results of Thayer-Calder & Randall imply that downdrafts are not the primary regulator of boundary-layer entropy.

slide-43
SLIDE 43

Relevance to BLQE

The results of Thayer-Calder & Randall imply that downdrafts are not the primary regulator of boundary-layer entropy. It might be possible to rescue BLQE by modifying it to allow both downdrafts and entrainment to regulate the boundary-layer entropy.

slide-44
SLIDE 44

Relevance to BLQE

The results of Thayer-Calder & Randall imply that downdrafts are not the primary regulator of boundary-layer entropy. It might be possible to rescue BLQE by modifying it to allow both downdrafts and entrainment to regulate the boundary-layer entropy. A problem with that idea is that cumulus convection does not strongly influence the rate of entrainment at the mixed-layer top (at least, as far as I can see), so it’s not clear how a negative feedback loop would work.

slide-45
SLIDE 45

The plot thickens…

slide-46
SLIDE 46

Free Tropospheric Quasi-Equilibrium
 (FTQE)

My first reaction:

Guang Zhang argues that the convective response is whatever is needed to cancel destabilization by the cloud-layer forcing. He says that the mixed-layer forcing should be ignored. FTQE is the antithesis of BLQE. It is, however, consistent with Wayne’s thesis.

Convective quasi-equilibrium in midlatitude continental environment and its effect on convective parameterization

Guang J. Zhang

Center for Atmospheric Sciences, Scripps Institution of Oceanography, La Jolla, California, USA Received 28 June 2001; revised 10 December 2001; accepted 17 December 2001; published 31 July 2002.

JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 107, NO. D14, 10.1029/2001JD001005, 2002

slide-47
SLIDE 47

But there is now a body of work…

Zhang, G. J., 2002: Convective quasi-equilibrium in midlatitude continental environment and its effect on convective parameterization. J. Geophys. Res.: Atmospheres, 107, 4220. Zhang, G. J., 2003a: Convective quasi-equilibrium in the tropical western Pacific: Comparison with midlatitude continental environment. J. Geophys. Res.: Atmospheres, 108, 4592. Zhang, G.J., 2003b: Roles of tropospheric and boundary layer forcing in the diurnal cycle of convection in the US southern Great Plains. Geophys. Res. Lett., 30. Zhang, G. J. and M. Mu., 2005a: Effects of modifications to the Zhang-McFarlane convection parameterization on the simulation of the tropical precipitation in the National Center for Atmospheric Research Community Climate Model, version 3. J. Geophys. Res: Atmospheres, 110. Zhang, G.J. and Mu, M., 2005b: Simulation of the Madden–Julian oscillation in the NCAR CCM3 using a revised Zhang–McFarlane convection parameterization

  • scheme. J. Climate, 18, 4046-4064.

Zhang, G. J. and H. Wang, 2006: Toward mitigating the double ITCZ problem in NCAR CCSM3. Geophysical Research Letters, 33. Bechtold, P., N. Semane, P. Lopez, J.-P. Chaboureau, A. Beljaars, and N. Bormann, 2014: Representing equilibrium and nonequilibrium convection in large-scale

  • models. J. Atmos. Sci., 71, 734-753.

Song, F. and Zhang, G.J., 2016. Effects of southeastern Pacific sea surface temperature on the double-ITCZ bias in NCAR CESM1. Journal of Climate, 29, 7417-7433. Song, X. and G. J. Zhang, 2018: The roles of convection parameterization in the formation of double ITCZ syndrome in the NCAR CESM: I. Atmospheric

  • processes. J. Adv. Modeling Earth Syst., 10, 842-866.

Song, X. and G. J. Zhang, 2009: Convection parameterization, tropical Pacific double ITCZ, and upper-ocean biases in the NCAR CCSM3. Part I: Climatology and atmospheric feedback. J.Climate, 22, 4299-4315. Song, X. and G. J. Zhang, 2018: The roles of convection parameterization in the formation of double ITCZ syndrome in the NCAR CESM: I. Atmospheric

  • processes. Journal of Advances in Modeling Earth Systems, 10, pp.842-866.
slide-48
SLIDE 48

RESEARCH ARTICLE

10.1002/2017MS001191

The Roles of Convection Parameterization in the Formation of Double ITCZ Syndrome in the NCAR CESM: I. Atmospheric Processes

Xiaoliang Song1 and Guang J. Zhang1

1Scripps Institution of Oceanography, University of California, San Diego, La Jolla, CA, USA

Abstract Several improvements are implemented in the Zhang-McFarlane (ZM) convection scheme to

investigate the roles of convection parameterization in the formation of double intertropical convergence zone (ITCZ) bias in the NCAR CESM1.2.1. It is shown that the prominent double ITCZ biases of precipitation, sea surface temperature (SST), and wind stress in the standard CESM1.2.1 are largely eliminated in all sea- sons with the use of these improvements in convection scheme. This study for the first time demonstrates that the modifications of convection scheme can eliminate the double ITCZ biases in all seasons, including boreal winter and spring. Further analysis shows that the elimination of the double ITCZ bias is achieved not by improving other possible contributors, such as stratus cloud bias off the west coast of South America and cloud/radiation biases over the Southern Ocean, but by modifying the convection scheme itself. This study demonstrates that convection scheme is the primary contributor to the double ITCZ bias in the CESM1.2.1, and provides a possible solution to the long-standing double ITCZ problem. The atmospheric model simulations forced by observed SST show that the original ZM convection scheme tends to produce double ITCZ bias in high SST scenario, while the modified convection scheme does not. The impact of changes in each core component of convection scheme on the double ITCZ bias in atmospheric model is identified and further investigated.

  • 1. Introduction

Key Points:

Improvements in convection scheme

largely eliminate the double intertropical convergence zone (ITCZ) bias in all seasons in Community Earth System Model

Analyses demonstrate that

convection scheme is the primary contributor to the double ITCZ syndrome

Impact of each modification to

convection scheme on ITCZ simulation in the atmospheric model is identified and investigated Supporting Information:

Supporting Information S1

Correspondence to:

  • X. Song,

xisong@ucsd.edu Citation: Song, X., & Zhang, G. J. (2018). The roles of convection parameterization in the formation of double ITCZ syndrome in the NCAR CESM: I. Atmospheric processes. Journal of Advances in Modeling Earth Systems, 10, 842–866. https://doi.org/10.1002/ 2017MS001191

Journal of Advances in Modeling Earth Systems

slide-49
SLIDE 49

From Song & Zhang 2018

Results from CESM 1.2.1 (CAM 5.3)

slide-50
SLIDE 50

We need to understand this.

slide-51
SLIDE 51

We need to understand this.

Guang Zhang’s tests show that FTQE works well in the CAM.

slide-52
SLIDE 52

We need to understand this.

Guang Zhang’s tests show that FTQE works well in the CAM. I am not aware of any physical argument that purports to explain why FTQE should work.

slide-53
SLIDE 53

We need to understand this.

Guang Zhang’s tests show that FTQE works well in the CAM. I am not aware of any physical argument that purports to explain why FTQE should work. It just works.

slide-54
SLIDE 54

We need to understand this.

Guang Zhang’s tests show that FTQE works well in the CAM. I am not aware of any physical argument that purports to explain why FTQE should work. It just works.

slide-55
SLIDE 55

Why does it work?

Possible reasons:

slide-56
SLIDE 56

Why does it work?

Possible reasons:

Maybe the mixed-layer forcing really is small. 
 It appears that some magic would be needed for this to be true.

slide-57
SLIDE 57

Why does it work?

Possible reasons:

Maybe the mixed-layer forcing really is small. 
 It appears that some magic would be needed for this to be true. If BLQE and ASQE were both true, FTQE would follow by subtraction. 
 But that’s impossible.

slide-58
SLIDE 58

Why does it work?

Possible reasons:

Maybe the mixed-layer forcing really is small. 
 It appears that some magic would be needed for this to be true. If BLQE and ASQE were both true, FTQE would follow by subtraction. 
 But that’s impossible. Maybe the mixed-layer forcing does not exist when deep convection is intense, because under those conditions the surface fluxes do not converge inside the boundary layer.

slide-59
SLIDE 59

DYNAMO simulation with SAM

Moistening due to large-scale vertical motion g kg-1 day-1 Moistening due to horizontal advection g kg-1 day-1

270 280 290 300 310 320 330 340 350 360 370 time, day of year 10 20 30 40 50 60 70 80 90 mm/day DYNAMO-1600m64L Precipitation rate

Precipitation rate mm day-1 height, km height, km

slide-60
SLIDE 60

DYNAMO simulation results binned by precipitation rate

DYNAMO-1600m64L Qv,Qc tendency due to resolved eddy flux plus SGS, g kg-1 day-1 5 10 15 20 25 30 35 40 45 precip rate, mm/day 2 4 6 8 10 12 height, km
  • 10
  • 8
  • 6
  • 4
  • 2
2 4 6 8 10 DYNAMO-1600m64L Qv,Qc tendency due to SGS flux, g kg-1 day-1 5 10 15 20 25 30 35 40 45 precip rate, mm/day 2 4 6 8 10 12 height, km
  • 10
  • 8
  • 6
  • 4
  • 2
2 4 6 8 10

Moistening due to SGS fluxes g kg-1 day-1 Moistening due to resolved fluxes g kg-1 day-1 Moistening due to resolved and SGS fluxes g kg-1 day-1 height, km

DYNAMO-1600m64L Qv+Qc tendency due to flux divergence, net, g kg -1 day-1 5 10 15 20 25 30 35 40 45 precip rate, mm/day 2 4 6 8 10 12 height, km
  • 10
  • 8
  • 6
  • 4
  • 2
2 4 6 8 10

height, km height, km Precipitation rate mm day-1 Precipitation rate mm day-1 Precipitation rate mm day-1

slide-61
SLIDE 61

DYNAMO simulation binned by precipitation rate

DYNAMO-1600m64L Qv+Qc tendency due to flux divergence, net, g kg -1 day-1 5 10 15 20 25 30 35 40 45 precip rate, mm/day 2 4 6 8 10 12 height, km
  • 10
  • 8
  • 6
  • 4
  • 2
2 4 6 8 10 DYNAMO-1600m64L Qv+Qc tendency due to flux divergence, upward, g kg -1 day-1 5 10 15 20 25 30 35 40 45 precip rate, mm/day 2 4 6 8 10 12 height, km
  • 10
  • 8
  • 6
  • 4
  • 2
2 4 6 8 10 DYNAMO-1600m64L Qv+Qc tendency due to flux divergence, downward, g kg -1 day-1 5 10 15 20 25 30 35 40 45 precip rate, mm/day 2 4 6 8 10 12 height, km
  • 10
  • 8
  • 6
  • 4
  • 2
2 4 6 8 10

Moistening due to all vertical fluxes g kg-1 day-1 Moistening due to updrafts g kg-1 day-1 Moistening due to downdrafts g kg-1 day-1 height, km height, km height, km Precipitation rate mm day-1 Precipitation rate mm day-1 Precipitation rate mm day-1

slide-62
SLIDE 62 DYNAMO-1600m64L Qv,Qc tendency due to microphyics, g kg -1 day-1 5 10 15 20 25 30 35 40 45 precip rate, mm/day 2 4 6 8 10 12 height, km
  • 10
  • 8
  • 6
  • 4
  • 2
2 4 6 8 10

DYNAMO simulation binned by precipitation rate

DYNAMO-1600m64L Qv+Qc tendency due to flux divergence, net, g kg -1 day-1 5 10 15 20 25 30 35 40 45 precip rate, mm/day 2 4 6 8 10 12 height, km
  • 10
  • 8
  • 6
  • 4
  • 2
2 4 6 8 10

Moistening due to all vertical fluxes g kg-1 day-1 Moistening due to microphysics g kg-1 day-1

DYNAMO-1600m64L Qv,Qc total tendency minus LS, g kg -1 day-1 5 10 15 20 25 30 35 40 45 precip rate, mm/day 2 4 6 8 10 12 height, km
  • 10
  • 8
  • 6
  • 4
  • 2
2 4 6 8 10

Moistening due to fluxes and microphysics g kg-1 day-1 height, km height, km height, km Precipitation rate mm day-1 Precipitation rate mm day-1 Precipitation rate mm day-1

slide-63
SLIDE 63 DYNAMO-1600m64L Qv,Qc tendency due to microphyics, g kg -1 day-1 5 10 15 20 25 30 35 40 45 precip rate, mm/day 2 4 6 8 10 12 height, km
  • 10
  • 8
  • 6
  • 4
  • 2
2 4 6 8 10

DYNAMO simulation binned by precipitation rate

DYNAMO-1600m64L Qv+Qc tendency due to flux divergence, net, g kg -1 day-1 5 10 15 20 25 30 35 40 45 precip rate, mm/day 2 4 6 8 10 12 height, km
  • 10
  • 8
  • 6
  • 4
  • 2
2 4 6 8 10

Moistening due to all vertical fluxes g kg-1 day-1 Moistening due to microphysics g kg-1 day-1

DYNAMO-1600m64L Qv,Qc total tendency minus LS, g kg -1 day-1 5 10 15 20 25 30 35 40 45 precip rate, mm/day 2 4 6 8 10 12 height, km
  • 10
  • 8
  • 6
  • 4
  • 2
2 4 6 8 10

Moistening due to fluxes and microphysics g kg-1 day-1 height, km height, km height, km Precipitation rate mm day-1 Precipitation rate mm day-1 Precipitation rate mm day-1

Mixed-layer forcing?

slide-64
SLIDE 64 DYNAMO-1600m64L Qv,Qc tendency due to microphyics, g kg -1 day-1 5 10 15 20 25 30 35 40 45 precip rate, mm/day 2 4 6 8 10 12 height, km
  • 10
  • 8
  • 6
  • 4
  • 2
2 4 6 8 10

DYNAMO simulation binned by precipitation rate

DYNAMO-1600m64L Qv+Qc tendency due to flux divergence, net, g kg -1 day-1 5 10 15 20 25 30 35 40 45 precip rate, mm/day 2 4 6 8 10 12 height, km
  • 10
  • 8
  • 6
  • 4
  • 2
2 4 6 8 10

Moistening due to all vertical fluxes g kg-1 day-1 Moistening due to microphysics g kg-1 day-1

DYNAMO-1600m64L Qv,Qc total tendency minus LS, g kg -1 day-1 5 10 15 20 25 30 35 40 45 precip rate, mm/day 2 4 6 8 10 12 height, km
  • 10
  • 8
  • 6
  • 4
  • 2
2 4 6 8 10

Moistening due to fluxes and microphysics g kg-1 day-1 height, km height, km height, km

The mixed-layer forcing does not exist when deep convection is intense, because the moisture flux does not converge inside the mixed layer.

Precipitation rate mm day-1 Precipitation rate mm day-1 Precipitation rate mm day-1

Mixed-layer forcing?

slide-65
SLIDE 65

DYNAMO-1600m64L Qv,Qc total tendency minus LS, g kg-1 day-1

5 10 15 20 25 30 35 40 45

precip rate, mm/day

2 4 6 8 10 12

height, km

  • 3
  • 2
  • 1

1 2 3

Moistening due to fluxes and microphysics g kg-1 day-1

Same plot with a different color bar

Precipitation rate mm day-1 Mixed-layer forcing?

slide-66
SLIDE 66

A more basic issue:

slide-67
SLIDE 67

A more basic issue:

Can we really separate the forcing from the response?

slide-68
SLIDE 68

A more basic issue:

Can we really separate the forcing from the response?

Surface fluxes are influenced by deep convection.

slide-69
SLIDE 69

A more basic issue:

Can we really separate the forcing from the response?

Surface fluxes are influenced by deep convection. Stratiform precipitation is influenced by deep convection.

slide-70
SLIDE 70

A more basic issue:

Can we really separate the forcing from the response?

Surface fluxes are influenced by deep convection. Stratiform precipitation is influenced by deep convection. Radiatively active stratiform clouds are influenced by deep convection.

slide-71
SLIDE 71

A more basic issue:

Can we really separate the forcing from the response?

Surface fluxes are influenced by deep convection. Stratiform precipitation is influenced by deep convection. Radiatively active stratiform clouds are influenced by deep convection. Randall and Pan (1993, p. 143): “… it is not always clear which processes are convective and which are not.”

Randall, D. A., and D.-M. Pan, 1993: Implementation of the Arakawa-Schubert cumulus parameterization with a prognostic closure. In The Representation of Cumulus Convection in Numerical Models, a Meteorological Monograph published by the American Meteorological Society, K. Emanuel and D. Raymond, Eds.,

  • pp. 137-144.
slide-72
SLIDE 72

An example: SP vs. SPX

SPX minus SP

20 40 60 80 100 120 140 Latent Heat Flux (W/m^2) 5 10 15 20 Precipitation Rate (mm/day)

  • 50 -40 -30 -20 -10 0
10 20 30 40 50

SP FV Non-Rotating

20 40 60 80 100 120 140 Latent Heat Flux (W/m^2) 5 10 15 20 Precipitation Rate (mm/day)

10 20 30 40 50 60 70 80 90 100 110

SP-X FV Non-Rotating

20 40 60 80 100 120 140 Latent Heat Flux (W/m^2) 5 10 15 20 Precipitation Rate (mm/day)

10 20 30 40 50 60 70 80 90 100 110

In SP , the surface fluxes are computed on the GCM grid and passed to the CRM, which uses the same fluxes for all grid columns of its fine grid. In SPX, the surface fluxes are computed on the CRM’s grid. 
 Averages over the CRM’s grid are sent back to the GCM for use as diagnostics. SP SPX SP minus SPX

slide-73
SLIDE 73

Is there an alternative to forcing and response?

slide-74
SLIDE 74

Is there an alternative to forcing and response?

Sure!

slide-75
SLIDE 75

Two alternatives

Q.

  • J. R. Meteorol. SOC.

(1998), 124, pp. 94%981

A cumulus parametrization with a prognostic closure

By DZONG-MING PAN* and DAVID. A. RANDALL

Colorado State University, USA

(Received 5 September 1996; revised 5 June 1997)

SUMMARY The paper describes the introduction of a prognostic cumulus kinetic energy (CKE) as a replacement for the quasi-equilibrium closure hypothesis of Arakawa and Schubert (AS). In the original version of the AS parametrization, the cloud work function, a measure of the convective available potential energy, is assumed to be maintained at ‘small’ values through a quasi-equilibrium between the cumulus convection and the ‘large-scale forcing’. It is argued here, however, that the distinction between the convective and large-scale processes is ambiguous and subjective. It is demonstrated that the need for such a distinction can be avoided by relaxing the quasi-equilibrium assumption, through the introduction of a prognostic

CKE; referred to as prognostic closure. A dimensional parameter, a,

is introduced to relate the CKE to the square

  • f the cloud-base convective mass flux. It is shown that ‘adjustment time’ defined by AS is related to a,

so that

when the adjustment time approaches zero the prognostic closure reduces to quasi-equilibrium closure. A second dimensional parameter, TD, is used to determine the rate at which the CKE is dissipated. In the limit of small CY and

?D,

the convective mass flux is formally independent of both CY and TD if the environmental sounding is assumed

to be given, but in reality the results of a prognostic model do depend on these two parameters because they affect

the time-dependent sounding. For simplicity, a single constant value of a is used for all cloud types in tests with a general-circulation model, and this gives reasonably good results. Larger values of a lead to more frequent shallow cumulus convection and a cooler and more humid troposphere, in which stratiform condensation is more active and more large-scale precipitation can reach the surface. A longer dissipation time-scale leads to a warmer tropical troposphere. The interactions between stratiform cloudiness and convection prove to be quite important, leading to the conclusion that the convection parametrization really cannot be evaluated independently of the stratiform cloud parametrization with which it interacts. KEYWORDS: Adjustment time Convections Cumulus kinetic energy Planetary boundary-layer

1.

INTRODUCTION

Cumulus convection plays an essential role in the atmospheric general circulation. Large-scale numerical models, such as general-circulation models (GCMs), have grid cells with horizontal dimensions on the order of a hundred kilometres or more, and so cannot resolve the cumuli, which have diameters on the order of 1-10 km. Representation of the cumulus effects in terms of grid-scale variables is called cumulus parametrization. The simulated atmospheric circulation produced by a GCM is extremely sensitive to the formulation of the cumulus parametrization. Variants of the cumulus parametrization proposed by Arakawa and Schubert (I 974; hereafter AS) are being used in many GCMs today. Key elements of the AS parametrization are as follows:

0 The use of a convective mass flux to parametrize the vertical transports by the

convective updraughts. This idea was first proposed by Arakawa (1969). It has now been

almost universally adopted, especially considering that moist convective adjustment can be formulated in terms of a mass flux (see Suarez et al. 1983).

a The introduction of

a simple but explicit model o

f

a cumulus cloud as a conceptual component o

f

the cumulus parametrization. The particular cloud model used by AS to

represent the life-cycle averaged properties of each cloud type consists of an entraining plume with a constant fractional entrainment rate, and detrainment only at the cloud- top level. Although this cloud model has been criticized by many authors (e.g. Warner 1970; Raymond and Blyth 1986), Lin (1994) has recently shown, using a cloud-resolving

* Corresponding author: Department of Atmospheric Science, Colorado State University, Ft. Collins, CO 80523,

USA. 949

slide-76
SLIDE 76

Reasons to use prognostic closure

There is no need to distinguish between forcing and response. The convection has memory. Prognostic closure is simpler and computationally faster.

slide-77
SLIDE 77

Reasons to use prognostic closure

There is no need to distinguish between forcing and response. The convection has memory. Prognostic closure is simpler and computationally faster.

*

slide-78
SLIDE 78

If we avoid defining forcing and response,
 can we still talk about QE?

slide-79
SLIDE 79

If we avoid defining forcing and response,
 can we still talk about QE?

Sure!

slide-80
SLIDE 80

Kinetic energy QE

280 290 300 310 320 330 340 350 360 day of year 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 day

V.Int.EddyKE/V.Int.BuoyFlux: DYNAMO SAM 1.6km V.Int.BF>0.5W/m-2

280 290 300 310 320 330 340 350 360 day of year 1 2 3 4 5 6 J m-2 104

DYNAMO SAM1.6km - Vertically integrated eddy KE

horizontal vertical total

Vertically integrated CKE 104 J m2 W m2 days Vertically integrated buoyancy flux The ratio, which is a time scale

See Lord & Arakawa (1980, “Part II”)

270 280 290 300 310 320 330 340 350 360 370 day of year
  • 1
1 2 3 4 5 6 W m-2
  • Vert. Int. Buoyancy Flux: DYNAMO SAM 1.6km
280 290 300 310 320 330 340 350 360 day of year
slide-81
SLIDE 81

We need to understand success of FTQE. In the process we will learn something. The mixed-layer forcing is not well defined when deep convection is intense. It’s best to avoid the forcing-and-response paradigm. Prognostic closure and super- parameterization do that. Even without the forcing-and-response paradigm, quasi-equilibrium can still be discussed in terms of the cumulus kinetic energy.

Conclusions

slide-82
SLIDE 82

Thanks