Engineering topics

Walker, G. K.

Publications and source records attributed to Walker, G. K..

At least 19 records

NASA NTRS2015

Interactive Vegetation Phenology, Soil Moisture, and Monthly Temperature Forecasts

The time scales that characterize the variations of vegetation phenology are generally much longer than those that characterize atmospheric processes. The explicit modeling of phenological processes in an atmospheric forecast system thus has the potential to provide skill to subseasonal or seasonal forecasts. We examine this possibility here using a forecast system fitted with a dynamic vegetation phenology model. We perform three experiments, each consisting of 128 independent warm-season monthly forecasts: 1) an experiment in which both soil moisture states and carbon states (e.g., those determining leaf area index) are initialized realistically, 2) an experiment in which the carbon states are prescribed to climatology throughout the forecasts, and 3) an experiment in which both the carbon and soil moisture states are prescribed to climatology throughout the forecasts. Evaluating the monthly forecasts of air temperature in each ensemble against observations, as well as quantifying the inherent predictability of temperature within each ensemble, shows that dynamic phenology can indeed contribute positively to subseasonal forecasts, though only to a small extent, with an impact dwarfed by that of soil moisture.

Seasonal Forecasting↗

NASA NTRS2009

A GCM Study of Responses of the Atmospheric Water Cycle of West Africa and the Atlantic to Saharan Dust Radiative Forcing

The responses of the atmospheric water cycle and climate of West Africa and the Atlantic to radiative forcing of Saharan dust are studied using the NASA finite volume general circulation model (fvGCM), coupled to a mixed layer ocean. We find evidence of an "elevated heat pump" (EHP) mechanism that underlines the responses of the atmospheric water cycle to dust forcing as follow. During the boreal summerr, as a result of large-scale atmospheric feedback triggered by absorbing dust aerosols, rainfall and cloudiness are ehanIed over the West Africa/Eastern Atlantic ITCZ, and suppressed over the West Atlantic and Caribbean region. Shortwave radiation absorption by dust warms the atmosphere and cools the surface, while longwave has the opposite response. The elevated dust layer warms the air over West Africa and the eastern Atlantic. As the warm air rises, it spawns a large-scale onshore flow carrying the moist air from the eastern Atlantic and the Gulf of Guinea. The onshore flow in turn enhances the deep convection over West Africa land, and the eastern Atlantic. The condensation heating associated with the ensuing deep convection drives and maintains an anomalous large-scale east-west overturning circulation with rising motion over West Africa/eastern Atlantic, and sinking motion over the Caribbean region. The response also includes a strengthening of the West African monsoon, manifested in a northward shift of the West Africa precipitation over land, increased low-level westerlies flow over West Africa at the southern edge of the dust layer, and a near surface westerly jet underneath the dust layer overr the Sahara. The dust radiative forcing also leads to significant changes in surface energy fluxes, resulting in cooling of the West African land and the eastern Atlantic, and warming in the West Atlantic and Caribbean. The EHP effect is most effective for moderate to highly absorbing dusts, and becomes minimized for reflecting dust with single scattering albedo at0.95 or higher.

Lau, K. M.↗

NASA NTRS2008

Effects of Doubled CO2 on Tropical Sea-Surface Temperature (SSTs) for Onset of Deep Convection and Maximum SST-GCM Simulations Based Inferences

A primary concern of CO2-induced warming is the associated rise of tropical (10S-10N) seasurface temperatures (SSTs). GISS Model-E was used to produce two sets of simulations-one with the present-day and one with doubled CO2 in the atmosphere. The intrinsic usefulness of model guidance in the tropics was confirmed when the model simulated realistic convective coupling between SSTs and atmospheric soundings and that the simulated-data correlations between SSTs and 300 hPa moiststatic energies were found to be similar to the observed. Model predicted SST limits: (i) one for the onset of deep convection and (ii) one for maximum SST, increased in the doubled C02 case. Changes in cloud heights, cloud frequencies, and cloud mass-fractions showed that convective-cloud changes increased the SSTs, while warmer mixed-layer of the doubled CO2 contained approximately 10% more water vapor; clearly that would be conducive to more intense storms and hurricanes.

Sud, Y. C.↗

NASA NTRS2007

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.

Sud, Y. C.↗

NASA NTRS2004

Effects of Cloud-Microphysics on Tropical Atmospheric Hydrologic Processes in the GEOS GCM

The sensitivity of tropical atmospheric hydrologic processes to cloud-microphysics is investigated using the NASA GEOS GCM. Results show that a faster autoconversion - rate produces more warm rain and less clouds at all levels. Fewer clouds enhances longwave cooling and reduces shortwave heating in the upper troposphere, while more warm rain produces increased condensation heating in the lower troposphere. This vertical heating differential destablizes the tropical atmosphere, producing a positive feedback resulting in more rain over the tropics. The feedback is maintained via a two-cell secondary circulation. The lower cell is capped by horizontal divergence and maximum cloud detrainment near the melting/freezing, with rising motion in the warm rain region connected to descending motion in the cold rain region. The upper cell is found above the freezing/melting level, with longwave-induced subsidence in the warm rain and dry regions, coupled to forced ascent in the deep convection region. The tropical large scale circulation is found to be very sensitive to the radiative-dynamic effects induced by changes in autoconversion rate. Reduced cloud-radiation processes feedback due to a faster autoconversion rate results in intermittent but more energetic eastward propagating Madden and Julian Oscillations (MJO). Conversely,-a slower autconversion rate, with increased cloud radiation produces MJO's with more realistic westward propagating transients, resembling a supercloud cluster structure. Results suggests that warm rain and associated low and mid level clouds, i.e., cumulus congestus, may play a critical role in regulating the time-intervals of deep convections and hence the fundamental time scales of the MJO.

Lau, K. M.↗

NASA NTRS2004

Influence of Sub-grid-Scale Isentropic Transports on McRAS Evaluations using ARM-CART SCM Datasets

In GCM-physics evaluations with the currently available ARM-CART SCM datasets, McRAS produced very similar character of near surface errors of simulated temperature and humidity containing typically warm and moist biases near the surface and cold and dry biases aloft. We argued it must have a common cause presumably rooted in the model physics. Lack of vertical adjustment of horizontal transport was thought to be a plausible source. Clearly, debarring such a freedom would force the incoming air to diffuse into the grid-cell which would naturally bias the surface air to become warm and moist while the upper air becomes cold and dry, a characteristic feature of McRAS biases. Since, the errors were significantly larger in the two winter cases that contain potentially more intense episodes of cold and warm advective transports, it further reaffirmed our argument and provided additional motivation to introduce the corrections. When the horizontal advective transports were suitably modified to allow rising and/or sinking following isentropic pathways of subgrid scale motions, the outcome was to cool and dry (or warm and moisten) the lower (or upper) levels. Ever, crude approximations invoking such a correction reduced the temperature and humidity biases considerably. The tests were performed on all the available ARM-CART SCM cases with consistent outcome. With the isentropic corrections implemented through two different numerical approximations, virtually similar benefits were derived further confirming the robustness of our inferences. These results suggest the need for insentropic advective transport adjustment in a GCM due to subgrid scale motions.

Sud, Y. C.↗

NASA NTRS2003

Influence of Ice-phase of Hydrometeors on Moist-Convection

Climate models often ignore the influence of ice-phase physics (IPP) of hydrometeors as a second order effect. This has also been true for McRAS (Microphysics of clouds with Relaxed Arakawa Schubert Scheme) developed by the authors. Recognizing that the temperature sounding is critical for moist-convection, and, that IPP would modify it, we investigated the influence of introducing IPP into McRAS coupled to FvGCM (finite volume General Circulation Model with NCAR physics). We analyzed three 3-yr long simulations; the first called Control Case, CC and had no IPP; the other two called Experiments El and E2 had IPP introduced with two different in-cloud freezing assumptions. Simulation El assumed that all hydrometeors remain liquid in the updraft and freeze upon detrainment. Simulation E2 invoked the in-cloud freezing of new condensate generated at subfreezing temperatures in the updraft while old cloud water continued to ascend as liquid. Upon detrainment, this cloud water also froze like in E1. With these assumptions, about 50% of hydrometeors froze in the tower and the rest froze in the anvil. However, in both El and E2, the frozen hydrometeors melted during fall at the first encounter of above freezing ambient temperature. Comparative analysis revealed that El simulated far more mid-level and far less deep clouds while E2 had modified deep and more mid-level clouds as compared to CC along with some major changes around the melt-level. We infer that IPP produced a more realistic response in E2. At the basic level, the results show that ice-phase processes influence convective detrainment at mid- and deep levels in accord with TOGAGOARE observations. The results suggest that IPP can help to mitigate less-than-observed mid-level and over-abundance of deep convective clouds in McRAS.

Sud, Y. C.↗

NASA NTRS2002

New Concepts for Refinement of Cumulus Parameterization in GCM's the Arakawa-Schubert Framework

Several state-of-the-art models including the one employed in this study use the Arakawa-Schubert framework for moist convection, and Sundqvist formulation of stratiform. clouds, for moist physics, in-cloud condensation, and precipitation. Despite a variety of cloud parameterization methodologies developed by several modelers including the authors, most of the parameterized cloud-models have similar deficiencies. These consist of: (a) not enough shallow clouds, (b) too many deep clouds; (c) several layers of clouds in a vertically demoralized model as opposed to only a few levels of observed clouds, and (d) higher than normal incidence of double ITCZ (Inter-tropical Convergence Zone). Even after several upgrades consisting of a sophisticated cloud-microphysics and sub-grid scale orographic precipitation into the Data Assimilation Office (DAO)'s atmospheric model (called GEOS-2 GCM) at two different resolutions, we found that the above deficiencies remained persistent. The two empirical solutions often used to counter the aforestated deficiencies consist of a) diffusion of moisture and heat within the lower troposphere to artificially force the shallow clouds; and b) arbitrarily invoke evaporation of in-cloud water for low-level clouds. Even though helpful, these implementations lack a strong physical rationale. Our research shows that two missing physical conditions can ameliorate the aforestated cloud-parameterization deficiencies. First, requiring an ascending cloud airmass to be saturated at its starting point will not only make the cloud instantly buoyant all through its ascent, but also provide the essential work function (buoyancy energy) that would promote more shallow clouds. Second, we argue that training clouds that are unstable to a finite vertical displacement, even if neutrally buoyant in their ambient environment, must continue to rise and entrain causing evaporation of in-cloud water. These concepts have not been invoked in any of the cloud parameterization schemes so far. We introduced them into the DAO-GEOS-2 GCM with McRAS (Microphysics of Clouds with Relaxed Arakawa-Schubert Scheme).

Sud, Y. C.↗

NASA NTRS2001

Simulation Study of Effects of Solar Irradiance and Sea Surface Temperature on Monsoons and Global Circulation

A recent version of the GEOS 2 GCM was used to isolate the roles of the annual cycles of solar irradiation and/or sea-surface temperatures (SSTs) on the simulated circulation and rainfall. Four 4-year long integrations were generated with the GCM. The first integration, called Control Case, used daily-interpolated SSTs from a 30 year monthly SST climatology that was obtained from the analyzed SST-data, while the solar irradiation at the top of the atmosphere was calculated normally at hourly intervals. The next two cases prescribed the SSTs or the incoming solar irradiance at the top of the atmosphere at their annual mean values, respectively while everything else was kept the same as in the Control Case. In this way the influence of the annual cycles of both external forcings was isolated.

Sud, Y. C.↗

NASA NTRS2001

Importance of the Annual Cycles of SST and Solar Irradiance for Circulation and Rainfall: A Climate Model Simulation Study

Annual cycle of climate and precipitation is related to annual cycle of sunshine and sea-surface temperatures. Understanding its behavior is important for the welfare of humans worldwide. For example, failure of Asian monsoons can cause widespread famine and grave economic disaster in the subtropical regions. For centuries meteorologists have struggled to understand the importance of the summer sunshine and associated heating and the annual cycle of sea-surface temperatures (SSTs) on rainfall in the subtropics. Because the solar income is pretty steady from year to year, while SSTs depict large interannual variability as consequence of the variability of ocean dynamics, the influence of SSTs on the monsoons are better understood through observational and modeling studies whereas the relationship of annual rainfall to sunshine remains elusive. However, using NASA's state of the art climate model(s) that can generate realistic climate in a computer simulation, one can answer such questions. We asked the question: if there was no annual cycle of the sunshine (and its associated land-heating) or the SST and its associated influence on global circulation, what will happen to the annual cycle of monsoon rains? By comparing the simulation of a 4-year integration of a baseline Control case with two parallel anomaly experiments: 1) with annual mean solar and 2) with annual mean sea-surface temperatures, we were able to draw the following conclusions: (1) Tropical convergence zone and rainfall which moves with the Sun into the northern and southern hemispheres, specifically over the Indian, African, South American and Australian regions, is strongly modulated by the annual cycles of SSTs as well as solar forcings. The influence of the annual cycle of solar heating over land, however, is much stronger than the corresponding SST influence for almost all regions, particularly the subtropics; (2) The seasonal circulation patterns over the vast land-masses of the Northern Hemisphere at mid and high latitudes also get strongly influenced by the annual cycles of solar heating. The SST influence is largely limited to the oceanic regions of these latitudes; (3) The annual mode of precipitation over Amazonia has an equatorial regime revealing a maxima in the month of March associated with SST, and another maxima in the month of January associated with the solar annual cycles, respectively. The baseline simulation, which has both annual cycles, depicts both annual modes and its rainfall is virtually equal to the sum of those two modes; (4) Rainfall over Sahelian-Africa is significantly reduced (increased) in simulations lacking (invoking) solar irradiation with (without) the annual cycle. In fact, the dominant influence of solar irradiation emerges in almost all monsoonal-land regions: India, Southeast Asia, as well as Australia. The only exception is the Continental United States, where solar annual cycle shows only a relatively minor influence on the annual mode of rainfall.

Sud, Yogesh C.↗

NASA NTRS2000

Influence of the Biosphere on Precipitation: July 1995 Studies with the ARM-CART Data

Ensemble sets of simulation experiments were conducted with a single column model (SCM) using the Goddard GEOS II GCM physics containing a recent version of the Cumulus Scheme (McRAS) and a biosphere based land-fluxes scheme (SSiB). The study used the 18 July to 5 August 1995 ARM-CART (Atmospheric Radiation Measurement-Cloud Atmospheric Radiation Test-bed) data, which was collected at the ARM-CART site in the mid-western United States and analyzed for single column modeling (SCM) studies. The new findings affirm the earlier findings that the vegetation, which increases the solar energy absorption at the surface together with soil and soil-moisture dependent processes, which modulate the surface, fluxes (particularly evapotranspiration) together help to increase the local rainfall. In addition, the results also show that for the particular study period roughly 50% of the increased evaporation over the ARM-CART site would be converted into rainfall with the Column, while the remainder would be advected out to the large-scale. Notwithstanding the limitations of only one-way interaction (i.e., the large-scale influencing the regional physics and not vice versa), the current SCM simulations show a very robust relationship. The evaporation-precipitation relationship turns out to be independent of the soil types, and soil moisture; however, it is weakly dependent on the vegetation cover because of its surface-albedo effect. Clearly, these inferences are prone to weaknesses of the SCM physics, the assumptions of the large-scale being unaffected by gridscale (SCM-scale) changes in moist processes, and other limitations of the evaluation procedures.

Sud, Y. C.↗

NASA NTRS1998

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.

Sud, Y. C.↗

NASA NTRS1998

Design, Evaluation and GCM-Performance of a New Parameterization for Microphysics of Clouds with Relaxed Arakawa-Schubert Scheme (McRas)

A prognostic cloud scheme named McRAS (Microphysics of clouds with Relaxed Arakawa-Schubert Scheme) was developed with the aim of improving cloud-microphysics, and cloud-radiation interactions in GCMs. McRAS distinguishes convective, stratiform, and boundary-layer clouds. The convective clouds merge into stratiform clouds on an hourly time-scale, while the boundary-layer clouds do so instantly. The cloud condensate transforms into precipitation following the auto-conversion relations of Sundqvist that contain a parametric adaptation for the Bergeron-Findeisen process of ice crystal growth and collection of cloud condensate by precipitation. All clouds convect, advect, and diffuse both horizontally and vertically with a fully active cloud-microphysics throughout its life-cycle, while the optical properties of clouds are derived from the statistical distribution of hydrometeors and idealized cloud geometry. An evaluation of McRAS in a single column model (SCM) with the GATE Phase III data has shown that McRAS can simulate the observed temperature, humidity, and precipitation without discernible systematic errors. An evaluation with the ARM-CART SCM data in a cloud model intercomparison exercise shows reasonable but not an outstanding accurate simulation. Such a discrepancy is common to almost all models and is related, in part, to the input data quality. McRAS was implemented in the GEOS II GCM. A 50 month integration that was initialized with the ECMWF analysis of observations for January 1, 1987 and forced with the observed sea-surface temperatures and sea-ice distribution and vegetation properties (biomes, and soils), with prognostic soil moisture, snow-cover, and hydrology showed a very realistic simulation of cloud process, incloud water and ice, and cloud-radiative forcing (CRF). The simulated ITCZ showed a realistic time-mean structure and seasonal cycle, while the simulated CRF showed sensitivity to vertical distribution of cloud water which can be easily altered by the choice of time constant and incloud critical cloud water amount regulators for auto-conversion. The CRF and its feedbacks also have a profound effect on the ITCZ. Even though somewhat weaker than observed, the McRAS-GCM simulation produces robust 30-60 day oscillations in the 200 hPa velocity potential. Two ensembles of 4-summer (July, August, September) simulations, one each for 1987 and 1988 show that the McRAS-GCM simulates realistic and statistically significant precipitation differences over India, Central America, and tropical Africa. Several seasonal simulations were performed with McRAS-GEOS II GCM for the summer (June-July- August) and winter (December-January-February) periods to determine how the simulated clouds and CRFs would be affected by: i) advection of clouds; ii) cloud top entrainment instability, iii) cloud water inhomogeneity correction, and (iv) cloud production and dissipation in different cloud-processes. The results show that each of these processes contributes to the simulated cloud-fraction and CRF.

Sud, Y. C.↗

NASA NTRS1995

Impact of the ongoing Amazonian deforestation on local precipitation: A GCM simulation study

Numerical simulation experiments were conducted to delineate the influence of in situ deforestation data on episodic rainfall by comparing two ensembles of five 5-day integrations performed with a recent version of the Goddard Laboratory for Atmospheres General Circulation Model (GCM) that has a simple biosphere model (SiB). The first set, called control cases, used the standard SiB vegetation cover (comprising 12 biomes) and assumed a fully forested Amazonia, while the second set, called deforestation cases, distinguished the partially deforested regions of Amazonia as savanna. Except for this difference, all other initial and prescribed boundary conditions were kept identical in both sets of integrations. The differential analyses of these five cases show the following local effects of deforestation. (1) A discernible decrease in evapotranspiration of about 0.80 mm/d (roughly 18%) that is quite robust in the averages for 1-, 2-, and 5-day forecasts. (2) A decrease in precipitation of about 1.18 mm/d (roughly 8%) that begins to emerge even in 1-2 day averages and exhibits complex evolution that extends downstream with the winds. (3) A significant decrease in the surface drag force (as a consequence of reduced surface roughness of deforested regions) that, in turn, affects the dynamical structure of moisture convergence and circulation. The surface winds increase significantly during the first day, and thereafter the increase is well maintained even in the 2- and 5-day averages.

Walker, G. K.↗

NASA NTRS1994

A class of the van Leer-type transport schemes and its application to the moisture transport in a general circulation model

A generalized form of the second-order van Leer transport scheme is derived. Several constraints to the implied subgrid linear distribution are discussed. A very simple positive-definite scheme can be derived directly from the generalized form. A monotonic version of the scheme is applied to the Goddard Laboratory for Atmospheres (GLA) general circulation model (GCM) for the moisture transport calculations, replacing the original fourth-order center-differencing scheme. Comparisons with the original scheme are made in idealized tests as well as in a summer climate simulation using the full GLA GCM. A distinct advantage of the monotonic transport scheme is its ability to transport sharp gradients without producing spurious oscillations and unphysical negative mixing ratio. Within the context of low-resolution climate simulations, the aforementioned characteristics are demonstrated to be very beneficial in regions where cumulus convection is active. The model-produced precipitation pattern using the new transport scheme is more coherently organized both in time and in space, and correlates better with observations. The side effect of the filling algorithm used in conjunction with the original scheme is also discussed, in the context of idealized tests. The major weakness of the proposed transport scheme with a local monotonic constraint is its substantial implicit diffusion at low resolution. Alternative constraints are discussed to counter this problem.

Lin, Shian-Jiann↗

NASA NTRS1994

Climatology and natural variability of the global hydrologic cycle in the GLA atmospheric general circulation model

Time average climatology and low-frequency variabilities of the global hydrologic cycle (GHC) in the Goddard Laboratory for Atmospheres (GLA) general circulation model (GCM) were investigated in the present work. A 730-day experiment was conducted with the GLA GCM forced by insolation, sea surface temperature, and ice-snow undergoing climatological annual cycles. Ifluences of interactive soil moisture on time average climatology and natural variability of the GHC were also investigated by conducting 365-day experiments with and without interactive soil moisture. Insolation, sea surface temperature, and ice-snow were fixed at their July levels in the latter two experiments. Results show that the model's time average hydrologic cycle variables for July in all three experiments agree reasonably well with observations. Except in the case of precipitable water, the zonal average climates of the annual cycle experiment and the two perpetual July experiments are alike, i.e., their differences are within limits of the natural variability of the model's climate. Statistics of various components of the GHC, i.e., water vapor, evaporation, and precipitation, are significantly affected by the presence of interactive soil moisture. A long-term trend is found in the principal empirical modes of variability of ground wetness, evaporation, and sensible heat. Dominant modes of variability of these quantities over land are physically consistent with one another and with land surface energy balance requirements. The dominant mode of precipitation variability is found to be closely related to organized convection over the tropical western Pacific Ocean. The precipitation variability has timescales in the range of 2 to 3 months and can be identified with the stationary component of the Madden-Julian Oscillation. The precipitation mode is not sensitive to the presence of interactive soil moisture but is closely linked to both the rotational and divergent components of atmospheric moisture transport. The present results indicate that globally coherent natural variability of the GHC in the GLA GCM has two basic timescales in the absence of annual cycles of external forcings: a long-term trend associated with atmosphere-soil moisture interaction which affects the model atmosphere mostly over midlatitude continental regions and a large-scale 2- to 3-month variability associated with atmospheric moist processes over the western Pacific Ocean.

Lau, K.-M.↗

NASA NTRS1993

Global fields of soil moisture and land surface evapotranspiration derived from observed precipitation and surface air temperature

The global fields of normal monthly soil moisture and land surface evapotranspiration are derived with a simple water budget model that has precipitation and potential evapotranspiration as inputs. The precipitation is observed and the potential evapotranspiration is derived from the observed surface air temperature with the empirical regression equation of Thornthwaite (1954). It is shown that at locations where the net surface radiation flux has been measured, the potential evapotranspiration given by the Thornthwaite equation is in good agreement with those obtained with the radiation-based formulations of Priestley and Taylor (1972), Penman (1948), and Budyko (1956-1974), and this provides the justification for the use of the Thornthwaite equation. After deriving the global fields of soil moisture and evapotranspiration, the assumption is made that the potential evapotranspiration given by the Thornthwaite equation and by the Priestley-Taylor equation will everywhere be about the same; the inverse of the Priestley-Taylor equation is used to obtain the normal monthly global fields of net surface radiation flux minus ground heat storage. This and the derived evapotranspiration are then used in the equation for energy conservation at the surface of the earth to obtain the global fields of normal monthly sensible heat flux from the land surface to the atmosphere.

Mintz, Y.↗

NASA NTRS1993

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.

Liston, G. E.↗