Assessment of a Mediterranean shallow lentic ecosystem (Lake Pamvotis, Greece) using benthic community diversity: Response to environmental parameters

Macroinvertebrates play a key role in freshwater lentic and lotic ecosystems. The macroinvertebrate benthic community of a shallow Mediterranean lake (Lake Pamvotis, NW Greece) was studied. The benthic assemblage was sampled monthly at five sites during a period of 1 year (Apr. 1998–Mar. 1999). In addition hypolimnetic water quality variables were monitored over the same period at each site.

The aim of the study was (a) to describe the intra-annual and spatial variability in benthic communities, (b) to relate possible community changes to environmental conditions and (c) to evaluate the responses of the lake's ecological status on community indices.

The benthic fauna of Lake Pamvotis was found to be very limited with a total of 10 species belonging to five taxonomic groups. The oligochaete community comprised 80% of the total benthic fauna with Potamothrix bavaricus as a new record for the Lake Pamvotis and Potamothrix hammoniensis, being the dominant benthic species represented more than 61% of the total benthic fauna. Chironomus plumosus was the most abundant chironomid species contributing with about 6% of the total benthic fauna, and Chaoborus flavicans with 19% was the important dipteran. Almost all benthic species showed the same intra-annual seasonal pattern, with peak population densities during spring and early summer except P. hammoniensis which predominated during the whole sampling period. Dissolved oxygen and temperature seemed to be the main environmental factors affecting community indices.

Benthic communities are affected by human disturbances in Lake Pamvotis shifting their composition to more tolerant taxa, reflecting also the eutrophic to hypertophic character of the lake.     

Applications of hybrid genetic algorithms in seismic tomography

Almost all earth sciences inverse problems are nonlinear and involve a large number of unknown parameters, making the application of analytical inversion methods quite restrictive. In practice, most analytical methods are local in nature and rely on a linearized form of the problem equations, adopting an iterative procedure which typically employs partial derivatives in order to optimize the starting (initial) model by minimizing a misfit (penalty) function. Unfortunately, especially for highly non-linear cases, the final model strongly depends on the initial model, hence it is prone to solution-entrapment in local minima of the misfit function, while the derivative calculation is often computationally inefficient and creates instabilities when numerical approximations are used. An alternative is to employ global techniques which do not rely on partial derivatives, are independent of the misfit form and are computationally robust. Such methods employ pseudo-randomly generated models (sampling an appropriately selected section of the model space) which are assessed in terms of their data-fit. A typical example is the class of methods known as genetic algorithms (GA), which achieves the aforementioned approximation through model representation and manipulations, and has attracted the attention of the earth sciences community during the last decade, with several applications already presented for several geophysical problems.In this paper, we examine the efficiency of the combination of the typical regularized least-squares and genetic methods for a typical seismic tomography problem. The proposed approach combines a local (LOM) and a global (GOM) optimization method, in an attempt to overcome the limitations of each individual approach, such as local minima and slow convergence, respectively. The potential of both optimization methods is tested and compared, both independently and jointly, using the several test models and synthetic refraction travel-time date sets that employ the same experimental geometry, wavelength and geometrical characteristics of the model anomalies. Moreover, real data from a crosswell tomographic project for the subsurface mapping of an ancient wall foundation are used for testing the efficiency of the proposed algorithm. The results show that the combined use of both methods can exploit the benefits of each approach, leading to improved final models and producing realistic velocity models, without significantly increasing the required computation time.     

Earthquake precursory signatures in electromagnetic radiation measurements in terms of day-to-day fractal spectral exponent variation: analysis of the eastern Aegean 13/04/2017–20/07/2017 seismic activity

The recent seismic activity in the eastern Aegean Sea, foregrounded by two major (M_L?=?6.1 and M_L?=?6.2) earthquakes occurring near Lesvos and Kos islands, respectively, is investigated in this work. Electromagnetic radiation measurements across different frequency bands and antenna orientations from a monitoring station in Agia Paraskevi, Lesvos Island, are analysed in order to reveal earthquake precursory signatures hidden in the electromagnetic data sequence. A straightforward, data-driven approach is employed in which day-to-day variations of the fractal characteristics of the measurements are adaptively monitored via a fractal spectral exponent similarity measure. The evolution of the fractal day-to-day variation in a 99-day period, which includes the two major earthquakes, shows a sustained, sudden increase lasting 1 to 4 days before every earthquake cluster of M_L?=?4.0 and above. Most importantly, this day-to-day variation subsides shortly after and remains relatively low during the absence of earthquakes, thus alleviating the emergence of false alarms.     

Dynamics of photosynthetic photon flux density (PPFD) and estimates in coastal northern California

Plants require solar radiation for photosynthesis and their growth is directly related to the amount received, assuming that other environmental parameters are not limiting. Therefore, precise estimation of photosynthetically active radiation (PAR) is necessary to enhance overall accuracies of plant growth models. This study aimed to explore the PAR radiant flux in the San Francisco Bay Area of northern California. During the growing season (March through August) for 2?years 2007?C2008, the on-site magnitudes of photosynthetic photon flux densities (PPFD) were investigated and then processed at both the hourly and daily time scales. Combined with global solar radiation (R _S) and simulated extraterrestrial solar radiation, five PAR-related values were developed, i.e., flux density-based PAR (PPFD), energy-based PAR (PARE), from-flux-to-energy conversion efficiency (fFEC), and the fraction of PAR energy in the global solar radiation (fE), and a new developed indicator??lost PARE percentages (LPR)??when solar radiation penetrates from the extraterrestrial system to the ground. These PAR-related values indicated significant diurnal variation, high values occurring at midday, with the low values occurring in the morning and afternoon hours. During the entire experimental season, the overall mean hourly value of fFEC was found to be 2.17???mol?J^?1, while the respective fE value was 0.49. The monthly averages of hourly fFEC and fE at the solar noon time ranged from 2.15 in March to 2.39???mol?J^?1 in August and from 0.47 in March to 0.52 in July, respectively. However, the monthly average daily values were relatively constant, and they exhibited a weak seasonal variation, ranging from 2.02?mol?MJ^?1 and 0.45 (March) to 2.19?mol?MJ^?1 and 0.48 (June). The mean daily values of fFEC and fE at the solar noon were 2.16?mol?MJ^?1 and 0.47 across the entire growing season, respectively. Both PPFD and the ever first reported LPR showed strong diurnal patterns. However, they had opposite trends. PPFD was high around noon, resulting in low values of LPR during the same time period. Both were found to be highly correlated with global solar radiation R _S, solar elevation angle h, and the clearness index K _t. Using the best subset selection of variables, two parametric models were developed for estimating PPFD and LPR, which can easily be applied in radiometric sites, by recording only global solar radiation measurements. These two models were found to be involved with the most commonly measured global solar radiation (R _S) and two large-scale geometric parameters, i.e., extraterrestrial solar radiation and solar elevation. The models were therefore insensitive to local weather conditions such as temperature. In particular, with two test data sets collected in USA and Greece, it was verified that the models could be extended across different geographical areas, where they performed well. Therefore, these two hourly based models can be used to provide precise PAR-related values, such as those required for developing precise vegetation growth models.     

Monodisperse and polydisperse colloid transport in water-saturated fractures with various orientations: Gravity effects

Numerical experiments are conducted to examine the effects of gravity on monodisperse and polydisperse colloid transport in water-saturated fractures with uniform aperture. Dense colloids travel in water-saturated fractures by advection and diffusion while subject to the influence of gravity. Colloids are assumed to neither attach onto the fracture walls nor penetrate the rock matrix based on the assumptions that they are inert and their size is larger than the pore size of the surrounding solid matrix. Both the size distribution of a colloid plume and colloid density are shown to be significant factors impacting their transport when gravitational forces are important. A constant-spatial-step particle-tracking code simulates colloid plumes with increasing densities transporting in water-saturated fractures while accounting for three forces acting on each particle: a deterministic advective force due to the Poiseuille flow field within the fracture, a random force caused by Brownian diffusion, and the gravitational force. Integer angles of fracture orientation with respect to the horizontal ranging from ±90° are considered: three lognormally distributed colloid plumes with mean particle size of 1 μm (averaged on a volumetric basis) and standard deviation of 0.6, 1.2 and 1.8 μm are examined. Colloid plumes are assigned densities of 1.25, 1.5, 1.75 and 2.0 g/cm³. The first four spatial moments and the first two temporal moments are estimated as functions of fracture orientation angle and colloid density. Several snapshots of colloid plumes in fractures of different orientations are presented. In all cases, larger particles tend to spread over wider sections of the fracture in the flow direction, but smaller particles can travel faster or slower than larger particles depending on fracture orientation angle.     

A cervelleite-like mineral and other Ag-Cu-Te-S minerals [Ag2CuTeS and (Ag,Cu)2TeS] in gold-bearing veins in metamorphic rocks of the Cycladic Blueschist Unit, Kallianou, Evia Island, Greece

A cervelleite-like mineral, two unnamed silver sulfotellurides in the system Ag-Cu-Te-S [Ag₂CuTeS, (Ag,Cu)₂TeS], Te-rich polybasite and cadmian tetrahedrite occur in gold-bearing quartz veins in metapelites and faults within brecciated marbles of the Cycladic Blueschist Unit in the Kallianou area (southern Evia Island, Greece). The quartz veins and faults are discordant to syn-metamorphic structures and formed during ductile to brittle deformation in the final stages of exhumation of the Styra Nappe extrusion wedge (~21?Ma). Te-rich polybasite (up to 7.4 wt. % Te), cadmian tetrahedrite (up to 12.4 wt. % Cd), together with electrum (23?C54 wt. % Ag) and the sulfotellurides, are the main silver carriers in the mineralization. The two unnamed sulfotellurides, Ag₂CuTeS and (Ag,Cu)₂TeS are believed to be new quaternary minerals in the system Ag-Cu-Te-S. These minerals and the cervelleite-like phase could have exsolved from galena during cooling (below 200°C). Initial temperatures for the formation of the sulfotellurides, in the form of hessite-intermediate solid solution, at Kallianou may be up to 300°C under logfS₂ values between?~ ?11.5 to ?8.3, and logfTe₂ from?~ ?14.8 to ?7.8. The values of logfTe₂ and logfS₂ during re-equilibration (at ~200°C) were constrained to ?19.5 to ?15.2 and to ?15.8 to ?11.5 respectively.     

Fast and reliable Markov chain Monte Carlo technique for cosmological parameter estimation

Markov chain Monte Carlo (MCMC) techniques are now widely used for cosmological parameter estimation. Chains are generated to sample the posterior probability distribution obtained following the Bayesian approach. An important issue is how to optimize the efficiency of such sampling and how to diagnose whether a finite-length chain has adequately sampled the underlying posterior probability distribution. We show how the power spectrum of a single such finite chain may be used as a convergence diagnostic by means of a fitting function, and discuss strategies for optimizing the distribution for the proposed steps. The methods developed are applied to current cosmic microwave background and large-scale structure data interpreted using both a pure adiabatic cosmological model and a mixed adiabatic/isocurvature cosmological model including possible correlations between modes. For the latter application, because of the increased dimensionality and the presence of degeneracies, the need for tuning MCMC methods for maximum efficiency becomes particularly acute.