SIS: a program to generate draft genome sequence scaffolds for prokaryotes

Background: Decreasing costs of DNA sequencing have made prokaryotic draft genome sequences increasingly common. A contig scaffold is an ordering of contigs in the correct orientation. A scaffold can help genome comparisons and guide gap closure efforts. One popular technique for obtaining contig scaffolds is to map contigs onto a reference genome. However, rearrangements that may exist between the query and reference genomes may result in incorrect scaffolds, if these rearrangements are not taken into account. Large-scale inversions are common rearrangement events in prokaryotic genomes. Even in draft genomes it is possible to detect the presence of inversions given sufficient sequencing coverage and a sufficiently close reference genome. Results: We present a linear-time algorithm that can generate a set of contig scaffolds for a draft genome sequence represented in contigs given a reference genome. The algorithm is aimed at prokaryotic genomes and relies on the presence of matching sequence patterns between the query and reference genomes that can be interpreted as the result of large-scale inversions; we call these patterns inversion signatures. Our algorithm is capable of correctly generating a scaffold if at least one member of every inversion signature pair is present in contigs and no inversion signatures have been overwritten in evolution. The algorithm is also capable of generating scaffolds in the presence of any kind of inversion, even though in this general case there is no guarantee that all scaffolds in the scaffold set will be correct. We compare the performance of SIS, the program that implements the algorithm, to seven other scaffold-generating programs. The results of our tests show that SIS has overall better performance. Conclusions: SIS is a new easy-to-use tool to generate contig scaffolds, available both as stand-alone and as a web server. The good performance of SIS in our tests adds evidence that large-scale inversions are widespread in prokaryotic genomes.

Calibration of a Species-Specific Spectral Vegetation Index for Leaf Area Index (LAI) Monitoring: Example with MODIS Reflectance Time-Series on Eucalyptus Plantations

The leaf area index (LAI) is a key characteristic of forest ecosystems. Estimations of LAI from satellite images generally rely on spectral vegetation indices (SVIs) or radiative transfer model (RTM) inversions. We have developed a new and precise method suitable for practical application, consisting of building a species-specific SVI that is best-suited to both sensor and vegetation characteristics. Such an SVI requires calibration on a large number of representative vegetation conditions. We developed a two-step approach: (1) estimation of LAI on a subset of satellite data through RTM inversion; and (2) the calibration of a vegetation index on these estimated LAI. We applied this methodology to Eucalyptus plantations which have highly variable LAI in time and space. Previous results showed that an RTM inversion of Moderate Resolution Imaging Spectroradiometer (MODIS) near-infrared and red reflectance allowed good retrieval performance (R-2 = 0.80, RMSE = 0.41), but was computationally difficult. Here, the RTM results were used to calibrate a dedicated vegetation index (called "EucVI") which gave similar LAI retrieval results but in a simpler way. The R-2 of the regression between measured and EucVI-simulated LAI values on a validation dataset was 0.68, and the RMSE was 0.49. The additional use of stand age and day of year in the SVI equation slightly increased the performance of the index (R-2 = 0.77 and RMSE = 0.41). This simple index opens the way to an easily applicable retrieval of Eucalyptus LAI from MODIS data, which could be used in an operational way.

Methane airborne measurements and comparison to global models during BARCA

Tropical regions, especially the Amazon region, account for large emissions of methane (CH4). Here, we present CH4 observations from two airborne campaigns conducted within the BARCA (Balanco Atmosferico Regional de Carbono na Amazonia) project in the Amazon basin in November 2008 (end of the dry season) and May 2009 (end of the wet season). We performed continuous measurements of CH4 onboard an aircraft for the first time in the Amazon region, covering the whole Amazon basin with over 150 vertical profiles between altitudes of 500 m and 4000 m. The observations support the finding of previous ground-based, airborne, and satellite measurements that the Amazon basin is a large source of atmospheric CH4. Isotope analysis verified that the majority of emissions can be attributed to CH4 emissions from wetlands, while urban CH4 emissions could be also traced back to biogenic origin. A comparison of five TM5 based global CH4 inversions with the observations clearly indicates that the inversions using SCIAMACHY observations represent the BARCA observations best. The calculated CH4 flux estimate obtained from the mismatch between observations and TM5-modeled CH4 fields ranges from 36 to 43 mg m(-2) d(-1) for the Amazon lowland region.

The rupture process of recent tsunamigenic earthquakes by geophysical data inversion

Subduction zones are the favorite places to generate tsunamigenic earthquakes, where friction between oceanic and continental plates causes the occurrence of a strong seismicity. The topics and the methodologies discussed in this thesis are focussed to the understanding of the rupture process of the seismic sources of great earthquakes that generate tsunamis. The tsunamigenesis is controlled by several kinematical characteristic of the parent earthquake, as the focal mechanism, the depth of the rupture, the slip distribution along the fault area and by the mechanical properties of the source zone. Each of these factors plays a fundamental role in the tsunami generation. Therefore, inferring the source parameters of tsunamigenic earthquakes is crucial to understand the generation of the consequent tsunami and so to mitigate the risk along the coasts. The typical way to proceed when we want to gather information regarding the source process is to have recourse to the inversion of geophysical data that are available. Tsunami data, moreover, are useful to constrain the portion of the fault area that extends offshore, generally close to the trench that, on the contrary, other kinds of data are not able to constrain. In this thesis I have discussed the rupture process of some recent tsunamigenic events, as inferred by means of an inverse method. I have presented the 2003 Tokachi-Oki (Japan) earthquake (Mw 8.1). In this study the slip distribution on the fault has been inferred by inverting tsunami waveform, GPS, and bottom-pressure data. The joint inversion of tsunami and geodetic data has revealed a much better constrain for the slip distribution on the fault rather than the separate inversions of single datasets. Then we have studied the earthquake occurred on 2007 in southern Sumatra (Mw 8.4). By inverting several tsunami waveforms, both in the near and in the far field, we have determined the slip distribution and the mean rupture velocity along the causative fault. Since the largest patch of slip was concentrated on the deepest part of the fault, this is the likely reason for the small tsunami waves that followed the earthquake, pointing out how much the depth of the rupture plays a crucial role in controlling the tsunamigenesis. Finally, we have presented a new rupture model for the great 2004 Sumatra earthquake (Mw 9.2). We have performed the joint inversion of tsunami waveform, GPS and satellite altimetry data, to infer the slip distribution, the slip direction, and the rupture velocity on the fault. Furthermore, in this work we have presented a novel method to estimate, in a self-consistent way, the average rigidity of the source zone. The estimation of the source zone rigidity is important since it may play a significant role in the tsunami generation and, particularly for slow earthquakes, a low rigidity value is sometimes necessary to explain how a relatively low seismic moment earthquake may generate significant tsunamis; this latter point may be relevant for explaining the mechanics of the tsunami earthquakes, one of the open issues in present day seismology. The investigation of these tsunamigenic earthquakes has underlined the importance to use a joint inversion of different geophysical data to determine the rupture characteristics. The results shown here have important implications for the implementation of new tsunami warning systems – particularly in the near-field – the improvement of the current ones, and furthermore for the planning of the inundation maps for tsunami-hazard assessment along the coastal area.

Investigations on maar-diatreme volcanoes by inversion of magnetic and gravity data from the Eifel area, Germany

A study of maar-diatreme volcanoes has been perfomed by inversion of gravity and magnetic data. The geophysical inverse problem has been solved by means of the damped nonlinear least-squares method. To ensure stability and convergence of the solution of the inverse problem, a mathematical tool, consisting in data weighting and model scaling, has been worked out. Theoretical gravity and magnetic modeling of maar-diatreme volcanoes has been conducted in order to get information, which is used for a simple rough qualitative and/or quantitative interpretation. The information also serves as a priori information to design models for the inversion and/or to assist the interpretation of inversion results. The results of theoretical modeling have been used to roughly estimate the heights and the dip angles of the walls of eight Eifel maar-diatremes — each taken as a whole. Inversemodeling has been conducted for the Schönfeld Maar (magnetics) and the Hausten-Morswiesen Maar (gravity and magnetics). The geometrical parameters of these maars, as well as the density and magnetic properties of the rocks filling them, have been estimated. For a reliable interpretation of the inversion results, beside the knowledge from theoretical modeling, it was resorted to other tools such like field transformations and spectral analysis for complementary information. Geologic models, based on thesynthesis of the respective interpretation results, are presented for the two maars mentioned above. The results gave more insight into the genesis, physics and posteruptive development of the maar-diatreme volcanoes. A classification of the maar-diatreme volcanoes into three main types has been elaborated. Relatively high magnetic anomalies are indicative of scoria cones embeded within maar-diatremes if they are not caused by a strong remanent component of the magnetization. Smaller (weaker) secondary gravity and magnetic anomalies on the background of the main anomaly of a maar-diatreme — especially in the boundary areas — are indicative for subsidence processes, which probably occurred in the late sedimentation phase of the posteruptive development. Contrary to postulates referring to kimberlite pipes, there exists no generalized systematics between diameter and height nor between geophysical anomaly and the dimensions of the maar-diatreme volcanoes. Although both maar-diatreme volcanoes and kimberlite pipes are products of phreatomagmatism, they probably formed in different thermodynamic and hydrogeological environments. In the case of kimberlite pipes, large amounts of magma and groundwater, certainly supplied by deep and large reservoirs, interacted under high pressure and temperature conditions. This led to a long period phreatomagmatic process and hence to the formation of large structures. Concerning the maar-diatreme and tuff-ring-diatreme volcanoes, the phreatomagmatic process takes place due to an interaction between magma from small and shallow magma chambers (probably segregated magmas) and small amounts of near-surface groundwater under low pressure and temperature conditions. This leads to shorter time eruptions and consequently to structures of smaller size in comparison with kimberlite pipes. Nevertheless, the results show that the diameter to height ratio for 50% of the studied maar-diatremes is around 1, whereby the dip angle of the diatreme walls is similar to that of the kimberlite pipes and lies between 70 and 85°. Note that these numerical characteristics, especially the dip angle, hold for the maars the diatremes of which — estimated by modeling — have the shape of a truncated cone. This indicates that the diatreme can not be completely resolved by inversion.

Advanced modelling of time domain electromagnetic data with updated hydrogeological interpretations

Several countries have acquired, over the past decades, large amounts of area covering Airborne Electromagnetic data. Contribution of airborne geophysics has dramatically increased for both groundwater resource mapping and management proving how those systems are appropriate for large-scale and efficient groundwater surveying. We start with processing and inversion of two AEM dataset from two different systems collected over the Spiritwood Valley Aquifer area, Manitoba, Canada respectively, the AeroTEM III (commissioned by the Geological Survey of Canada in 2010) and the “Full waveform VTEM” dataset, collected and tested over the same survey area, during the fall 2011. We demonstrate that in the presence of multiple datasets, either AEM and ground data, due processing, inversion, post-processing, data integration and data calibration is the proper approach capable of providing reliable and consistent resistivity models. Our approach can be of interest to many end users, ranging from Geological Surveys, Universities to Private Companies, which are often proprietary of large geophysical databases to be interpreted for geological and\or hydrogeological purposes. In this study we deeply investigate the role of integration of several complimentary types of geophysical data collected over the same survey area. We show that data integration can improve inversions, reduce ambiguity and deliver high resolution results. We further attempt to use the final, most reliable output resistivity models as a solid basis for building a knowledge-driven 3D geological voxel-based model. A voxel approach allows a quantitative understanding of the hydrogeological setting of the area, and it can be further used to estimate the aquifers volumes (i.e. potential amount of groundwater resources) as well as hydrogeological flow model prediction. In addition, we investigated the impact of an AEM dataset towards hydrogeological mapping and 3D hydrogeological modeling, comparing it to having only a ground based TEM dataset and\or to having only boreholes data.

Zur Karbonatsedimentologie, Biofazies und sequenzstratigraphischen Architektur eines fossilen Hochenergie-Schelfs aus dem Neogen der Algarve (Miozän, Südportugal)

Die neogene Lagos-Portimão Formation (Unter- bis Mittelmiozän) bildet einen Teil der Steilküste der Algarve (S-Portugal) und besteht aus einer zyklischen Wechsellagerung von Karbonaten und Sand-steinen. Die vorliegende Arbeit bietet ein Modell zur sedimentologischen, faziellen und stratigraphischen Entwicklung dieser Einheit an. Basierend auf Profilen entlang der gesamten lateralen Erstreckung der Einheit wurden verschiedene Gelände- und Labormethoden angewandt, um ein Modell entwickeln zu können. Messungen des Sr87/86-Isotopenverhältnisses sollten Klarheit bezüglich der stratigraphischen Position bringen. Die laterale Korrelation der Profile erfolgte über lithologische und fazielle Ansprachen. Unterstützend wurden einzelne Profile mit einem tragbaren Gammaray-Spektrometer gemessen. Es wurden vier Leithorizonte etabliert, die sich durch fazielle Merkmale und spezielle Fossilführung defi-nieren lassen. Die Mikrofazies wurde qualitativ und quantitativ analysiert. Als statistisches Verfahren wurde unter anderem eine hierarchische Clusteranalyse durchgeführt, über welche fünf Biofaziestypen des warm-temperierten Klimabereichs unterschieden werden. Die Fossilführung wird von Mollusken, Bryozoen und Rotalgen dominiert. Ausnahmen bilden stratigraphisch isolierte Vorkommen von kolo-nialen Korallen, die jedoch keine Riffkörper aufbauen. Die Ergebnisse aller zuvor erwähnten Untersuchungen deuten auf Ablagerungen eines nicht-tropischen Hochenergie-Schelfs hin. Sedimentäre Zyklen sind oftmals unvollständig, es treten Hartgründe und Auf-arbeitungs- bzw- Kondesationshorizonte auf. Die geochemische Altersdatierung weist Alterssprünge und -inversionen auf. Ein Vergleich mit dem SW-australischen Schelf und dem von James et al. (1994) eingeführten Modell des shaved shelf bietet sich aufgrund der Ähnlichkeit der Sedimentgesteine und des ozeanographischen Settings an. Weiterhin werden zeitgleiche bzw. faziell ähnliche Becken vergleichend diskutiert. Das Sedimentationsgeschehen der Lagos-Portimão Formation wird maßgeblich durch eine halokinetisch bedingte unregelmäßige Subsidenz und Hebung beeinflußt. Der Salzdom von Albufeira war während der Sedimentation der Einheit mehrfach in Bewegung. Rutschungspakete, Entlastungsspalten und Sanddikes zeugen davon. Die sequenzstratigraphische Interpretation bietet einen neuen Ansatz, in dem sie von Hochstand-Sandsteinen und Tiefstand-Karbonaten ausgeht.

Constraining global methane emissions and uptake by ecosystems

Natural methane (CH4) emissions from wet ecosystems are an important part of today's global CH4 budget. Climate affects the exchange of CH4 between ecosystems and the atmosphere by influencing CH4 production, oxidation, and transport in the soil. The net CH4 exchange depends on ecosystem hydrology, soil and vegetation characteristics. Here, the LPJ-WHyMe global dynamical vegetation model is used to simulate global net CH4 emissions for different ecosystems: northern peatlands (45°–90° N), naturally inundated wetlands (60° S–45° N), rice agriculture and wet mineral soils. Mineral soils are a potential CH4 sink, but can also be a source with the direction of the net exchange depending on soil moisture content. The geographical and seasonal distributions are evaluated against multi-dimensional atmospheric inversions for 2003–2005, using two independent four-dimensional variational assimilation systems. The atmospheric inversions are constrained by the atmospheric CH4 observations of the SCIAMACHY satellite instrument and global surface networks. Compared to LPJ-WHyMe the inversions result in a~significant reduction in the emissions from northern peatlands and suggest that LPJ-WHyMe maximum annual emissions peak about one month late. The inversions do not put strong constraints on the division of sources between inundated wetlands and wet mineral soils in the tropics. Based on the inversion results we diagnose model parameters in LPJ-WHyMe and simulate the surface exchange of CH4 over the period 1990–2008. Over the whole period we infer an increase of global ecosystem CH4 emissions of +1.11 Tg CH4 yr−1, not considering potential additional changes in wetland extent. The increase in simulated CH4 emissions is attributed to enhanced soil respiration resulting from the observed rise in land temperature and in atmospheric carbon dioxide that were used as input. The long-term decline of the atmospheric CH4 growth rate from 1990 to 2006 cannot be fully explained with the simulated ecosystem emissions. However, these emissions show an increasing trend of +3.62 Tg CH4 yr−1 over 2005–2008 which can partly explain the renewed increase in atmospheric CH4 concentration during recent years.

FISH mapping of 10 canine BAC clones harbouring genes and microsatellites in the arctic fox and the Chinese raccoon dog genomes

Cytogenetic mapping of the arctic fox and the Chinese raccoon dog were performed using a set of canine probes derived from the Bacterial Artificial Chromosome (BAC) library. Altogether, 10 BAC clones containing sequences of selected genes (PAX3, HBB, ATP2A2, TECTA, PIT1, ABCA4, ESR2, TPH1, HTR2A, MAOA) and microsatellites were mapped by fluorescence in situ hybridization (FISH) experiments to chromosomes of the canids studied. At present, the cytogenetic map on the arctic fox and Chinese raccoon dog consists of 45 loci each. Chromosomal localization of the BAC clones was in agreement with data obtained by earlier independent comparative chromosome painting. However, two events of telomere-to-centromere inversions were tentatively identified while compared with assignments in the dog karyotype.

THE EFFECTS OF IRRADIATION ON HOT JOVIAN ATMOSPHERES: HEAT REDISTRIBUTION AND ENERGY DISSIPATION

Hot Jupiters, due to the proximity to their parent stars, are subjected to a strong irradiating flux that governs their radiative and dynamical properties. We compute a suite of three-dimensional circulation models with dual-band radiative transfer, exploring a relevant range of irradiation temperatures, both with and without temperature inversions. We find that, for irradiation temperatures T irr lsim 2000 K, heat redistribution is very efficient, producing comparable dayside and nightside fluxes. For T irr ≈ 2200-2400 K, the redistribution starts to break down, resulting in a high day-night flux contrast. Our simulations indicate that the efficiency of redistribution is primarily governed by the ratio of advective to radiative timescales. Models with temperature inversions display a higher day-night contrast due to the deposition of starlight at higher altitudes, but we find this opacity-driven effect to be secondary compared to the effects of irradiation. The hotspot offset from the substellar point is large when insolation is weak and redistribution is efficient, and decreases as redistribution breaks down. The atmospheric flow can be potentially subjected to the Kelvin-Helmholtz instability (as indicated by the Richardson number) only in the uppermost layers, with a depth that penetrates down to pressures of a few millibars at most. Shocks penetrate deeper, down to several bars in the hottest model. Ohmic dissipation generally occurs down to deeper levels than shock dissipation (to tens of bars), but the penetration depth varies with the atmospheric opacity. The total dissipated Ohmic power increases steeply with the strength of the irradiating flux and the dissipation depth recedes into the atmosphere, favoring radius inflation in the most irradiated objects. A survey of the existing data, as well as the inferences made from them, reveals that our results are broadly consistent with the observational trends.

Usefulness of electrocardiographic parameters for risk prediction in arrhythmogenic right ventricular dysplasia

The value of electrocardiographic findings predicting adverse outcome in patients with arrhythmogenic right ventricular dysplasia (ARVD) is not well known. We hypothesized that ventricular depolarization and repolarization abnormalities on the 12-lead surface electrocardiogram (ECG) predict adverse outcome in patients with ARVD. ECGs of 111 patients screened for the 2010 ARVD Task Force Criteria from 3 Swiss tertiary care centers were digitized and analyzed with a digital caliper by 2 independent observers blinded to the outcome. ECGs were compared in 2 patient groups: (1) patients with major adverse cardiovascular events (MACE: a composite of cardiac death, heart transplantation, survived sudden cardiac death, ventricular fibrillation, sustained ventricular tachycardia, or arrhythmic syncope) and (2) all remaining patients. A total of 51 patients (46%) experienced MACE during a follow-up period with median of 4.6 years (interquartile range 1.8 to 10.0). Kaplan-Meier analysis revealed reduced times to MACE for patients with repolarization abnormalities according to Task Force Criteria (p = 0.009), a precordial QRS amplitude ratio (∑QRS mV V1 to V3/∑QRS mV V1 to V6) of ≤ 0.48 (p = 0.019), and QRS fragmentation (p = 0.045). In multivariable Cox regression, a precordial QRS amplitude ratio of ≤ 0.48 (hazard ratio [HR] 2.92, 95% confidence interval [CI] 1.39 to 6.15, p = 0.005), inferior leads T-wave inversions (HR 2.44, 95% CI 1.15 to 5.18, p = 0.020), and QRS fragmentation (HR 2.65, 95% CI 1.1 to 6.34, p = 0.029) remained as independent predictors of MACE. In conclusion, in this multicenter, observational, long-term study, electrocardiographic findings were useful for risk stratification in patients with ARVD, with repolarization criteria, inferior leads TWI, a precordial QRS amplitude ratio of ≤ 0.48, and QRS fragmentation constituting valuable variables to predict adverse outcome.

Stable oxygen isotope record from Greenland ice cores

Twenty ice cores drilled in medium to high accumulation areas of the Greenland ice sheet have been used to extract seasonally resolved stable isotope records. Relationships between the seasonal stable isotope data and Greenland and Icelandic temperatures as well as atmospheric flow are investigated for the past 150-200 years. The winter season stable isotope data are found to be influenced by the North Atlantic Oscillation (NAO) and very closely related to SW Greenland temperatures. The linear correlation between the first principal component of the winter season stable isotope data and Greenland winter temperatures is 0.71 for seasonally resolved data and 0.83 for decadally filtered data. The summer season stable isotope data display higher correlations with Stykkisholmur summer temperatures and North Atlantic SST conditions than with SW Greenland temperatures. The linear correlation between Stykkisholmur summer temperatures and the first principal component of the summer season stable isotope data is 0.56, increasing to 0.66 for decadally filtered data. Winter season stable isotope data from ice core records that reach more than 1400 years back in time suggest that the warm period that began in the 1920s raised southern Greenland temperatures to the same level as those that prevailed during the warmest intervals of the Medieval Warm Period some 900-1300 years ago. This observation is supported by a southern Greenland ice core borehole temperature inversion. As Greenland borehole temperature inversions are found to correspond better with winter stable isotope data than with summer or annual average stable isotope data it is suggested that a strong local Greenland temperature signal can be extracted from the winter stable isotope data even on centennial to millennial time scales.

Magnetic petrology of ODP Leg 121 holes

Given the importance of the inversion of seamount magnetic anomalies, particularly to the motion of the Pacific plate, it is important to gain a better understanding of the nature of the magnetic source of these features. Although different in detail, Ninetyeast Ridge is composed of submarine and subaerial igneous rocks that are similar to those found at many seamounts, making it a suitable proxy. We report here on the magnetic petrology of a collection of samples from Ninetyeast Ridge in the Indian Ocean. Our purpose is to determine the relationship between primary petrology, subsequent alteration, and magnetic properties of the recovered rocks. Such information will eventually lead to a more complete understanding of the magnetization of seamounts and presumably improvements in the accuracy of anomaly inversions. Three basement sites were drilled on Ninetyeast Ridge, with recovery of subaerial basalt flows at the first two (Sites 756 and 757) and submarine massive and pillow flows at the final one (Site 758). The three sites were distinctly different. Site 756 was dominated by ilmenite. What titanomagnetite was present had undergone deuteric alteration and secondary hematite was present in many samples. The magnetization was moderate and stable although it yielded a paleolatitude somewhat lower than expected. Site 757 was highly oxidized, presumably while above sea level. It was dominated by primary titanomagnetite, which was deuterically altered. Secondary hematite was common. Magnetization was relatively weak but quite stable. The paleolatitude for all but the lowermost flows was approximately 40° lower than expected. Site 758 was also dominated by primary titanomagnetite. There was relatively little oxidation with most primary titanomagnetite showing no evidence of high-temperature alteration. No secondary hematite was in evidence. This site had the highest magnetization of the three (although somewhat low relative to other seamounts) but was relatively unstable with significant viscous remanence in many samples. Paleolatitude was close to the expected value. It is not possible, at present, to confidently associate these rocks with specific locations in a seamount structure. A possible and highly speculative model would place rocks similar to Site 757 near the top of the edifice, Site 756 lower down but still erupted above sea level, and Site 758 underlying these units, erupted while the seamount was still below sea level.

Turbiditic trench deposits at the South-Chilean active margin

The active plate margin of South America is characterized by a frequent occurrence of large and devastating subduction earthquakes. Here we focus on marine sedimentary records off Southern Chile that are archiving the regional paleoseismic history over the Holocene and Late Pleistocene. The investigated records - Ocean Drilling Program (ODP) Site 1232 and SONNE core 50SL - are located at ~40°S and ~38°S, within the Perú-Chile trench, and are characterized by frequent interbedded strata of turbiditic and hemipelagic origin. On the basis of the sedimentological characteristics and the association with the active margin of Southern Chile, we assume that the turbidites are mainly seismically triggered, and may be considered as paleo-megaearthquake indicators. However, the long-term changes in turbidite recurrence times appear to be strongly influenced by climate and sea level changes as well. During sea level highstands in the Holocene and Marine Isotope Stage (MIS) 5, recurrence times of turbiditic layers are substantially higher, primarily reflecting a climate-induced reduction of sediment availability and enhanced slope stability. In addition, segmented tectonic uplift changes and related drainage inversions likely influenced the postglacial decrease in turbidite frequencies. Glacial turbidite recurrence times (including MIS 2, MIS 3, cold substages of MIS 5, and MIS 6), on the other hand, are within the same order of magnitude as earthquake recurrence times derived from the historical record and other terrestrial paleoseismic archives of the region. Only during these cold stages sediment availability and slope instability were high enough to enable recording of the complete sequence of large earthquakes in Southern Chile. Our data thus suggest that earthquake recurrence times on the order of 100 to 200 years are a persistent feature at least during the last glacial period.

IBAMar DATABASE: 4 decades of oceanographic sampling on the Western Mediterranean Sea

IBAMar (http://www.ba.ieo.es/ibamar) is a regional database that puts together all physical and biochemical data obtained by multiparametric probes (CTDs equipped with different sensors), during the cruises managed by the Balearic Center of the Spanish Institute of Oceanography (COB-IEO). It has been recently extended to include data obtained with classical hydro casts using oceanographic Niskin or Nansen bottles. The result is a database that includes a main core of hydrographic data: temperature (T), salinity (S), dissolved oxygen (DO), fluorescence and turbidity; complemented by bio-chemical data: dissolved inorganic nutrients (phosphate, nitrate, nitrite and silicate) and chlorophyll-a. In IBAMar Database, different technologies and methodologies were used by different teams along the four decades of data sampling in the COB-IEO. Despite of this fact, data have been reprocessed using the same protocols, and a standard QC has been applied to each variable. Therefore it provides a regional database of homogeneous, good quality data. Data acquisition and quality control (QC): 94% of the data are CTDs Sbe911 and Sbe25. S and DO were calibrated on board using water samples, whenever a Rossetta was available (70% of the cases). All CTD data from Seabird CTDs were reviewed and post processed with the software provided by Sea-Bird Electronics. Data were averaged to get 1 dbar vertical resolution. General sampling methodology and pre processing are described in https://ibamardatabase.wordpress.com/home/). Manual QC include visual checks of metadata, duplicate data and outliers. Automatic QC include range check of variables by area (north of Balearic Islands, south of BI and Alboran Sea) and depth (27 standard levels), check for spikes and check for density inversions. Nutrients QC includes a preliminary control and a range check on the observed level of the data to detect outliers around objectively analyzed data fields. A quality flag is assigned as an integer number, depending on the result of the QC check.