Assessment of the Land Surface and Boundary Layer Models in Two Operational Versions of the NCEP Eta Model Using FIFE Data

Alan K. Betts Pittsford, Vermont

Search for other papers by Alan K. Betts in
Current site
Google Scholar
PubMed
Close
,
Fei Chen NCEP/EMC, Washington, DC

Search for other papers by Fei Chen in
Current site
Google Scholar
PubMed
Close
,
Kenneth E. Mitchell NCEP/EMC, Washington, DC

Search for other papers by Kenneth E. Mitchell in
Current site
Google Scholar
PubMed
Close
, and
Zaviša I. Janjić NCEP/EMC, Washington, DC

Search for other papers by Zaviša I. Janjić in
Current site
Google Scholar
PubMed
Close
Restricted access

Abstract

Data from the 1987 summer FIFE experiment for four pairs of days are compared with corresponding 48-h forecasts from two different versions of the Eta Model, both initialized from the NCEP–NCAR (National Centers for Environmental Prediction–National Center for Atmospheric Research) global reanalysis. One used the late 1995 operational Eta Model physics, the second, with a new soil and land surface scheme and revisions to the surface layer and boundary layer schemes, used the physics package that became operational on 31 January 1996. Improvements in the land surface parameterization and its interaction with the atmosphere are one key to improved summer precipitation forecasts. The new soil thermal model is an improvement over the earlier slab soil model, although the new skin temperature generally now has too large a diurnal cycle (whereas the old surface temperature had too small a diurnal cycle) and is more sensitive to net radiation errors. The nighttime temperature minima are often too low, because of a model underestimate of the downwelling radiation, despite improvements in the coupling of the surface and boundary layer at night. The daytime incoming solar radiation has a substantial high bias in both models, because of some coding errors (which have now been corrected), insufficient atmospheric shortwave absorption, and underestimates of cloud.

The authors explore evaporation before and after a midsummer heavy rain event with the two models. The late 1995 operational model uses a soil moisture bucket physics, with a specified annual-mean fixed field soil moisture climatology, so the surface evaporation responds primarily to the atmospheric forcing. While the surface fluxes in the new model show this strong rain event more dramatically, because its soil moisture comes from the global reanalysis rather than climatology, there remain problems with soil moisture initialization. It appears that a fully continuous Eta data assimilation system (which is under development), likely with more than two soil layers and assimilation of observed hourly precipitation, will be needed to get an adequate soil moisture initialization. Evaporation in the new two-layer soil model falls too much from forecast day 1 to day 2, as the first shallow 10-cm layer dries out (as it also does in the 1995 model with the bucket physics). This appears to be related to the specified low vegetation fraction and the bare soil evaporation model. Although the new boundary layer scheme is better coupled to the surface at night, both versions underestimate entrainment at the top of the mixed layer. The improvement in the surface evaporation resulting from using a climatological green vegetation fraction (derived from satellite data) and a revised bare soil evaporation formulation are shown. These changes were incorporated in a model physics revision in February 1997. An encouraging result from one case study, when it rained in the model, shows that the interaction between the surface, boundary layer, and convection schemes during precipitation is satisfactory, although the model underestimates the impact of cloud cover on the incoming solar radiation.

Corresponding author address: Dr. Alan K. Betts, RR 3, Box 3125, Pittsford, VT 05763.

Abstract

Data from the 1987 summer FIFE experiment for four pairs of days are compared with corresponding 48-h forecasts from two different versions of the Eta Model, both initialized from the NCEP–NCAR (National Centers for Environmental Prediction–National Center for Atmospheric Research) global reanalysis. One used the late 1995 operational Eta Model physics, the second, with a new soil and land surface scheme and revisions to the surface layer and boundary layer schemes, used the physics package that became operational on 31 January 1996. Improvements in the land surface parameterization and its interaction with the atmosphere are one key to improved summer precipitation forecasts. The new soil thermal model is an improvement over the earlier slab soil model, although the new skin temperature generally now has too large a diurnal cycle (whereas the old surface temperature had too small a diurnal cycle) and is more sensitive to net radiation errors. The nighttime temperature minima are often too low, because of a model underestimate of the downwelling radiation, despite improvements in the coupling of the surface and boundary layer at night. The daytime incoming solar radiation has a substantial high bias in both models, because of some coding errors (which have now been corrected), insufficient atmospheric shortwave absorption, and underestimates of cloud.

The authors explore evaporation before and after a midsummer heavy rain event with the two models. The late 1995 operational model uses a soil moisture bucket physics, with a specified annual-mean fixed field soil moisture climatology, so the surface evaporation responds primarily to the atmospheric forcing. While the surface fluxes in the new model show this strong rain event more dramatically, because its soil moisture comes from the global reanalysis rather than climatology, there remain problems with soil moisture initialization. It appears that a fully continuous Eta data assimilation system (which is under development), likely with more than two soil layers and assimilation of observed hourly precipitation, will be needed to get an adequate soil moisture initialization. Evaporation in the new two-layer soil model falls too much from forecast day 1 to day 2, as the first shallow 10-cm layer dries out (as it also does in the 1995 model with the bucket physics). This appears to be related to the specified low vegetation fraction and the bare soil evaporation model. Although the new boundary layer scheme is better coupled to the surface at night, both versions underestimate entrainment at the top of the mixed layer. The improvement in the surface evaporation resulting from using a climatological green vegetation fraction (derived from satellite data) and a revised bare soil evaporation formulation are shown. These changes were incorporated in a model physics revision in February 1997. An encouraging result from one case study, when it rained in the model, shows that the interaction between the surface, boundary layer, and convection schemes during precipitation is satisfactory, although the model underestimates the impact of cloud cover on the incoming solar radiation.

Corresponding author address: Dr. Alan K. Betts, RR 3, Box 3125, Pittsford, VT 05763.

Save
  • Beljaars, A. C. M., 1995: The parameterization of surface fluxes in large-scale models under free convection. Quart. J. Roy. Meteor. Soc.,121, 255–270.

  • ——, P. Viterbo, M. J. Miller, and A. K. Betts, 1996: The anomalous rainfall over the United States during July 1993: Sensitivity to land surface parameterization. Mon. Wea. Rev.,124, 362–383.

  • Berbery, E. H., E. M. Rasmusson, and K. E. Mitchell, 1996: Studies of North American continental-scale hydrology using Eta model forecast products. J. Geophys. Res.,101, 7305–7321.

  • Betts, A. K., 1974: Reply to comment on the paper “Non-precipitating cumulus convection and its parameterization.” Quart. J. Roy. Meteor. Soc.,100, 464–471.

  • ——, 1992: FIFE atmospheric boundary layer budget methods. J. Geophys. Res.,97, 18523–18532.

  • ——, and M. J. Miller, 1986: A new convective adjustment scheme. Part II: Single column tests using GATE-wave, BOMEX, ATEX, and Arctic Airmass data sets. Quart. J. Roy. Meteor. Soc., 112, 693–710.

  • ——, and J. H. Ball, 1994: Budget analysis of FIFE-1987 sonde data. J. Geophys. Res.,99, 3655–3666.

  • ——, and ——, 1995: The FIFE surface diurnal cycle climate. J. Geophys. Res.,100, 25674–25693.

  • ——, and A. G. Barr, 1996: FIFE 1987 Sonde budget revisited. J. Geophys. Res.,101, 23285–23288.

  • ——, and J. H. Ball, 1997: FIFE surface climate and site-average dataset 1987–89. J. Atmos. Sci., in press.

  • ——, ——, and A. C. M. Beljaars, 1993: Comparison between the land surface response of the European Centre model and the FIFE-1987 data. Quart. J. Roy. Meteor. Soc.,119, 975–1001.

  • ——, S.-Y. Hong, and H.-L. Pan, 1996a: Comparison of NCEP/NCAR reanalysis with 1987 FIFE data. Mon. Wea. Rev.,124, 1480–1498.

  • ——, J. H. Ball, A. C. M. Beljaars, M. J. Miller, and P. Viterbo, 1996b: The land surface–atmosphere interaction: A review based on observational and global modeling perspectives. J. Geophys. Res.,101, 7209–7225.

  • Black, T. L., 1994: The new NMC Mesoscale Eta Model: Description and forecast examples. Wea. Forecasting,9, 265–278.

  • ——, J. H. Ward, and Z. I. Janjić, 1989: Tropical storm forecasts using Betts–Miller convection in NMC’s Eta coordinate regional model. Research Activities in Atmospheric and Oceanic Modelling, WCRP 13, 5.38–5.39.[Available from World Meteorological Organization, Case Postale 2300, CH-1211 Geneva, Switzerland.].

  • Cess, R. D., and Coauthors, 1995: Absorption of solar radiation by clouds: Observations versus models. Science,267, 496–499.

  • Chen, F., K. Mitchell, J. Schaake, Y. Xue, H.-L. Pan, V. Koren, Q. Duan, and A. Betts, 1996: Modeling of land-surface evaporation by four schemes and comparison with FIFE observations. J. Geophys. Res.,101, 7251–7268.

  • ——, Z. I. Janjić, and K. Mitchell, 1997: Impact of atmospheric surface layer parameterization in the new land-surface scheme of the NCEP Mesoscale Eta Numerical Model. Bound.-Layer Meteor., in press.

  • Clapp, R. B., and G. M. Hornberger, 1978: Empirical equations for some soil hydraulic properties. Water Resour. Res.,14, 601–604.

  • Cosby, B. J., G. M. Hornberger, R. B. Clapp, and T. R. Ginn, 1984: A statistical exploration of the relationships of soil moisture characteristics to the physical properties of soils. Water Resour. Res.,20, 682–690.

  • Gerrity, J. P., Jr., T. L. Black, and R. E. Treadon, 1994: The numerical solution of the Mellor–Yamada level-2.5 turbulent kinetic energy equation in the Eta model. Mon. Wea. Rev.,122, 1640–1646.

  • Gutman, G., and A. Ignatov, 1996: Derivation of green vegetation fraction from NOAA/AVHRR for use in numerical weather prediction models. [Available from NOAA/NESDIS, Office of Research and Applications, Washington, DC 20233.].

  • Janjić, Z. I., 1990: The step-mountain coordinate: Physical package. Mon. Wea. Rev.,118, 1429–1443.

  • ——, 1994: The step-mountain eta coordinate model: Further developments of the convection, viscous sublayer, and turbulence closure schemes. Mon. Wea. Rev.,122, 927–945.

  • ——, 1996a: The Mellor–Yamada level 2.5 turbulence closure scheme in the NCEP Eta Model. Research Activities in Atmospheric and Oceanic Modeling, WMO, Geneva, CAS/JSC WGNE, 4.14–4.15. [Available from World Meteorological Organization, Case Postale 2300, CH-1211 Geneva, Switzerland.].

  • ——, 1996b: The surface layer parameterization in the NCEP Eta Model. Research Activities in Atmospheric and Oceanic Modeling, WMO, Geneva, CAS/JSC WGNE, 4.16–4.17. [Available from World Meteorological Organization, Case Postale 2300, CH-1211 Geneva, Switzerland.].

  • ——, 1996c: The surface layer in the NCEP Eta Model. Preprints, 11th Conf. on Numerical Weather Prediction, Norfolk, VA, Amer. Meteor. Soc., 354–356.

  • Kalnay, E., and Coauthors, 1996: The NCEP/NCAR 40-year reanalysis project. Bull. Amer. Meteor. Soc.,77, 437–471.

  • Kim, J. and S. B. Verma, 1990: Components of surface energy balance in a temperate grassland ecosystem. Bound.-Layer Meteor.,51, 401–417.

  • Lanzić, L., 1990: Forecasts of AMEX tropical cyclones with step-mountain model. Aust. Meteor. Mag.,38, 207–216.

  • ——, 1993a: Eta model forecasts of tropical cyclones from Australian Monsoon Experiment: Dynamical adjustment of initial conditions. Meteor. Atmos. Phys.,52, 101–111.

  • ——, 1993b: Eta model forecasts of tropical cyclones from Australian Monsoon Experiment: The model sensitivity. Meteor. Atmos. Phys.,52, 113–127.

  • Leese, J. A., 1993: Implementation plan for the GEWEX Continental-Scale International Project (GCIP). Int. GEWEX Project Office 6, 148 pp. [Available from GPO, 1100 Wayne Ave, Suite 1225, Silver Springs, MD 20910.].

  • ——, 1994: Major activities plan for 1995, 1996 and outlook for 1997 for the GCIP. IGP0 12, 166 pp. [Available from IGPO, 1100 Wayne Ave, Suite 1225, Silver Springs, MD 20910.].

  • Lobocki, L., 1993: A procedure for the derivation of surface-layer bulk relationships from simplified second-order closure models. J. Appl. Meteor.,32, 126–138.

  • Mahfouf, J. F., and J. Noilhan, 1991: Comparative study of various formulations of evaporation from bare soil using in situ data. J. Appl. Meteor.,30, 1354–1365.

  • Mahrt, L., and H.-L. Pan, 1984: A two layer model for soil hydrology. Bound.-Layer Meteor.,29, 1–20.

  • Manabe, S., 1969: Climate and the ocean circulation: I. The atmospheric circulation and the hydrology of the earth’s surface. Mon. Wea. Rev.,97, 739–774.

  • Mellor, G. L., and T. Yamada, 1982: Development of a turbulence closure model for geophysical fluid problems. Rev. Geophys. Space Phys.,20, 851–875.

  • Mesinger, F., 1993: Sensitivity of the definition of a cold front to the parameterization of turbulent fluxes in the NMC’s Eta model. Research Activities in Atmospheric and Oceanic Modeling 18, 4.25–4.38. [Available from World Meteorological Organization, Case Postale 2300, CH-1211 Geneva, Switzerland.].

  • ——, 1996: Improvements in quantitative precipitation forecasts with the Eta regional model at the U.S. National Centers for Environmental Prediction. The 48-km upgrade. Bull. Amer. Meteor. Soc.,77, 2637–2649.

  • ——, and T. L. Black, 1992: On the impact on forecast accuracy of the step mountain (Eta) vs. sigma coordinate. Meteor. Atmos. Phys.,50, 47–60.

  • ——, Z. I. Janjić, S. Nicković, D. Gavrilov, and D. G. Deaven, 1988: The step-mountain coordinate: Model description and performance for cases of Alpine lee cyclogenesis and for a case of an Appalechian redevelopment. Mon. Wea. Rev.,116, 1493–1518.

  • ——, T. L. Black, and M. E. Baldwin, 1997: Impact of resolution and of the eta coordinate on skill of the Eta Model precipitation forecasts. Numerical Methods in Atmospheric Modelling. The Andre J. Robert Memorial Volume, C. A. Lin, R. Laprise, and H. Ritchie, Eds., Canadian Meteorological and Oceanic Society, 399–423.

  • Pan, H. L., 1990: A simple parameterization scheme of evapotranspiration over land for the NMC medium-range forecast model. Mon. Wea. Rev.,118, 2500–2512.

  • ——, and L. Mahrt, 1987: Interaction between soil hydrology and boundary-layer development. Bound.-Layer Meteor.,38, 185–202.

  • Paulson, C. A., 1970: The mathematical representation of wind speed and temperature profiles in the unstable atmospheric surface layer. J. Appl. Meteor., 9, 857–861.

  • Rogers, E., D. G. Deaven, and G. J. DiMego, 1995: The regional analysis system for the operational “early” eta model: Original 80-km configuration and recent changes. Wea. Forecasting,10, 810–825.

  • Smith, E. A., W. L. Crosson, and B. D. Tanner, 1992a: Estimation of surface heat and moisture fluxes over a prairie grassland. 1: In situ energy budget measurements incorporating a cooled mirror dew point hygrometer. J. Geophys. Res.,97, 18557–18582.

  • ——, and Coauthors, 1992b: Area-averaged surface fluxes and their time-space variability over the FIFE experimental domain. J. Geophys. Res.,97, 18599–18622.

  • Strebel, D. E., D. R. Landis, K. F. Huemmrich, and B. W. Meeson, 1994: Collected data of the First ISLSCP Field Experiment. Surface Observations and Nonimage Data Sets, Vol. 1, NASA Goddard Space Flight Center, CD-ROM.

  • Sugita, M., and W. Brutsaert, 1990a: Wind velocity measurements in the neutral boundary layer above hilly prairie. J. Geophys. Res.,95, 7617–7624.

  • ——, and ——, 1990b: How similar are temperature and humidity profiles in the unstable boundary layer? J. Appl. Meteor.,29, 489–497.

  • Viterbo, P., and A. C. M. Beljaars, 1995: A new land surface parameterization scheme in the ECMWF model and its validation. J. Climate,8, 2716–2748.

  • Ward, D. M., 1995: Comparison of the surface solar radiation budget derived from satellite data with that simulated by the NCAR CCM2. J. Climate,8, 2824–2842.

  • Zilitinkevich, S. S., 1995: Non-local turbulent transport: Pollution dispersion aspects of coherent structure of convective flows. Air Pollution III, H. Power, N. Moussiopoulos, and C. A. Brebbia, Eds., Air Pollution Theory and Simulation, Vol. 1, Computational Mechanics Publications, 53–60.

  • Zobler, L., 1986: A world soil file for global climate modeling. NASA Tech. Memo. 87802, 33 pp. [Available from NASA Goddard Space Flight Center, Institute for Space Studies, 2800 Broadway, New York, NY 10025.].

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 3199 2649 121
PDF Downloads 480 114 8