Search NASASearch

Engineering topics

Dickinson, R. E.

Publications and source records attributed to Dickinson, R. E..

Cabauw Experimental Results from the Project for Intercomparison of Land-Surface Parameterization Schemes

In the Project for Intercomparison of Land-Surface Parameterization Schemes phase 2a experiment, meteorological data for the year 1987 from Cabauw, the Netherlands, were used as inputs to 23 land-surface flux schemes designed for use in climate and weather models. Schemes were evaluated by comparing their outputs with long-term measurements of surface sensible heat fluxes into the atmosphere and the ground, and of upward longwave radiation and total net radiative fluxes, and also comparing them with latent heat fluxes derived from a surface energy balance. Tuning of schemes by use of the observed flux data was not permitted. On an annual basis, the predicted surface radiative temperature exhibits a range of 2 K across schemes, consistent with the range of about 10 W/m in predicted surface net radiation. Most modeled values of monthly net radiation differ from the observations by less than the estimated maximum monthly observational error (+/- 10 W/sq m). However, modeled radiative surface temperature appears to have a systematic positive bias in most schemes; this might be explained by an error in assumed emissivity and by models' neglect of canopy thermal heterogeneity. Annual means of sensible and latent heat fluxes, into which net radiation is partitioned, have ranges across schemes of 30 W/sq m and 25 W/sq m, respectively. Annual totals of evapotranspiration and runoff, into which the precipitation is partitioned, both have ranges of 315 mm. These ranges in annual heat and water fluxes were approximately halved upon exclusion of the three schemes that have no stomatal resistance under non-water-stressed conditions. Many schemes tend to underestimate latent heat flux and overestimate sensible heat flux in summer, with a reverse tendency in winter. For six schemes, root-mean-square deviations of predictions from monthly observations are less than the estimated upper bounds on observation errors (5 W/m for sensible beat flux and 10 W/m for latent heat flux). Actual runoff at the site is believed to be dominated by vertical drainage to ground water, but several schemes produced significant amounts of runoff as overland flow or interflow. There is a range across schemes of 184 mm (40% of total pore volume) in the simulated annual mean root-zone soil moisture. Unfortunately, no measurements of soil moisture were available for model evaluation. A theoretical analysis suggested that differences in boundary conditions used in various schemes are not sufficient to explain the large variance in soil moisture. However, many of the extreme values of soil moisture could be explained in terms of the particulars of experimental setup or excessive evapotranspiration.

Chen, Tian Hong

Modeling the bidirectional reflectance distribution function of mixed finite plant canopies and soil

An analytical model of the bidirectional reflectance for optically semi-infinite plant canopies has been extended to describe the reflectance of finite depth canopies contributions from the underlying soil. The model depends on 10 independent parameters describing vegetation and soil optical and structural properties. The model is inverted with a nonlinear minimization routine using directional reflectance data for lawn (leaf area index (LAI) is equal to 9.9), soybeans (LAI, 2.9) and simulated reflectance data (LAI, 1.0) from a numerical bidirectional reflectance distribution function (BRDF) model (Myneni et al., 1988). While the ten-parameter model results in relatively low rms differences for the BRDF, most of the retrieved parameters exhibit poor stability. The most stable parameter was the single-scattering albedo of the vegetation. Canopy albedo could be derived with an accuracy of less than 5% relative error in the visible and less than 1% in the near-infrared. Sensitivity were performed to determine which of the 10 parameters were most important and to assess the effects of Gaussian noise on the parameter retrievals. Out of the 10 parameters, three were identified which described most of the BRDF variability. At low LAI values the most influential parameters were the single-scattering albedos (both soil and vegetation) and LAI, while at higher LAI values (greater than 2.5) these shifted to the two scattering phase function parameters for vegetation and the single-scattering albedo of the vegetation. The three-parameter model, formed by fixing the seven least significant parameters, gave higher rms values but was less sensitive to noise in the BRDF than the full ten-parameter model. A full hemispherical reflectance data set for lawn was then interpolated to yield BRDF values corresponding to advanced very high resolution radiometer (AVHRR) scan geometries collected over a period of nine days. The resulting parameters and BRDFs are similar to those for the full sampling geometry, suggesting that the limited geometry of AVHRR measurements might be used to reliably retrieve BRDF and canopy albedo with this model.

Schluessel, G.

The Project for Intercomparison of Land-surface Parameterization Schemes

The Project for Intercomparison of Land-surface Parameterization Schemes (PILPS) is described and the first stage science plan outlined. PILPS is a project designed to improve the parameterization of the continental surface, especially the hydrological, energy, momentum, and carbon exchanges with the atmosphere. The PILPS Science Plan incorporates enhanced documentation, comparison, and validation of continental surface parameterization schemes by community participation. Potential participants include code developers, code users, and those who can provide datasets for validation and who have expertise of value in this exercise. PILPS is an important activity because existing intercomparisons, although piecemeal, demonstrate that there are significant differences in the formulation of individual processes in the available land surface schemes. These differences are comparable to other recognized differences among current global climate models such as cloud and convection parameterizations. It is also clear that too few sensitivity studies have been undertaken with the result that there is not yet enough information to indicate which simplifications or omissions are important for the near-surface continental climate, hydrology, and biogeochemistry. PILPS emphasizes sensitivity studies with and intercomparisons of existing land surface codes and the development of areally extensive datasets for their testing and validation.

Henderson-Sellers, A.

How will changes in carbon dioxide and methane modify the mean structure of the mesosphere and thermosphere?

A global average model of the coupled mesosphere, thermosphere, and ionosphere is used to examine the effect of trace gas variations on the overall structure of these regions. In particular, the variations caused by CO2 and CH4 doublings and halvings from present day mixing ratios are presented. The results indicate that the mesosphere and thermosphere temperatures will cool by about 10 K and 50 K, respectively, as the CO2 and CH4 mixing ratios are doubled. These regions are heated by similar amounts when the trace gas mixing ratios are halved. Compositional redistributions also occur in association with changes in the temperature profile. The results show that global change will occur in the upper atmosphere and ionosphere as well as in the lower atmosphere during the 21st century.

Roble, R. G.

Taylor instability in rhyolite lava flows

A refined Taylor instability model is developed to describe the surface morphology of rhyolite lava flows. The effect of the downslope flow of the lava on the structures resulting from the Taylor instability mechanism is considered. Squire's (1933) transformation is developed for this flow in order to extend the results to three-dimensional modes. This permits assessing why ridges thought to arise from the Taylor instability mechanism are preferentially oriented transverse to the direction of lava flow. Measured diapir and ridge spacings for the Little and Big Glass Mountain rhyolite flows in northern California are used in conjunction with the model in order to explore the implications of the Taylor instability for flow emplacement. The model suggests additional lava flow features that can be measured in order to test whether the Taylor instability mechanism has influenced the flows surface morphology.

Baum, B. A.

Venus mesosphere and thermosphere. III - Three-dimensional general circulation with coupled dynamics and composition

The three-dimensional structure and circulation of Venus' upper mesosphere and thermosphere is examined by means of a modification of NCAR's earth thermosphere general circulation model, using the parameterizations from an earlier two-dimensional Venus model that included eddy diffusion and wave drag. Many of the observed Venus thermosphere features are found to be reproduced by the model, which thereby serves as a benchmark on which to incorporate additional minor constituents and test new self-consistent parameterizations for wave drag and superrotation.

Bougher, S. W.

On the global mean structure of the thermosphere

A self-consistent model of the global mean structure of the thermosphere is used to examine the ionospheric and thermospheric processes that control the structure of the thermosphere for solar minimum and maximum conditions during geomagnetic quiet periods. The model includes the physical and chemical processes believed to be dominant in the thermosphere and ionosphere. The total solar energy input available in the extreme ultraviolet and in the Schumann-Runge continuum portion of the solar spectrum is 14.1 x 10 to the 11th W and 21.5 x 10 to the 11th W for solar minumum and maximum, respectively. The calculated exospheric temperature is too small if maintained only by radiative processes, and an auroral heat source is necessary to bring the calculated global mean structure into agreement with the MSIS-83 empirical model for both solar minimum and maximum conditions. These results support the consensus that there is high-latitude heat source present at all times.

Roble, R. G.

Venus mesosphere and thermosphere. II - Global circulation, temperature, and density variations

The Dickinson and Ridley (1977) symmetric, two-dimensional hydrodynamical model framework is presently used as the basis of a reexamination of the circulation and structure of the Venus thermosphere recently revealed by Pioneer Venus observations. The observed day-to-night variation of composition and temperatures can largely be derived by a wave-drag parameterization yielding a weaker circulation system than that predicted prior to Pioneer Venus. It is also suggested that eddy diffusion is a minor contributor to the maintenance of observed day and nightside densities, and that eddy coefficients are smaller than than those of one-dimensional composition models previously employed.

Bougher, S. W.

Applications of a Venus thermospheric circulation model

A variety of Pioneer Venus observations suggest a global scale, day-to-night Venus thermospheric circulation. Model studies of the dynamics and energetics of the Venus thermosphere are presented in order to address new driving, mixing and cooling mechanisms for an improved model simulation. The adopted approach was to reexamine the circulation by first using a previous two dimensional code to quantify those physical processes which can be inferred from the Pioneer Venus observations. Specifically, the model was used to perform sensitivity studies to determine the degree to which eddy cooling, eddy or wave drag, eddy diffusion and 15 micrometer radiational cooling are necessary to bring the model temperature and composition fields into agreement with observations. Three EUV heating cases were isolated for study. Global temperature and composition fields in good agreement with Pioneer data were obtained. Large scale horizontal winds 220 m/s were found to be consistent with the observed cold nightside temperatures and dayside bulges of O, CO and CO2. Observed dayside temperatures were obtained by using a 7 to 19% EUV heating efficiency profile. The enhanced 15 micrometer cooling needed for thermal balance is obtained using the best rate coefficient available for atomic O collisional excitation of CO2(0,1,0). Eddy conduction was not found to be a viable cooling mechanism due to the weakened global circulation. The strong 15 micrometer damping and low EUV efficiency imply a very weak dependence of the general circulation to solar cycle variability. The NCAR terrestrial thermospheric general circulation model was adapted for Venus inputs using the above two dimensional model parameters, to give a three dimensional benchmark for future Venus modelling work.

Bougher, S. W.

Simulation of the thermospheric tides by use of the NCAR thermospheric general circulation model

Numerical calculations of the thermospheric tidal winds and temperatures at equinox are presented. The calculations were made using the National Center for Atmospheric Research Thermospheric General Circulation Model (TGCM), which includes the effects of viscosity, conductivity, diffusion, ion drag, winds, and temperature gradients. The thermospheric diurnal and semidiurnal tides are excited in situ by solar heating and by ion-neutral momentum coupling. The semidiurnal tidal calculations also include the effects of upward propagating waves generated by heating in the lower atmosphere. This semidiurnal propagating component is modeled by use of the classical tidal perturbations as lower boundary conditions. The model is tuned by adjusting the propagating tidal forcing term until calculated semidiurnal wind and temperature fields best approximate incoherent scatter observations. The tidal TGCM results are consistent with previous theoretical work and successfully reproduce high altitude temperature and meridional velocity data, but they give significantly lower magnitudes for velocities and temperatures near 160 km than are seen by observations.

Fensen, C. G.

Venus mesosphere and thermosphere. I - Heat budget and thermal structure

The calculation of a reasonable global average temperature structure is the first step in the simulation of the dynamics and observed characteristics of Venus's upper mesosphere and thermosphere. Various physical processes that influence Venus temperatures are examined, using updated inputs and new parameterizations. It is demonstrated, using a one-dimensional NLTE radiative transfer code, that temperatures observed during the Pioneer Venus mission can be reproduced using an O-CO2 collisional excitation rate coefficient of 4 x 10 to the -13th cu cm/s for strong 15-micron cooling, balanced by 9.5 percent-efficient EUV heating. Cooling by eddy mixing is at best a minor contribution to the total cooling required. Exospheric temperatures are calculated to vary by 60 K or less over a solar cycle, in rough agreement with observations. It is concluded that CO2 cooling effectively buffers against such solar perturbations, due to its nonlinear temperature dependence.

Dickinson, R. E.

Radiative processes

Solar radiation and the processes that control its deposition in the Earth atmosphere are considered. The published data obtained since 1978 define a reference solar spectral irradiance for use in atmospheric chemical and dynamical studies, while long term satellite measurements are now providing information on variations in the output of the Sun over a range of time scales. As concerns absorption of solar radiation in the atmosphere, new cross section data for molecular oxygen and ozone are now available. Line-by-line calculations used to predict infrared flux divergences, both as regards assumptions made in radiative transfer calculations and in the spectroscopic parameters used as inputs are examined. Also examined are the influence of radiative processes on planetary scale wave activity, photochemical acceleration of radiative damping, and the breakdown of local thermodynamic equilibrium at mesospheric altitudes.

Frederick, J. E.

Models of Venus neutral upper atmosphere - Structure and composition

Models of the Venus neutral upper atmosphere, based on both in situ and remote sensing measurements, are provided for the height interval from 100 to 3500 km. The general approach in model formulation was to divide the atmosphere into three regions: 100-150 km, 150-250 km, and 250-3500 km. Boundary conditions at 150 km are consistent with both drag and mass spectrometer measurements. A paramount consideration was to keep the models simple enough to be used conveniently. Available observations are reviewed.

Keating, G. M.

Thermospheric circulation, temperature, and compositional structure of the southern hemisphere polar cap during October-November 1981

Results of NCAR thermospheric general circulation model (TGCM) simulations were compared with Dynamics Explorer (DE2) data on the neutral wind, temperature and O and N2 number variations over the south polar cap in the F region in October-November 1981. Attention was focused on variations in the thermospheric characteristics. All variations reached a maximum over the magnetic polar cap, with the largest O and N2 changes being present in the dawn sector of the auroral oval, where Joule heating was maximized. The TGCM simulations for conditions similar to the DE2 passes agreed well with the satellite data, indicating that the variations in the wind, temperature and compositional characteristics are dependent on the magnitude of magnetospheric convection and influenced by the separation between the geomagnetic and geographic poles.

Roble, R. G.

The high latitude circulation and temperature structure of the thermosphere near solstice

NCAR thermospheric-general-circulation-model (TGCM) computations of solar-maximum thermospheric neutral-gas temperature and circulation around the December solstice are presented and discussed. The TGCM uses a 5 x 5-deg grid and 24 constant-pressure layers, corresponding to altitudes of about 97-500 km. The results are mapped as electron-density contours, polar plots, cylindrical equidistant projections, meridional cross sections, and F-region polar plots comparing the TGCM predictions with DE-2 satellite observations. The significant differences between summer and winter high-latitude F-region winds are attributed to the ion drag momentum associated with magnetospheric convection. The TGCM wind predictions follow the same pattern as the satellite measurements but are too small; possible model corrections are considered.

Roble, R. G.

Effect of chlorofluoromethane infrared radiation on zonal atmospheric temperatures

Estimates are made of changes in the atmospheric climate due to the radiative effects of 10 ppb of chlorofluoromethanes (CFM's). The estimates are derived on the basis of a 12-layer stratospheric general circulation model with a specified change of ocean temperature. Two tropical maxima in zonal average temperature change were observed: one in the upper troposphere and one centered at the tropopause. The temperature change exceeds the surface temperature change by a factor of at least two. If the 1975 CFM emission rate were to continue indefinitely, stratospheric water-vapor concentrations would increase by up to 60% due to CFM radiative effects. This would reduce ozone concentrations by an additional 4% of the natural ozone column.

Dickinson, R. E.

Study of the relationship between solar activity and terrestrial weather

Evidence for some connection between weather and solar related phenomena is presented. Historical data of world wide temperature variations with relationship to change in solar luminosity are examined. Several test methods for estimating the statistical significance of such phenomena are discussed in detail.

Sturrock, P. A.