DOCSLIB.ORG

Potential Vorticity in a Moist Atmosphere

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Potential Vorticity in a Moist Atmosphere 3148 JOURNAL OF THE ATMOSPHERIC SCIENCES VOLUME 58 Potential Vorticity in a Moist Atmosphere WAYNE H. SCHUBERT,SCOTT A. HAUSMAN, AND MATTHEW GARCIA Department of Atmospheric Science, Colorado State University, Fort Collins, Colorado KATSUYUKI V. O OYAMA Hurricane Research Division, Atlantic Oceanographic and Meteorological Laboratory, Miami, Florida HUNG-CHI KUO Department of Atmospheric Sciences, National Taiwan University, Taipei, Taiwan (Manuscript received 11 October 2000, in ®nal form 27 April 2001) ABSTRACT The potential vorticity principle for a nonhydrostatic, moist, precipitating atmosphere is derived. An appropriate 21 generalization of the well-known (dry) Ertel potential vorticity is found to be P 5 r (2V1= 3 u)´=ur, where r is the total density, consisting of the sum of the densities of dry air, airborne moisture (vapor and cloud condensate), RaPa/c and precipitation; u is the velocity of the dry air and airborne moisture; and ur 5 Tr(p0/p) is the virtual potential temperature, with Tr 5 p/(rRa) the virtual temperature, p the total pressure (the sum of the partial pressures of dry air and water vapor), p0 the constant reference pressure, Ra the gas constant for dry air, and cPa the speci®c heat at constant pressure for dry air. Since ur is a function of total density and total pressure only, its use as the thermodynamic variable in P leads to the annihilation of the solenoidal term, that is, =ur ´(=r 3 =p) 5 0. In the special case of an absolutely dry atmosphere, P reduces to the usual (dry) Ertel potential vorticity. For balanced ¯ows, there exists an invertibility principle that determines the balanced mass and wind ®elds from the spatial distribution of P. It is the existence of this invertibility principle that makes P such a fundamentally important dynamical variable. In other words, P (in conjunction with the boundary conditions associated with the invertibility principle) carries all the essential dynamical information about the slowly evolving balanced part of the ¯ow. 1. Introduction to the parameterization of the microphysics of the pre- cipitation process. At present, the dynamical basis for global numerical The purpose of the present paper is to extend the weather prediction and climate models is the set of qua- potential vorticity conservation principle to nonhydro- si-static primitive equations. Since the nonhydrostatic static models with Ooyama's (1990, 2001) form of moist motions typical of dry and moist atmospheric convec- dynamics and thermodynamics. We begin by reviewing tion have such small horizontal scales, and since they the exact, nonhydrostatic primitive equations for a moist cannot be accurately simulated with the quasi-static atmosphere in section 2. In section 3 we derive the primitive equations, the collective effects of these mo- generalized potential vorticity principle. Equations (20) tions must be parameterized in quasi-static primitive and (21) are our main results, the latter form being useful equation models. However, in the not-too-distant future, for physical interpretation. In section 4 we discuss one it will be possible to construct and run global numerical of the many possible invertibility principles (depending weather prediction and climate models based on the on the particular balance conditions) associated with the exact primitive equations, with much more accurate generalized potential vorticity. Section 5 contains a der- treatments of the moist thermodynamics, and with ivation of the ``equivalent potential vorticity'' and a cloud-resolving spatial discretization. In such nonhy- discussion of why this form is unacceptable from the drostatic models, cumulus cloud ®elds will be explicitly standpoint of possessing an invertibility principle. Con- simulated, eliminating the need for cumulus parame- clusions are given in section 6. terization. This pushes the frontier of empiricism back 2. Nonhydrostatic primitive equations for a moist Corresponding author address: Wayne H. Schubert, Department atmosphere of Atmospheric Science, Colorado State University, Fort Collins, CO 80523. Consider atmospheric matter to consist of dry air, air- E-mail: [email protected] borne moisture (vapor and cloud condensate), and pre- q 2001 American Meteorological Society Unauthenticated | Downloaded 09/25/21 05:25 AM UTC 1NOVEMBER 2001 SCHUBERT ET AL. 3149 1 cipitation. Let ra denote the mass density of dry air, rm total pressure p. In the determination of U, sr, Qr, Qs, 5 ry 1 rc the mass density of airborne moisture (con- and p, several other diagnostic variables are introduced. sisting of the sum of the mass densities of water vapor The additional diagnostic variables introduced in (6)± ry and airborne condensed water rc), and rr the mass (13) are the mass density of dry air ra, the thermody- density of precipitating water substance. The total mass namically possible temperatures T1 and T 2, the actual density r is given by r 5 ra 1 ry 1 rc 1 rr 5 ra 1 temperature T, the partial pressure of dry air pa, and the rm 1 rr. The ¯ux forms of the prognostic equations for partial pressure of water vapor py , the formulas for r, rm, and rr are given in (1)±(3), in which u denotes which are given below. To summarize, the diagnostic the velocity (relative to the rotating earth) of dry air and relations are airborne moisture, and u 1 U denotes the velocity (rel- r 5 r 2 r 2 r , (6) ative to the rotating earth) of the precipitating water sub- amr stance, so that U is the velocity of the precipitating water S2(ram, r 1 r r, T2) 5 s, (7) substance relative to the dry air and airborne moisture. The term Qr, on the right-hand sides of (2) and (3), is srr5 r C(T2), (8) the rate of conversion from airborne moisture to precip- S ( , , T ) , (9) itation; this term can be positive (e.g., the collection of 1 ramr 1 5 s 2 sr cloud droplets by rain) or negative (e.g., the evaporation T 5 max(T12, T ), (10) of precipitation falling through unsaturated air). p 5 r RT, (11) The total entropy density is s 5 sa 1 sm 1 sr, aaa consisting of the sum of the entropy densities of dry r 5 r , r 5 0, p 5 r RT, air, airborne moisture, and precipitation. Since the en- y mc yyy tropy ¯ux is given by s u 1 s u 1 s (u 1 U) 5 su if T 5 T12. T (absence of condensate), a m r (12) 1 srU, we can write the ¯ux form of the entropy con- ryy5 r*(T), rcm5 r 2 ryy, p 5 E(T), servation principle as (4), where Qs denotes diabatic processes such as radiation. if T 5 T21. T (saturated vapor), Next, we can write the equation of motion as (5), where p 5 p 1 p . (13) z 5 2V 1 = 3 u is the absolute vorticity; F the potential a y for the sum of the Newtonian gravitational force and the Note that, although there are four types of matter (with 21 centrifugal force; r =p the pressure gradient force, with densities ra, ry , rc, rr), there are only three prognostic p 5 pa 1 py the sum of the partial pressures of dry air mass continuity equations [(1)±(3)]. The separation of pa and water vapor py ; and F the frictional force per unit the predicted total airborne moisture density rm into the mass. The derivation of (5) is given in Ooyama (2001) vapor density ry and the cloud condensate density rc is and is reviewed in appendix B. accomplished diagnostically in the two alternatives of Consolidating our results so far, the prognostic equa- (12). In the absence of condensate, all the predicted tions for the total mass density r, the mass density of airborne moisture rm is in vapor form so that ry 5 rm airborne moisture rm, the mass density of precipitation and rc 5 0, while if the vapor is saturated, ry 5 r*y (T) content rr, the entropy density s, and the three-dimen- and rc 5 rm 2 ry . The formulas for the entropy density sional velocity vector u are functions S1 and S 2, from which the temperatures T1 and T are diagnosed, are given in appendix A. ]r 2 1 = ´(ru 1 rrU) 5 0, (1) In the context of numerical modeling, the procedure ]t for advancing from one time level to the next consists ]rm of computing new values of the prognostic variables r, 1 = ´(rmru) 52Q , (2) ]t rm, rr, s, and u from (1)±(5). The diagnostic variables required for the prognostic stage are determined by se- ]r r 1 = ´[r (u 1 U)] 5 Q , (3) quential evaluation of (6)±(13), namely, the diagnosis ]t rr of the dry air density ra from (6), the thermodynamically ]s possible (wet bulb) temperature T 2 from (7), the entropy 1 = ´(su 1 srU) 5 Qs, (4) density of precipitation from (8), the thermodynamically ]t possible temperature T1 from (9), the choice of the actual ]u 11 temperature from (10), the dry air partial pressure pa 1 z 3 u 1 = u ´ u 1F 1 =p 5 F. (5) from (11), the water vapor mass density r , the airborne ]t 122 r y condensate mass density r , and the water vapor partial The diagnostic variables appearing explicitly in (1)±(5) c pressure py from the appropriate condition in (12), and are the terminal fall velocity U, the entropy density of the total pressure from (13). Since they are not essential precipitation sr, the source terms Qr and Qs, and the to our discussion here, we have omitted the parameter- ization formulas for the terminal fall velocity U and the 1 The notation used here follows Ooyama (2001).
Load more

Recommended publications
  • Chapter 8 Atmospheric Statics and Stability
    Chapter 8 Atmospheric Statics and Stability 1. The Hydrostatic Equation • HydroSTATIC – dw/dt = 0! • Represents the balance between the upward directed pressure gradient force and downward directed gravity. ρ = const within this slab dp A=1 dz Force balance p-dp ρ p g d z upward pressure gradient force = downward force by gravity • p=F/A. A=1 m2, so upward force on bottom of slab is p, downward force on top is p-dp, so net upward force is dp. • Weight due to gravity is F=mg=ρgdz • Force balance: dp/dz = -ρg 2. Geopotential • Like potential energy. It is the work done on a parcel of air (per unit mass, to raise that parcel from the ground to a height z. • dφ ≡ gdz, so • Geopotential height – used as vertical coordinate often in synoptic meteorology. ≡ φ( 2 • Z z)/go (where go is 9.81 m/s ). • Note: Since gravity decreases with height (only slightly in troposphere), geopotential height Z will be a little less than actual height z. 3. The Hypsometric Equation and Thickness • Combining the equation for geopotential height with the ρ hydrostatic equation and the equation of state p = Rd Tv, • Integrating and assuming a mean virtual temp (so it can be a constant and pulled outside the integral), we get the hypsometric equation: • For a given mean virtual temperature, this equation allows for calculation of the thickness of the layer between 2 given pressure levels. • For two given pressure levels, the thickness is lower when the virtual temperature is lower, (ie., denser air). • Since thickness is readily calculated from radiosonde measurements, it provides an excellent forecasting tool.
    [Show full text]
  • Comparison Between Observed Convective Cloud-Base Heights and Lifting Condensation Level for Two Different Lifted Parcels
    AUGUST 2002 NOTES AND CORRESPONDENCE 885 Comparison between Observed Convective Cloud-Base Heights and Lifting Condensation Level for Two Different Lifted Parcels JEFFREY P. C RAVEN AND RYAN E. JEWELL NOAA/NWS/Storm Prediction Center, Norman, Oklahoma HAROLD E. BROOKS NOAA/National Severe Storms Laboratory, Norman, Oklahoma 6 January 2002 and 16 April 2002 ABSTRACT Approximately 400 Automated Surface Observing System (ASOS) observations of convective cloud-base heights at 2300 UTC were collected from April through August of 2001. These observations were compared with lifting condensation level (LCL) heights above ground level determined by 0000 UTC rawinsonde soundings from collocated upper-air sites. The LCL heights were calculated using both surface-based parcels (SBLCL) and mean-layer parcels (MLLCLÐusing mean temperature and dewpoint in lowest 100 hPa). The results show that the mean error for the MLLCL heights was substantially less than for SBLCL heights, with SBLCL heights consistently lower than observed cloud bases. These ®ndings suggest that the mean-layer parcel is likely more representative of the actual parcel associated with convective cloud development, which has implications for calculations of thermodynamic parameters such as convective available potential energy (CAPE) and convective inhibition. In addition, the median value of surface-based CAPE (SBCAPE) was more than 2 times that of the mean-layer CAPE (MLCAPE). Thus, caution is advised when considering surface-based thermodynamic indices, despite the assumed presence of a well-mixed afternoon boundary layer. 1. Introduction dry-adiabatic temperature pro®le (constant potential The lifting condensation level (LCL) has long been temperature in the mixed layer) and a moisture pro®le used to estimate boundary layer cloud heights (e.g., described by a constant mixing ratio.
    [Show full text]
  • MSE3 Ch14 Thunderstorms
    Chapter 14 Copyright © 2011, 2015 by Roland Stull. Meteorology for Scientists and Engineers, 3rd Ed. thunderstorms Contents Thunderstorms are among the most violent and difficult-to-predict weath- Thunderstorm Characteristics 481 er elements. Yet, thunderstorms can be Appearance 482 14 studied. They can be probed with radar and air- Clouds Associated with Thunderstorms 482 craft, and simulated in a laboratory or by computer. Cells & Evolution 484 They form in the air, and must obey the laws of fluid Thunderstorm Types & Organization 486 mechanics and thermodynamics. Basic Storms 486 Thunderstorms are also beautiful and majestic. Mesoscale Convective Systems 488 Supercell Thunderstorms 492 In thunderstorms, aesthetics and science merge, making them fascinating to study and chase. Thunderstorm Formation 496 Convective Conditions 496 Thunderstorm characteristics, formation, and Key Altitudes 496 forecasting are covered in this chapter. The next chapter covers thunderstorm hazards including High Humidity in the ABL 499 hail, gust fronts, lightning, and tornadoes. Instability, CAPE & Updrafts 503 CAPE 503 Updraft Velocity 508 Wind Shear in the Environment 509 Hodograph Basics 510 thunderstorm CharaCteristiCs Using Hodographs 514 Shear Across a Single Layer 514 Thunderstorms are convective clouds Mean Wind Shear Vector 514 with large vertical extent, often with tops near the Total Shear Magnitude 515 tropopause and bases near the top of the boundary Mean Environmental Wind (Normal Storm Mo- layer. Their official name is cumulonimbus (see tion) 516 the Clouds Chapter), for which the abbreviation is Supercell Storm Motion 518 Bulk Richardson Number 521 Cb. On weather maps the symbol represents thunderstorms, with a dot •, asterisk , or triangle Triggering vs. Convective Inhibition 522 * ∆ drawn just above the top of the symbol to indicate Convective Inhibition (CIN) 523 Trigger Mechanisms 525 rain, snow, or hail, respectively.
    [Show full text]
  • Thunderstorm Predictors and Their Forecast Skill for the Netherlands
    Atmospheric Research 67–68 (2003) 273–299 www.elsevier.com/locate/atmos Thunderstorm predictors and their forecast skill for the Netherlands Alwin J. Haklander, Aarnout Van Delden* Institute for Marine and Atmospheric Sciences, Utrecht University, Princetonplein 5, 3584 CC Utrecht, The Netherlands Accepted 28 March 2003 Abstract Thirty-two different thunderstorm predictors, derived from rawinsonde observations, have been evaluated specifically for the Netherlands. For each of the 32 thunderstorm predictors, forecast skill as a function of the chosen threshold was determined, based on at least 10280 six-hourly rawinsonde observations at De Bilt. Thunderstorm activity was monitored by the Arrival Time Difference (ATD) lightning detection and location system from the UK Met Office. Confidence was gained in the ATD data by comparing them with hourly surface observations (thunder heard) for 4015 six-hour time intervals and six different detection radii around De Bilt. As an aside, we found that a detection radius of 20 km (the distance up to which thunder can usually be heard) yielded an optimum in the correlation between the observation and the detection of lightning activity. The dichotomous predictand was chosen to be any detected lightning activity within 100 km from De Bilt during the 6 h following a rawinsonde observation. According to the comparison of ATD data with present weather data, 95.5% of the observed thunderstorms at De Bilt were also detected within 100 km. By using verification parameters such as the True Skill Statistic (TSS) and the Heidke Skill Score (Heidke), optimal thresholds and relative forecast skill for all thunderstorm predictors have been evaluated.
    [Show full text]
  • 1 Module 4 Water Vapour in the Atmosphere 4.1 Statement of The
    Module 4 Water Vapour in the Atmosphere 4.1 Statement of the General Meteorological Problem D. Brunt (1941) in his book Physical and Dynamical Meteorology has stated, “The main problem to be discussed in connection to the thermodynamics of the moist air is the variation of temperature produced by changes of pressure, which in the atmosphere are associated with vertical motion. When damp air ascends, it must eventually attain saturation, and further ascent produces condensation, at first in the form of water drops, and as snow in the later stages”. This statement of the problem emphasizes the role of vertical ascends in producing condensation of water vapour. However, several text books and papers discuss this problem on the assumption that products of condensation are carried with the ascending air current and the process is strictly reversible; meaning that if the damp air and water drops or snow are again brought downwards, the evaporation of water drops or snow uses up the same amount of latent heat as it was liberated by condensation on the upward path of the air. Another assumption is that the drops fall out as the damp air ascends but then the process is not reversible, and Von Bezold (1883) termed it as a pseudo-adiabatic process. It must be pointed out that if the products of condensation are retained in the ascending current, the mathematical treatment is easier in comparison to the pseudo-adiabatic case. There are four stages that can be discussed in connection to the ascent of moist air. (a) The air is saturated; (b) The air is saturated and contains water drops at a temperature above the freezing-point; (c) All the water drops freeze into ice at 0°C; (d) Saturated air and ice at temperatures below 0°C.
    [Show full text]
  • Convective Cores in Continental and Oceanic Thunderstorms: Strength, Width, And
    Convective Cores in Continental and Oceanic Thunderstorms: Strength, Width, and Dynamics THESIS Presented in Partial Fulfillment of the Requirements for the Degree Master of Science in the Graduate School of The Ohio State University By Alexander Michael McCarthy, B.S. Graduate Program in Atmospheric Sciences The Ohio State University 2017 Master's Examination Committee: Jialin Lin, Advisor Jay Hobgood Bryan Mark Copyrighted by Alexander Michael McCarthy 2017 Abstract While aircraft penetration data of oceanic thunderstorms has been collected in several field campaigns between the 1970’s and 1990’s, the continental data they were compared to all stem from penetrations collected from the Thunderstorm Project in 1947. An analysis of an updated dataset where modern instruments allowed for more accurate measurements was used to make comparisons to prior oceanic field campaigns that used similar instrumentation and methodology. Comparisons of diameter and vertical velocity were found to be similar to previous findings. Cloud liquid water content magnitudes were found to be significantly greater over the oceans than over the continents, though the vertical profile of oceanic liquid water content showed a much more marked decrease with height above 4000 m than the continental profile, lending evidence that entrainment has a greater impact on oceanic convection than continental convection. The difference in buoyancy between vigorous continental and oceanic convection was investigated using a reversible CAPE calculation. It was found that for the strongest thunderstorms over the continents and oceans, oceanic CAPE values tended to be significantly higher than their continental counterparts. Alternatively, an irreversible CAPE calculation was used to investigate the role of entrainment in reducing buoyancy for continental and oceanic convection where it was found that entrainment played a greater role in diluting oceanic buoyancy than continental buoyancy.
    [Show full text]
  • Stability Analysis, Page 1 Synoptic Meteorology I
    Synoptic Meteorology I: Stability Analysis For Further Reading Most information contained within these lecture notes is drawn from Chapters 4 and 5 of “The Use of the Skew T, Log P Diagram in Analysis and Forecasting” by the Air Force Weather Agency, a PDF copy of which is available from the course website. Chapter 5 of Weather Analysis by D. Djurić provides further details about how stability may be assessed utilizing skew-T/ln-p diagrams. Why Do We Care About Stability? Simply put, we care about stability because it exerts a strong control on vertical motion – namely, ascent – and thus cloud and precipitation formation on the synoptic-scale or otherwise. We care about stability because rarely is the atmosphere ever absolutely stable or absolutely unstable, and thus we need to understand under what conditions the atmosphere is stable or unstable. We care about stability because the purpose of instability is to restore stability; atmospheric processes such as latent heat release act to consume the energy provided in the presence of instability. Before we can assess stability, however, we must first introduce a few additional concepts that we will later find beneficial, particularly when evaluating stability using skew-T diagrams. Stability-Related Concepts Convection Condensation Level The convection condensation level, or CCL, is the height or isobaric level to which an air parcel, if sufficiently heated from below, will rise adiabatically until it becomes saturated. The attribute “if sufficiently heated from below” gives rise to the convection portion of the CCL. Sufficiently strong heating of the Earth’s surface results in dry convection, which generates localized thermals that act to vertically transport energy.
    [Show full text]
  • The Static Atmsophere Equation of State, Hydrostatic Balance, and the Hypsometric Equation
    The Static Atmsophere Equation of State, Hydrostatic Balance, and the Hypsometric Equation Equation of State The atmosphere at any point can be described by its temperature, pressure, and density. These variables are related to each other through the Equation of State for dry air: pα = Rd T v (1) Where α is the specific volume (1/ρ), p = pressure, Tv = virtual temperature (K), and R is the gas constant for dry air (R = 287 J kg-1 K-1). The virtual temperature is the temperature air would have if we account for the moisture content of it. Above the surface, where moisture is scarce, Tv and T are nearly equal. Near the surface you can use the following equation to calculate Tv: TTv = (1+ 0.6078q) where q = specific humidity (T must be in K, q must be dimensionless). More on virtual temperature a little later….. We can also write the equation of state as: p= ρ Rd T v (2) The equation of state is also known as the “Ideal Gas Law”. The atmosphere is treated as an “Ideal Gas”, meaning that the molecules are considered to be of negligible size so that they exert no intermolecular forces. Hydrostatic Equation Except in smaller-scale systems such as an intense squall line, the vertical pressure gradient in the atmosphere is balanced by the gravity force. We call this the hydrostatic balance. dp = −ρg (3) dz Integrating (3) from a height z to the top of the atmosphere yields: ∞ p() z= ∫ ρ gdz (4) z This says that the pressure at any point in the atmosphere is equal to the weight of the air above it.
    [Show full text]
  • M.Sc. in Meteorology Physical Meteorology Part 2 Atmospheric
    M.Sc. in Meteorology Physical Meteorology Part 2 Prof Peter Lynch Atmospheric Thermodynamics Mathematical Computation Laboratory Dept. of Maths. Physics, UCD, Belfield. 2 Atmospheric Thermodynamics Outline of Material Thermodynamics plays an important role in our quanti- • 1 The Gas Laws tative understanding of atmospheric phenomena, ranging from the smallest cloud microphysical processes to the gen- • 2 The Hydrostatic Equation eral circulation of the atmosphere. • 3 The First Law of Thermodynamics The purpose of this section of the course is to introduce • 4 Adiabatic Processes some fundamental ideas and relationships in thermodynam- • 5 Water Vapour in Air ics and to apply them to a number of simple, but important, atmospheric situations. • 6 Static Stability The course is based closely on the text of Wallace & Hobbs • 7 The Second Law of Thermodynamics 3 4 We resort therefore to a statistical approach, and consider The Kinetic Theory of Gases the average behaviour of the gas. This is the approach called The atmosphere is a gaseous envelope surrounding the Earth. the kinetic theory of gases. The laws governing the bulk The basic source of its motion is incoming solar radiation, behaviour are at the heart of thermodynamics. We will which drives the general circulation. not consider the kinetic theory explicitly, but will take the thermodynamic principles as our starting point. To begin to understand atmospheric dynamics, we must first understand the way in which a gas behaves, especially when heat is added are removed. Thus, we begin by studying thermodynamics and its application in simple atmospheric contexts. Fundamentally, a gas is an agglomeration of molecules.
    [Show full text]
  • An R Package for Atmospheric Thermodynamics
    An R package for atmospheric thermodynamics Jon S´aenza,b,∗, Santos J. Gonz´alez-Roj´ıa, Sheila Carreno-Madinabeitiac,a, Gabriel Ibarra-Berastegid,b aDept. Applied Physics II, Universidad del Pa´ısVasco-Euskal Herriko Unibertsitatea (UPV/EHU), Barrio Sarriena s/n, 48940-Leioa, Spain bBEGIK Joint Unit IEO-UPV/EHU, Plentziako Itsas Estazioa (PIE, UPV/EHU), Areatza Pasealekua, 48620 Plentzia, Spain cMeteorology Area, Energy and Environment Division, TECNALIA R&I, Basque Country, Spain. dDept. Nuc. Eng. and Fluid Mech., Universidad del Pa´ısVasco-Euskal Herriko Unibertsitatea (UPV/EHU), Alameda Urquijo s/n, 48013-Bilbao, Spain. Abstract This is a non-peer reviewed preprint submitted to EarthArxiv. It is the authors' version initially submitted to Computers and Geo- sciences corresponding to final paper "Analysis of atmospheric ther- modynamics using the R package aiRthermo", published in Com- puters and Goesciences after substantial changes during the review process with doi: http://doi.org/10.1016/j.cageo.2018.10.007, which the readers are encouraged to cite. In this paper the publicly available R package aiRthermo is presented. It allows the user to process information relative to atmospheric thermodynam- ics that range from calculating the density of dry or moist air and converting between moisture indices to processing a full sounding, obtaining quantities such as Convective Available Potential Energy, additional instability indices or adiabatic evolutions of particles. It also provides the possibility to present the information using customizable St¨uve diagrams. Many of the functions are writ- ten inside a C extension so that the computations are fast. The results of an application to five years of real soundings over the Iberian Peninsula are also shown as an example.
    [Show full text]
  • Moist Potential Vorticity Generation in Extratropical Cyclones Zuohao
    Moist Potential Vorticity Generation in Extratropical Cyclones bv Zuohao Cao A thesis submitted in conformity with the requirements for the Degree o f Doctor o f Philosophy in the Department o f Physics o f the University o f Toronto © Copyright by Zuohao Cao 1995 Reproduced with permission of the copyright owner. Further reproduction prohibited without permission For my parents Mr. Jie Cao and Mrs. Sujuan Vang Reproduced with permission of the copyright owner. Further reproduction prohibited without permission. Moist Potential Vorticity Generation in Extratropical Cyclones Zuohao Cao, Ph.D. 1995 Department of Physics, University of Toronto Abstract The mechanism o f moist potential vorticity (M PV) generation in a three- dimensional moist adiabatic and frictionless flow is investigated. It is found that MPV generation is governed by baroclinic vectors and moisture gradients. Negative (positive) M PV can be generated in the region where baroclinic vectors have a component along (against) the direction o f moisture gradients. Numerical simulations of extratropical cyclones with different moisture distributions show that at the different stages of cvclogenesis, negative MPV usually appears in the warm sector near the north part o f the cold front, the bent-back warm front, the warm core and the cold front. The effects o f the Boussinesq approximation on the distribution of vorticity and M PV are examined. The Boussinesq approximation neglects several components o f the ^olenoidal term in the vorticity equation. As a result, it underestimates the thermally direct circulations in the cold front and the bent-back warm front by 25% to 30%. This effect is more pronounced when latent heat release is taken into account.
    [Show full text]
  • Estimating Air Density Using Observations and Re-Analysis Outputs for Wind Energy Purposes
    energies Article Estimating Air Density Using Observations and Re-Analysis Outputs for Wind Energy Purposes Rogier Floors * and Morten Nielsen DTU Wind Energy, Technical University of Denmark, Risø Campus, 4000 Roskilde, Denmark; [email protected] * Correspondence: rofl@dtu.dk; Tel.: +45-4480-5090 Received: 26 April 2019; Accepted: 23 May 2019; Published: 28 May 2019 Abstract: A method to estimate air density as a function of elevation for wind energy resource assessments is presented. The current practice of using nearby measurements of pressure and temperature is compared with a method that uses re-analysis data. It is found that using re-analysis data to estimate air density gives similar or smaller mean absolute errors compared to using measurements that were on average located 40 km away. A method to interpolate power curves that are valid for different air densities is presented. The new model is implemented in the industry-standard model for wind resource assessment and compared with the current version of that model and shown to lead to more accurate assessment of the air density at different elevations. Keywords: air density; wind energy; power curves; WAsP; ERA5; CFSR 1. Introduction Wind power capacity has been growing rapidly over the past decades and reached approximately 539 GW globally in 2017 and continued growth is expected over the next five years [1]. The success of new wind projects is largely dependent on the estimated annual energy production (AEP) and its uncertainty at a proposed wind farm site. The wind power density at a turbine site depends on the wind speed cubed, so there is usually a strong focus on measuring and modelling the wind speed frequency distribution with an accuracy that is as high as possible [2–4].
    [Show full text]