The Simulation of the Opposing Fluxes of Latent Heat and CO<Subscript>2</Subscript> over Various Land-Use Types: Coupling a Gas Exchange Model to a Mesoscale Atmospheric Model

A mesoscale meteorological model (FOOT3DK) is coupled with a gas exchange model to simulate surface fluxes of CO₂ and H₂O under field conditions. The gas exchange model consists of a C3 single leaf photosynthesis sub-model and an extended big leaf (sun/shade) sub-model that divides the canopy into sunlit and shaded fractions. Simulated CO₂ fluxes of the stand-alone version of the gas exchange model correspond well to eddy-covariance measurements at a test site in a rural area in the west of Germany. The coupled FOOT3DK/gas exchange model is validated for the diurnal cycle at singular grid points, and delivers realistic fluxes with respect to their order of magnitude and to the general daily course. Compared to the Jarvis-based big leaf scheme, simulations of latent heat fluxes with a photosynthesis-based scheme for stomatal conductance are more realistic. As expected, flux averages are strongly influenced by the underlying land cover. While the simulated net ecosystem exchange is highly correlated with leaf area index, this correlation is much weaker for the latent heat flux. Photosynthetic CO₂ uptake is associated with transpirational water loss via the stomata, and the resulting opposing surface fluxes of CO₂ and H₂O are reproduced with the model approach. Over vegetated surfaces it is shown that the coupling of a photosynthesis-based gas exchange model with the land-surface scheme of a mesoscale model results in more realistic simulated latent heat fluxes.     

Minor- and trace-element zoning in plagioclase: implications for magma chamber processes at Parinacota volcano, northern Chile

Textural and compositional zoning in plagioclase phenocrysts in a sample from Parinacota volcano (Chile) was investigated using backscattered electron images and electron microprobe analysis of major and trace elements. Large (2 mm) oscillatory zoned crystals (type I) with resorption surfaces of moderate An discontinuities (Ⲓ% An) and decreasing trace-element contents (Sr, Mg, Ti) towards the rim reflect melt differentiation and turbulent convection in the main magma body. Early recharge with a low-Sr mafic magma is seen in the core. Small-scale Sr variations in the core indicate limited diffusion and thus residence and differentiation times of the magma shorter than a few thousand years. Smaller crystals (type II) with low trace-element/An ratio reflect the influence of an H₂O-rich melt probably from a differentiated boundary layer. Closed-system in-situ crystallisation, mafic magma recharge and the role of a water-rich differentiated boundary layer can be distinguished from the An-trace element relationships. Crystals apparently move relatively freely between different parts and regimes in the magma chamber, evidence for "convective crystal dispersion". High-Sr type II crystals indicate an earlier input of Sr-rich mafic magma. Recharge of two distinct mafic magma types is thus identified (high-Sr and low-Sr), which must have been present - at increasing recharge rates with time - in the plumbing system throughout the volcano's history.     

The Signature of Sea Spray in the Hexos Turbulent Heat Flux Data

The role of sea spray intransferring heat and moisture across the air-sea interface has remained elusive. Some studies have reported that sea spray does not affect the turbulent air-sea heat fluxes for 10-m wind speeds up to at least 25 m s^-1, while others have reported important spray contributions for wind speeds as low as 12 m s^-1. One goal of the HEXOS (Humidity Exchange over the Sea) program was to quantify spray's contribution to the turbulent air-sea heat fluxes, but original analyses of the HEXOS flux data found the spray signal to be too small to be reliably identified amid the scatter in the data. We look at the HEXOS data again in the context of the TOGA-COARE bulk flux algorithm and a sophisticated microphysical spray model. This combination of quality data andstate-of-the-art modelling reveals a distinct spray signature in virtually all HEXOS turbulent heat flux data collected in winds of 15 m s^-1 and higher. Spray effects are most evident in the latent heat flux data, where spray contributes roughly 10% of the total turbulent flux in winds of 10 m s^-1 and between 10 and 40% in winds of 15–18 m s^-1. The spray contribution to the total sensible heat flux is also at least 10% in winds above 15 m s^-1. These results lead to a new, unified parameterization for the turbulent air-sea heat fluxes that should be especially useful in high winds because it acknowledges both the interfacial and spray routes by which the sea exchanges heat and moisture with the atmosphere.     

Internal sedimentary structure of linear dunes modelled with a cellular automaton

Linear dunes are the most common type of dune found on Earth and exist on several extra-terrestrial bodies, but despite this abundance their internal stratigraphy has not been commonly agreed. A cellular automaton is deployed to simulate the development of linear dunes, starting from a flat bed, under bi-modal oblique wind regimes of varying degrees of asymmetry. The internal stratigraphy of the linear dunes is monitored by keeping track of (buried) erosion surfaces, avalanche deposits and vertical accumulation, as well as the age of last subaerial exposure of the sediments. The simulations show the initial pattern-coarsening of a network of small dunes into fewer larger longitudinal ridges via bedform interactions and Y-junction dynamics. Three newly recognized types of bedform interaction are identified that relate to initial Y-junction dynamics: longitudinal crest-splitting, which creates free dune tips that can interact with adjacent dunes, and laterally oscillating interactions that lead to ephemeral Y-junctions (normal or reverse). The results show that these three bedform interactions leave no persistent signatures in the stratigraphic record. However, a further three bedform interactions involving the superposition of one dune onto another – merging, cannibalizing and repulsion (known from transverse dune field dynamics) – do leave specific evidence in the internal stratigraphy of the remaining dune, a buried interaction surface at a specific inclination. The preservation potential of this interaction surface varies between the three types. After the initial pattern-coarsening phase, the linear dunes become larger and more independent and their crest orientation follows the net resultant transport direction. The stratigraphies of mature dunes under wind regimes of differing asymmetry show that under (nearly) symmetrical winds the dune accumulates mainly vertically, with strata dipping parallel to the flanks, while under more asymmetrical wind regimes the internal stratigraphy resembles that of transverse dunes.     

Monazite as a monitor for melt‐rock interaction during cooling and exhumation

Granulite facies cordierite–garnet–biotite gneisses from the southeastern Reynolds Range, central Australia, contain both orthopyroxene‐bearing and orthopyroxene‐free quartzofeldspathic leucosomes. Mineral reaction microstructures at the interface of gneiss and leucosome observed in outcrop and petrographically, reflect melt‐rock interaction during crystallization. Accessory monazite, susceptible to fluid alteration, dissolution and recrystallization at high temperature, is tested for its applicability to constrain the chemical and P–T–time evolution of melt‐rock reactions during crystallization upon cooling. Bulk rock geochemistry and phase equilibria modelling constrain peak pressure and temperature conditions to 6.5–7.5 kbar and ~850°C, and U–Pb geochronology constrains the timing of monazite crystallization to 1.55 Ga, coeval with the Chewings Orogeny. Modelling predicts the presence of up to 15 vol.% melt at peak metamorphic conditions. Upon cooling below 800°C, melt extraction and in situ crystallization of melt decrease the melt volume to less than 7%, at which time it becomes entrapped and melt pockets induce replacement reactions in the adjacent host rock. Replacement reactions of garnet, orthopyroxene and K‐feldspar liberate Y, REE, Eu and U in addition to Mg, Fe, Al, Si and K. We demonstrate that distinguishing between monazite varieties solely on the basis of U–Pb ages cannot solve the chronological order of events in this study, nor does it tie monazite to the evolution of melt or stability of rock‐forming minerals. Rather, we argue that analyses of various internal monazite textures, their composition and overprinting relations allow us to identify the chronology of events following the metamorphic peak. We infer that retrograde reactions involving garnet, orthopyroxene and K‐feldspar can be attributed to melt‐rock interaction subsequent to partial melting, which is reflected in the development of compositionally distinct monazite textural domains. Internal monazite textures and their composition are consistent with dissolution and precipitation reactions induced by a high‐T melt. Monazite rims enriched in Y, HREE, Eu and U indicate an increased availability of these elements, consistent with the breakdown of orthopyroxene, garnet and K‐feldspar observed petrographically. Our study indicates that compositional and textural analysis of monazite in relation to major rock‐forming minerals can be used to infer the post‐peak chemical evolution of partial melts during high‐ to ultrahigh‐temperature metamorphism.     

Mobility and H2O loss from fluid inclusions in natural quartz crystals

Petrographic observations on quartz crystals from the Mole Granite (Australia) and other localities shed new light on the mechanisms of post-entrapment modification of fluid inclusions. These modifications include migration away from pseudosecondary trails, changes in fluid salinity and density, shape distortion and the formation of “sweat-haloes” around strongly deformed inclusions. Increases in fluid salinity, which usually are associated with inclusion migration, indicate water-losses of up to 50%. However, LA-ICP-MS-analysis of unmobilized and mobilized inclusions of the same trail reveals basically unchanged ratios of major – and trace element cations, with the exception of Li, which seems to be incorporated into the crystal lattice during migration. Despite the fact that all these modifications are closely related to deformation processes, they occur not only in mechanically deformed quartz, but also in free-standing crystals. In the latter samples, stress has been generated internally as a result of brazil-twinned growth and compositional zonation. These observations and their interpretation leads to a list of practical criteria that should help in differentiating between reliable and suspect fluid inclusions in other samples. Received: 17 November 1998 / Accepted: 16 April 1999     

Interpretation of tracer tests performed in fractured rock of the Lange Bramke basin, Germany

 Two multitracer tests performed in one of the major cross-fault zones of the Lange Bramke basin (Harz Mountains, Germany) confirm the dominant role of the fault zone in groundwater flow and solute transport. Tracers having different coefficients of molecular diffusion (deuterium, bromide, uranine, and eosine) yielded breakthrough curves that can only be explained by a model that couples the advective–dispersive transport in the fractures with the molecular diffusion exchange in the matrix. For the scale of the tests (maximum distance of 225 m), an approximation was used in which the influence of adjacent fractures is neglected. That model yielded nearly the same rock and transport parameters for each tracer, which means that the single-fracture approximation is acceptable and that matrix diffusion plays an important role. The hydraulic conductivity of the fault zone obtained from the tracer tests is about 1.5×10^–2 m/s, whereas the regional hydraulic conductivity of the fractured rock mass is about 3×10^–7 m/s, as estimated from the tritium age and the matrix porosity of about 2%. These values show that the hydraulic conductivity along the fault is several orders of magnitude larger than that of the remaining fractured part of the aquifer, which confirms the dominant role of the fault zones as collectors of water and conductors of fast flow. Received, April 1997 Revised, January 1998, August 1998 Accepted, August 1998     

A Density Correction for Lagrangian Particle Dispersion Models

Current Lagrangian particle dispersion models, used to simulate the dispersion of passive tracers in the turbulent planetary boundary layer (PBL), assume that the density is constant within the PBL. In deep PBLs, where the density at the boundary-layer top may be lower by more than 20% than at the surface, this assumption leads to errors in the tracer concentrations on the order of 10%. In the presence of a vertical wind shear, this also leads to inaccurate calculations of the horizontal tracer transport. To remove this deficiency, a Langevin equation is presented that contains a density correction term. The effect of the density correction is studied using data from a large-scale tracer experiment. It is found that for this experiment, the main effect of the density correction is an increase in the surface tracer concentrations, whereas the horizontal tracer transport patterns remain largely unaffected.     

The Transport Time pdf for a Simple Turbulent Flow

In a turbulent fluid, the time a particle needs to travel from a point source to the observation point, can be considered as a random variable. It is shown that the probability density function (pdf) for this random variable is determined by the Lagrangian particle position pdf. The characteristics of the transport time pdf are discussed for the simple case of a turbulent fluid moving with a constant mean velocity.