Engineering PapersSearch

Engineering topics

Ghil, M.

Publications and source records attributed to Ghil, M..

53 records · Page 3

Some experiments on the effect of remote sounding temperatures upon weather forecasting

A time-continuous statistical method is presented for the four dimensional assimilation of remote sounding temperatures based on radiance measurements from polar orbiting satellites. This method is applied to DST 6 data from the NOAA 4 and Nimbus 6 satellites. The state of the atmosphere throughout the test period was determined using a varying amount of satellite data from the NOAA 4 satellite only, from Nimbus 6 only, and from both satellites together. The methods tested included different variations of the statistical method, as well as more traditional methods. It is concluded that satellite derived temperature data can have a modest, but statistically significant positive impact on numerical weather prediction in the two to three day range, and that this impact is highly sensitive to the quantity of data available and to the assimilation method used.

Halem, M.

A stochastic-dynamic model for global atmospheric mass field statistics

A model that yields the spatial correlation structure of atmospheric mass field forecast errors was developed. The model is governed by the potential vorticity equation forced by random noise. Expansion in spherical harmonics and correlation function was computed analytically using the expansion coefficients. The finite difference equivalent was solved using a fast Poisson solver and the correlation function was computed using stratified sampling of the individual realization of F(omega) and hence of phi(omega). A higher order equation for gamma was derived and solved directly in finite differences by two successive applications of the fast Poisson solver. The methods were compared for accuracy and efficiency and the third method was chosen as clearly superior. The results agree well with the latitude dependence of observed atmospheric correlation data. The value of the parameter c sub o which gives the best fit to the data is close to the value expected from dynamical considerations.

Ghil, M.

Optimal interpolation and the Kalman filter

The estimation theory of stochastic-dynamic systems is described and used in a numerical study of optimal interpolation. The general form of data assimilation methods is reviewed. The Kalman-Bucy, KB filter, and optimal interpolation (OI) filters are examined for effectiveness in performance as gain matrices using a one-dimensional form of the shallow-water equations. Control runs in the numerical analyses were performed for a ten-day forecast in concert with the OI method. The effects of optimality, initialization, and assimilation were studied. It was found that correct initialization is necessary in order to localize errors, especially near boundary points. Also, the use of small forecast error growth rates over data-sparse areas was determined to offset inaccurate modeling of correlation functions near boundaries.

Cohn, S.

The compatible balancing approach to initialization, and four-dimensional data assimilation

A unified theoretical approach to both the four-dimensional assimilation of asynoptic data and the initialization problem is attempted. This approach relies on the derivation of certain relationships between geopotential tendencies and tendencies of the horizontal velocity field in primitive-equation models of atmospheric flow. The approach is worked out and analyzed in detail for some simple barotropic models. Certain independent results of numerical experiments for the time-continuous assimilation of real asynoptic meteorological data into a complex, baroclinic weather prediction model are discussed in the context of the present approach. Tentative inferences are drawn for practical assimilation procedures.

Ghil, M.

An energy-balance model of glaciation cycles

A one dimensional energy balance model is presented which contains a time lagged albedo to account for the delayed dependence of continental ice sheets on temperature; it also includes a smoothing of temperature gradients in the tropics to account for the effect of the Hadley circulation on the strong flattening of meridional temperature profiles there. The model exhibits finite amplitude, sustained free oscillations; these oscillations are triggered by a change in the insulation parameter and occur in the absence of any external forcing. This model behavior strongly suggests the possibility of an almost-intransitive mechanism playing a role in glaciation cycles. This behavior also occurs on shorter time scales which might be relevant to the interannual variability of the atmosphere.

Ghil, M.

Free oscillations in a climate model with ice-sheet dynamics

A study of stable periodic solutions to a simple nonlinear model of the ocean-atmosphere-ice system is presented. The model has two dependent variables: ocean-atmosphere temperature and latitudinal extent of the ice cover. No explicit dependence on latitude is considered in the model. Hence all variables depend only on time and the model consists of a coupled set of nonlinear ordinary differential equations. The globally averaged ocean-atmosphere temperature in the model is governed by the radiation balance. The reflectivity to incoming solar radiation, i.e., the planetary albedo, includes separate contributions from sea ice and from continental ice sheets. The major physical mechanisms active in the model are (1) albedo-temperature feedback, (2) continental ice-sheet dynamics and (3) precipitation-rate variations. The model has three-equilibrium solutions, two of which are linearly unstable, while one is linearly stable. For some choices of parameters, the stability picture changes and sustained, finite-amplitude oscillations obtain around the previously stable equilibrium solution. The physical interpretation of these oscillations points to the possibility of internal mechanisms playing a role in glaciation cycles.

Kallen, E.

A fast Cauchy-Riemann solver

The inhomogeneous Cauchy-Riemann equations in a rectangle are discretized by a finite difference approximation. Several different boundary conditions are treated explicitly, leading to algorithms which have overall second-order accuracy. All boundary conditions with either u or v prescribed along a side of the rectangle can be treated by similar methods. The algorithms presented here have nearly minimal time and storage requirements and seem suitable for development into a general-purpose direct Cauchy-Riemann solver for arbitrary boundary conditions.

Ghil, M.

Time-continuous assimilation of remote-sounding data and its effect on weather forecasting

Methods were developed for the time-continuous assimilation of satellite-sounding temperature data: direct insertion method, asynoptic successive correction method, and local linear regression method. These methods were applied to DST-6 data from the operational and experimental temperature sounders aboard the NOAA 4 and Nimbus 6 satellites. Attention is given to a comparison of these methods and their effect on the accuracy of the initial states obtained, as well as the resulting forecasts generated from these initial states. The results suggest that (1) satellite-derived temperature data can have a modest though statistically significant positive impact on numerical weather prediction in the 2-3 day range; (2) the impact is highly sensitive to the quantity of data available; and (3) the method of satellite-data assimilation can substantially influence the magnitude of the impact obtained for the same data.

Ghil, M.

An energy-balance model with multiply-periodic and quasi-chaotic free oscillations

A slightly modified version of the one-dimensional time-dependent energy-balance climate model of Ghil and Bhattacharya (1978) is presented. The albedo-temperature parameterization has been reformulated and the smoothing of the temperature distribution in the tropics has been eliminated. The model albedo depends on time-lagged temperature in order to account for finite growth and decay time of continental ice sheets. Two distinct regimes of oscillatory behavior which depend on the value of the albedo-temperature time lag are considered.

Bhattacharya, K.

Effects of sounding temperature assimilation on weather forecasting - Model dependence studies

In comparing various methods for the assimilation of remote sounding information into numerical weather prediction (NWP) models, the problem of model dependence for the different results obtained becomes important. The paper investigates two aspects of the model dependence question: (1) the effect of increasing horizontal resolution within a given model on the assimilation of sounding data, and (2) the effect of using two entirely different models with the same assimilation method and sounding data. Tentative conclusions reached are: first, that model improvement as exemplified by increased resolution, can act in the same direction as judicious 4-D assimilation of remote sounding information, to improve 2-3 day numerical weather forecasts. Second, that the time continuous 4-D methods developed at GLAS have similar beneficial effects when used in the assimilation of remote sounding information into NWP models with very different numerical and physical characteristics.

Ghil, M.

Subjective evaluation of the combined influence of satellite temperature sounding data and increased model resolution on numerical weather forecasting

The present evaluation is concerned with (1) the significance of prognostic differences resulting from the inclusion of satellite-derived temperature soundings, (2) how specific differences between the SAT and NOSAT prognoses evolve, and (3) comparison of two experiments using the Goddard Laboratory for Atmospheric Sciences general circulation model. The subjective evaluation indicates that the beneficial impact of sounding data is enhanced with increased resolution. It is suggested that satellite sounding data posses valuable information content which at times can correct gross analysis errors in data sparse regions.

Atlas, R.

The GISS sounding temperature impact test

The impact of DST 5 and DST 6 satellite sounding data on mid-range forecasting was studied. The GISS temperature sounding technique, the GISS time-continuous four-dimensional assimilation procedure based on optimal statistical analysis, the GISS forecast model, and the verification techniques developed, including impact on local precipitation forecasts are described. It is found that the impact of sounding data was substantial and beneficial for the winter test period, Jan. 29 - Feb. 21. 1976. Forecasts started from initial state obtained with the aid of satellite data showed a mean improvement of about 4 points in the 48 and 772 hours Sub 1 scores as verified over North America and Europe. This corresponds to an 8 to 12 hour forecast improvement in the forecast range at 48 hours. An automated local precipitation forecast model applied to 128 cities in the United States showed on an average 15% improvement when satellite data was used for numerical forecasts. The improvement was 75% in the midwest.

Halem, M.

A balanced diagnostic system compatible with a barotropic prognostic model

A system of diagnostic equations for the velocity field, or wind laws, for a barotropic primitive-equation model of large-scale atmospheric flow is derived. Attention is given to the classical balance equation and its ellipticity condition. Numerical solutions of the diagnostic system are presented, including examples of cases of the mixed elliptic-hyperbolic type and cases with non-zero divergence. Procedures for implementing such a system are outlined, along with a review of factors in using the technique for operational numerical weather prediction.

Ghil, M.

A statistical method for time-continuous assimilation of remote sounding temperatures

A time-continuous statistical method is presented for the four-dimensional assimilation of remote sounding temperatures based on radiometric measurements from polar-orbiting satellites. This method is applied to DST-6 data from the NOAA-4 and Nimbus-6 satellites. Experiments are reported in which the state of the atmosphere throughout the test period was determined using a varying amount of satellite data and in which different methods were used for their assimilation. Data from the NOAA-4 satellite only, from Nimbus-6 only, and from both satellites together were used; the methods tested include different variations of the statistical method as well as more traditional methods. The conclusions are that: (1) satellite-derived temperature data can have a modest, but statistically significant positive impact on numerical weather prediction in the two-to-three day range; (2) this impact is highly sensitive to the quantity of data available; and (3) the assimilation method plays a major role in the magnitude of the impact for the same data.

Ghil, M.

Climate stability for a Sellers-type model

We study a diffusive energy-balance climate model governed by a nonlinear parabolic partial differential equation. Three positive steady-state solutions of this equation are found; they correspond to three possible climates of our planet: an interglacial (nearly identical to the present climate), a glacial, and a completely ice-covered earth. We consider also models similar to the main one studied, and determine the number of their steady states. All the models have albedo continuously varying with latitude and temperature, and entirely diffusive horizontal heat transfer. The diffusion is taken to be nonlinear as well as linear. We investigate the stability under small perturbations of the main model's climates. A stability criterion is derived, and its application shows that the 'present climate' and the 'deep freeze' are stable, whereas the model's glacial is unstable. A variational principle is introduced to confirm the results of this stability analysis. For a sufficient decrease in solar radiation (about 2%) the glacial and interglacial solutions disappear, leaving the ice-covered earth as the only possible climate.

Ghil, M.

Wind laws for shockless initialization

A system of diagnostic equations for the velocity field, or wind laws, was derived for each of a number of models of large-scale atmospheric flow. The derivation in each case is mathematically exact and does not involve any physical assumptions not already present in the prognostic equations, such as nondivergence or vanishing of derivatives of the divergence. Therefore, initial states computed by solving these diagnostic equations should be compatible with the type of motion described by the prognostic equations of the model and should not generate initialization shocks when inserted into the model. Numerical solutions of the diagnostic system corresponding to a barotropic model are exhibited. Some problems concerning the possibility of implementing such a system in operational numerical weather prediction are discussed.

Ghil, M.

Steady-state solutions of a diffusive energy-balance climate model and their stability

A diffusive energy-balance climate model, governed by a nonlinear parabolic partial differential equation, was studied. Three positive steady-state solutions of this equation are found; they correspond to three possible climates of our planet: an interglacial (nearly identical to the present climate), a glacial, and a completely ice-covered earth. Models similar to the main one are considered, and the number of their steady states was determined. All the models have albedo continuously varying with latitude and temperature, and entirely diffusive horizontal heat transfer. The stability under small perturbations of the main model's climates was investigated. A stability criterion is derived, and its application shows that the present climate and the deep freeze are stable, whereas the model's glacial is unstable. The dependence was examined of the number of steady states and of their stability on the average solar radiation.

Ghil, M.