Identification of tsunami-induced deposits using numerical modeling and rock magnetism techniques: A study case of the 1755 Lisbon tsunami in Algarve, Portugal

Storm- and tsunami-deposits are generated by similar depositional mechanisms making their discrimination hard to establish using classic sedimentologic methods. Here we propose an original approach to identify tsunami-induced deposits by combining numerical simulation and rock magnetism. To test our method, we investigate the tsunami deposit of the Boca do Rio estuary generated by the 1755 earthquake in Lisbon which is well described in the literature. We first test the 1755 tsunami scenario using a numerical inundation model to provide physical parameters for the tsunami wave. Then we use concentration (MS. SIRM) and grain size (chi(ARM), ARM, B1/2, ARM/SIRM) sensitive magnetic proxies coupled with SEM microscopy to unravel the magnetic mineralogy of the tsunami-induced deposit and its associated depositional mechanisms. In order to study the connection between the tsunami deposit and the different sedimentologic units present in the estuary, magnetic data were processed by multivariate statistical analyses. Our numerical simulation show a large inundation of the estuary with flow depths varying from 0.5 to 6 m and run up of similar to 7 m. Magnetic data show a dominance of paramagnetic minerals (quartz) mixed with lesser amount of ferromagnetic minerals, namely titanomagnetite and titanohematite both of a detrital origin and reworked from the underlying units. Multivariate statistical analyses indicate a better connection between the tsunami-induced deposit and a mixture of Units C and D. All these results point to a scenario where the energy released by the tsunami wave was strong enough to overtop and erode important amount of sand from the littoral dune and mixed it with reworked materials from underlying layers at least 1 m in depth. The method tested here represents an original and promising tool to identify tsunami-induced deposits in similar embayed beach environments.

Identification of tsunami-induced deposits using numerical modeling and rock magnetism techniques: A study case of the 1755 Lisbon tsunami in Algarve, Portugal

Storm- and tsunami-deposits are generated by similar depositional mechanisms making their discrimination hard to establish using classic sedimentologic methods. Here we propose an original approach to identify tsunami-induced deposits by combining numerical simulation and rock magnetism. To test our method, we investigate the tsunami deposit of the Boca do Rio estuary generated by the 1755 earthquake in Lisbon which is well described in the literature. We first test the 1755 tsunami scenario using a numerical inundation model to provide physical parameters for the tsunami wave. Then we use concentration (MS. SIRM) and grain size (chi(ARM), ARM, B1/2, ARM/SIRM) sensitive magnetic proxies coupled with SEM microscopy to unravel the magnetic mineralogy of the tsunami-induced deposit and its associated depositional mechanisms. In order to study the connection between the tsunami deposit and the different sedimentologic units present in the estuary, magnetic data were processed by multivariate statistical analyses. Our numerical simulation show a large inundation of the estuary with flow depths varying from 0.5 to 6 m and run up of similar to 7 m. Magnetic data show a dominance of paramagnetic minerals (quartz) mixed with lesser amount of ferromagnetic minerals, namely titanomagnetite and titanohematite both of a detrital origin and reworked from the underlying units. Multivariate statistical analyses indicate a better connection between the tsunami-induced deposit and a mixture of Units C and D. All these results point to a scenario where the energy released by the tsunami wave was strong enough to overtop and erode important amount of sand from the littoral dune and mixed it with reworked materials from underlying layers at least 1 m in depth. The method tested here represents an original and promising tool to identify tsunami-induced deposits in similar embayed beach environments.

Identification of tsunami-induced deposits using numerical modeling and rock magnetism techniques: Astudy case of the 1755 Lisbon tsunami in Algarve, Portugal (vol 182, pg 187, 2010)

Storm- and tsunami-deposits are generated by similar depositional mechanisms making their discrimination hard to establish using classic sedimentologic methods. Here we propose an original approach to identify tsunami-induced deposits by combining numerical simulation and rock magnetism. To test our method, we investigate the tsunami deposit of the Boca do Rio estuary generated by the 1755 earthquake in Lisbon which is well described in the literature. We first test the 1755 tsunami scenario using a numerical inundation model to provide physical parameters for the tsunami wave. Then we use concentration (MS. SIRM) and grain size (chi(ARM), ARM, B1/2, ARM/SIRM) sensitive magnetic proxies coupled with SEM microscopy to unravel the magnetic mineralogy of the tsunami-induced deposit and its associated depositional mechanisms. In order to study the connection between the tsunami deposit and the different sedimentologic units present in the estuary, magnetic data were processed by multivariate statistical analyses. Our numerical simulation show a large inundation of the estuary with flow depths varying from 0.5 to 6 m and run up of similar to 7 m. Magnetic data show a dominance of paramagnetic minerals (quartz) mixed with lesser amount of ferromagnetic minerals, namely titanomagnetite and titanohematite both of a detrital origin and reworked from the underlying units. Multivariate statistical analyses indicate a better connection between the tsunami-induced deposit and a mixture of Units C and D. All these results point to a scenario where the energy released by the tsunami wave was strong enough to overtop and erode important amount of sand from the littoral dune and mixed it with reworked materials from underlying layers at least 1 m in depth. The method tested here represents an original and promising tool to identify tsunami-induced deposits in similar embayed beach environments.

Metallogenic model of the Trepca Pb-Zn-Ag skarn deposit, Kosovo: Evidence from fluid inclusions, Rare Earth Elements, and stable isotope data

The Trepca Pb-Zn-Ag skarn deposit (29 Mt of ore at 3.45% Pb, 2.30% Zn, and 80 g/t Ag) is located in the Kopaonik block of the western Vardar zone, Kosovo. The mineralization, hosted by recrystallized limestone of Upper Triassic age, was structurally and lithologically controlled. Ore deposition is spatially and temporally related with the postcollisional magmatism of Oligocene age (23-26 Ma). The deposit was formed during two distinct mineralization stages: an early prograde closed-system and a later retrograde open-system stage. The prograde mineralization consisting mainly of pyroxenes (Hd(54-100)Jo(0-45)Di(0-45)) resulted from the interaction of magmatic fluids associated with Oligocene (23-26 Ma) postcollisional magmatism. Whereas there is no direct contact between magmatic rocks and the mineralization, the deposit is classified as a distal Pb-Zn-Ag skarn. Abundant pyroxene reflects low oxygen fugacity (<10(-31) bar) and anhydrous environment. Fluid inclusion data and mineral assemblage limit the prograde stage within a temperature range between 390 degrees and 475 degrees C. Formation pressure is estimated below 900 bars. Isotopic composition of aqueous fluid, inclusions hosted by hedenbergite (delta D = -108 to -130 parts per thousand; delta O-18 = 7.5-8.0 parts per thousand), Mn-enriched mineralogy and high REE content of the host carbonates at the contact with the skarn mineralization suggest that a magmatic fluid was modified during its infiltration through the country rocks. The retrograde mineral assemblage comprises ilvaite, magnetite, arsenopyrite, pyrrhotite, marcasite, pyrite, quartz, and various carbonates. Increases in oxygen and sulfur fugacities, as well as a hydrous character of mineralization, require an open-system model. The opening of the system is related to phreatomagmatic explosion and formation of the breccia. Arsenopyrite geothermometer limits the retrograde stage within the temperature range between 350 degrees and 380 degrees C and sulfur fugacity between 10(-8.8) and 10(-7.2) bars. The principal ore minerals, galena, sphalerite, pyrite, and minor chalcopyrite, were deposited from a moderately saline Ca-Na chloride fluid at around 350 degrees C. According to the isotopic composition of fluid inclusions hosted by sphalerite (delta D = -55 to -74 parts per thousand; delta O-18 = -9.6 to -13.6 parts per thousand), the fluid responsible for ore deposition was dominantly meteoric in origin. The delta S-31 values of the sulfides spanning between -5.5 and +10 parts per thousand point to a magmatic origin of sulfur. Ore deposition appears to have been largely contemporaneous with the retrograde stage of the skarn development. Postore stage accompanied the precipitation of significant amount of carbonates including the travertine deposits at the deposit surface. Mineralogical composition of travertine varies from calcite to siderite and all carbonates contain significant amounts of Mn. Decreased formation temperature and depletion in the REE content point to an influence of pH-neutralized cold ground water and dying magmatic system.

Volcanic stratigraphy of Hannah Point, Livingston Island, South Shetland Islands, Antarctica

The Upper Cretaceous volcanic succession of Hannah Point is the best exposure of the Antarctic Peninsula Volcanic Group on L ivingston Island. The aim of the present paper is to contribute to the characterisation of the stratigr a p hy and petrogr a p hy of this little studied succession, and briefly discuss some aspects of the eru p t ive style of its volcanism. The succession is about 470 m thick and is here subdivided into five lithostratigraphic units (A to E from base to top). Unit A, approximately 120 m thick, is mainly composed of polymict clast-supported volcaniclastic breccias and also includes a dacitic lava laye r. Interstratified in the breccias of this unit, there is a thin laminated devitrified layer which shows some degree of welding. Unit B, approx imately 70 m thick, is almost entirely composed of volcaniclastic breccias, and includes a volcaniclastic conglomerate laye r. Breccias in this unit can be subdivided into two distinct types; polymict clast-supported breccias, and monomict matrix-supported breccias rich in juvenile components and displaying incipient welding. Unit C, about 65 m thick, is mainly composed of basaltic lavas, which are interlayered with minor vo lcaniclastic breccias. Unit D, approximately 65 m thick, is lithologically similar to unit B, composed of an alternation of polymict clasts upported breccias and matrix-supported breccias, and includes a volcaniclastic conglomerate laye r. Unit E, about 150 m thick, is mainly formed of thick andesitic lava layers. Minor basaltic dykes and a few normal faults cut the succession, and the contact betwe e n units A and B can be interpreted both as an unconformity or a fault. The matrix-supported breccias included in the succession of Hannah Point have high contents of juvenile components and incipient welding, which suggest that part of the succession is the result of pyroclastic fragmentation and emplacement from pyroclastic flows. In contrast, the polymict clast-supported breccias suggest reworking of previous deposits and deposition from cool mass flows. The lavas indicate eff u s ive volcanic eruptions, and the absence of features indicative of subaqueous volcanism suggests that at least these portions of the succession were emplaced in a subaerial environment .

Macular dystrophy associated with the mitochondrial DNA A3243G mutation: pericentral pigment deposits or atrophy? Report of two cases and review of the literature.

BACKGROUND: The A3243G point mutation in mitochondrial DNA (mtDNA) is associated with MELAS (mitochondrial encephalomyopathy with lactic acidosis and stroke-like episodes) and MIDD syndromes (maternally inherited diabetes and deafness). Both MELAS and MIDD patients can present with visual symptoms due to a retinopathy, sometimes before the genetic diagnosis is made. CASE PRESENTATION: Patient 1: 46 year-old woman with diabetes mellitus and hearing loss was referred for an unspecified maculopathy detected during screening evaluation for diabetic retinopathy. Visual acuity was 20/20 in both eyes. Fundus examination showed bilateral macular and peripapillary hyperpigmented/depigmented areas.Patient 2: 45 year-old woman was referred for recent vision loss in her left eye. History was remarkable for chronic fatigue, migraine and diffuse muscular pain. Visual acuity was 20/20 in her right eye and 20/30 in her left eye. Fundus exhibited several nummular perifoveal islands of retinal pigment epithelium atrophy and adjacent pale deposits in both eyes.Retinal anatomy was investigated with autofluorescence, retinal angiography and optical coherence tomography. Retinal function was assessed with automated static perimetry, full-field and multifocal electroretinography and electro-oculography. Genetic testing of mtDNA identified a point mutation at the locus 3243. CONCLUSION: Observation of RPE abnormalities in the context of suggestive systemic findings should prompt mtDNA testing.

Deposits of the sandy braided South Saskatchewan River: Implications for the use of modern analogs in reconstructing channel dimensions in reservoir characterization

Estimation of the dimensions of fluvial geobodies from core data is a notoriously difficult problem in reservoir modeling. To try and improve such estimates and, hence, reduce uncertainty in geomodels, data on dunes, unit bars, cross-bar channels, and compound bars and their associated deposits are presented herein from the sand-bed braided South Saskatchewan River, Canada. These data are used to test models that relate the scale of the formative bed forms to the dimensions of the preserved deposits and, therefore, provide an insight as to how such deposits may be preserved over geologic time. The preservation of bed-form geometry is quantified by comparing the Alluvial architecture above and below the maximum erosion depth of the modem channel deposits. This comparison shows that there is no significant difference in the mean set thickness of dune cross-strata above and below the basal erosion surface of the contemporary channel, thus suggesting that dimensional relationships between dune deposits and the formative bed-form dimensions are likely to be valid from both recent and older deposits. The data show that estimates of mean bankfull flow depth derived from dune, unit bar, and cross-bar channel deposits are all very similar. Thus, the use of all these metrics together can provide a useful check that all components and scales of the alluvial architecture have been identified correctly when building reservoir models. The data also highlight several practical issues with identifying and applying data relating to cross-strata. For example, the deposits of unit bars were found to be severely truncated in length and width, with only approximately 10% of the mean bar-form length remaining, and thus making identification in section difficult. For similar reasons, the deposits of compound bars were found to be especially difficult to recognize, and hence, estimates of channel depth based on this method may be problematic. Where only core data are available (i.e., no outcrop data exist), formative flow depths are suggested to be best reconstructed using cross-strata formed by dunes. However, theoretical relationships between the distribution of set thicknesses and formative dune height are found to result in slight overestimates of the latter and, hence, mean bankfull flow depths derived from these measurements. This article illustrates that the preservation of fluvial cross-strata and, thus, the paleohydraulic inferences that can be drawn from them, are a function of the ratio of the size and migration rate of bed forms and the time scale of aggradation and channel migration. These factors must thus be considered when deciding on appropriate length:thickness ratios for the purposes of object-based modeling in reservoir characterization.

Perfect Simulation of a Coupling Achieving the (d)over-bar-distance Between Ordered Pairs of Binary Chains of Infinite Order

We explicitly construct a stationary coupling attaining Ornstein`s (d) over bar -distance between ordered pairs of binary chains of infinite order. Our main tool is a representation of the transition probabilities of the coupled bivariate chain of infinite order as a countable mixture of Markov transition probabilities of increasing order. Under suitable conditions on the loss of memory of the chains, this representation implies that the coupled chain can be represented as a concatenation of i.i.d. sequences of bivariate finite random strings of symbols. The perfect simulation algorithm is based on the fact that we can identify the first regeneration point to the left of the origin almost surely.

Method for determination of vertical bar U(e)3 vertical bar in neutrino oscillation appearance experiments

We point out that determination of the MNS matrix element \U-e3\ = s(13) in long-baseline nu(mu) --> nu(e) neutrino oscillation experiments suffers from large intrinsic uncertainty due to the unknown CP violating phase delta and sign of Deltam(13)(2). We propose a new strategy for accurate determination of theta(13); tune the beam energy at the oscillation maximum and do the measurement both in neutrino and antineutrino channels. We show that it automatically resolves the problem of parameter ambiguities which involves delta, theta(13), and the sign of Deltam(13)(2). (C) 2002 Elsevier B.V. B.V. All rights reserved.

Measurement of the CP-violating phase phi sJ/psi phi using the flavor-tagged decay Bs(0) -> J/psi phi in 8 fb(-1) of p(p)over-bar collisions

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

(X)OVER-BAR CHARTS WITH VARIABLE SAMPLE-SIZE

The usual practice in using a control chart to monitor a process is to take samples of size n from the process every h hours This article considers the properties of the XBAR chart when the size of each sample depends on what is observed in the preceding sample. The idea is that the sample should be large if the sample point of the preceding sample is close to but not actually outside the control limits and small if the sample point is close to the target. The properties of the variable sample size (VSS) XBAR chart are obtained using Markov chains. The VSS XBAR chart is substantially quicker than the traditional XBAR chart in detecting moderate shifts in the process.

Perturbative bosonization from two-point correlation functions

Here we address the problem of bosonizing massive fermions without making expansions in the fermion masses in both massive QED(2) and QED(3) with N fermion flavors including also a Thirring coupling. We start from two-point correlators involving the U(1) fermionic current and the gauge field. From the tensor structure of those correlators we prove that the U(1) current must be identically conserved (topological) in the corresponding bosonized theory in both D=2 and D=3 dimensions. We find an effective generating functional in terms of bosonic fields which reproduces these two-point correlators and from that we obtain a map of the Lagrangian density (ψ) over bar (r)(ipartial derivative-m)psi(r) into a bosonic one in both dimensions. This map is nonlocal but it is independent of the electromagnetic and Thirring couplings, at least in the quadratic approximation for the fermionic determinant.

Some Comments Regarding the Synthetic (X)over-bar Chart

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

The steady-state behavior of the synthetic and side-sensitive synthetic double sampling (X)over-bar charts

The steady-state average run length is used to measure the performance of the recently proposed synthetic double sampling (X) over bar chart (synthetic DS chart). The overall performance of the DS X chart in signaling process mean shifts of different magnitudes does not improve when it is integrated with the conforming run length chart, except when the integrated charts are designed to offer very high protection against false alarms, and the use of large samples is prohibitive. The synthetic chart signals when a second point falls beyond the control limits, no matter whether one of them falls above the centerline and the other falls below it; with the side-sensitive feature, the synthetic chart does not signal when they fall on opposite sides of the centerline. We also investigated the steady-state average run length of the side-sensitive synthetic DS X chart. With the side-sensitive feature, the overall performance of the synthetic DS X chart improves, but not enough to outperform the non-synthetic DS X chart. Copyright (C) 2014 John Wiley &Sons, Ltd.

A new method for calculating the convergent point of a moving group

Context. Convergent point (CP) search methods are important tools for studying the kinematic properties of open clusters and young associations whose members share the same spatial motion. Aims. We present a new CP search strategy based on proper motion data. We test the new algorithm on synthetic data and compare it with previous versions of the CP search method. As an illustration and validation of the new method we also present an application to the Hyades open cluster and a comparison with independent results. Methods. The new algorithm rests on the idea of representing the stellar proper motions by great circles over the celestial sphere and visualizing their intersections as the CP of the moving group. The new strategy combines a maximum-likelihood analysis for simultaneously determining the CP and selecting the most likely group members and a minimization procedure that returns a refined CP position and its uncertainties. The method allows one to correct for internal motions within the group and takes into account that the stars in the group lie at different distances. Results. Based on Monte Carlo simulations, we find that the new CP search method in many cases returns a more precise solution than its previous versions. The new method is able to find and eliminate more field stars in the sample and is not biased towards distant stars. The CP solution for the Hyades open cluster is in excellent agreement with previous determinations.