SEARCH · Search NASA

Results for “Precipitation parameterization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

NASA NTRS2019

An Integrated Approach to Weather Radar Calibration and Monitoring Using Ground Clutter and Satellite Comparisons

The stability and accuracy of weather radar reflectivity calibration are imperative for quantitative applications, such as rainfall estimation, severe weather monitoring and nowcasting, and assimilation in numerical weather prediction models. Various radar calibration and monitoring techniques have been developed, but only recently have integrated approaches been proposed, that is, using different calibration techniques in combination. In this paper the following three techniques are used: 1) ground clutter monitoring, 2) comparisons with spaceborne radars, and 3) the self-consistency of polarimetric variables. These techniques are applied to a C-band polarimetric radar (CPOL) located in the Australian tropics since 1998. The ground clutter monitoring technique is applied to each radar volumetric scan and provides a means to reliably detect changes in calibration, relative to a baseline. It is remarkably stable to within a standard deviation of 0.1 dB (decibels). To obtain an absolute calibration value, CPOL observations are compared to spaceborne radars on board TRMM (Tropical Rainfall Measuring Mission) and GPM (Global Precipitation Measurement) using a volume-matching technique. Using an iterative procedure and stable calibration periods identified by the ground echoes technique, we improve the accuracy of this technique to about 1 dB. Finally, we review the self-consistency technique and constrain its assumptions using results from the hybrid TRMM-GPM and ground echo technique. Small changes in the self-consistency parameterization can lead to 5 dB of variation in the reflectivity calibration. We find that the drop-shape model of Brandes et al. with a standard deviation of the canting angle of 12 degrees best matches our dataset.

Louf, Valentin↗

NASA NTRS2019

DYAMOND: The DYnamics of the Atmospheric General Circulation Modeled on Non-Hydrostatic Domains

A review of the experimental protocol and motivation for DYAMOND, the first intercomparison project of global storm-resolving models, is presented. Nine models submitted simulation output for a 40-day (1 August–10 September 2016) intercomparison period. Eight of these employed a tiling of the sphere that was uniformly less than 5 km. By resolving the transient dynamics of convective storms in the tropics, global storm-resolving models remove the need to parameterize tropical deep convection, providing a fundamentally more sound representation of the climate system and a more natural link to commensurately high-resolution data from satellite-borne sensors. The models and some basic characteristics of their output are described in more detail, as is the availability and planned use of this output for future scientific study. Tropically and zonally averaged energy budgets, precipitable water distributions, and precipitation from the model ensemble are evaluated, as is their representation of tropical cyclones and the predictability of column water vapor, the latter being important for tropical weather.

Bjorn Stevens↗

NASA NTRS2021

Cold Season Performance of the NU-WRF Regional Climate Model in the Great Lakes Region

As Earth’s largest collection of freshwater, the Laurentian Great Lakes have enormous ecological and socio-economic value. Their basin has become a regional hotspot of climatic and limnological change, potentially threatening its vital natural resources. Consequentially, there is a need to assess the current state of climate models regarding their performance across the Great Lakes region and develop the next generation of high-resolution regional climate models to address complex limnological processes and lake-atmosphere interactions. In response to this need, the current paper focuses on the generation and analysis of a 20-member ensemble of 3-km National Aeronautics and Space Administration (NASA)-Unified Weather Research and Forecasting (NU-WRF) simulations for the 2014-2015 cold season. The study aims to identify the model’s strengths and weaknesses; optimal configuration for the region; and the impacts of different physics parameterizations, coupling to a 1D lake model, time-variant lake-surface temperatures, and spectral nudging. Several key biases are identified in the cold-season simulations for the Great Lakes region, including an atmospheric cold bias that is amplified by coupling to a 1D lake model but diminished by applying the Community Atmosphere Model radiation scheme and Morrison microphysics scheme; an excess precipitation bias; anomalously early initiation of fall lake turnover and subsequent cold lake bias; excessive and overly persistent lake ice cover; and insufficient evaporation over Lakes Superior and Huron. The research team is currently addressing these key limitations by coupling NU-WRF to a 3D lake model in support of the next generation of regional climate models for the critical Great Lakes Basin.

Michael Notaro↗

NASA NTRS2022

Future Climate Change Under SSP Emission Scenarios with GISS-E2.1

This paper presents the response to anthropogenic forcing in the GISS-E2.1 climate models for the 21st century Shared Socioeconomic Pathways emission scenarios within the Coupled Model Intercomparison Project Phase 6 (CMIP6). The experiments were performed using an updated and improved version of the NASA Goddard Institute for Space Studies (GISS) coupled general circulation model that includes two different versions for atmospheric composition: A non-interactive version (NINT) with prescribed composition and a tuned aerosol indirect effect and the One-Moment Aerosol model (OMA) version with fully interactive aerosols which includes a parameterized first indirect aerosol effect on clouds. The effective climate sensitivities are 3.0°C and 2.9°C for the NINT and OMA models, respectively. Each atmospheric version is coupled to two different ocean general circulation models: The GISS ocean model (E2.1-G) and HYCOM (E2.1-H). We describe the global mean responses for all future scenarios and spatial patterns of change for surface air temperature and precipitation for four of the marker scenarios: SSP1-2.6, SSP2-4.5, SSP4-6.0, and SSP5-8.5. By 2100, global mean warming ranges from 1.5°C to 5.2°C relative to 1,850–1,880 mean temperature. Two high-mitigation scenarios SSP1-1.9 and SSP1-2.6 limit the surface warming to below 2°C by the end of the 21st century, except for the NINT E2.1-H model that simulates 2.2°C of surface warming. For the high emission scenario SSP5-8.5, the range is 4.6–5.2°C at 2100. Due to about 15% larger effective climate sensitivity and stronger transient climate response in both NINT and OMA CMIP6 models compared to CMIP5 versions, there is a stronger warming by 2100 in the SSP emission scenarios than in the comparable Representative Concentration Pathway (RCP) scenarios in CMIP5. Changes in sea ice area are highly correlated to global mean surface air temperature anomalies and show steep declines in both hemispheres, with the largest sea ice area decreases occurring during September in the Northern Hemisphere in both E2.1-G (−1.21 × 106 km2/°C) and E2.1-H models (−0.94 × 106 km2/°C). Both coupled models project decreases in the Atlantic overturning stream function by 2100. The largest decrease of 56%–65% in the 21st century overturning stream function is produced in the warmest scenario SSP5-8.5 in the E2.1-G model, comparable to the reduction in the corresponding CMIP5 GISS-E2 RCP8.5 simulation. Both low-end scenarios SSP1-1.9 and SSP1-2.6 also simulate substantial reductions of the overturning (9%–37%) with slow recovery of about 10% by the end of the 21st century (relative to the maximum decrease at the middle of the 21st century).

climate models↗

NASA NTRS1993

The annual pressure cycle on Mars: Results from the LMD Martian atmospheric general circulation model

We have been developing a General Circulation Model (GCM) of the martian atmosphere since 1989. The model has been described rather extensively elsewhere and only the main characteristics are given here. The dynamical part of the model, adapted from the LMD terrestrial climate model, is based on a finite-difference formulation of the classical 'primitive equations of meteorology.' The radiative transfer code includes absorption and emission by CO2 (carefully validated by comparison to line-by-line calculations) and dust in the thermal range and absorption and scattering by dust in the visible range. Other physical parameterizations are included: modeling of vertical turbulent mixing, dry convective adjustment (in order to prevent vertical unstable temperature profiles), and a multilayer model of the thermal conduction in the soil. Finally, the condensation-sublimation of CO2 is introduced through specification of a pressure-dependent condensation temperature. The atmospheric and surface temperatures are prevented from falling below this critical temperature by condensation and direct precipitation onto the surface of atmospheric CO2. The only prespecified spatial fields are the surface thermal inertia, albedo, and topography.

Hourdin, Frederic↗

NASA NTRS1984

The effect of physical parameterizations and initial data on the numerical prediction of the President's Day cyclone

Results are presented from a series of forecast experiments which were conducted to assess the importance of large-scale dynamical processes, diabatic heating, and initial data to the prediction of the President's Day cyclone. The synoptic situation and NMC model forecasts for this case are summarized, and the analysis/forecast system and experiments are described. The GLAS Model forecast from the GLAS analysis at 0000 GMT 18 February is found to have correctly predicted intense coastal cyclogenesis and heavy precipitation. A forecast with surface heat and moisture fluxes eliminated failed to predict any cyclogenesis while a similar forecast with only the surface moisture flux excluded showed weak development. Diabatic heating resulting from oceanic fluxes significantly contributed to the generation of low-level cyclonic vorticity and the intensification and slow rate of movement of an upper level ridge over the western Atlantic.

Atlas, R.↗

NASA NTRS2015

Cloud Feedbacks on Greenhouse Warming in a Multi-Scale Modeling Framework with a Higher-Order Turbulence Closure

Five-year simulation experiments with a multi-scale modeling Framework (MMF) with a advanced intermediately prognostic higher-order turbulence closure (IPHOC) in its cloud resolving model (CRM) component, also known as SPCAM-IPHOC (super parameterized Community Atmospheric Model), are performed to understand the fast tropical (30S-30N) cloud response to an instantaneous doubling of CO2 concentration with SST held fixed at present-day values. SPCAM-IPHOC has substantially improved the low-level representation compared with SPCAM. It is expected that the cloud responses to greenhouse warming in SPCAM-IPHOC is more realistic. The change of rising motion, surface precipitation, cloud cover, and shortwave and longwave cloud radiative forcing in SPCAM-IPHOC from the greenhouse warming will be presented in the presentation.

Cheng, Anning↗

NASA NTRS2013

Estimating the Concentration of Large Raindrops from Polarimetric Radar and Disdrometer Observations

Estimation of rainfall integral parameters, including radar observables, and empirical relations between them are sensitive to the truncation of the drop size distribution (DSD), particularly at the large drop end. The sensitivity of rainfall integral parameters to the maximum drop diameter (D(sub max)) is exacerbated at C‐band since resonance effects are pronounced for large drops in excess of 5 mm diameter (D). Due to sampling limitations, it is often difficult to reliably estimate D(sub max) with disdrometers. The resulting uncertainties in D(sub max0 potentially increase errors in radar retrieval methods, particularly at C‐band, that rely on disdrometer observations for DSD input to radar models. In fact, D(sub max) is typically an assumed DSD parameter in the development of radar retrieval methods. Because of these very uncertainties, it is difficult to independently confirm disdrometer estimates of D(sub max) with polarimetric radar observations. A couple of approaches can be taken to reduce uncertainty in large drop measurement. Longer integration times can be used for the collection of larger disdrometer samples. However, integration periods must be consistent with a radar resolution volume (RRV) and the temporal and spatial scales of the physical processes affecting the DSD therein. Multiple co‐located disdrometers can be combined into a network to increase the sample size within a RRV. However, over a reasonable integration period, a single disdrometer sample volume is many orders of magnitudes less than a RRV so it is not practical to devise a network of disdrometers that has an equivalent volume to a typical RRV. Since knowledge of DSD heterogeneity and large drop occurrence in time and space is lacking, the specific accuracy or even general representativeness of disdrometer based D(sub max) and large drop concentration estimates within a RRV are currently unknown. To address this complex issue, we begin with a simpler question. Is the frequency of occurrence of large rain drops (D > 5 mm) in disdrometer observations, either stand alone or networked, generally representative and consistent with polarimetric radar observations? We first show from simulations that the concentration of large (D > 5 mm) rain drops (N(sub T5)) can be estimated from polarimetric observations of specific differential phase (K(sub dp)) and differential reflectivity (Z(sub dr)), N(sub T5)=F(K(sub dp),Z(sub dr)), or horizontal reflectivity (Z(sub h)) and Z(sub dr), N(sub T5)=(Z(sub h),Z(sub dr)). We assess the error associated with polarimetric retrieval of N(sub T5), including sensitivity to D(sub max) parameterization assumptions and measurement error in the radar simulations. Polarimetric measurements at S‐band and C‐band will then be used to retrieve estimates of N(sub T5) and compared to disdrometer estimates of N(sub T5). After careful consideration of retrieval error, we will check consistency between disdrometer and polarimetric radar estimates of N(sub T5) and the frequency of occurrence of large rain drops in a variety of precipitating regimes using data from NASA's Global Precipitation Measurement (GPM) Ground Validation (GV) program, including field campaigns such as MC3E (Oklahoma) and IFloodS (Iowa) and extended measurements over Huntsville, Alabama and NASA Wallops Flight Facility in coastal Virginia.

Carey, Lawrence D.↗

NASA NTRS2005

The Effect of Environmental Conditions on Tropical Deep Convective Systems Observed from the TRMM Satellite

This study uses measurements of radiation and cloud properties taken between January and August 1998 by three Tropical Rainfall Measuring Mission (TRMM) instruments, the Clouds and the Earth's Radiant Energy System (CERES) scanner, the TRMM Microwave Imager (TMI), and the Visible and InfraRed Scanner (VIRS), to evaluate the variations of tropical deep convective systems (DCS) with sea surface temperature (SST) and precipitation. This study finds that DCS precipitation efficiency increases with SST at a rate of approx. 2%/K. Despite increasing rainfall efficiency, the cloud areal coverage rises with SST at a rate of about 7%/K in the warm tropical seas. There, the boundary layer moisture supply for deep convection and the moisture transported to the upper troposphere for cirrus-anvil cloud formation increase by approx. 6.3%/K and approx. 4.0%/K, respectively. The changes in cloud formation efficiency, along with the increased transport of moisture available for cloud formation, likely contribute to the large rate of increasing DCS areal coverage. Although no direct observations are available, the increase of cloud formation efficiency with rising SST is deduced indirectly from measurements of changes in the ratio of DCS ice water path and boundary layer water vapor amount with SST. Besides the cloud areal coverage, DCS cluster effective sizes also increase with precipitation. Furthermore, other cloud properties, such as cloud total water and ice water paths, increase with SST. These changes in DCS properties will produce a negative radiative feedback for the earth's climate system due to strong reflection of shortwave radiation by the DCS. These results significantly differ from some previous hypothesized dehydration scenarios for warmer climates, and have great potential in testing current cloud-system resolving models and convective parameterizations of general circulation models.

Lin, Bing↗

NASA NTRS2014

Assessment of NASA GISS CMIP5 and Post-CMIP5 Simulated Clouds and TOA Radiation Budgets Using Satellite Observations. Part I: Cloud Fraction and Properties

Although many improvements have been made in phase 5 of the Coupled Model Intercomparison Project (CMIP5), clouds remain a significant source of uncertainty in general circulation models (GCMs) because their structural and optical properties are strongly dependent upon interactions between aerosol/cloud microphysics and dynamics that are unresolved in such models. Recent changes to the planetary boundary layer (PBL) turbulence and moist convection parameterizations in the NASA GISS Model E2 atmospheric GCM(post-CMIP5, hereafter P5) have improved cloud simulations significantly compared to its CMIP5 (hereafter C5) predecessor. A study has been performed to evaluate these changes between the P5 and C5 versions of the GCM, both of which used prescribed sea surface temperatures. P5 and C5 simulated cloud fraction (CF), liquid water path (LWP), ice water path (IWP), cloud water path (CWP), precipitable water vapor (PWV), and relative humidity (RH) have been compared to multiple satellite observations including the Clouds and the Earth's Radiant Energy System-Moderate Resolution Imaging Spectroradiometer (CERES-MODIS, hereafter CM), CloudSat- Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO; hereafter CC), Atmospheric Infrared Sounder (AIRS), and Advanced Microwave Scanning Radiometer for Earth Observing System (AMSR-E). Although some improvements are observed in the P5 simulation on a global scale, large improvements have been found over the southern midlatitudes (SMLs), where correlations increased and both bias and root-mean-square error (RMSE) significantly decreased, in relation to the previous C5 simulation, when compared to observations. Changes to the PBL scheme have resulted in improved total column CFs, particularly over the SMLs where marine boundary layer (MBL) CFs have increased by nearly 20% relative to the previous C5 simulation. Globally, the P5 simulated CWPs are 25 gm22 lower than the previous C5 results. The P5 version of the GCM simulates PWV and RH higher than its C5 counterpart and agrees well with the AMSR-E and AIRS observations. The moister atmospheric conditions simulated by P5 are consistent with the CF comparison and provide a strong support for the increase in MBL clouds over the SMLs. Over the tropics, the P5 version of the GCM simulated total column CFs and CWPs are slightly lower than the previous C5 results, primarily as a result of the shallower tropical boundary layer in P5 relative to C5 in regions outside the marine stratocumulus decks.

General circulation models↗

NASA NTRS2003

Optical Emissions from Proton Aurora

Hydrogen emissions are the signature of proton aurora. The Doppler-shifted hydrogen emission lines can be interpreted in terms of the mean energy of the precipitating protons. A red shifted component of the line profiles observed from the ground indicates upward going hydrogen atoms due to angular redistribution of the precipitation. Secondary electrons from ionization and stripping collisions also contribute to the auroral emissions. Since the energy distribution of these secondaries has a lower mean energy than secondary electrons in electron aurora, the relative brightness of eniission features differs from that in electron aurora. The secondaries contribute little to additional ionization. These differences between proton and electron aurora can lead to misinterpretation when brightness ratios are used to derive ionospheric conductances with parameterizations that are based on electron aurora.

Lummerzheim, D.↗

NASA NTRS2014

The Impacts of Microphysics and Planetary Boundary Layer Physics on Model Simulations of U. S. Deep South Summer Convection

Inspection of output from various configurations of high-resolution, explicit convection forecast models such as the Weather Research and Forecasting (WRF) model indicates significant sensitivity to the choices of model physics parameterizations employed. Some of the largest apparent sensitivities are related to the specifications of the cloud microphysics and planetary boundary layer physics packages. In addition, these sensitivities appear to be especially pronounced for the weakly-sheared, multicell modes of deep convection characteristic of the Deep South of the United States during the boreal summer. Possible ocean-land sensitivities also argue for further examination of the impacts of using unique ocean-land surface initialization datasets provided by the NASA Short-term Prediction Research and Transition (SPoRT Center to select NOAA/NWS weather forecast offices. To obtain better quantitative understanding of these sensitivities and also to determine the utility of the ocean-land initialization data, we have executed matrices of regional WRF forecasts for selected convective events near Mobile, AL (MOB), and Houston, TX (HGX). The matrices consist of identically initialized WRF 24-h forecasts using any of eight microphysics choices and any of three planetary boundary layer choices. The resulting 24 simulations performed for each event within either the MOB or HGX regions are then compared to identify the sensitivities of various convective storm metrics to the physics choices. Particular emphasis is placed on sensitivities of precipitation timing, intensity, and coverage, as well as amount and coverage of lightning activity diagnosed from storm kinematics and graupel in the mixed phase layer. The results confirm impressions gleaned from study of the behavior of variously configured WRF runs contained in the ensembles produced each spring at the Center for the Analysis and Prediction of Storms, but with the benefit of more straightforward control of the physics package choices. The design of the experiments thus allows for more direct interpretation of the sensitivities to each possible physics combination. The results should assist forecasters in their efforts to anticipate and correct for possible biases in simulated WRF convection patterns, and help the modeling community refine their model parameterizations.

McCaul, E. W., Jr.↗

NASA NTRS1989

Effects of cumulus entrainment and multiple cloud types on a January global climate model simulation

An improved version of the GISS Model II cumulus parameterization designed for long-term climate integrations is used to study the effects of entrainment and multiple cloud types on the January climate simulation. Instead of prescribing convective mass as a fixed fraction of the cloud base grid-box mass, it is calculated based on the closure assumption that the cumulus convection restores the atmosphere to a neutral moist convective state at cloud base. This change alone significantly improves the distribution of precipitation, convective mass exchanges, and frequencies in the January climate. The vertical structure of the tropical atmosphere exhibits quasi-equilibrium behavior when this closure is used, even though there is no explicit constraint applied above cloud base.

Yao, Mao-Sung↗

NASA NTRS2003

Precipitation Processes developed during ARM (1997), TOGA COARE (1992), GATE (1974), SCSMEX (1998), and KWAJEX (1999), Consistent 2D, semi-3D and 3D Cloud Resolving Model Simulations

Real clouds and cloud systems are inherently three-dimensional (3D). Because of the limitations in computer resources, however, most cloud-resolving models (CRMs) today are still two-dimensional (2D). A few 3D CRMs have been used to study the response of clouds to large-scale forcing. In these 3D simulations, the model domain was small, and the integration time was 6 hours. The major objectives of this paper are: (1) to assess the performance of the super-parameterization technique (i.e. is 2D or semi-3D CRM appropriate for the super-parameterization?); (2) calculate and examine the surface energy (especially radiation) and water budgets; (3) identify the differences and similarities in the organization and entrainment rates of convection between simulated 2D and 3D cloud systems.

Tao, Wei-Kuo↗

NASA NTRS1985

Sensitivity of an energy balance climate model with predicted snowfall rates

A snowfall parameterization and a polar-ice-sheet model are developed and applied to the two-level zonally averaged seasonal energy-balance climate model of Held and Suarez (1979), and sensitivity experiments involving changes in insolation are performed both with and without ice sheets. The results are presented in tables and graphs, and the hydrological-cycle response to insolation changes is found to be similar to that predicted by global-circulation models employing prescribed precipitation levels, with a somewhat higher sensitivity in the snow line. The area covered by ice sheets in the ice-sheet models is shown to be greater than that covered by permanent snow in the models without ice sheets, an effect attributed to lower surface temperatures over the ice. It is inferred that an increase in the solar constant can cause increased high-latitude precipitation but not an ice age.

Bowman, K. P.↗

NASA NTRS2007

Simulations of Tropospheric NO2 by the Global Modeling Initiative (GMI) Model Utilizing Assimilated and Forecast Meteorological Fields: Comparison to Ozone Monitoring Instrument (OMI) Measurements

We present simulations of the tropospheric composition for the years 2004 and 2005, carried out by the GMI Combined Stratosphere-Troposphere (Combo) model, at a resolution of 2degx2.5deg. The model includes a new parameterization of lightning sources of NO(x) which is coupled to the cloud mass fluxes in the adopted meteorological fields. These simulations use two different sets of input meteorological fields: a)late-look assimilated fields from the Global Modeling and Assimilation Office (GMAO), GEOS-4 system and b) 12-hour forecast fields initialized with the assimilated data. Comparison of the forecast to the assimilated fields indicates that the forecast fields exhibit less vigorous convection, and yield tropical precipitation fields in better agreement with observations. Since these simulations include a complete representation of the stratosphere, they provide realistic stratosphere-tropospheric fluxes of O3 and NO(y). Furthermore, the stratospheric contribution to total columns of different troposheric species can be subtracted in a consistent fashion, and the lightning production of NO(y) will depend on the adopted meteorological field. We concentrate here on the simulated tropospheric columns of NO2, and compare them to observations by the OM1 instrument for the years 2004 and 2005. The comparison is used to address these questions: a) is there a significant difference in the agreement/disagreement between simulations for these two different meteorological fields, and if so, what causes these differences?; b) how do the simulations compare to OMI observations, and does this comparison indicate an improvement in simulations with the forecast fields? c) what are the implications of these simulations for our understanding of the NO2 emissions over continental polluted regions?

Rodriquez, J. M.↗

NASA NTRS2016

A Model-Model and Data-Model Comparison for the Early Eocene Hydrological Cycle

A range of proxy observations have recently provided constraints on how Earth's hydrological cycle responded to early Eocene climatic changes. However, comparisons of proxy data to general circulation model (GCM) simulated hydrology are limited and inter-model variability remains poorly characterised. In this work, we undertake an intercomparison of GCM-derived precipitation and P - E distributions within the extended EoMIP ensemble (Eocene Modelling Intercomparison Project; Lunt et al., 2012), which includes previously published early Eocene simulations performed using five GCMs differing in boundary conditions, model structure, and precipitation-relevant parameterisation schemes. We show that an intensified hydrological cycle, manifested in enhanced global precipitation and evaporation rates, is simulated for all Eocene simulations relative to the preindustrial conditions. This is primarily due to elevated atmospheric paleo-CO2, resulting in elevated temperatures, although the effects of differences in paleogeography and ice sheets are also important in some models. For a given CO2 level, globally averaged precipitation rates vary widely between models, largely arising from different simulated surface air temperatures. Models with a similar global sensitivity of precipitation rate to temperature (dP=dT ) display different regional precipitation responses for a given temperature change. Regions that are particularly sensitive to model choice include the South Pacific, tropical Africa, and the Peri-Tethys, which may represent targets for future proxy acquisition. A comparison of early and middle Eocene leaf-fossil-derived precipitation estimates with the GCM output illustrates that GCMs generally underestimate precipitation rates at high latitudes, although a possible seasonal bias of the proxies cannot be excluded. Models which warm these regions, either via elevated CO2 or by varying poorly constrained model parameter values, are most successful in simulating a match with geologic data. Further data from low-latitude regions and better constraints on early Eocene CO2 are now required to discriminate between these model simulations given the large error bars on paleoprecipitation estimates. Given the clear differences between simulated precipitation distributions within the ensemble, our results suggest that paleohydrological data offer an independent means by which to evaluate model skill for warm climates.

Boundary conditions↗

NASA NTRS1999

Land-Atmosphere Interactions: Successes, Problems and Prospects

After two decades of active research, a much better understanding of the broader role of biospheric processes on the local climate has emerged. A surface-albedo increase, particularly in desert border regions of the subtropics (as well as the deforested tropical regions), leads to a net surface energy deficit, which in turn leads to a relative sinking and reduced rainfall. On the other hand, studies of the influence of altered ratios of evapotranspiration and sensible fluxes, in situations where the net solar income is unchanged, show that evapotranspiration is a more desirable flux for increased precipitation and vitality of the biosphere. Besides providing water vapor and convective available potential energy (CAPE) to the lower troposphere, evapotranspiration helps in building larger CAPE before "turning on" the moist-convection. Larger CAPE in the lower troposphere enables convection to reach into the deeper atmosphere thereby heating the upper troposphere; indeed, moist-convection is also accompanied by the evaporation of falling precipitation that cools and moistens the lower atmosphere. While convective, as opposed to stratiform, precipitation reduces the fractional cloud cover; it also allows more solar radiation to reach the surface thereby invigorating surface fluxes. These, together with moist convection and associated downdrafts help to maintain the characteristic upper temperature limit(s) of the moist-land as well as oceanic regions. Regardless of the above understanding, several important problems continue to hinder the accurate simulation of a realistic land atmosphere interaction in a numerical model (both GCM and/or Meso-scale models). Among the unsolved problems are parameterization of sub-grid scale land processes that include small-scale variability of soil moisture, snow-cover and snow-physics, the biodiversity of the biosphere, orography, local drainage characteristics under natural conditions, and surface flow over the natural terrain. A well-known non-linear response of surface fluxes to these variations makes the problem of parameterizing land-atmosphere interaction processes hard-to-address and simulate, particularly in a GCM. In our presentation, we will discuss how orographic, snow-cover, and water table interactions can be included into a Simple Biosphere Model such as SiB/SSiB. Figure I shows how, in the Russian region, spring snowmelt affects the soil moisture profile. Corresponding figure 2 shows how interaction with the water table decreases the natural evapotranspiration in the Sahel region simulation. While these simulations need better validation with data, the simulations reveal that surface processes are sensitive to these parameterizations. With these developments, we continue to advance our understanding of the interaction of land with the atmosphere aloft, but the intrinsic variability of the newer parameters, e. g., hydraulic properties of the soil, diminish the positive influences of these advances on the improved climate simulation with GCMs.

Sud, Y. C.↗