1 / 60

Zong-Liang Yang

Introduction to the NOAH Land Surface Model. N : N ational Centers for Environmental Prediction (NCEP) O : O regon State University (Dept of Atmospheric Sciences) A : A ir Force (both AFWA and AFRL - formerly AFGL, PL)

mikel
Download Presentation

Zong-Liang Yang

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Introduction to the NOAH Land Surface Model N: National Centers for Environmental Prediction (NCEP) O: Oregon State University (Dept of Atmospheric Sciences) A: Air Force (both AFWA and AFRL - formerly AFGL, PL) H: Hydrologic Research Lab - NWS (now Office of Hydrologic Dev – OHD) • Outline • What is an LSM? • What is the NOAH LSM? • NOAH @ UT Zong-Liang Yang Prepared for Surface Water Hydrology February 24, 2005

  2. Outline • What is an LSM? (29 slides) • What is the NOAH LSM? (13 slides) • NOAH @ UT (14 slides)

  3. Computer Code to Model Land Surface Processes • Exchanges with the atmosphere • Momentum • Energy • Vapor • Trace gases/dust/aerosols • Exchanges with the oceans • Runoff • Sediments / nitrogen • Land-memory processes • Topography • Snow cover • Soil moisture • Vegetation • Human activities • Land use (agriculture, deforestation, urbanization…) • Inland water pollution

  4. Earth’s Radiation Budget Shortwave:Earth reflects 30% directly back to space, absorbs about 20% in the atmosphere, and absorbs about 50% at the surface. Longwave: Greenhouse effects keep the surface warm. Turbulent heat fluxes: sensible and latent heat keeps the surface cool.

  5. Land-atmosphere coupling strength (JJA) Koster et al., 2004 (Science) Hot spots: initial soil moisture impacts rainfall forecasts Locations for the routine monitoring of soil moisture

  6. L Qa Out Sa Qa In Pa Pm E Ql Interpretations • Rainfall in a region depends on atmospheric moisture transport from upwind areas and local evapotranspiration. • Moisture transport depends on large-scale circulations • Evapotranspiration depends on vegetation, soil, and weather conditions

  7. Do Land Surface Processes Matter to Climate Prediction?

  8. What Are Land Surface Processes • Land surface consists of • urban areas, soil, vegetation, snow, topography, inland water (lake, river) … • Land surface processes describe • exchanges of momentum, energy, water vapor, and other trace gases between land surface and the overlying atmosphere • states of land surface (e.g., soil moisture, soil temperature, canopy temperature, snow water equivalent) • characteristics of land surface (e.g., roughness, albedo, emissivity, soil texture, vegetation type, cover extent, leaf area index, and seasonality)

  9. What Are Land Surface Processes • Land surface processes function as • lower boundary condition in Atmospheric Models • Atmospheric Boundary Layer Simulation Climate Simulation • Numerical Weather Prediction 4-D Data Assimilation • upper boundary condition in Hydrological Models • Water Resources Estimation; Crop Water Use; Runoff Simulation • interface for coupled Atmospheric/Hydrological/Ecological Models

  10. Land Surface Models (LSMs) • Computer code describing land surface processes (also called LSSs, LSPs, SVATs) • FORTRAN, C, C++, PASCAL, … ... • Tens to thousands of lines • There are a huge number of LSMs (100+ examples in literature) • many are just “research models’’, local-scale oriented, with specific process emphasis • up to ~100 canopy, ~100 soil, ~100 snow, even ~100 atmosphere layers! • LSMs in GCMs and Hydrological Models are less diverse • one dimensional, with 1-2 canopy, 1-10 soil, 1-10 snow layers • three general classes • “Bucket” Models (no vegetation canopy) • “Micrometeorological” Models (detailed soil/snow/canopy processes) • “Intermediate” Models (some soil/snow/canopy features)

  11. Land Surface Models (LSMs) • Four basic requirements • frequently-sampled (hourly or sub-hourly) weather “data” to “drive” LSMs • precipitation (rate; coverage, large-scale/convective in GCMs) • radiation (shortwave, longwave) • temperature • wind components (u, v) • specific humidity • surface pressure • initialization of state variables • soil moisture (liquid, frozen) deep soil temperature • specification of surface characteristics • albedo roughness vegetation cover • validation of simulations of state variables and fluxes • soil moisture sensible/latent heat fluxes

  12. Land Surface Models (LSMs) • Minimum Requirement is to describe: • bulk momentum exchange (roughness length, zero displacement height) • exchange of radiation for: • solar radiation (0-3 micrometers): albedo • long-wave (3-100 micrometers): surface emissivity • partition of radiant energy between soil heat, latent heat and sensible heat, along with its relationship to water availability • Micrometeorological Models typically • have 20-50+ parameters • require to prescribe or predefine them before running the models • Research Issues • Obtaining and applying relevant “pure biome” data to test or calibrate LSMs • Dealing with spatial/temporal heterogeneity • defining area-average parameters; scaling up surface parameters from local to regional • defining space-time structure of atmospheric inputs • Making best use of remote sensing data for initialization, specification and validation • Improving key processes • Snow/Frozen soil • Runoff generation/routing • “Greening” of LSMs (carbon balance and vegetation dynamics)

  13. Land Surface Models (LSMs) • Project for Intercomparison of Land-surface Parameterization Schemes (PILPS) • is helping define needs, understanding processes, from sensitivity tests and field data to coupled modeling • ongoing project (4 phases, starting from 1992, 20-30 LSMs participating) • Phase 1: GCM-generated forcing data for three “biomes”: • models’ response time to initialization of soil moisture • sensitivity experiments with albedo, canopy/soil water storage, aerodynamic resistances • Phase 2 : observed data • 2(a) Cabauw, Netherlands (grass; 1-year) • 2(b) HAPEX-MOBILHY, France (crop; 1-year) • 2(c ) Red-Arkansas River Basin, USA (mixed land covers; 10-year) • 2(d) Valdai, Russia (snow/frozen soil; 18-year) • 2(e) Sweden Torne/Kalix Basins, 58000 sq. km, 218 grids, 1979-98 • Phase 3: joint with Atmospheric Model Intercomparison Project -AMIP • Phase 4: selected LSMs coupled to one host model • Regional Model Global Model Coupling Strategy

  14. Simple “Bucket” LSMs • “Bucket” Models” first in hydrological models, later in GCMs • “Bucket” is filled by precipitation, and emptied by evaporation • Runoff occurs when “Bucket” is full (soil water level exceeds the field capacity) • Evaporation is based on soil-water limited “Potential Evaporation” • The limiting factor is a linear function of soil moisture • In original “bare soil” form, e.g., Manabe (1969) include no vegetation so specify (fixed) parameter values for • albedo (Solar radiation) • surface emissivity (Longwave radiation) • roughness length (Momentum exchange) • soil thermal properties (Soil heat storage) • field capacity (15 cm) for all land points (Soil hydrology) • snow areal coverage (covers 100% gridbox) (Snow energy/water budgets) • Some variants of “Bucket” Models • allowing leakage by implicit transpiration or gravitational drainage • including surface resistance in evaporation formula • extending to 2 connected “Buckets”: Deardorff (1977) includes a thin surface layer and a deep bulk layer

  15. Micrometeorological Model-Based LSMs • Micrometeorological Model-Based LSMs are • more complex than bucket models • less complex than “multi-layer” or “second order closure” canopy models • Best known examples: • “Biosphere-Atmosphere Transfer Scheme (BATS)” • “Simple Biosphere Model (SiB)” • SiB more complex (not at run time) some whole-canopy parameters from “pre-processing”, one calculation, prior to actual model run • Five components of radiation exchange considered in BATS and SiB • Direct beam photosynthetically active radiation (PAR) (<0.72 micrometers) • Diffuse solar photosynthetically active radiation (PAR) (<0.72 micrometers) • Direct beam near infrared (NIR) (0.72-4.0 micrometers) • Diffuse solar near infrared (NIR) (0.72-4.0 micrometers) • Thermal infrared (>4.0 micrometers)

  16. Micrometeorological Model-Based LSMs Biosphere-Atmosphere Transfer Scheme BATS

  17. Micrometeorological Model-Based LSMs • “Two-Stream Approximation” are used to calculate solar radiation exchange in SiB for both PAR and NIR, diffuse and direct beam components from • Can obtain single-scattering albedo, diffuse fluxes within canopy (for stomatal resistances) and at soil surface, and total albedo of the vegetated surface • Vegetation albedo depends on specifications of • leaf area index scattering coefficients of leaves and soil • leaf angle distribution function angle of the direct-beam incident radiation • proportion of different components of the incident radiation • Not suitable for application over heterogeneous or ‘clumpy’ canopies (coniferous forests); require three-dimensional radiation models • BATS prescribes vegetation albedo; 4-layer model for stomata

  18. Micrometeorological Model-Based LSMs • Thermal Radiation from surface is computed from “Stefan-Boltzmann” law: • Whole-canopy turbulent transfer is important for exchange of momentum, heat and mass (water vapor, trace gases) • is a complex process • requires higher-order closure models, Markov-chain numerical simulations, or eddy simulation techniques • simpler treatments in LSMs because of • computational expense and storage requirements • aerodynamic resistance roughly an order of magnitude smaller than aerodynamic resistance • prescribed for BATS (and most other models); aerodynamic resistances from curve-fitting functions, iteration for nonneutral corrections • calculated in SiB: pre-processed from canopy morphology using a “K-Theory” equivalent of Shaw-Pereira (1982); aerodynamic resistances from canopy morphology using “K-Theory”, iteration for nonneutral corrections

  19. Micrometeorological Model-Based LSMs • Canopy interactions described by a “big leaf” model in BATS and SiB and others [SiB allows an “understorey”; SiB2 allows only one layer] • Dry canopy parameterized through surface resistance based on the model of Jarvis (1976). For a single leaf: • Canopy surface resistance is obtained by integrating above equation over leaf angle and canopy depth. This integration and those dependencies are “coded” in different forms in BATS and SiB and many other models.

  20. Micrometeorological Model-Based LSMs • Different forms of coding for the integration: • In SiB, an analytical integration is performed by using the Goudriaan formulation for some regular leaf angle distribution • In BATS, the integration is done by using a 4-layer numerical scheme for canopy • Different forms of coding for the stress factors: • primarily determined from curves fits to data (e.g. Jarvis, 1976), but • in SiB was a linear function of the leaf water potential, which, in turn, was determined by a catenary model of water flow from root zone to leaf • in BATS, was introduced only when the atmospheric demand exceeds the water supply from the soil through roots • in SiB taken as a smooth curve varying from zero at some minimum temperature around freezing to unity at an optimum temperature between 25 degree C and 35 degree C, depending on species, and dropping sharply down to zero at 45 degree C to 55 degree C • in BATS parameterized as a quadratic function of canopy temperature, varying from 0 at freezing point to unity at 25 degree C, and remaining at unity when temperature is greater than 25 degree C

  21. Micrometeorological Model-Based LSMs • Calibration of stomatal resistance largely speculative: guessed from single leaf data, but occasional “Calibration Data Sets” available e.g. for forest, Shuttleworth (1989); Gash et al. (1996) • Wet Canopy Description is a simplified “Rutter” Model: S, proportional to Leaf Area Index ~[(0.1-0.25) x LAI] (mm) • Both BATS and SiB allow partially wet canopy

  22. Micrometeorological Model-Based LSMs • “Plot-scale” calibration of interception models adequate, given tuning of S against LAI • BATS/SiB have simple “vertical hydrology”: soil water diffuses vertically through 2-3 layers, but is ultimately inaccessible by draining from lowest layer • Global Distributions of different land classes required, BATS has 15, SiB has 12, with parameters defined for each class: Calibration of these Parameters Remains Speculative

  23. “Intermediate” LSMs • “Intermediate” LSMs trade process detail for simplicity, justified by still speculative need for more complete process models, and their poor calibration • Often have distributed vegetation class, but simpler physics: e.g. • UKMO simple canopy model • French ISBA simple soil model • Simple SiB Simplified empirical parameterization fitted to time consuming, SiB-calculated functions • “Super Bucket” basic bucket, with “reference crop” estimates instead of “potential” rate(+Several others) • Permissible, practical alternative, pending evidence of need for more complex micrometeorological models

  24. LSMs Calibration Issues • Primary LSM short-coming is lack of reliable calibration: 3 needs • Obtaining and applying relevant “pure biome” data • Understanding how to define area-average parameters • Making best use of Remote Sensing Data • Single Biome Calibration • Needs: • sets of continuous “driving variables” for whole years (frequently-sampled, near-surface, equivalent to GCM) • some substantial periods with measured exchanged fluxes (preferably for several seasons, all exchanges, time scales) • Some biomes done or in-hand • Rain forest (ARME, ABRACOS) • Prairie grass (FIFE) • Boreal forest (BOREAS) • Sonora Desert, Tucson, Arizona (Example of multi-objective calibration)

  25. LSMs Calibration Issues: Area-average • Defining Area-average LSM parameters is active research which involves data collection, analysis and coupled models for mixed land covers, and for selected regions • mixed vegetation cover and/or soils uneven precipitation sloping terrain • Examples: • Hydrologic-Atmospheric Pilot Experiment (HAPEX) - S.W. France, 1986, (mixed temperate; agriculture, forest) • First ISLSCP Field Experiment (FIFE) - Kansas, 1987,9 (tall prairie grass) • HAPEX-Sahel - Niger, 1992, (semi-arid Sahelian savannah-across rainfall gradient) • BOREAS - Canada, 1993,4, Boreal Forest • AMAZONIA (LBA) - S. America, 1997,8 - Amazonian land covers • GEWEX

  26. LSMs Calibration Issues: Atmospheric Data • Calibration or model development depends on atmospheric input data • Runoff estimation error is determined by the accuracy of the precipitation forcing as measured by the density of rain gauge

  27. Micrometeorological Model-Based LSMs • New developments: The Greening of LSMs (improved canopy physiology) • enzyme kinetics-electron transport model (Farquhar et al., 1980), relating leaf photosynthetic rate to PAR, leaf temperature, leaf internal carbon dioxide concentration, and leaf carboxylase concentration • photosynthesis modeled as the minimum of three potential capacities to fix carbon

  28. Micrometeorological Model-Based LSMs • New developments: The Greening of LSMs (improved canopy physiology) • Ball-Berry model of stomatal conductance (Ball et al., 1986): robust semi-empirical method, relating leaf stomatal conductance to photosynthetic rate and environmental parameters • Other forms: Leuning (1995), Foley et al. (1996), Dickinson et al. (1998): using vapor pressure deficit instead of relative humidity

  29. Micrometeorological Model-Based LSMs • New developments: The Greening of LSMs (improved canopy physiology) • C3 (all trees and many herbaceous) and C4 (warm grasses) plants are similarly modeled as the minimum of three potential capacities to fix carbon, but • The three potential capacities are defined differently for C3 and C4 plants • The parameters are different for both types of plants • The Ball-Berry model of stomatal conductance replaces all the terms in the Jarvis model except for the soil water stress term

  30. Micrometeorological Model-Based LSMs • New developments: The Greening of LSMs (improved vegetation phenology) • Traditionally, annual cycle in areal vegetation coverage and leaf area index prescribed as a quadratic function of deep soil temperature • Some recent attempts: relating vegetation to climate using empirical rules involving monthly or annual mean temperature and precipitation • New developments: • simple rule-based formulation to describe behavior of winter-deciduous and drought-deciduous plants (Foley et al., 1996), using • daily average temperature, yearly carbon balance • competition between trees and grasses (through allocation/shading) • interactive canopy model (Dickinson et al., 1998) • short time-scale leaf dynamics • consistent treatment of stomatal conductance and assimilation • assimilated carbon allocated into leaf, fine root and wood • soil carbon model • fast pool: leaf/root turnover (senescence, herbivory or mechanical), leaf loss (cold and drought stress), soil respiration • slow pool: inert for seasonal variations

  31. Recent LSM Developments • Improved understanding in physical hydrology • Snow hydrology • Topographic controls on soil moisture/runoff • Enhanced linkage with biogeochemistry/ecology • Carbon and nitrogen • Vegetation growth/competition • Closed water cycle through river routing • Advanced techniques for using remotely-sensed land parameters • Extended field datasets and model calibration

  32. Outline • What is an LSM? (29 slides) • What is the NOAH LSM? (13 slides) • NOAH @ UT (14 slides)

  33. Land-Atmosphere Interaction

  34. NOAH LSM Development Team NCEP/EMC: NCEP Environmental Modeling Center (EMC) (Mitchell, Chen, Ek, Lin, Marshall, Janjic, Manikin, Lohmann, Grunmann, Pan) OSU: Oregon State University (Mahrt, Pan, Ek, Kim, Rusher) HL: NWS Hydrology Lab - formerly Office of Hydrology (Schaake, Koren, Duan) AFWA: Air Force Weather Agency - formerly AFGWC (Moore, Mitchell, Gayno) AFRL: Air Force Research Lab - formerly AFGL and PL (Mitchell, Hahn, Chang, Yang)

  35. NOAH LSM Evaluation Team E.H. Berbery and Rasmusson: U. Maryland (ARM/CART) C. Marshall and Crawford U. Oklahoma (OU Mesonet) I. Yucel and Shuttleworth: U. Arizona (ARM/CART, AZNET) A.K. Betts: Atmospheric Res Inc (ISLSCP/FIFE) C.D. Peters-Lidard, Wood Princeton U. (TOPLATS extensions) L. Hinkelman and Ackerman: Penn State U. (ARM/CART) T.H. Chen, W. Qu, Henderson-Sellers, et al. RMIT (PILPS-2a) E. Wood, Lettenmaier, Liang, Lohmann: Princeton U. (PILPS-2c) A. Schlosser, A.G. Slater, A. Robock, et al. U. Maryland (PILPS-2d) R. Angevine NOAA/AL (Flatland Exp) 

  36. References 1. User’s Guide, Public Release Version 2.2, Last Updated 25 June 2001, filename README_2.2.doc at ftp://ftp.ncep.noaa.gov/pub/gcp/ldas/noahlsm/ver_2.2 2. Author / Point of Contact: Kenneth.Mitchell@noaa.gov (NCEP/EMC), phone 301-763-8000 3. PROGRAM HISTORY LOG 01 Mar 99: Ver_1.0 08 Mar 99: Ver_1.1 27 Jul 99: Ver_2.0 23 Oct 00: Ver_2.1 25 Mar 01: Ver_2.2

  37. Timeline of NOAH LSM Evolution

  38. Atmospheric Forcing Variables Air temperature at height Z above ground Air humidity at height Z above ground Surface pressure at height Z above ground Wind speed at height Z above ground Surface downward longwave radiation Surface downward solar radiation Precipitation

  39. Initial Conditions • SMC: total volumetric soil moisture (liquid & frozen) in each soil layer • SH2O: liquid volumetric soil moisture in each soil layer • STC: temperature in each soil layer • T1: skin temperature • CMC: canopy water content • SNOWH: snow depth • SNEQV: water-equivalent snow depth

  40. Model Parameters • There are five kinds of land-surface parameters, reviewed in order below. • single universal values • values dependent on the soil class index (default categories are 1–9) • values dependent on the vegetation class index (default categories are 1–13) • values dependent on the surface slope index (default categories are 1–7) • parameters specifying the numbers of vegetation, soil, and slope classes

  41. Model Soil Parameters SMCMAX: maximum volumetric soil moisture (porosity) SMCREF: soil moisture threshold for onset of some transpiration stress SMCWLT: soil moisture wilting point at which transpiration ceases SMCDRY: top layer soil moisture threshold at which direct evaporation from soil ceases DKSAT: saturated soil hydraulic conductivity PSISAT: saturated soil matric potential B: the “b” parameter in hydraulic functions DWSAT: saturated soil water diffusivity QUARTZ: quartz content, used to compute soil thermal diffusivity FRZFACT: a parameter used with FRZK to compute the value of parameter FRZX

  42. Model Vegetation Parameters Z0 (m): roughness length RCMIN (s/m) : minimal stomatal resistance used in canopy resistance of routine CANRES RGL: radiation stress parameter used in F1 term in canopy resistance of routine CANRES HS: coefficient used in vapor pressure deficit term F2 in canopy resistance of routine CANRES LAI: presently set to universal value of 3.0 across all vegetation classes Note: seasonality of vegetation greenness carried by fraction of green vegetation (SHDFAC) NROOT: number of soil layers from top down reached by roots SNUP: the water-equivalent snowdepth upper threshold at which 1)    100 percent snow cover is achieved for given veg class 2) maximum snow albedo is achieved for given veg class

  43. Forecasts of Dew Point Temperature Old version New version

  44. Forecasts of Dew Point Temperature

  45. Outline • What is an LSM? (29 slides) • What is the NOAH LSM? (13 slides) • NOAH @ UT (14 slides)

  46. NOAH @ UT • To put a dynamic vegetation component in the NOAH LSM • understand the feedback mechanisms between precipitation, soil moisture, and vegetation dynamics • To improve weather and climate forecasts on time scales from hourly to seasonal • To interpret and transfer weather and climate forecasts for water resources management applications • To put a biogenic VOC emission component in the NOAH LSM • understand the impacts of climate change and land use/land cover change on BVOC emissions and the formation of near ground ozone

  47. The Modeling Domains: Multiple interactive nesting grids: Weather Research and Forecast (WRF) Model, version 2.02 (http://www.wrf-model.org/) The model is developed by NCAR, NCEP and Universities. 30 km grid; 31 Levels Temperature, wind, humidity, pressure predicted 3-10 km grid; 31 Levels Temperature, wind, humidity, pressure predicted 90 km grid; 31 Levels Temperature, wind, humidity, pressure prescribed from large-scale analyses

  48. Vegetation Distribution in the Model Other Land Data: Vegetation LAI Albedo Roughness length Stomatal parameters Soil Topography Validation: Weather stations Hydrologic stations Remote sensing Field data

  49. Stomatal Conductance in NOAH • Transpiration parameterized through surface conductance based on the model of Jarvis (1976). For a single leaf:

More Related