The application of cellular automata to the erosion of landforms

Two approaches to modelling surface erosion are described which can be used to explain the development of geomorphological features. The first approach uses non-linear wave theory and shows how surface gradient discontinuities can arise as the erosion of an initially smooth surface proceeds. For isotropic erosion, surface shapes may be predicted by a simple application of the Huyghens wavefront construction and Snell's law at media boundaries. A second approach uses a computational cellular method. The general shapes predicted by the computer simulations are shown to agree with those determined from the wave approach. The application of the cellular method to problems of geomorphological interest are discussed.     

Rainfall estimation using raingages and radar — A Bayesian approach: 1. Derivation of estimators

Procedures for estimating rainfall from radar and raingage observations are constructed in a Bayesian framework. Given that the number of raingage measurements is typically very small, mean and variance of gage rainfall are treated as uncertain parameters. Under the assumption that log gage rainfall and log radar rainfall are jointly multivariate normal, the estimation problem is equivalent to lognormal co-kriging with uncertain mean and variance of the gage rainfall field.The posterior distribution is obtained under the assumption that the prior for the mean and inverse of the variance of log gage rainfall is normal-gamma 2. Estimate and estimation variance do not have closed-form expressions, but can be easily evaluated by numerically integrating two single integrals. To reduce computational burden associated with evaluating sufficient statistics for the likelihood function, an approximate form of parameter updating is given. Also, as a further approximation, the parameters are updated using raingage measurements only, yielding closed-form expressions for estimate and estimation variance in the Gaussian domain.With a reduction in the number of radar rainfall data in constructing covariance matrices, computational requirements for the estimation procedures are not significantly greater than those for simple co-kriging. Given their generality, the estimation procedures constructed in this work are considered to be applicable in various estimation problems involving an undersampled main variable and a densely sampled auxiliary variable.     

地面对晴空大气热辐射结构的影响

本文利用1987年美国FIFE资料及1986年6月TIPMEX资料,并利用大气长、短波一维辐射传输数值模式,计算了晴空大气状态下地面辐射参数(地面热辐射温度、地表反照率等)对大气长波冷却率、短波加热率廓线的影响。结果表明,地球表面对大气边界层内的热辐射结构影响很大,特别是在高原地区这种影响尤其显著。     

Kinematic interpretation of remnant en-echelon rift segments of the southwestern Japan Sea

The southwestern Japan Sea and the margin of southwestern Japan feature parallel elongate extensional sub-basins and grabens. These structures are oblique to the trend of southwestern Japan and are interpreted as the remnant of an en-echelon array. Such a right-stepping remnant en-echelon pattern indicates sinistral rifting of the southwestern Japan block from Asia. The major extensional fault trends are parallel over a distance of 600 km. Infinitesimal kinematic relations show that it is unlikely that these faults initiated in response to rotation about the nearby pole determined from paleomagnetic studies. Rather, an earlier stage of rifting about a distant pole is indicated. The rotational pole of this early rifting was probably located at least 6000 km away from the rift and the pole of rotation was offset from the rift trend by a similar distance. Rifting about this distant rotational pole loosened the southwestern Japan continental block from the Asian mainland making it susceptible to a large rotation about a nearby pole.     

Glacial meltwater cooling of the Gulf of Mexico: GCM implications for Holocene and present-day climates

An atmospheric general circulation model, the NCAR CCM, has been used to investigate the possible effects that reduced Gulf of Mexico sea surface temperatures (SST) could have on regional and hemispheric climates. ¹⁸O records and terrestrial evidence indicate at least two major glacial meltwater discharges into the Gulf of Mexico subsequent to the last glacial maximum. It is probable that these discharges reduced Gulf of Mexico SST. We have conducted three numerical experiments, with imposed gulf-wide SST coolings of 3°C, 6°C, and 12°C, and find in all three experiments significant reductions in the North Atlantic storm-track intensity, along with a strong decrease in transient eddy water vapor transport out of the Gulf of Mexico. Surface pressures are higher over the North Atlantic, indicating a reduction of the climatological Icelandic low. The region is generally cooler and drier, with a reduction in precipitation that agrees well with evidence from Greenland ice cores. Other statistically significant changes occur across the Northern Hemisphere, but vary between the three experiments. In particular, warmer, wetter conditions are found over Europe for both the 6°C and 12°C SST reductions, but cooler conditions are found for the 3°C reduction. This indicates a dependence, in both the sign and magnitude of the model response, on the magnitude of the imposed SST anomaly. The results suggest that the present-day North Atlantic storm track is dependent on warm Gulf of Mexico SST for much of its intensity. They also suggest that meltwater-induced coolings may help account, in part, for some of the climatic oscillations that occurred during the last glacial/interglacial transition.     

A note on the characteristics and possible origins of desert varnishes from Southeast Morocco. A reply to comments by R. I. Dorn

This reply to Dorn (1989) argues that a biotic origin is not the only way in which desert varnishes may be formed. It is suggested that a variety of initial conditions may result in similar end products. It is agreed that the dissimilar views of the authors may result from their different philosophies concerning the classification of Mn/Fe varnishes.     

An intercomparison of seawater filtration procedures

Filtration systems, used for the preparation of coastal seawater samples for trace analysis, have been intercompared to determine the extent of contamination arising during the filtration process. The results show that all the systems intercompared are able to provide uncontaminated filtrates for copper, nickel, cobalt and manganese and, with two exceptions, for cadmium and iron, at levels that would be encountered in uncontaminated coastal waters.     

Convective instabilities in a variable viscosity fluid cooled from above

Whether in the mantle or in magma chambers, convective flows are characterized by large variations of viscosity. We study the influence of the viscosity structure on the development of convective instabilities in a viscous fluid which is cooled from above. The upper and lower boundaries of the fluid are stress-free. A viscosity dependence with depth of the form ν₀ + ν₁ exp(?γ.z) is assumed. After the temperature of the top boundary is lowered, velocity and temperature perturbations are followed numerically until convective breakdown occurs. Viscosity contrasts of up to 10⁷ and Rayleigh numbers of up to 10⁸ are studied.For intermediate viscosity contrasts (around 10³), convective breakdown is characterized by the almost simultaneous appearance of two modes of instability. One involves the whole fluid layer, has a large horizontal wavelength (several times the layer depth) and exhibits plate-like behaviour. The other mode has a much smaller wavelength and develops below a rigid lid. The “whole layer” mode dominates for small viscosity contrasts but is suppressed by viscous dissipation at large viscosity contrasts.For the “rigid lid” mode, we emphasize that it is the form of the viscosity variation which determines the instability. For steep viscosity profiles, convective flow does not penetrate deeply in the viscous region and only weak convection develops. We propose a simple method to define the rigid lid thickness. We are thus able to compute the true depth extent and the effective driving temperature difference of convective flow. Because viscosity contrasts in the convecting region do not exceed 100, simple scaling arguments are sufficient to describe the instability. The critical wavelength is proportional to the thickness of the thermal boundary layer below the rigid lid. Convection occurs when a Rayleigh number defined locally exceeds a critical value of 160–200. Finally, we show that a local Rayleigh number can be computed at any depth in the fluid and that convection develops below depth z_r (the rigid lid thickness) such that this number is maximum.The simple similarity laws are applied to the upper mantle beneath oceans and yield estimates of 5 × 10¹⁵?5 × 10¹⁶ m² s^?1 for viscosity in the thermal boundary layer below the plate.