Calculation of surface temperature and surface fluxes in the GLAS GOM

Because the GLAS model's surface fluxes of sensible and latent heat exhibit strong 2 delta t oscillations at the individual grid points as well as in the zonal hemispheric averages and because a basic weakness of the GLAS model lower evaporation over oceans and higher evaporation over land in a typical monthly simulation, the GLAS model PBL parameterization was changed to calculate the mixed layer temperature gradient by solution of a quadratic equation for a stable PBL and by a curve fit relation for an unstable PBL. The new fluxes without any 2 delta t oscillation. Also, the geographical distributions of the surface fluxes are improved. The parameterization presented is incorporated into the new GLAS climate model. Some results which compare the evaporation over land and ocean between old and new calculations are appended

A review of recent research on improvement of physical parameterizations in the GLA GCM

A systematic assessment of the effect of a series of improvements in physical parameterizations of the Goddard Laboratory for Atmospheres (GLA) general circulation model (GCM) are summarized. The implementation of the Simple Biosphere Model (SiB) in the GCM is followed by a comparison of SiB GCM simulations with that of the earlier slab soil hydrology GCM (SSH-GCM) simulations. In the Sahelian context, the biogeophysical component of desertification was analyzed for SiB-GCM simulations. Cumulus parameterization is found to be the primary determinant of the organization of the simulated tropical rainfall of the GLA GCM using Arakawa-Schubert cumulus parameterization. A comparison of model simulations with station data revealed excessive shortwave radiation accompanied by excessive drying and heating to the land. The perpetual July simulations with and without interactive soil moisture shows that 30 to 40 day oscillations may be a natural mode of the simulated earth atmosphere system

A study of the influence of soil moisture on future precipitation

Forty years of precipitation and surface temperature data observed over 261 Local Climatic Data (LCD) stations in the Continental United States was utilized in a ground hydrology model to yield soil moisture time series at each station. A month-by-month soil moisture dataset was constructed for each year. The monthly precipitation was correlated with antecedent monthly precipitation, soil moisture and vapotranspiration separately. The maximum positive correlation is found to be in the drought prone western Great Plains region during the latter part of summer. There is also some negative correlation in coastal regions. The correlations between soil moisture and precipitation particularly in the latter part of summer, suggest that large scale droughts over extended periods may be partially maintained by the feedback influence of soil moisture on rainfall. In many other regions the lack of positive correlation shows that there is no simple answer such as higher land-surface evapotranspiration leads to more precipitation, and points out the complexity of the influence of soil moisture on the ensuring precipitation

Soil Surface Runoff Scheme for Improving Land-Hydrology and Surface Fluxes in Simple SiB (SSiB)

Evapotranspiration on land is hard to measure and difficult to simulate. On the scale of a GCM grid, there is large subgrid-scale variability of orography, soil moisture, and vegetation. Our hope is to be able to tune the biophysical constants of vegetation and soil parameters to get the most realistic space-averaged diurnal cycle of evaporation and its climatology. Field experiments such as First ISLSCP Field Experiment (FIFE), Boreal Ecosystem-Atmosphere Study (BOREAS), and LBA help a great deal in improving our evapotranspiration schemes. However, these improvements have to be matched with, and coupled to, consistent improvement in land-hydrology; otherwise, the runoff problems will intrinsically reflect on the soil moisture and evapotranspiration errors. Indeed, a realistic runoff simulation also ensures a reasonable evapotranspiration simulation provided the precipitation forcing is reliable. We have been working on all of the above problems to improve the simulated hydrologic cycle. Through our participation in the evaluation and intercomparison of land-models under the behest of Global Soil Wetness Project (GSWP), we identified a few problems with Simple SiB (SSIB; Xue et al., 1991) hydrology in regions of significant snowmelt. Sud and Mocko (1999) show that inclusion of a separate snowpack model, with its own energy budget and fluxes with the atmosphere aloft and soil beneath, helps to ameliorate some of the deficiencies of delayed snowmelt and excessive spring season runoff. Thus, much more realistic timing of melt water generation was simulated with the new snowpack model in the subsequent GSWP re-evaluations using 2 years of ISLSCP Initiative I forcing data for 1987 and 1988. However, we noted an overcorrection of the low meltwater infiltration of SSiB. While the improvement in snowmelt timing was found everywhere, the snowmelt infiltration has became excessive in some regions, e.g., Lena river basin. This leads to much reduced runoff in many basins as compared to observations. We believe this is a consequence of neglect of the influence of subgrid-scale variations in orography that affects the production of surface runoff

Snowmelt and Infiltration Deficiencies of SSiB and Their Resolution with a New Snow-Physics Scheme

A two-year 1987-1988 integration of SSiB forced with ISLSCP Initiative I surface data (as part of the Global Soil Wetness Project, GSWP, evaluation and intercomparison) produced generally realistic land surface fluxes and hydrology. Nevertheless, the evaluation also helped to identify some of the deficiencies of the current version of the Simplified Simple Biosphere (SSiB) model. The simulated snowmelt was delayed in most regions, along with excessive runoff and lack of an spring soil moisture recharge. The SSIB model had previously been noted to have a problem producing accurate soil moisture as compared to observations in the Russian snowmelt region. Similarly, various GSWP implementations of SSIB found deficiencies in this region of the simulated soil moisture and runoff as compared to other non-SSiB land-surface models (LSMs). The origin of these deficiencies was: 1) excessive cooling of the snow and ground, and 2) deep frozen soil disallowing snowmelt infiltration. The problem was most severe in regions that experience very cold winters. In SSiB, snow was treated as a unified layer with the first soil layer, causing soil and snow to cool together in the winter months, as opposed to snow cover acting as an insulator. In the spring season, a large amount of heat was required to thaw a hard frozen snow plus deep soil layers, delaying snowmelt and causing meltwater to become runoff over the frozen soil rather than infiltrate into it

Sensitivity of Boreal-Summer Circulation and Precipitation to Atmospheric Aerosols in Selected Regions

Aerosol perturbations over selected land regions are imposed in Version-4 of the Goddard Earth Observing System (GEOS-4) general circulation model (GCM) to assess the influence of increasing aerosol concentrations on regional circulation patterns and precipitation in four selected regions: India, Africa, and North and South America. Part 1 of this paper addresses the responses to aerosol perturbations in India and Africa. This paper presents the same for aerosol perturbations over the Americas. GEOS-4 is forced with prescribed aerosols based on climatological data, which interact with clouds using a prognostic scheme for cloud microphysics including aerosol nucleation of water and ice cloud hydrometeors. In clear-sky conditions the aerosols interact with radiation. Thus the model includes comprehensive physics describing the aerosol direct and indirect effects on climate (hereafter ADE and AIE respectively). Each simulation is started from analyzed initial conditions for 1 May and was integrated through June-July-August of each of the six years: 1982 1987 to provide a 6-ensemble set. Results are presented for the difference between simulations with double the climatological aerosol concentration and one-half the climatological aerosol concentration for three experiments: two where the ADE and AIE are applied separately and one in which both the ADE and AIE are applied. The ADE and AIE both yield reductions in net radiation at the top of the atmosphere and surface while the direct absorption of shortwave radiation contributes a net radiative heating in the atmosphere. A large net heating of the atmosphere is also apparent over the subtropical North Atlantic Ocean that is attributable to the large aerosol perturbation imposed over Africa. This atmospheric warming and the depression of the surface pressure over North America contribute to a northward shift of the inter-Tropical Convergence Zone over northern America, an increase in precipitation over Central America and the Caribbean, and an enhancement of convergence in the North American monsoon region

Mechanisms Regulating Deep Moist Convection and Sea-Surface Temperatures of the Tropics

Despite numerous previous studies, two relationships between deep convection and the sea-surface temperature (SST) of the tropics remain unclear. The first is the cause for the sudden emergence of deep convection at about 28 deg SST, and the second is its proximity to the highest observed SST of about 30 C. Our analysis provides a rational explanation for both by utilizing the Improved Meteorological (IMET) buoy data together with radar rainfall retrievals and atmospheric soundings provided by the Tropical Ocean Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA-COARE). The explanation relies on the basic principles of moist convection as enunciated in the Arakawa-Schubert cumulus parameterization. Our analysis shows that an SST range of 28-29 C is necessary for "charging" the atmospheric boundary layer with sufficient moist static energy that can enable the towering convection to reach up to the 200 hPa level. In the IMET buoy data, the changes in surface energy fluxes associated with different rainfall amounts show that the deep convection not only reduces the solar flux into the ocean with a thick cloud cover, but it also generates downdrafts which bring significantly cooler and drier air into the boundary-layer thereby augmenting oceanic cooling by increased sensible and latent heat fluxes. In this way, the ocean seasaws between a net energy absorber for non-raining and a net energy supplier for deep-convective raining conditions. These processes produce a thermostat-like control of the SST. The data also shows that convection over the warm pool is modulated by dynamical influences of large-scale circulation embodying tropical easterly waves (with a 5-day period) and MJOs (with 40-day period); however, the quasi-permanent feature of the vertical profile of moist static energy, which is primarily maintained by the large-scale circulation and thermodynamical forcings, is vital for both the 28 C SST for deep convection and its upper limit at about 30 C

Influences of Local Sea-Surface Temperatures and Large-scale Dynamics on Monthly Precipitation Inferred from Two 10-year GCM-Simulations

Two parallel sets of 10-year long: January 1, 1982 to December 31, 1991, simulations were made with the finite volume General Circulation Model (fvGCM) in which the model integrations were forced with prescribed sea-surface temperature fields (SSTs) available as two separate SST-datasets. One dataset contained naturally varying monthly SSTs for the chosen period, and the oth& had the 12-monthly mean SSTs for the same period. Plots of evaporation, precipitation, and atmosphere-column moisture convergence, binned by l C SST intervals show that except for the tropics, the precipitation is more strongly constrained by large-scale dynamics as opposed to local SST. Binning data by SST naturally provided an ensemble average of data contributed from disparate locations with same SST; such averages could be expected to mitigate all location related influences. However, the plots revealed: i) evaporation, vertical velocity, and precipitation are very robust and remarkably similar for each of the two simulations and even for the data from 1987-ENSO-year simulation; ii) while the evaporation increased monotonically with SST up to about 27 C, the precipitation did not; iii) precipitation correlated much better with the column vertical velocity as opposed to SST suggesting that the influence of dynamical circulation including non-local SSTs is stronger than local-SSTs. The precipitation fields were doubly binned with respect to SST and boundary-layer mass and/or moisture convergence. The analysis discerned the rate of change of precipitation with local SST as a sum of partial derivative of precipitation with local SST plus partial derivative of precipitation with boundary layer moisture convergence multiplied by the rate of change of boundary-layer moisture convergence with SST (see Eqn. 3 of Section 4.5). This analysis is mathematically rigorous as well as provides a quantitative measure of the influence of local SST on the local precipitation. The results were recast to examine the dependence of local rainfall on local SSTs; it was discernible only in the tropics. Our methodology can be used for computing relationship between any forcing function and its effect(s) on a chosen field

Design of a global soil moisture initialization procedure for the simple biosphere model

Global soil moisture and land-surface evapotranspiration fields are computed using an analysis scheme based on the Simple Biosphere (SiB) soil-vegetation-atmosphere interaction model. The scheme is driven with observed precipitation, and potential evapotranspiration, where the potential evapotranspiration is computed following the surface air temperature-potential evapotranspiration regression of Thomthwaite (1948). The observed surface air temperature is corrected to reflect potential (zero soil moisture stress) conditions by letting the ratio of actual transpiration to potential transpiration be a function of normalized difference vegetation index (NDVI). Soil moisture, evapotranspiration, and runoff data are generated on a daily basis for a 10-year period, January 1979 through December 1988, using observed precipitation gridded at a 4 deg by 5 deg resolution