CE 394 K 2 Hydrology Atmospheric Water and

  • Slides: 72
Download presentation
CE 394 K. 2 Hydrology Atmospheric Water and Precipitation • Literary quote for today:

CE 394 K. 2 Hydrology Atmospheric Water and Precipitation • Literary quote for today: “In Köhln, a town of monks and bones, And pavements fang'd with murderous stones And rags, and hideous wenches; I counted two and seventy stenches, All well defined, and several stinks! Ye nymphs that reign o'er sewers and sinks, The river Rhine, it is well known, Doth wash your city of Cologne; But tell me, nymphs, what power devine Shall henceforth wash the river Rhine? ” Samuel Taylor Coleridge, “The City of Cologne”, 1800 Contributed by Eric Hersh

Questions for today (1) How is net radiation to the earth’s surface partitioned into

Questions for today (1) How is net radiation to the earth’s surface partitioned into latent heat, sensible heat and ground heat flux and how does this partitioning vary with location on the earth? (2) What are the factors that govern the patterns of atmospheric circulation over the earth? (3) What are the key variables that describe atmospheric water vapor and how are they connected? (4) What causes precipitation to form and what are the factors that govern the rate of precipitation? (5) How is precipitation measured and described? (Some slides in this presentation were prepared by Venkatesh Merwade)

Questions for today (1) How is net radiation to the earth’s surface partitioned into

Questions for today (1) How is net radiation to the earth’s surface partitioned into latent heat, sensible heat and ground heat flux and how does this partitioning vary with location on the earth? (2) What are the factors that govern the patterns of atmospheric circulation over the earth? (3) What are the key variables that describe atmospheric water vapor and how are they connected? (4) What causes precipitation to form and what are the factors that govern the rate of precipitation? (5) How is precipitation measured and described? (Some slides in this presentation were prepared by Venkatesh Merwade)

Heat energy • Energy – Potential, Kinetic, Internal (Eu) • Internal energy – Sensible

Heat energy • Energy – Potential, Kinetic, Internal (Eu) • Internal energy – Sensible heat – heat content that can be measured and is proportional to temperature – Latent heat – “hidden” heat content that is related to phase changes

Energy Units • In SI units, the basic unit of energy is Joule (J),

Energy Units • In SI units, the basic unit of energy is Joule (J), where 1 J = 1 kg x 1 m/s 2 • Energy can also be measured in calories where 1 calorie = heat required to raise 1 gm of water by 1°C and 1 kilocalorie (C) = 1000 calories (1 calorie = 4. 19 Joules) • We will use the SI system of units

Energy fluxes and flows • Water Volume [L 3] (acre-ft, m 3) • Water

Energy fluxes and flows • Water Volume [L 3] (acre-ft, m 3) • Water flow [L 3/T] (cfs or m 3/s) • Water flux [L/T] (in/day, mm/day) • Energy amount [E] (Joules) • Energy “flow” in Watts [E/T] (1 W = 1 J/s) • Energy flux [E/L 2 T] in Watts/m 2 Energy flow of 1 Joule/sec Area = 1 m 2

Mega. Joules • When working with evaporation, its more convenient to use Mega. Joules,

Mega. Joules • When working with evaporation, its more convenient to use Mega. Joules, MJ (J x 106) • So units are – Energy amount (MJ) – Energy flow (MJ/day, MJ/month) – Energy flux (MJ/m 2 -day, MJ/m 2 -month)

Internal Energy of Water vapor Water Ice Water Heat Capacity (J/kg-K) 2220 4190 Latent

Internal Energy of Water vapor Water Ice Water Heat Capacity (J/kg-K) 2220 4190 Latent Heat (MJ/kg) 0. 33 2. 5/0. 33 = 7. 6 2. 5 Water may evaporate at any temperature in range 0 – 100°C Latent heat of vaporization consumes 7. 6 times the latent heat of fusion (melting)

Latent heat flux • Water flux – Evaporation rate, E (mm/day) • Energy flux

Latent heat flux • Water flux – Evaporation rate, E (mm/day) • Energy flux – Latent heat flux (W/m 2), Hl r = 1000 kg/m 3 lv = 2. 5 MJ/kg 28. 94 W/m 2 = 1 mm/day Area = 1 m 2

Radiation • Two basic laws – Stefan-Boltzman Law • R = emitted radiation (W/m

Radiation • Two basic laws – Stefan-Boltzman Law • R = emitted radiation (W/m 2) • e = emissivity (0 -1) • s = 5. 67 x 10 -8 W/m 2 -K 4 • T = absolute temperature (K) All bodies emit radiation – Wiens Law • l = wavelength of emitted radiation (m) Hot bodies (sun) emit short wave radiation Cool bodies (earth) emit long wave radiation

Net Radiation, Rn Ri Incoming Radiation Re Ro =a. Ri Reflected radiation a= albedo

Net Radiation, Rn Ri Incoming Radiation Re Ro =a. Ri Reflected radiation a= albedo (0 – 1) Rn Net Radiation Average value of Rn over the earth and over the year is 105 W/m 2

Net Radiation, Rn H – Sensible Heat LE – Evaporation G – Ground Heat

Net Radiation, Rn H – Sensible Heat LE – Evaporation G – Ground Heat Flux Rn Net Radiation Average value of Rn over the earth and over the year is 105 W/m 2

Energy Balance of Earth 6 70 20 100 6 26 4 38 15 19

Energy Balance of Earth 6 70 20 100 6 26 4 38 15 19 21 51 Sensible heat flux 7 Latent heat flux 23 http: //www. uwsp. edu/geo/faculty/ritter/geog 101/textbook/energy/radiation_balance. html

Energy balance at earth’s surface Downward short-wave radiation, Jan 2003 600 Z

Energy balance at earth’s surface Downward short-wave radiation, Jan 2003 600 Z

Energy balance at earth’s surface Downward short-wave radiation, Jan 2003 900 Z

Energy balance at earth’s surface Downward short-wave radiation, Jan 2003 900 Z

Energy balance at earth’s surface Downward short-wave radiation, Jan 2003 1200 Z

Energy balance at earth’s surface Downward short-wave radiation, Jan 2003 1200 Z

Energy balance at earth’s surface Downward short-wave radiation, Jan 2003 1500 Z

Energy balance at earth’s surface Downward short-wave radiation, Jan 2003 1500 Z

Energy balance at earth’s surface Downward short-wave radiation, Jan 2003 1800 Z

Energy balance at earth’s surface Downward short-wave radiation, Jan 2003 1800 Z

Energy balance at earth’s surface Downward short-wave radiation, Jan 2003 2100 Z

Energy balance at earth’s surface Downward short-wave radiation, Jan 2003 2100 Z

Latent heat flux, Jan 2003, 1500 z

Latent heat flux, Jan 2003, 1500 z

Questions for today (1) How is net radiation to the earth’s surface partitioned into

Questions for today (1) How is net radiation to the earth’s surface partitioned into latent heat, sensible heat and ground heat flux and how does this partitioning vary with location on the earth? (2) What are the factors that govern the patterns of atmospheric circulation over the earth? (3) What are the key variables that describe atmospheric water vapor and how are they connected? (4) What causes precipitation to form and what are the factors that govern the rate of precipitation? (5) How is precipitation measured and described? (Some slides in this presentation were prepared by Venkatesh Merwade)

Heating of earth surface • Heating of earth surface is uneven – Solar radiation

Heating of earth surface • Heating of earth surface is uneven – Solar radiation strikes perpendicularly near the equator (270 W/m 2) – Solar radiation strikes at an oblique angle near the poles (90 Amount of energy transferred from W/m 2) equator to the poles is • Emitted radiation is approximately 4 x 109 MW more uniform than incoming radiation

Hadley circulation Warm air rises, cool air descends creating two huge convective cells.

Hadley circulation Warm air rises, cool air descends creating two huge convective cells.

Coriolis Force Cone is moving southward towards the pole Camera fixed in the outer

Coriolis Force Cone is moving southward towards the pole Camera fixed in the outer space Camera fixed on to the globe (cone appears moving straight) (looking southward, cone appears deflecting to the right) the force that deflects the path of the wind on account of earth rotation is called Coriolis force. The path of the wind is deflected to the right in the Northern Hemisphere and the to left in the Southern Hemisphere.

Atmospheric circulation Circulation cells Polar Cell 1. Hadley cell 2. Ferrel Cell 3. Polar

Atmospheric circulation Circulation cells Polar Cell 1. Hadley cell 2. Ferrel Cell 3. Polar cell Winds 1. Tropical Easterlies/Trades 2. Westerlies 3. Polar easterlies Latitudes 1. Intertropical convergence zone (ITCZ)/Doldrums 2. Horse latitudes 3. Subpolar low 4. Polar high

Effect of land mass distribution Uneven distribution of land ocean, coupled with different thermal

Effect of land mass distribution Uneven distribution of land ocean, coupled with different thermal properties creates spatial variation in atmospheric circulation A) Idealized winds generated by pressure gradient and Coriolis Force. B) Actual wind patterns owing to land mass distribution

Shifting in Intertropical Convergence Zone (ITCZ) Owing to the tilt of the Earth's axis

Shifting in Intertropical Convergence Zone (ITCZ) Owing to the tilt of the Earth's axis in orbit, the ITCZ shifts north and south. Southward shift in January Creates wet Summers (Monsoons) and dry winters, especially in India and SE Asia Northward shift in July

ITCZ movement http: //iri. ldeo. columbia. edu/%7 Ebgordon/ITCZ. html

ITCZ movement http: //iri. ldeo. columbia. edu/%7 Ebgordon/ITCZ. html

Questions for today (1) How is net radiation to the earth’s surface partitioned into

Questions for today (1) How is net radiation to the earth’s surface partitioned into latent heat, sensible heat and ground heat flux and how does this partitioning vary with location on the earth? (2) What are the factors that govern the patterns of atmospheric circulation over the earth? (3) What are the key variables that describe atmospheric water vapor and how are they connected? (4) What causes precipitation to form and what are the factors that govern the rate of precipitation? (5) How is precipitation measured and described? (Some slides in this presentation were prepared by Venkatesh Merwade)

Structure of atmosphere

Structure of atmosphere

Atmospheric water • Atmospheric water exists – Mostly as gas or water vapor –

Atmospheric water • Atmospheric water exists – Mostly as gas or water vapor – Liquid in rainfall and water droplets in clouds – Solid in snowfall and in hail storms • Accounts for less than 1/100, 000 part of total water, but plays a major role in the hydrologic cycle

Water vapor Suppose we have an elementary volume of atmosphere d. V and we

Water vapor Suppose we have an elementary volume of atmosphere d. V and we want quantify how much water vapor it contains Water vapor density Air density d. V ma = mass of moist air mv = mass of water vapor Atmospheric gases: Nitrogen – 78. 1% Oxygen – 20. 9% Other gases ~ 1% http: //www. bambooweb. com/articles/e/a/Earth's_atmosphere. html

Specific Humidity, qv • Specific humidity measures the mass of water vapor per unit

Specific Humidity, qv • Specific humidity measures the mass of water vapor per unit mass of moist air • It is dimensionless

Vapor pressure, e • Vapor pressure, e, is the pressure that water vapor exerts

Vapor pressure, e • Vapor pressure, e, is the pressure that water vapor exerts on a surface • Air pressure, p, is the total pressure that air makes on a surface • Ideal gas law relates pressure to absolute temperature T, Rv is the gas constant for water vapor • 0. 622 is ratio of mol. wt. of water vapor to avg mol. wt. of dry air

Dalton’s Law of Partial Pressures John Dalton studied the effect of gases in a

Dalton’s Law of Partial Pressures John Dalton studied the effect of gases in a mixture. He observed that the Total Pressure of a gas mixture was the sum of the Partial Pressure of each gas. P total = P 1 + P 2 + P 3 +. . . . Pn The Partial Pressure is defined as the pressure of a single gas in the mixture as if that gas alone occupied the container. In other words, Dalton maintained that since there was an enormous amount of space between the gas molecules within the mixture that the gas molecules did not have any influence on the motion of other gas molecules, therefore the pressure of a gas sample would be the same whether it was the only gas in the container or if it were among other gases. http: //members. aol. com/profchm/dalton. html

Avogadro’s law Equal volumes of gases at the same temperature and pressure contain the

Avogadro’s law Equal volumes of gases at the same temperature and pressure contain the same number of molecules regardless of their chemical nature and physical properties. This number (Avogadro's number) is 6. 023 X 1023 in 22. 41 L for all gases. Dry air ( z = x+y molecules) Moist air (x dry and y water vapor) Dry air Water vapor rd = (x+y) * Md/Volume rm = (x* Md + y*Mv)/Volume rm < rd, which means moist air is lighter than dry air!

Saturation vapor pressure, es Saturation vapor pressure occurs when air is holding all the

Saturation vapor pressure, es Saturation vapor pressure occurs when air is holding all the water vapor that it can at a given air temperature Vapor pressure is measured in Pascals (Pa), where 1 Pa = 1 N/m 2 1 k. Pa = 1000 Pa

Relative humidity, Rh es e Relative humidity measures the percent of the saturation water

Relative humidity, Rh es e Relative humidity measures the percent of the saturation water content of the air that it currently holds (0 – 100%)

Dewpoint Temperature, Td e Td T Dewpoint temperature is the air temperature at which

Dewpoint Temperature, Td e Td T Dewpoint temperature is the air temperature at which the air would be saturated with its current vapor content

Water vapor in an air column • We have three equations describing column: 2

Water vapor in an air column • We have three equations describing column: 2 – Hydrostatic air pressure, dp/dz = -rag – Lapse rate of temperature, d. T/dz = - a – Ideal gas law, p = ra. Ra. T • Combine them and integrate over column to get pressure variation elevation Column Element, dz 1

Precipitable Water • In an element dz, the mass of water vapor is dmp

Precipitable Water • In an element dz, the mass of water vapor is dmp • Integrate over the whole atmospheric column to get precipitable water, mp • mp/A gives precipitable water per unit area in kg/m 2 2 Column Element, dz 1 Area = A

Precipitable Water, Jan 2003

Precipitable Water, Jan 2003

Precipitable Water, July 2003

Precipitable Water, July 2003

January July

January July

Questions for today (1) How is net radiation to the earth’s surface partitioned into

Questions for today (1) How is net radiation to the earth’s surface partitioned into latent heat, sensible heat and ground heat flux and how does this partitioning vary with location on the earth? (2) What are the factors that govern the patterns of atmospheric circulation over the earth? (3) What are the key variables that describe atmospheric water vapor and how are they connected? (4) What causes precipitation to form and what are the factors that govern the rate of precipitation? (5) How is precipitation measured and described? (Some slides in this presentation were prepared by Venkatesh Merwade)

Precipitation • Precipitation: water falling from the atmosphere to the earth. – Rainfall –

Precipitation • Precipitation: water falling from the atmosphere to the earth. – Rainfall – Snowfall – Hail, sleet • Requires lifting of air mass so that it cools and condenses.

Mechanisms for air lifting 1. Frontal lifting 2. Orographic lifting 3. Convective lifting

Mechanisms for air lifting 1. Frontal lifting 2. Orographic lifting 3. Convective lifting

Definitions • Air mass : A large body of air with similar temperature and

Definitions • Air mass : A large body of air with similar temperature and moisture characteristics over its horizontal extent. • Front: Boundary between contrasting air masses. • Cold front: Leading edge of the cold air when it is advancing towards warm air. • Warm front: leading edge of the warm air when advancing towards cold air.

Frontal Lifting • Boundary between air masses with different properties is called a front

Frontal Lifting • Boundary between air masses with different properties is called a front • Cold front occurs when cold air advances towards warm air • Warm front occurs when warm air overrides cold air Cold front (produces cumulus cloud) Cold front (produces stratus cloud)

Orographic lifting Orographic uplift occurs when air is forced to rise because of the

Orographic lifting Orographic uplift occurs when air is forced to rise because of the physical presence of elevated land.

Convective lifting Convective precipitation occurs when the air near the ground is heated by

Convective lifting Convective precipitation occurs when the air near the ground is heated by the earth’s warm surface. This warm air rises, cools and creates precipitation. Hot earth surface

Condensation • Condensation is the change of water vapor into a liquid. For condensation

Condensation • Condensation is the change of water vapor into a liquid. For condensation to occur, the air must be at or near saturation in the presence of condensation nuclei. • Condensation nuclei are small particles or aerosol upon which water vapor attaches to initiate condensation. Dust particulates, sea salt, sulfur and nitrogen oxide aerosols serve as common condensation nuclei. • Size of aerosols range from 10 -3 to 10 mm.

Precipitation formation • Lifting cools air masses so moisture condenses • Condensation nuclei –

Precipitation formation • Lifting cools air masses so moisture condenses • Condensation nuclei – Aerosols – water molecules attach • Rising & growing – 0. 5 cm/s sufficient to carry 10 mm droplet – Critical size (~0. 1 mm) – Gravity overcomes and drop falls

Forces acting on rain drop • Three forces acting on rain drop – Gravity

Forces acting on rain drop • Three forces acting on rain drop – Gravity force due to weight – Buoyancy force due to displacement of air – Drag force due to friction with surrounding air D Fb Fd Fd Fg

Terminal Velocity • Terminal velocity: velocity at which the forces acting on the raindrop

Terminal Velocity • Terminal velocity: velocity at which the forces acting on the raindrop are in equilibrium. • If released from rest, the raindrop will accelerate until it reaches its terminal velocity D Fb Fd Fd Fg At standard atmospheric pressure (101. 3 kpa) and temperature (20 o. C), rw = 998 kg/m 3 and ra = 1. 20 kg/m 3 V • Raindrops are spherical up to a diameter of 1 mm • For tiny drops up to 0. 1 mm diameter, the drag force is specified by Stokes law

Precipitation Variation • Influenced by – Atmospheric circulation and local factors • Higher near

Precipitation Variation • Influenced by – Atmospheric circulation and local factors • Higher near coastlines • Seasonal variation – annual oscillations in some places • Variables in mountainous areas • Increases in plains areas • More uniform in Eastern US than in West

Rainfall patterns in the US

Rainfall patterns in the US

Global precipitation pattern

Global precipitation pattern

Spatial Representation • Isohyet – contour of constant rainfall • Isohyetal maps are prepared

Spatial Representation • Isohyet – contour of constant rainfall • Isohyetal maps are prepared by interpolating rainfall data at gaged points. Austin, May 1981 Wellsboro, PA 1889

Texas Rainfall Maps

Texas Rainfall Maps

Temporal Representation • Rainfall hyetograph – plot of rainfall depth or intensity as a

Temporal Representation • Rainfall hyetograph – plot of rainfall depth or intensity as a function of time • Cumulative rainfall hyetograph or rainfall mass curve – plot of summation of rainfall increments as a function of time • Rainfall intensity – depth of rainfall per unit time

Rainfall Depth and Intensity

Rainfall Depth and Intensity

Incremental Rainfall Hyetograph

Incremental Rainfall Hyetograph

Cumulative Rainfall Mass Curve

Cumulative Rainfall Mass Curve

Arithmetic Mean Method • Simplest method for determining areal average P 1 = 10

Arithmetic Mean Method • Simplest method for determining areal average P 1 = 10 mm P 1 P 2 = 20 mm P 3 = 30 mm P 2 P 3 • Gages must be uniformly distributed • Gage measurements should not vary greatly about the mean

Thiessen polygon method • • • Any point in the watershed receives the same

Thiessen polygon method • • • Any point in the watershed receives the same amount of rainfall as that at the nearest gage Rainfall recorded at a gage can be applied to any point at a distance halfway to the next station in any direction Steps in Thiessen polygon method 1. Draw lines joining adjacent gages 2. Draw perpendicular bisectors to the lines created in step 1 3. Extend the lines created in step 2 in both directions to form representative areas for gages 4. Compute representative area for each gage 5. Compute the areal average using the following formula P 1 A 1 P 2 A 2 P 3 A 3 P 1 = 10 mm, A 1 = 12 Km 2 P 2 = 20 mm, A 2 = 15 Km 2 P 3 = 30 mm, A 3 = 20 km 2

Isohyetal method • Steps – Construct isohyets (rainfall contours) – Compute area between each

Isohyetal method • Steps – Construct isohyets (rainfall contours) – Compute area between each pair of adjacent isohyets (Ai) – Compute average precipitation for each pair of adjacent isohyets (pi) – Compute areal average using the following formula 10 20 P 1 A 1=5 , p 1 = 5 A 2=18 , p 2 = 15 P 2 A 3=12 , p 3 = 25 30 P 3 A 4=12 , p 3 = 35

Inverse distance weighting • Prediction at a point is more influenced by nearby measurements

Inverse distance weighting • Prediction at a point is more influenced by nearby measurements than that by distant measurements • The prediction at an ungaged point is inversely proportional to the distance to the measurement points • Steps – Compute distance (di) from ungaged point to all measurement points. – Compute the precipitation at the ungaged point using the following formula P 1=10 P 2= 20 d 1=25 d 2=15 P 3=30 d 3=10 p

Rainfall interpolation in GIS • Data are generally available as points with precipitation stored

Rainfall interpolation in GIS • Data are generally available as points with precipitation stored in attribute table.

Rainfall maps in GIS Nearest Neighbor “Thiessen” Polygon Interpolation Spline Interpolation

Rainfall maps in GIS Nearest Neighbor “Thiessen” Polygon Interpolation Spline Interpolation

NEXRAD • NEXt generation RADar: is a doppler radar used for obtaining weather information

NEXRAD • NEXt generation RADar: is a doppler radar used for obtaining weather information • A signal is emitted from the radar which returns after striking a rainfall drop • Returned signals from the radar are analyzed to compute the rainfall intensity and integrated over time to get the precipitation NEXRAD Tower Working of NEXRAD

NEXRAD data • NCDC data (JAVA viewer) – http: //www. ncdc. noaa. gov/oa/radar/jnx/ •

NEXRAD data • NCDC data (JAVA viewer) – http: //www. ncdc. noaa. gov/oa/radar/jnx/ • West Gulf River Forecast Center – http: //www. srh. noaa. gov/wgrfc/ • National Weather Service Animation – http: //weather. noaa. gov/radar/mosaic. loop/DS. p 19 r 0/ar. us. conus. shtml