Effects of crab–halophytic plant interactions on creek growth in a S.W. Atlantic salt marsh: A Cellular Automata model

The Bahía Blanca Estuary (38° 50′ S, and 62° 30′ W) presents salt marshes where interactions between the local flora (Sarcocornia perennis) and fauna (Chasmagnathus granulatus) generate some kind of salt pans that alter the normal water circulation and condition its flow and course towards tidal creeks. The crab–vegetation dynamics in the salt marsh presents variations that cannot be quantified in a reasonable period of time. The interaction between S. perennis plant and C. granulatus crab is based on simple laws, but its result is a complex biological mechanism that causes an erosive process on the salt marsh and favors the formation of tidal creeks. To study it, a Cellular Automata model is proposed, based on the laws deduced from the observation of these phenomena in the field, and then verified with measurable data within macroscale time units. Therefore, the objective of this article is to model how the interaction between C. granulatus and S. perennis modifies the landscape of the salt marsh and influences the path of tidal creeks. The model copies the basic laws that rule the problem based on purely biological factors.The Cellular Automata model proved capable of reproducing the effects of the interaction between plants and crabs in the salt marsh. A study of the water drainage of the basins showed that this interaction does indeed modify the development of tidal creeks. Model dynamics would likewise follow different laws, which would provide a different formula for the probability of patch dilation. The patch shape can be obtained changing the pattern that dilates.     

On the longitudinal distribution of the solar surges

After re-examination of data presented by Verma (1984), some misuse of solar parameteres published in Solar Geophysical Data is presented. Consequently some conclusions made by Verma (1984, 1986) on about the longitudinal distribution of cool solar surges based on these misuses are discussed.     

Mirdita Zone ophiolites and associated sediments in Albania reveal Neotethys Ocean origin

The Mirdita Ophiolite Zone in Albania is associated with widespread mélanges containing components of up to nappe-size. We dated matrix and components of the mélange by radiolarians, conodonts, and other taxa. The components consist of radiolarites, pelagic limestones and shallow-water limestones, all of Triassic age, as well as ophiolites. Triassic radiolarite as a primary cover of ophiolite material proves Middle Triassic onset of Mirdita ocean-floor formation. The mélange contains a turbiditic radiolarite-rich matrix (“radiolaritic flysch”), dated as Late Bajocian to Early Oxfordian. It formed as a synorogenic sediment during west-directed thrusting of ophiolite and sediment-cover nappes representing ocean floor and underplated fragments of the western continental margin. The tectonic structures formed during these orogenic events (“Younger Kimmeridian or Eohellenic Orogeny”) are sealed by Late Jurassic platform carbonates. The geological history conforms with that of the Inner Dinarides and adjoining areas; we therefore correlate the Mirdita-Pindos Ophiolite Zone with the Vardar Zone and explain its present position by far-distance west-directed thrusting.     

Simulations of 3D DC borehole resistivity measurements with a goal-oriented <Emphasis Type="Italic">hp</Emphasis> finite-element method. Part II: through-casing resistivity instruments

We simulate direct current (DC) borehole resistivity measurements acquired in steel-cased deviated wells for the assessment of rock formation properties. The assumed data acquisition configuration considers one current (emitter) and three voltage (collector) electrodes that are utilized to measure the second difference of the electric potential along the well trajectory. We assume a homogeneous, 1.27-cm-thick steel casing with resistivity equal to 10^− 5 Ω· m. Simulations are performed with two different numerical methodologies. The first one is based on transferring two-dimensional (2D) axisymmetric optimal grids to a three-dimensional (3D) simulation software. The second one automatically produces optimal 3D grids yielded by a 3D self-adaptive goal-oriented algorithm. Both methodologies utilize high-order finite elements (FE) that are specially well-suited for problems with high-contrast coefficients and rapid spatial variations of the electric field, as it occurs in simulations that involve steel-cased wells. The method based on transferring 2D-optimal grids is efficient in terms of CPU time (few seconds per logging position). Unfortunately, it may produce inaccurate 3D simulations in deviated wells, even though the error remains below 1% for the axisymmetric (vertical) well. The method based on optimal 3D grids, although less efficient in terms of CPU time (few hours per logging position), produces more accurate results that are validated by a built-in a posteriori error estimator. This paper provides the first existing simulations of through-casing resistivity measurements in deviated wells. Simulated resistivity measurements indicate that, for a 30° deviated well, measurements in conductive layers 0.01 Ω· m) are similar to those obtained in vertical wells. However, in resistive layers (10,000 Ω· m), we observe 100% larger readings in the 30° deviated well. This difference becomes 3,000% for the case of a 60° deviated well. For this highly-deviated well, readings corresponding to the conductive formation layer are about 30% smaller in magnitude than those in a vertical well. Shoulder effects significantly vary in deviated wells.     

Application of FLATModel, a 2D finite volume code, to debris flows in the northeastern part of the Iberian Peninsula

FLATModel is a two-dimensional shallow-water approximation code with corrections and modifications that create a simulation tool adapted to debris-flows behaviour. FLATModel uses the finite volume method with the numerical implementation of the Godunov scheme and includes correction terms regarding the effect of flow over high slopes and curvature. Additionally, the stop-and-go phenomenon, the basal entrainment and a correction regarding the front inclination of the final deposit are incorporated into FLATModel. In addition, different flow resistance laws were integrated in the numerical code including Bingham, Herschel–Bulkley and Voellmy fluid model. Firstly, our numerical model was validated using analytical solutions of a dam-break scenario and published data on a laboratory experiment. Secondly, three real events, which occurred in the northeastern part of the Iberian Peninsula, were back-calculated. Although field observations of the three events are not very detailed, the back-analyses revealed interesting patterns on the flow dynamics, and the numerical results generally showed good agreement with field data. Comparing the different flow resistance laws, the Voellmy fluid model presents the best behaviour regarding both the flow behaviour and the deposit characteristics. Preliminary simulation runs incorporating the effect of basal entrainment offered satisfactory results, although the final volume is rather sensitive on the selected friction angle of channel-bed material. The outcomes regarding the correction of the calculated front inclination of the final deposit showed that this implementation strongly improves the simulation results and better represents steep fronts of final deposits.     

Giant ancient landslide in the Alma water gap (Crimean Mountains,Ukraine): notes to the predisposition,structure, and chronology

Large-scale ancient landslides of the area of more than 5 km² and volume exceeding 200 × 10⁶ m³ are characteristic features of the valleys incised in the northern periphery of the Crimean Mountains (Ukraine). The largely affected area is located in the outermost cuesta range of the Crimean Mountains which consists of rigid Sarmatian limestones overlying weak Middle Miocene and Upper Palaeogene deposits. A giant landslide arose in the Alma water gap as a reflection of several coincident preparatory factors such as suitable bedrock stratification, smectite-rich bedrock exposed to swelling activity, presence of faults parallel to the valley trend, and river capture event which preceded the landslide event. The occurrence of such ancient megaslides is particularly interesting in the area which is characterized by low precipitation (<500 mm/year) and weak contemporary seismicity. It probably reflects a more dynamic environment in humid phases of the Holocene; however, seismic triggering along the Mesozoic suture zone cannot be rejected. Compressional features such as gravitational folds in the central and distal parts of the landslide, which probably correlate with the whole landslide genesis or its significant reactivation, arose, according to the radiocarbon dating, during the Holocene climatic optimum in the Atlantic period. The slope deformation has been relatively quiescent since that time, except minor historic reactivization which took place in the frontal part of the landslide. We suppose that the studied landslide could be classified as a transitional type of slope deformation with some signs of spreading and translational block slides.     

The qualitative zoning record of minerals. A method for determining the duration of metamorphic events?

Summary The efficacy of ion exchange in petrological systems, for example the Fe-Mg exchange between garnet and biotite, is a function of grainsize, temperature, rates of temperature change and diffusion parameters. The combination of these variables determines the final zoning profile of minerals. Therefore, zoning profiles may be used to derive one of these variables if the others are known. For example, a mineral grain that experienced a short thermal event may still preserve a zoning profile characteristic of the heating path. In contrast, if that grain experienced a long thermal event it may develop a zoning profile typical of the cooling path. Conversely, for the same temperature-time cycle, large grains may not completely equilibrate at the metamorphic peak, and in smaller grains any previous record may have been erased. This is commonly observed in natural rocks where different grains within one thin section preserve often qualitatively different zoning profiles. Thus, a critical grainsize,l _crit, may exist that separates grains with qualitatively different zoning profiles so that grains of sizel > l _crit still retain zoning information about the heating path and grains of sizel <l _crit contain only information about the cooling path. If the critical grainsize can be measured and an independent estimate for the peak metamorphic temperature exists, the duration of the thermal event may, in principle, be estimated.The applicability of this method to natural garnets is hampered by a range of uncertainties. However, a parameterisation of the critical grain size may be used to illustrate the relative importance of grain size, temperature and event duration to the equilibration of minerals. In this paper, the critical grainsize is parameterized for various temperature-time cycles and its dependence on diffusion parameters is discussed. It is shown that, for Barrovian conditions,l _crit between 0.1 and 1 mm separates garnets retaining prograde information from garnets retaining cooling path information. For illustration, we compare the results with the critical grainsize of garnets from two metamorphic terranes, the Prydz Bay region (Antarctica) and the Koralm complex (Eastern Alps). Despite the large range of uncertainties attached to the method, it is shown that the critical grainsize of garnets in both terranes is consistent with a very short duration of the last thermal event that affected the two regions.

---

Die qualitative Form von Zonierungsprofilen in Mineralen. Information über die Dauer metamorpher Ereignisse?

Zusammenfassung Zonierungsprofile von Paragenesen im Ionenaustausch, zum Beispiel das Fe-Mg Austauschsgleichgewicht zwischen Granat und Biotit, sind eine Funktion von (i) Korngröße, (ii) Temperatur, (iii) Kühl- und Heizgeschwindigkeit sowie (iv) den Diffusionskonstanen. Es sollte daher möglich sein, Zonierungsprofile dazu zu verwenden, einen dieser Parameter zu bestimmen, wenn die anderen bekannt sind. Zonierungsprofile von Kristallen die ein kurzes thermisches Ereignis erfahren haben, mögen daher prograde Information beinhalten, wogegen derselbe Kristall ein retrogrades Profil aufweisen mag, wenn er einem langlebigen thermischen Ereignis unterlag (prograd und retrograd ist hier als: bei ansteigender Temperatur und bei abfallender Temperatur definiert). In Kristallen die dem gleichen Temperatur-Zeitpfad unterworfen waren, könnten kleine Körner nur die Kühlgeschichte dokumentieren, wogegen große noch Zonierungsinformation vom Heizpfad aufweisen. In natürlichen Gesteinen wird das oft dadurch beobachtet, daß verschiedene Körner auseinem Dünnschliff qualitativ verschiedene Zonierungsprofile aufweisen. Es ist daher möglich eine kritische Korngöße,l _crit, zu definieren, die Korngrößen mit qualitativ verschiedenen Zonierungsprofilen voneinander trennt. Körner mit einem Durchmesserl > l _crit haben, zumindest teilweise, noch prograde Profile, wogegen Körner mitl <l _crit nur retrograde Information dokumentieren. Wenn man diese kritische Korngröße messen kann, sollte sie dazu benutzt werden können, etwas über die Dauer des thermischen Ereignisses auszusagen.Die Anwendbarkeit dieser Methode ist durch eine Reihe von Fehlern limitiert. Nichtsdestotrotz, ist eine Parameterisierung der kritischen Korngröße nützlich um die relative Wichtigkeit von Korngröße, Temperatur und Zeit, für die Entwicklung von Zonierungsprofilen, zu illustrieren. Unser Modell kann dazu benutzt werden, um die Größenordnung von Granaten abzuschätzen, die noch prograde Information dokumentieren können. Die Abhängigkeit der kritischen Korngröße von Diffusionsparametern und der Form des Temperatur-Zeit Pfades wird ebenfalls diskutiert. Es wird gezeigt, daß für Metamorphose in der mittleren Grünschiefer und Amphibolit fazies,l _crit zwischen 0.1 and 1 mm prograd zonierte von retrograd zonierten Granaten trennt. Um das Modell zu illustrieren, werden die Ergebnisse auf die Prydz Bay Region (Antarktis) und den Koralm Komplex (Ostalpen) angewendet. Trotz der großen Fehlergrenzen der Methode glauben wir zeigen zu können, daß die kritische Korngröße beider Terrains andeutet, daß das letzte thermische Ereignis in beiden Gebieten kurzlebig war.

---

---

With 5 Figures     

Bayesian estimation of seismic hazard for two sites in Switzerland

Seismic hazard analysis is based on data and models, which both are imprecise and uncertain. Especially the interpretation of historical information into earthquake parameters, e.g. earthquake size and location, yields ambiguous and imprecise data. Models based on probability distributions have been developed in order to quantify and represent these uncertainties. Nevertheless, the majority of the procedures applied in seismic hazard assessment do not take into account these uncertainties, nor do they show the variance of the results. Therefore, a procedure based on Bayesian statistics was developed to estimate return periods for different ground motion intensities (MSK scale).Bayesian techniques provide a mathematical model to estimate the distribution of random variables in presence of uncertainties. The developed method estimates the probability distribution of the number of occurrences in a Poisson process described by the parameter . The input data are the historical occurrences of intensities for a particular site, represented by a discrete probability distribution for each earthquake. The calculation of these historical occurrences requires a careful preparation of all input parameters, i.e. a modelling of their uncertainties. The obtained results show that the variance of the recurrence rate is smaller in regions with higher seismic activity than in less active regions. It can also be demonstrated that long return periods cannot be estimated with confidence, because the time period of observation is too short. This indicates that the long return periods obtained by seismic source methods only reflects the delineated seismic sources and the chosen earthquake size distribution law.     

The contact binary system YY Eri

Relatively new and unanalysed photometric data-sets of the contact binary system YY Eri are presented. The light curves have been analysed using information limit optimization techniques, and the Binary Maker program of Bradstreet (1992).Comparison of the results allows insight into determinacy questions affecting the W UMa type of light curve, which relates to our general understanding of contact binaries. A limiting contact configuration cannot be ruled out on the basis of empirical data-analysis alone.Very careful analysis of the photospheric flux distribution over the surface of the Sun may help establish appropriate values of the gravity brightening parameter for cool dwarf stars. At present, however, empirical photometric information on contact binaries is non-discriminatory: i.e. alternative cosmogonies can find alternative support from available evidence. The implication is then for more and better observational data to allow better real independent parameter determination.