Tarawera 1886, New Zealand — A basaltic plinian fissure eruption

New Zealand's biggest and most destructive volcanic eruption of historical times was that of Tarawera in 1886. The resulting scoria fall has a dispersal very similar in extent to that of the Vesuvius A.D. 79 pumice fall and is one of the few known examples of a basaltic deposit of plinian type. A new estimate of the volume (2 km³) is significantly greater than previous estimates. The basalt came mainly from a 7-km length of fissure, and emission and exit velocity were fairly uniform along at least 4 km of it, this is one of the few documented examples of a plinian eruption from a fissure vent. Primary welding of the scoria fall resulted where the accumulation rate exceeded about 250 mm min⁻¹. A model of the eruption dynamics is proposed which leads to an estimate of 28 km for the height of the eruption cloud and implies a magma volatile fraction of 1.5–3%. Violent phreatic explosions occurred in the southwestern extension of the fissure across the Rotomahana geothermal field, and it is thought that some of the water responsible for the power of the plinian eruption came from this source, though its amount was not sufficient to turn the eruption into a phreatoplinian one.     

Hydrogen speciation in synthetic quartz

The dominant hydrogen impurity in synthetic quartz is molecular H₂O. H-OH groups also occur, but there is no direct evidence for the hydrolysis of Si-O-Si bonds to yield Si-OH HO-Si groups. Molecular H₂O concentrations in the synthetic quartz crystals studied range from less than 10 to 3,300 ppm (H/Si), and decrease smoothly by up to an order of magnitude with distance away from the seed. OH^? concentrations range from 96 to 715 ppm, and rise smoothly with distance away from the seed by up to a factor of three. The observed OH^? is probably all associated with cationic impurities, as in natural quartz. Molecular H₂O is the dominant initial hydrogen impurity in weak quartz. The hydrolytic weakening of quartz may be caused by the transformation H₂O + Si-O-Si → 2SiOH, but this may be a transitory change with the SiOH groups recombining to form H₂O, and the average SiOH concentration remaining very low. Synthetic quartz is strengthened when the H₂O is accumulated into fluid inclusions and cannot react with the quartz framework.     

Mixed valence of iron in minerals with cation clusters

The valence and distribution of iron in vivianite, lazulite, babingtonite, rockbridgeite, acmite, aegirine-augite, hedenbergite, and ilvaite were studied with optical and Mössbauer spectroscopy. Optically activated intervalence charge transfer between Fe²⁺ and Fe³⁺ in neighboring sites through common edges or faces is observed in all these minerals irrespective of the polymerization of the iron-oxygen polyhedra ranging from finite clusters to infinite structural units. However, a distinct decrease occurs in the energy of the corresponding optical absorption band with increasing number of Fe²⁺ and Fe³⁺ ions involved in the charge transfer process. Thermally activated electron delocalization between Fe²⁺ and Fe³⁺ occurs only if Fe²⁺ and Fe³⁺ occupy crystallographically equivalent or geometrically very similar neighboring sites which share common edges to form extended structural units such as the ribbon in ilvaite. If the Fe-O polyhedra form finite clusters of two, three, or four polyhedra (e.g., in vivianite, lazulite, and babingtonite, respectively) no thermally-activated mixed-valence states of iron are observed. In aegirine, extended regions of the M1 chain are statistically occupied by Fe²⁺ and Fe³⁺ giving rise to thermally-activated electron delocalization in addition to the intervalence band in the optical absorption spectrum. The intensity of the optical intervalence absorption has been measured in a number of systems: ? values range from 60 to 210.     

Numerical simulation and analysis of the mean coastal circulation off California

A two-dimensional numerical model is applied to a coastal ocean wherein alongshore elevation and density gradients, normally calculated by a three-dimensional model, are instead supplied by climatologically averaged data for the California Current System between 25 and 40°N. Surface wind stress is also obtained from climatological data. Both surface and bottom boundary layers are resolved in the model calculations; a second moment turbulence closure submodel supplies vertical diffusivities. Near steady state solutions are possible when surface buoyancy flux is imposed at the surface.Model results are as follows: Southward wind stress produces a broad equatorward current with an embedded coastal jet in accordance with previous studies. Positive wind stress curl reduces the jet current and produces a poleward undercurrent which then surfaces as the curl is increased. The jet currents are reduced and poleward flow increases as bottom steepness increases; to a lesser extent, inclusion of the beta effect has a similar effect. The existence of near bottom, poleward or equatorward flow is explained rather simply in terms of the bottom stress resulting from the alongshore balance of surface wind stress and vertically integrated pressure gradient, the latter involving the alongshore surface elevation and density gradient. A further finding is that the upwelling circulation associated with wind stress is confined to the top 200 to 300 m of the ocean along the California coast.     

Efficiency of kriging estimation for square, triangular, and hexagonal grids

Although several researchers have pointed out some advantages and disadvantages of various soil sampling designs in the presence of spatial autocorrelation, a more detailed study is presented herein which examines the geometrical relationship of three sampling designs, namely the square, the equilateral triangle, and the regular hexagon. Both advantages and disadvantages exist in the use of these designs with respect to estimation of the semivariogram and their effect on the mean square error or variance of error. This research could be used to design optimal sampling strategies; it is based on the theory of regionalized variables, in which the intrinsic hypothesis is satisfied. Among alternative designs, an equilateral triangle design gives the most reliable estimate of the semivariogram. It also gives the minimum maximum mean square error of point estimation of the concentration over the other two designs for the same number of measurements when the nugget effect is small relative to the variance. If the nugget effect is large (.90 ² or more), and the linear sampling density is >0.85r where r is the range, the hexagonal design is best. This study computes and compares the maximum mean square error for each of these designs.     

Stable-channel design in alluvial rivers

A theory of stable-channel design is deduced from the condition that, for any given steady flow, total sediment concentration maintains a balanced average in space and time within a stable channel. This condition for stability applies during changes in channel position but does not constrain them. A stability index is defined in terms of concentration of bedmaterial load: a function of steady, uniform or quasi-uniform water discharge, channel bedslope, water surface width and a typical diameter of cohesionless, granular bedmaterial. The index is constant only within a stable channel, and its mean magnitude may be determined from observations at cross-sections of more or less straight, wide reaches. Two distinct design methods, threshold and mobile-bed, are deduced from the theory. Design formula are presented for both methods as functions of the channel stability index, and are applicable to cross-sections or reaches, independent of planform, of both sand-silt and gravel-bed rivers. A relation describing variation along a channel in the ratio of bedload to suspended-sediment discharge, for a specified steady flow, is also inferred from the stability condition.     

Effect of a kelp forest on coastal currents

Ocean currents supply a kelp ecosystem with nutrients, planktonic food, and larvae. We have found that these currents in a kelp forest (Macrocystis pyrifera) are slower than currents outside. At the Pt. Loma, San Diego, California, site that we studied, current velocities were about a third of those outside. A comparison of frequency spectra shows that semi-diurnal frequencies are preferentially passed by the kelp. This effect of a kelp forest on the currents that nurture it is similar to that of a terrestrial forest on local winds.     

A world model of soil carbon dioxide

Mean growing season soil PCO₂ data were obtained for 19 regions of the world in nine countries. Bivariate and multiple linear regression analysis with soil log(PCO₂) as the dependent variable and TEMP, PRECIP, log(AET), and log(PET) as the four climatic independent variables demonstrated that AET was the best independent predictor of soil PCO₂. An improved soil PCO₂-AET model was developed by assuming (1) that as AET approaches zero, soil PCO₂ approaches the atmospheric value and (2) that there is an upper limit to soil PCO₂ at very high AET. This model has the form log(PCO₂) = ?3·47 + 2·09 (1 ?e^{?0·0172 AET}) where AET is in mm. It explains 67 per cent of the initial variation in the soil PCO₂ data, predicts a soil log(PCO₂) of ? 3·47 at AET = 0, and an upper limit of 3·5 per cent (log(PCO₂) = ? 1·45) for mean growing season soil PCO₂ at AET values of 2000 mm and above. The results of this study suggest that soil PCO₂ levels in tropical areas are, on average, higher than those in temperate, alpine, and arctic regions.