Dynamics of lava flow fronts, Arenal Volcano, Costa Rica

Arenal Volcano has effused basaltic andesite lava flows nearly continuously since September, 1968. The two different kinds of material in flows, lava and lava debris, have different rheologic properties and dynamic behavior. Flow morphology depends on the relationship between the amount and distribution of the lava and the debris, and to a lesser extent the ground morphology.Two main units characterize the flows: the channel zone and the frontal zone. The channel zone consists of two different units, the levées and the channel proper. A velocity profile in the channel shows a maximum value at the plug where the rate of shear is zero, and a velocity gradient increasing outward until, at the levées, the velocity becomes zero. Cooling produces a marked temperature gradient in the flow, leading to the formation of debris by brittle fracture when a critical value of shear rate to viscosity is reached. When the lava supply ceases, much of this debris and part of the lava is left behind after the flow nucleus drains out, forming a collapsed channel.Processes at the frontal zone include levée formation, debris formation, the change in shape of the front, and the choice of the flow path. These processes are controlled primarily by the rheological properties of the lava.Frontal zone dynamics can be understood by fixing the flow front as the point of reference. The lava flows through the channel into the front where it flows out into the levées, thereby increasing the length of the channel and permitting the front to advance. The front shows a relationship of critical height to the yield strength (τ₀) surface tension, and slope; its continued movement is activated by the pressure of the advancing lava in the channel behind. For an ideal flow (isothermal, homogeneous, and isotropic) the ratio of the section of channel proper to the section of levées is calculated and the distance the front will have moved at any time t_x can be determined once the amount of lava available to the front is known. Assuming that the velocity function of the front {G(t)} during the collapsing stage is proportional to the entrance pressure of the lava at the channel-front boundary, an exponential decrease of velocity through time is predicted, which shows good agreement with actual frontal velocity measurements taken on two flows. Local variations in slope have a secondary effect on frontal velocities.Under conditions of constant volume the frontal zone can be considered as a machine that consumes energy brought in by the lava to perform work (front advancement). While the front will use its potential energy to run the process, the velocity at which it occurs is controlled by the activation energy that enters the system as the kinetic energy of the lava flowing into the front. A relation for the energy contribution due to frontal acceleration is also derived. Finally the entrance pressure, that permits the front to deform, is calculated. Its small value confirms that the lava behaves very much like a Bingham plastic.     

Petrographic and geochemical evidence for the formation of primary, bacterially induced lacustrine dolomite: La Roda 'white earth' (Pliocene, central Spain)

Upper Pliocene dolomites (‘white earth’) from La Roda, Spain, offer a good opportunity to evaluate the process of dolomite formation in lakes. The relatively young nature of the deposits could allow a link between dolomites precipitated in modern lake systems and those present in older lacustrine formations. The La Roda Mg‐carbonates (dolomite unit) occur as a 3·5‐ to 4‐m‐thick package of poorly indurated, white, massive dolomite beds with interbedded thin deposits of porous carbonate displaying root and desiccation traces as well as local lenticular gypsum moulds. The massive dolomite beds consist mainly of loosely packed 1‐ to 2‐μm‐sized aggregates of dolomite crystals exhibiting poorly developed faces, which usually results in a subrounded morphology of the crystals. Minute rhombs of dolomite are sparse within the aggregates. Both knobbly textures and clumps of spherical bodies covering the crystal surfaces indicate that bacteria were involved in the formation of the dolomites. In addition, aggregates of euhedral dolomite crystals are usually present in some more clayey (sepiolite) interbeds. The thin porous carbonate (mostly dolomite) beds exhibit both euhedral and subrounded, bacterially induced dolomite crystals. The carbonate is mainly Ca‐dolomite (51–54 mol% CaCO₃), showing a low degree of ordering (degree of ordering ranges from 0·27 to 0·48). Calcite is present as a subordinate mineral in some samples. Sr, Mn and Fe contents show very low correlation coefficients with Mg/Ca ratios, whereas SiO₂ and K contents are highly correlated. δ¹⁸O‐ and δ¹³C‐values in dolomites range from ?3·07‰ to 5·40‰ PDB (mean=0·06, σ=1·75) and from ?6·34‰ to ?0·39‰ PDB (mean=?3·55, σ=1·33) respectively. Samples containing significant amounts of both dolomite and calcite do not in general show significant enrichment or depletion in ¹⁸O and ¹³C between the two minerals. The correlation coefficient between δ¹⁸O and δ¹³C for dolomite is extremely low and negative (r=?0·05), whereas it is higher and positive (r=0·47) for calcite. The lacustrine dolomite deposit from La Roda is interpreted mainly as a result of primary precipitation of dolomite in a shallow, hydrologically closed perennial lake. The lake was supplied by highly saturated HCO₃^?/CO₃^2? groundwater that leached dolomitic Mesozoic formations. Precipitation of dolomite from alkaline lake waters took place under a semi‐arid to arid climate. However, according to our isotopic data, strong evaporative conditions were not required for the formation of the La Roda dolomite. A significant contribution by bacteria to the formation of the dolomites is assumed in view of both petrographic and geochemical evidence.     

Optimum damping in linear isolation systems

Optimum isolation damping for minimum acceleration response of base-isolated structures subjected to stationary random excitation is investigated. Three linear models are considered to account for the energy dissipation mechanism of the isolation system: a Kelvin element, a linear hysteretic element and a standard solid linear element, commonly used viscoelastic models for isolation systems comprising natural rubber bearings and viscous dampers. The criterion selected for optimality is the minimization of the mean-square floor acceleration response. The effects of the frequency content of the excitation and superstructure properties on the optimum damping and on the mean-square acceleration response are addressed. The study basically shows that the attainable reduction in the floor acceleration largely depends on the energy dissipation mechanism assumed for the isolation system as well as on the frequency content of the ground acceleration process. Special care should be taken in accurately modelling the mechanical behaviour of the energy dissipation devices.     

CFD simulation of airflow over a regular array of cubes. Part II: analysis of spatial average properties

In the first part of this study, results of a computational fluid dynamics simulation over an array of cubes have been validated against a set of wind-tunnel measurements. In Part II, such numerical results are used to investigate spatially-averaged properties of the flow and passive tracer dispersion that are of interest for high resolution urban mesoscale modelling (e.g. non resolved obstacle approaches). The results show that vertical profiles of mean horizontal wind are linear within the canopy and logarithmic above. The drag coefficient, derived from the numerical results using the classical formula for the drag force, is height dependent (it decreases with height). However, a modification of the formula is proposed (accounting for subgrid velocity scales) that makes the drag coefficient constant with height. Results also show that the dispersive fluxes are similar in magnitude to the turbulent fluxes, and that they play a very important role within the canopy. Vertical profiles of turbulent length scales (to be used in k–l closure schemes, where k is the turbulent kinetic energy and l a turbulent length scale) are also derived. Finally the distribution of the values around the mean over the reference volumes are analysed for wind and tracer concentrations.     

On the Intensity Thresholds of the Network and Plage Regions

A period of minimum solar activity in April and May 1996 was used to analyze full-disk CaK spectroheliograms with the aim of deriving the intensity thresholds of the quiet network and the plages by applying the so-called inflexion point method. The average network intensity threshold is found to be influenced mainly by the seeing, whereas the average plage threshold shows an increasing trend from disk center towards the solar limb. Both parameters are compared with the results of other authors.     

Relocation of shallow earthquakes in southern Alaska using Joint Hypocenter Determination method

A subset of 2660 shallow earthquakes (0–50 km) that occurred from 1988 to 1996 in south central Alaska between 155°W and 145°W and 59°N and 63°N was relocated using the joint hypocenter determination (JHD) method. Both P- and S-wave observations recorded by the regional seismic network were used. Events were relocated in twenty different groups based on their geographic location and depth using two velocity models. As a result of the relocation, the majority of the hypocenters shifted downward, while the epicenter locations did not change significantly. The distribution of the shallow subduction zone earthquakes indicates the existence of two seismically independent blocks, with one block occupying the northeastern part and the other occupying the central and western parts of the study area. The boundary between the blocks is marked by a 15 to 20 km wide seismicity gap to the southeast of 149.5°W and 62°N. The analysis of the fault plane solutions for shallow subduction zone earthquakes shows that an overwhelming majority of the solutions represent normal, oblique-normal or strike-slip faulting with predominant WNW-ESE orientation of T-axes. This indicates a down-dip extensional regime for the subducting slab at shallow depths. Very few earthquakes yielded fault plane solutions consistent with thrusting on a contact zone between the overriding and subducting plates. This result may be an indication that currently either the strain energy is not released at the contact zone or it is associated with aseismic motion.     

The radioracemization of isovaline. Cosmochemical implications

The optically pure d- and l-enantiomers of isovaline (I), which cannot be racemized by ordinary chemical mechanisms involving α-hydrogen removal, and which has been isolated in apparently racemic form from the Murchison meteorite, have been subjected to partial radiolysis by the ionizing radiation from a 3000 Ci ⁶⁰Co γ-ray source. Both in the anhydrous and hydrated solid states and as solid sodium or hydrochloride salts each enantiomer suffered significant radioracemization of the undestroyed residue during its partial radiolysis. The sodium salt of isovaline in 0.1 M aqueous solution suffered extensive radiolysis with relatively small radiation doses, but showed no detectable radioracemization. The significance of these observations with respect to the primordial enantiomeric composition of the isovaline (and other amino acids) indigenous to meteroties is discussed.     

Levels of chlorinated hydrocarbons and trace metals in bivalves and nearshore sediments from the Dominican Republic

Bivalves (Crassostrea rizhophorae, Isognomon alatus, Codakia orbicularis Linné, and Tellina fausta) and sediment samples collected from 12 locations along the Dominican Republic coastline were used to assess the degree of nearshore contamination by selected chlorinated hydrocarbons and heavy metals. Concentrations of total chlordanes, DDTs and PCBs ranged from 0.51 to 7.47 ng/g, < MDL to 30.9 ng/g and 11.3 to 82.3 ng/g, and from < MDL to 1.73 ng/g, 0.21 to 12.5 ng/g and 0.46 to 41.9 ng/g in bivalves and surface sediments, respectively. The concentration ranges of predominantly anthropogenic Cu, Pb and Zn were 3.08–866 μg/g, 0.08–1.46 μg/g and 22.9–4380 μg/g in bivalve soft parts and 1.01–111 mg/g, 0.420–134 μg/g and 1.96–277 μg/g in surface sediments, respectively. Mercury concentrations ranged from 0.290 to 7.02 μg/g and from 0.096 to 0.565 μg/g in both matrices. The relatively high Hg concentrations found in bivalves from two stations in Bahia Samana, near to a known source of past Hg contamination in the sediments, require a more in-depth survey of the area. In general the concentrations of selected chlorinated hydrocarbons and trace metals, compared with those reported for the same matrices in the Caribbean and west Atlantic coastal areas, show at present the existence of only very localized contamination problems in the Dominican Republic.     

Geodetic constraints on active tectonics of the Western Mediterranean: Implications for the kinematics and dynamics of the Nubia-Eurasia plate boundary zone

We present GPS observations in Morocco and adjacent areas of Spain from 15 continuous (CGPS) and 31 survey-mode (SGPS) sites extending from the stable part of the Nubian plate to central Spain. We determine a robust velocity field for the W Mediterranean that we use to constrain models for the Iberia-Nubia plate boundary. South of the High Atlas Mountain system, GPS motions are consistent with Nubia plate motions from prior geodetic studies. We constrain shortening in the Atlas system to <1.5 mm/yr, 95% confidence level. North of the Atlas Mountains, the GPS velocities indicate Nubia motion with respect to Eurasia, but also a component of motion normal to the direction of Nubia-Eurasia motion, consisting of southward translation of the Rif Mountains in N Morocco at rates exceeding 5 mm/yr. This southward motion appears to be directly related to Miocene opening of the Alboran Sea. The Betic Mountain system north of the Alboran Sea is characterized by WNW motion with respect to Eurasia at ～1–2 mm/yr, paralleling Nubia-Eurasia relative motion. In addition, sites located in the Betics north of the southerly moving Rif Mountains also indicate a component of southerly motion with respect to Eurasia. We interpret this as indicating that deformation associated with Nubia-Eurasia plate motion extends into the southern Betics, but also that the Betic system may be affected by the same processes that are causing southward motion of the Rif Mountains south of the Alboran Sea. Kinematic modeling indicates that plate boundary geometries that include a boundary through the Straits of Gibraltar are most compatible with the component of motion in the direction of relative plate motion, but that two additional blocks (Alboran-Rif block, Betic Mountain block), independent of both Nubia and Eurasia are needed to account for the motions of the Rif and Betic Mountains normal to the direction of relative plate motion. We speculate that the southward motions of the Alboran-Rif and Betic blocks may be related to mantle flow, possibly induced by southward rollback of the subducted Nubian plate beneath the Alboran Sea and Rif Mountains.