Finite-element analysis of the transport of water and solutes in title-drained soils

The finite-element method based on a Galerkin technique was used to formulate the problem of simulating the two-dimensional (cross-sectional) transient movement of water and solute in saturated or partially saturated nonuniform porous media. The numerical model utilizes linear triangular elements. Nonreactive, as well as reactive solutes whose behaviour can be described by a distribution coefficient or first-order reaction term were considered. The flow portion of the model was tested by comparison of the model results with experimental and finite-difference results for transient flow in an unsaturated sand column and the solute transport portion of the model was tested by comparison with analytical solution results. The model was applied to a hypothetical case involving movement of water and solutes in tile-drained soils. The simulation results showed the development of distinct solute leaching patterns in the soil as drainage proceeded. Although applied to a tile drainage problem in this study, the model should be equally useful in the study of a wide range of two-dimensional water and solute migration problems.     

Melanoidins—kerogen precursors and geochemical lipid sinks: a study using pyrolysis gas chromatography (PGC)

Melanoidins, acidic polymeric products of amino acid/sugar condensation reactions, are shown to be capable of reacting with and binding functionalized lipid molecules. Pyrolysis gas chromatography is used to monitor the inclusion of lipids (alcohols and acids) into, and to monitor the changes in composition caused by artificial diagenesis of, the melanoidins. Artificial diagenesis of the lipidcontaining melanoidins produces a material resembling kerogen. The possible role of melanoidins in the formation of kerogens is discussed as are internal transformations within kerogens during diagenesis.     

The influence of minerals on the pyrolysis of kerogens

The effect of mineral matter on the laboratory pyrolysis of sediments, kerogens and coals, and of coal macerals mixed with either alumina, bentonite, kaolinite, or calcium carbonate has been investigated. Some minerals are more active than others in effecting changes in the composition of the pyrolysate. The relative content of low molecular weight pyrolysis products is higher for kerogens pyrolyzed in a mineral matrix than for isolated kerogens. Our limited data suggests that this is due to both condensation and gasification of higher molecular weight constituents in the ‘primary’ pyrolysate. Differences in the content of aromatic versus aliphatic compounds have been noted when pyrograms from coal macerals are compared with pyrograms from coal maceral-mineral mixtures. We conclude that mineralogy is important in controlling the composition of kerogen pyrolysis products.     

Geometric similarity for the temperature field in the unstable atmospheric surface layer

Davenport's geometric similarity is applied to the temperature field in the unstable atmospheric surface layer for displacements in all three orthogonal directions. The decay constants for lateral displacements are over 50 times as large as for vertical displacements. Such large lateral decay constants mean that the measurement of longitudinal constants will be very difficult.Some of the assumptions inherent in geometric similarity were not satisfied. Two regimes, plume scales and larger scales, having different decay constants, were observed. The possibility of a similar scale dependency of the decay constants for the wind field is discussed.     

Hemispheric simulation of the Asian summer monsoon

A three-level, -plane, filtered model is used to simulate the Northern Hemisphere summer monsoon. A time-averaged initial state, devoid of sub-planetary scale waves, is integrated through 30 days on a 5° latitude-longitude grid. Day 25 through day 30 integrations are then repeated on a 2.5° grid. The planetary-scale waves are forced by time-independent, spatially varying diabatic heating. Energy is extracted via internal and surface frictional processes. Orography is excluded to simplify synoptic-scale energy sources.During integration the model energy first increases, but stabilizes near day 10. Subsequent flow patterns closely resemble the hemisphere summer monsoon. Climatological features remain quasi-stationary. At 200 mb high pressure dominates the land area, large-scale troughs are found over the Atlantic and Pacific Oceans, the easterly jet forms south of Asia, and subtropical jets develop in the westerlies. At 800 mb subtropical highs dominate the oceans and the monsoon trough develops over the Asian land mass. The planetary scales at all levels develop a realistic cellular structure from the passage of transient synoptic-scale features, e.g., a baroclinic cyclone track develops near 55°N and westward propagating waves form in the easterlies.Barotropic redistribution of kinetic energy is examined over a low-latitude zonal strip using a Fourier wave-space. In contrast to higher latitudes where the zonal flow and both longer and shorter waves are fed by barotropic energy redistribution from the baroclinically unstable wavelengths, the low-latitude waves have a planetary-scale kinetic energy source. Wave numbers 1 and 2 maintain both the zonal flow and all shorter scales via barotropic transfers. Transient and standing wave processes are examined individually and in combination.Wave energy accumulates at wave numbers 7 and 8 at 200 mb and at wave number 11 in the lower troposphere. The 800-mb waves are thermally indirect and in the mean they give energy to the zonal flow. These characteristics agree with atmospheric observation. The energy source for these waves is the three wave barotropic transfer. The implications of examining barotropic processes in a Fourier wave-space, vice the more common approach of separating the flow into a mean plus a deviation are discussed.     

Stress in the lithosphere: Inferences from steady state flow of rocks

Mechanical data and flow processes from steady state deformation experiments may be used to infer the state of stress in the lithosphere and asthenosphere. Extrapolations of flow equations to a representative geologic strain rate of 10^–14/sec. for halite, marble, quartzite, dolomite, dunite and enstatolite are now warranted because the steady state flow processes in the experiments are identical to those in rocks and because the geotherms are reasonably well established. More direct estimates are obtained from free dislocation densities, subgrain sizes and recrystallized grain sizes all of which are functions only of stress. Using the last of these techniques, we have estimated stress profiles as a function of depth from xenoliths in basalts and kimberlites, whose depths of equilibration were determined by pyroxene techniques, from four different areas of subcontinental and suboceanic upper mantle. The results are similar and indicate stress differences of about 200 to 300 bars at 40 to 50 km, decaying to a few tens of bars at depths betow 100 km. These stresses are reasonable and are in accord with extrapolations of the mechanical data provided that allowance is made for a general increase in strain rate and decrease in viscosity with depth.