On the determination of the closure parameters in higher-order closure models

In higher-order closure models at least the pressure redistribution and the dissipation of the turbulent kinetic energy and the temperature variance have to be parameterized. Due to this, the introduction of proportionality coefficients — the so-called closure parameters — is forced, which have to be determined before the model is used. We compare a group of models which use the return-to-isotropy hypothesis (Rotta, 1951) to describe the pressure redistribution and assume local isotropy for the smallest eddies in order to parameterize the dissipation. Special concern is given to the method of Mellor and Yamada (1982). Some of the closure parameters are re-derived on the basis of sensitivity studies requiring that both shear production and buoyancy behave in a realistic way if pressure redistribution or dissipation is changed by varying the closure parameters. This set of parameters is compared with those obtained by fitting to experimental data, by use of the Monin-Obukhov similarity theory and by considering ratios of variances, covariances and mean flow gradients, respectively. It is shown that the various sets of closure parameters are at least of the same order. The differences give some insight into the advantages and disadvantages of the various determination procedures. However, the general accordance of the different parameter sets supports the assumption of universal constants.     

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.     

Coastal processes and environmental hazards: the Buenos Aires (Argentina) and Venetian (Italy) littorals

The Buenos Aires (Argentina) and Venice (Italy) coastlands have experienced significant saltwater contamination of the phreatic aquifer, coastal erosion, hydrodynamic changes and relative sea level rise processes due to natural and man-induced factors. These factors expose coastal areas to morpho-hydro-geological hazards, such as soil desertification, frequency and degree of flooding, littoral erosion, and the silting of river mouths and channels. Man-made interventions and actions, such as beach mining, construction of coastal structures and exploitation of aquifers without an adequate knowledge of the hydrology setting and an adequate management program, worsen these natural hazards. Uncontrolled human activity induces environmental damage to the overall coastal plains. The coastal plains play an important role in the social/economic development of the two regions based on land use, such as agriculture, horticulture, breeding, and tourism, as well as industry. Results of investigations on saltwater contamination, sea level rise and morphological changes recently performed in these two coastal areas are presented here.     

Ozone Production Efficiency in Savanna and Forested Areas during the EXPRESSO Experiment

NO, NO_x, NO_y and O₃ have been measuredduring the airborne EXPRESSO experiment, November 96,which took place near Bangui (Central Africa) at thebeginning of the dry season. This period correspondsto an intense burning activity. Chemical andphotochemical characteristics of the planetaryboundary layer, which corresponds most of the time tothe monsoon layer, and the Harmattan layer, which issituated above the latter, have been studied forsavanna as well as rain forest areas. These two layersare very different when considering the chemicalcomposition (especially for ozone and NO_z) andthe photochemical age.The relationship of photochemical ozone productionversus photochemical NO_x oxidation products hasbeen investigated. Results indicate an ozoneproduction efficiency (OPE) ranging from 6.3 to 14.8in the planetary boundary layer. Thus, this layer ischaracteristic of a photochemically young and activeair mass. In this layer, the ozone potentialproduction increases with the air mass photochemicalage. On the other hand, the Harmattan layer shows anOPE ranging from 4.6 to 7.4. These values arecharacteristic of photochemically well-aged airmasses. In this layer, the ozone potential productionseems to be exhausted with values around 4.0 (i.e., 4ozone molecules produced for each NO_x moleculeemitted).     

The Cloud Ice Mountain Experiment (CIME) 1998: experiment overview and modelling of the microphysical processes during the seeding by isentropic gas expansion

The second field campaign of the Cloud Ice Mountain Experiment (CIME) project took place in February 1998 on the mountain Puy de Dôme in the centre of France. The content of residual aerosol particles, of H₂O₂ and NH₃ in cloud droplets was evaluated by evaporating the drops larger than 5 μm in a Counterflow Virtual Impactor (CVI) and by measuring the residual particle concentration and the released gas content. The same trace species were studied behind a round jet impactor for the complementary interstitial aerosol particles smaller than 5 μm diameter. In a second step of experiments, the ambient supercooled cloud was converted to a mixed phase cloud by seeding the cloud with ice particles by the gas release from pressurised gas bottles. A comparison between the physical and chemical characteristics of liquid drops and ice particles allows a study of the fate of the trace constituents during the presence of ice crystals in the cloud.In the present paper, an overview is given of the CIME 98 experiment and the instrumentation deployed. The meteorological situation during the experiment was analysed with the help of a cloud scale model. The microphysics processes and the behaviour of the scavenged aerosol particles before and during seeding are analysed with the detailed microphysical model ExMix. The simulation results agreed well with the observations and confirmed the assumption that the Bergeron–Findeisen process was dominating during seeding and was influencing the partitioning of aerosol particles between drops and ice crystals. The results of the CIME 98 experiment give an insight on microphysical changes, redistribution of aerosol particles and cloud chemistry during the Bergeron–Findeisen process when acting also in natural clouds.     

Behaviour of H2O2, NH3, and black carbon in mixed-phase clouds during CIME

We investigated the partitioning of trace substances during the phase transition from supercooled to mixed-phase cloud induced by artificial seeding. Simultaneous determination of the concentrations of H₂O₂, NH₃ and black carbon (BC) in both condensed and interstitial phases with high time resolution showed that the three species undergo different behaviour in the presence of a mixture of ice crystals and supercooled droplets. Both H₂O₂ and NH₃ are efficiently scavenged by growing ice crystals, whereas BC stayed predominantly in the interstitial phase. In addition, the scavenging of H₂O₂ is driven by co-condensation with water vapour onto ice crystals while NH₃ uptake into the ice phase is more efficient than co-condensation alone. The high solubility of NH₄⁺ in the ice could explain this result. Finally, it appears that the H₂O₂–SO₂ reaction is very slow in the ice phase with respect to the liquid phase. Our results are directly applicable for clouds undergoing limited riming.     

Free cupric ion concentrations and Cu complexation in selected Swiss lakes and rivers

[Cu²⁺] and Cu complexation parameters in some selected freshwater systems in Switzerland were determined by the technique of ligand-exchange and DPCSV. Results from the water columns of some eutrophic and oligotrophic lakes are presented and compared to small acid lakes. Cu is strongly complexed by organic ligands which with very high stability constants at low concentrations are probably biologically produced, as indicated by the seasonal variations in the eutrophic lakes and by the relationship between Cu complexation and algal activity in the eutrophic (pCu=15–16), oligotrophic (pCu=13–14) and acidic (pCu=9–10) lakes. The extent of Cu complexation in river waters was generally lower than in the eutrophic lakes, at similar DOC levels. No obvious correlation between Cu complexation and DOC was observed, indicating that Cu complexing ligands are specific organic compounds.     

A viscoelastic shear zone model of compressional and extensional plate boundaries

A model is proposed describing the mechanical evolution of a shear zone along compressional and extensional plate boundaries, subject to constant strain rate. The shear zones are assumed as viscoelastic with Maxwell rheology and with elastic and rheological parameters depending on temperature and petrology. Stress and strain are computed as functions of time and depth. For both kinds of boundaries the model reproduces the existence of a shallow seismogenic zone, characterized by a stress concentration. The thickness of the seismogenic layer is evaluated considering the variations of shear stress and frictional strength on faults embedded in the shear zone. Assuming that a fault dislocation takes place, the brittle-ductile transition is assumed to occur at the depth at which the time derivative of total shear stress changes from positive to negative values. The effects of different strain rates and geothermal gradients on the depth of the brittle-ductile transition are studied. The model predictions are consistent with values inferred from seismicity data of different boundary zones.