LES Study of the Energy Imbalance Problem with Eddy Covariance Fluxes

The spatial representativeness of heat fluxes on the basis of single-tower measurements, and the mechanism of the so-called energy imbalance problem, are investigated through numerical experiments using large-eddy simulation (LES). LES experiments are done for the daytime atmospheric boundary layer heated over a flat surface, as a best-case scenario completely free of sensor errors and the uncertainties of field conditions. Imbalance is defined as the deviation of the `turbulent' heat flux at a grid point from the horizontally averaged `total' heat flux. Both the theoretical and numerical results of the present study suggest the limitation of single-tower measurements and the necessity of horizontally-distributed observation networks.The temporally averaged `turbulent' flux based on a point measurement systematically underestimates the `total' flux (negative imbalance). This is attributed to local advection effects caused by the existence of turbulent organized structures (TOS), whose time scale is much longer than that of thermal plumes. The temporal and spatial change of TOS patterns causes low-frequency trends in the velocity and temperature data resulting in large scatter of the flux estimates. The influences of geostrophic wind speed, averaging time, observation height, computational domain size and resolution on tower-measured fluxes are also discussed. Finally, it is suggested that a weak inhomogenity in surface heating may reduce the negative bias of flux estimates.     

Use of Plant Bioassays for the Detection of Genotoxins in the Aquatic Environment

The present article describes the use of plant bioassays for the detection of genotoxins in the aquatic environment and gives an overview of test methods, their detection spectrum for environmental mutagens and carcinogens and their limitations and pitfalls. The most widely used test systems are micronucleus assays with meiotic cells of Tradescantia and with meristematic root tip cells of Allium and Vicia. In the last years, protocols for single cell gel electrophoresis assays have been developed, which can be conducted with a variety of species. Also various gene mutation test procedures have been developed with plants but they have hardly ever been used in environmental studies. Plants detect a broad variety of environmentally relevant genotoxins, in particular directly acting compounds. Many pesticides and industrial chemicals caused positive results; plant assays are also a unique tool for the detection of DNA‐reactive carcinogenic heavy metals in the environment. In many studies with complex mixtures, positive results were obtained which indicates that plants are sufficiently sensitive to detect effects without concentration of water samples. One of the shortcomings of the use of plants as indicators is their lack of sensitivity towards certain classes of promutagens such as nitrosamines, heterocyclic amines, and polycyclic aromatic hydrocarbons (PAHs). However, the former compounds are hardly encountered in the environment and PAHs can easily be detected chemically and in other mutagenicity tests. Taken together, the currently available data show that plant bioassays are a useful component of test batteries for environmental monitoring.     

Population dynamics of spot,leiostomus xanthurus,in polyhaline tidal creeks of the York River estuary,Virginia

Most populations of estuarine-dependent, early life stages of marine fishes are open. As a result, it has been difficult to apply conventional population models to most systems. In this study, a marked population of young-of-year spot (Leiostomus xanthurus) was released into a polyhaline tidal creek within the Guinea Marshes of the York River estuary, Virginia. Over a 90-day study period, 221 marked fishes were recaptured. Plots of the ratio of marked to unmarked individuals (m_i/n_i) in subsequent samples indicated that the population was resident in the creek for up to 162 days with the average individual present for 81 days. When this population turnover rate was compared to the total population decay rate (marked plus unmarked fish), it was determined that exchange between habitats (immigration/emigration) accounted for about 36.4% of the total decay rate, with the remainder attributed to natural mortality. By correcting the overall disappearance rate for population turnover due to emigration and using this adjusted value as a measure of instantaneous mortality (z), the estimated production (over 90 days) in this population was 23,630 cal (98,870 J) per m². This figure agrees with a previously derived estimate for spot in the Guinea marshes and is nearly two orders of magnitude higher than other reported values for this species for all size classes over the entire growing season.     

Spatial representativeness of single tower measurements and the imbalance problem with eddy-covariance fluxes: results of a large-eddy simulation study

A large-eddy simulation (LES) study is presented that investigates the spatial variability of temporal eddy covariance fluxes and the systematic underestimation of representative fluxes linked to them. It extends a prior numerical study by performing high resolution simulations that allow for virtual measurements down to 20 m in a convective boundary layer, so that conditions for small tower measurement sites can be analysed. It accounts for different convective regimes as the wind speed and the near-surface heat flux are varied. Moreover, it is the first LES imbalance study that extends to the stable boundary layer. It reveals shortcomings of single site measurements and the necessity of using horizontally-distributed observation networks. The imbalances in the convective case are attributed to a locally non-vanishing mean vertical advection due to turbulent organised structures (TOS). The strength of the TOS and thus the imbalance magnitude depends on height, the horizontal mean wind and the convection type. Contrary to the results of a prior study, TOS cannot generally be responsible for large energy imbalances: at low observation heights (corresponding to small towers and near-surface energy balance stations) the TOS related imbalances are generally about one order of magnitude smaller than those in field experiments. However, TOS may cause large imbalances at large towers not only in the case of cellular convection and low wind speeds, as found in the previous study, but also in the case of roll convection at large wind speeds. In the stably stratified boundary layer for all observation heights neither TOS nor significant imbalances are observed. Attempting to reduce imbalances in convective situations by applying the conventional linear detrending method increases the systematic flux underestimation. Thus, a new filter method is proposed.     

Non-stationarity of the signal and noise characteristics of seasonal precipitation anomalies

In order to improve seasonal-to-interannual precipitation forecasts and their application by decision makers, there is a clear need to understand when, where, and to what extent seasonal precipitation anomalies are driven by potentially predictable surface–atmosphere interactions versus to chaotic interannual atmospheric dynamics. Using a simple Monte Carlo approach, interannual variability and linear trends in the SST-forced signal and potential predictability of boreal winter precipitation anomalies is examined in an ensemble of twentieth century AGCM simulations. Signal and potential predictability are shown to be non-stationary over more than 80% of the globe, while chaotic noise is shown to be stationary over most of the globe. Correlation analysis with respect to magnitudes of the four leading modes of global SST variability suggests that interannual variability and trends in signal and potential predictability over 35% of the globe is associated with ENSO-related SST variability; signal and potential predictability are not significantly associated with SST modes characterized by a global SST trend, North Atlantic SST variability, and North Pacific SST variability, respectively. Results suggest that mechanisms other than SST variability contribute to the non-stationarity of signal and noise characteristics of hydroclimatic variability over mid- and high-latitude regions.     

A New Aerodynamic Parametrization for Real Urban Surfaces

This study conducted large-eddy simulations (LES) of fully developed turbulent flow within and above explicitly resolved buildings in Tokyo and Nagoya, Japan. The more than 100 LES results, each covering a 1,000 $\times $ 1,000 m $^{2}$ area with 2-m resolution, provide a database of the horizontally-averaged turbulent statistics and surface drag corresponding to various urban morphologies. The vertical profiles of horizontally-averaged wind velocity mostly follow a logarithmic law even for districts with high-rise buildings, allowing estimates of aerodynamic parameters such as displacement height and roughness length using the von Karman constant $=$ 0.4. As an alternative derivation of the aerodynamic parameters, a regression of roughness length and variable Karman constant was also attempted, using a displacement height physically determined as the central height of drag action. Although both the regression methods worked, the former gives larger (smaller) values of displacement height (roughness length) by 20–25 % than the latter. The LES database clearly illustrates the essential difference in bulk flow properties between real urban surfaces and simplified arrays. The vertical profiles of horizontally-averaged momentum flux were influenced by the maximum building height and the standard deviation of building height, as well as conventional geometric parameters such as the average building height, frontal area index, and plane area index. On the basis of these investigations, a new aerodynamic parametrization of roughness length and displacement height in terms of the five geometric parameters described above was empirically proposed. The new parametrizations work well for both real urban morphologies and simplified model geometries.     

Investigation of droplet dynamics in a convective cloud using a Lagrangian cloud model

A precipitating convective cloud is simulated successfully using the Lagrangian cloud model, in which the flow field is simulated by large eddy simulation and the droplets are treated as Lagrangian particles, and the results are analyzed to investigate precipitation initiation and to examine the parameterization of cloud microphysics. It is found that raindrops appear initially near the cloud top, in which strong turbulence and broadened droplet spectrum are induced by the entrainment of dry air, but high liquid–water mixing ratio is maintained within cloud parts because of insufficient mixing. Statistical analysis of the downward vertical velocity of a droplet W reveals that the transition from cloud droplets to raindrops occurs in the range 20 μm < r < 100 μm, while the variation of W depends on turbulence as well as the droplet radius r. The general pattern of the raindrop size distribution is found to be consistent with the Marshall–Palmer distribution. The precipitation flux can be underestimated substantially, if the terminal velocity $w_{\text{s}}$ is used instead of W, but it is not sensitive to the choice of the critical droplet radius dividing cloud drops and raindrops. It is also found that precipitation starts earlier and becomes stronger if the effect of turbulence is included in the collection kernel.     

Geochronology of the Birim Supergroup of the West African craton in the Wa-Bolé region of west-central Ghana: Implications for the stratigraphic framework

The Birim rocks of the West African craton comprise belts of greenschist- to amphibolite-grade gneiss and schist, and subparallel basins of greenschist-grade phyllite of volcaniclastic and epiclastic origin, which were intruded by igneous rocks. The granitoids intruded between 2213 and 2060 Ma and overlap with the volcaniclastic units dated between 2211 and 2064 Ma. The simultaneous occurrence of the magmatic events and irregular distribution of the rock ages hamper the formulation of a stratigraphic succession. SHRIMP spot analyses were done on older cores, crystals and rims from 23 rocks from the Bolé-Wa region in west-central Ghana. The crystallization ages range from 2195 to 2118 Ma, the inherited ages from 2876 to 2130 Ma, and metamorphic ages from 2114 to 2090 Ma. Aided by metamorphic, structural and chemical studies an older geotectonic cycle (2195–2150 Ma), containing the Dole and Guropie Suite and Bolé Group, was established. These units were subjected to several orthogonal and shear deformation events. These events were followed by the contemporaneous Sawla calc-alkaline monzonitic plutonism (2132–2126 Ma) and deposition of the epiclastic Maluwe Group (2137–2125 Ma) of calc-alkaline felsic to tholeiitic volcanic origin. Deformation of the basin beds was succeeded by the intrusion of the Tanina Suite granitoids of 2122–2120 Ma, which, themselves, were deformed prior to 2119 Ma. At 2118 Ma syenite and gabbro intruded along conjugate extension fractures. The gabbro and syenite of the Wakawaka Suite were only affected by three events of brittle strike-slip faulting. The first had significant displacement along NNE- to NE-directed shear zones, while the latter only formed conjugate joint systems with limited transport. Palaeo- to Neoarchaean cores, the oldest yet reported in the Baoulé Mossi domain, are restricted to the gneissic Dole Suite biotite granites. The presence of Dole-, Guropie-, Sawla-, and Tanina-aged older cores and grains in younger rocks reflects continuous reworking of the developing crust during successive magmatic episodes. Zircon rim growth between 2105 and 2090 Ma indicates posttectonic crustal thickening. The low Rb/Sr R_o of ～0.7032 of gabbro and monzonite, and the recycling of the Birim-age crust confirm the primary and juvenile nature of the West African craton after ～2195 Ma. With the various ages obtained, it was possible to link deposition, magmatism and deformation to crustal processes, and establish the cyclic geotectonic evolution in the West Africa craton (individual basin opening and closure) over time as part of an intraoceanic arc–back-arc basin system.     

Urban and rural limestone weathering; the contribution of dust to black crust formation

A novel approach of studying dust, black crust and host limestones by using various techniques is presented in this paper. Samples were collected from limestone monuments located in rural, urban and industrial areas in Germany and Hungary. The paper focuses on differences in the mineralogical composition, major and trace element distribution of materials and their total polyaromatic hydrocarbons (PAH) concentration having different exposure to air pollutants. High concentration of gypsum was found in laminar and dendritic black crusts and even in the dust of both urban and rural areas, despite the low concentration of SO₂ in rural atmospheric environments. The black crusts show a typical microfabric with distinct layers of various gypsum–calcite mineral associations, silt-sized mineral fragments and black particle content. In dust, newly formed gypsum crystals were found along with mineral fragments, and siliceous and, less frequently, carbonaceous fly ash particles. High concentration of lead was found both in the samples of urban and rural areas of Germany and Hungary. Lead mostly accumulates in dust and also in the black crust. Highest lead concentrations were found not on the surface of the black crust, but close to the crust/limestone boundary according to LA-ICPMS analyses. It indicates that despite the ban of leaded petrol, lead is still present in the dust and the in the soiled zones of the built environment. Similarly to lead, no significant difference in the total PAH content of dust and crust samples were found in the two countries. Dusts are equally enriched in PAH in both countries, whereas black crusts contain fewer amounts. Most of the aqueous extracts of the samples were saturated with respect to sulphate, which also indicate the presence of gypsum, whilst in some German samples unsaturated conditions were detected. In summary, by using a combination of these analytical techniques, it was possible to distinguish rural and urban samples and to outline the role of dust in bringing pollutants to the stone surface. It is suggested that the analysed dust and stone samples indicate not only present, but also the past, pollution fluxes and as a consequence stones similar to sediments or surface waters can be used as environment indicators.     

Weathering of Fruchtschiefer building stones: mineral dissolution or rock disaggregation?

Building stones manufactured from contact metamorphic slates (Fruchtschiefer slate) from Theuma (Sachsen, Germany) were investigated for mineralogical alterations as well as for changes in porosity and surface roughness due to weathering. After weathering periods of several years to decades, the originally dark gray-colored slates show pale spots of several centimeters in size at the surface of building stones. The dark-colored and light-colored sections of the slate show no differences in mineralogy. Surface weathering did not result in newly precipitated minerals. It was also found that the observed differences in color are not caused by variations in sedimentary organic carbon concentration or in sulfide/sulfate concentrations. Obtained results instead indicate that dark surface sections may show a thin cover of recent organic matter (OM), e.g., living OM, soot, dirt, etc. Microscopic investigations suggested that this cover was exfoliated at light-colored surface sections. The observed disaggregation of the upper 2 mm of the building block material results in an increase in porosity. Porosity of black (unweathered) slate is <2 vol.%. Due to weathering, the slate’s pores with diameters >1 μm show a significant increase in frequency compared to the original pore size distribution. Porosity of weathered rock volumes increased to approx. 8 vol.%. Discolored surface sections show a higher surface roughness (root-mean-square roughness, Rq ~ 1 μm) compared to dark-colored slate surfaces (Rq ~ 200 nm), both data are for cleavage planes. Preferentially, the discolored surface sections are located close to the edges of cut stones. This and the alteration in porosity, pore size, and surface roughness indicate that color changes of the slate are largely influenced by rock disaggregation proceeding from the edges into the center rather than by mineral dissolution/precipitation processes.