Genetic diversity in old populations of sessile oak from Calabria assessed by nuclear and chloroplast SSR

In the framework of forest resources conservation, this study aims to understand the dynamic and the genetic structure of sessile oak forests in Calabria, Italy. Two old populations of sessile oak(Quercus petraea(Mattuschka) Liebl.) from two areas of Sila and Aspromonte massifs in Calabria were analyzed for genetic diversity and population structure based on 6 nuclear simple sequence repeat(nSSR) and 4 chloroplastic SSR(cpSSR) loci. The populations displayed high amount of genetic diversity, which was toughly structured according to their geographical origins. Number of alleles at SSR loci ranged from 11 to 20 with an average of 13.5 per locus. Gene diversity(expected heterozygosity, He) estimates ranged from 0.575 to 0.834 with a mean of 0.749. The observed heterozygosity(Ho) was on average 0.458 ranging from 0.150 to 0.682. Polymorphism information content(PIC) values ranged from 0.625 to 0.865 with an average of 0.787. The analysis of molecular variance(AMOVA) highlighted a significant higher estimated variance within populations compared to among populations. Finally, the analysis of haplotypes by using cpSSR suggested a higher diversification in the population from Sila. Hierarchical clustering analysis grouped the genotypes into two major clusters, which agreed with the geographic origin of populations, and was confirmed by the Discriminant Analysis of Principal Components(DAPC). The first cluster included plants/population from Sila massif, while the second encompassed mostly plants/population sampled in Aspromonte massif. Finally, model-based clustering by STRUCTURE analysis also supported the presence of clear genetic structuring in the collection with two major populations(K=2) supported to PCoA analysis as well. Finally, our data indicated the Aspromonte population as a marginal forest with fragmented distribution suggesting different strategies of preservation than in Sila massif.     

EXOTIME: searching for planets around pulsating subdwarf B stars

In 2007, a companion with planetary mass was found around the pulsating subdwarf B star V391 Pegasi with the timing method, indicating that a previously undiscovered population of substellar companions to apparently single subdwarf B stars might exist. Following this serendipitous discovery, the EXOTIME (http://www.na.astro.it/~silvotti/exotime/) monitoring program has been set up to follow the pulsations of a number of selected rapidly pulsating subdwarf B stars on time scales of several years with two immediate observational goals:

1. determine $\dot{P}$ of the pulsational periods P
2. search for signatures of substellar companions in O–C residuals due to periodic light travel time variations, which would be tracking the central star’s companion-induced wobble around the centre of mass

These sets of data should therefore, at the same time, on the one hand be useful to provide extra constraints for classical asteroseismological exercises from the $\dot{P}$ (comparison with “local” evolutionary models), and on the other hand allow one to investigate the preceding evolution of a target in terms of possible “binary” evolution by extending the otherwise unsuccessful search for companions to potentially very low masses. While timing pulsations may be an observationally expensive method to search for companions, it samples a different range of orbital parameters, inaccessible through orbital photometric effects or the radial velocity method: the latter favours massive close-in companions, whereas the timing method becomes increasingly more sensitive toward wider separations. In this paper we report on the status of the on-going observations and coherence analysis for two of the currently five targets, revealing very well-behaved pulsational characteristics in HS?0444+0458, while showing HS?0702+6043 to be more complex than previously thought.     

Upper mantle anisotropy beneath central and southwest Japan: An insight into subduction-induced mantle flow

Analysis of seismic anisotropy in the crust and mantle wedge above subduction zones gives much information about the dynamic processes inside the Earth. For this reason, we measure shear wave polarization anisotropy in the crust and upper mantle beneath central and southwestern Japan from local shallow, intermediate, and deep earthquakes occurring in the subducting Pacific slab. We analyze S phases from 198 earthquakes recorded at 42 Japanese F-net broadband seismic stations. This data set yields a total of 980 splitting parameter pairs for central and southwestern Japan. Dominant fast polarization directions of shear waves obtained at most stations in the Kanto–Izu–Tokai areas are oriented WNW–ESE, which are sub-parallel to the subduction direction of the Pacific plate. However, minor fast polarization directions are oriented in NNE–SSW directions being parallel to the strike of the Japan Trench, especially in the north of Izu Peninsula and the northern Tokai district. Generally, fast directions obtained at stations located in Kii Peninsula and the Chubu district are oriented ENE–WSW, almost parallel to the Nankai Trough, although some fast directions have NW–SE trends. The fast directions obtained at stations in northern central Honshu are oriented N–S. Delay times vary considerably and range from 0.1 to 1.25 s depending on the source depth and the degree of anisotropy along the ray path. These lateral variations in splitting character suggest that the nature of anisotropy is quite different between the studied areas. Beneath Kanto–Tokai, the observed WNW–ESE fast directions are probably caused by the olivine A-fabric induced by the corner flow. However, the slab morphology in this region is relatively complicated as the Philippine Sea slab is overriding the Pacific slab. This complex tectonic setting may induce lateral heterogeneity in the flow and stress state of the mantle wedge, and may have produced NNE–SSW orientations of fast directions. The ENE–WSW fast directions in Kii Peninsula and the Chubu district are more coherent and may be partly induced by the subduction of the Philippine Sea plate. The N–S fast directions in northern central Honshu might be produced by the trench-parallel stretching of the wedge due to the curved slab at the arc–arc junction.     

Rupture process of the Miyagi-Oki,Japan, earthquake of June 12, 1978

The faulting mechanism and multiple rupture process of the M = 7.4 Miyagi-Oki earthquake are studied using surface and body wave data from local and worldwide stations. The main results are as follows. (1) P-wave first motion data and radiation patterns of long-period surface waves indicate a predominantly thrust mechanism with strike N10° E, dip 20°W, and slip angle 76°. The seismic moment is 3.1 × 10²⁷ dyne-cm. (2) Farfield SH waveforms and local seismograms suggest that the rupture occurred in two stages, being concordant with the two zones of aftershock activity revealed by the microearthquake network of Tohoku University. The upper and lower zones, located along the westward-dipping plate interface, are separated by a gap at a depth of 35 km and have dimensions of 37 × 34 and 24 × 34 km², respectively. Rupture initiated at the southern end of the upper aftershock zone and propagated at N20°W subparallel to the trench axis. About 11 s later, the second shock, which was located 30 km landward (westward) of the first, initiated at the upper corner of the lower aftershock zone and propagated down-dip N80°W. Using Haskell modelling for this rupture process, synthetic seismograms were computed for teleseismic SH waves and nearfield body waves. Other parameters determined are: seismic moment $\text{M}{}_0\text{= 1.7 × 10}{}^{27}\text{dyne-cm, slip dislocation}\text{u}\text{= 1.9}\text{m}\text{, Δσ = 95}\text{bar, rupture velocity}\text{ν = 3.2}\text{km s}{}^{\mathrm{?1}}\text{,}\text{rise time}\text{τ = 2}\text{s}$, for the first event; $\text{M}{}_0\text{= 1.4 × 10}{}^{27}\text{dyne-cm}\text{,}\text{u}\text{= 2.4}\text{m}\text{, Δσ = 145}\text{bar}$, for the second event; and time separation between the two shocks ΔT = 11 s. The above two-segment model does not explain well the sharp onsets of the body waves at near-source stations. An initial break of a small subsegment on the upper zone, which propagated down-dip, was hypothesized to explain the observed near-source seismograms. (3) The multiple rupture of the event and the absence of aftershocks between the two fault zones suggests that the frictional and/or sliding characteristics along the plate interface are not uniform. The rupture of the first event was arrested, presumably by a region of high fracture strength between the two zones. The fracture energy of the barrier was estimated to be 10¹⁰ erg cm^?2. (4) The possible occurrence of a large earthquake has been noted for the region adjacent to and seaward of the area that ruptured during the 1978 event. The 1978 event does not appear to reduce the likelihood of occurrence of this expected earthquake.     

The role of pre-existing faults in the structural evolution of thrust systems: Insights from the Ligurian Alps (Italy)

New structural and stratigraphic data for a selected area of the Ligurian Alps are combined in order to assess and discuss the role played by extensional structures in the southernmost segment of the Western Alps during thrusting. Restored cross-sections and field data suggest that the structural style in the external sector of the chain may depend upon the presence of pre-orogenic normal faults ascribed to three extensional events linked to different geodynamic contexts: (i) Permian post-Variscan plate reorganisation, (ii) Mesozoic rifting–drifting phases leading to the opening of the Alpine Tethys, and (iii) Eocenic development of the European foreland basins. During positive inversion in Eocene times, a thin-skinned thrust system developed in this area, followed by a thick-skinned phase. In both situations the inherited extensional structures played fundamental roles: during the thin-skinned phase they conditioned the thrusting sequence, also producing large-scale buckle folds and partial reactivations; during the thick-skinned phase the strain was compartmentalized and partitioned by pre-existing faults.The kinematic model of the external sectors of the Ligurian chain also allows the re-assessment of the Alpine evolution of the front-foreland transition, including: (i) indirect confirmation that in the Eocene the Ligurian Briançonnais and Dauphinois domains were not separated by the Valais-Pyrenean oceanic basin; (ii) that the thin-skinned phase progressively changed into thick-skinned; (iii) the assertion that there were no significant deformations from the Oligocene to the present-day, and the Corsica–Sardinia block rotation only produced a change in orientation of previously formed structures and normal fault system development.     

MHYDAS‐Erosion: a distributed single‐storm water erosion model for agricultural catchments

In this paper, we present MHYDAS‐Erosion, a dynamic and distributed single‐storm water erosion model developed as a module of the existing hydrological MHYDAS model. As with many catchment erosion models, MHYDAS‐Erosion is able to simulate sediment transport, erosion and deposition by rill and interrill processes. Its originality stems from its capacity to integrate the impact of land management practices (LMP) as key elements controlling the sedimentological connectivity in agricultural catchments. To this end, the water‐sediment pathways are first determined by a specific process‐oriented procedure defined and controlled by the user, which makes the integration of LMP easier. The LMP dynamic behaviours are then integrated into the model as a time‐dependent function of hydrological variables and LMP characteristics. The first version of the model was implemented for vegetative filters and tested using water and sediment discharge measurements at three nested scales of a densely instrumented catchment (Roujan, OMERE Observatory, southern France). The results of discharge and soil loss for simulated rainfall events have been found to acceptably compare with available data. The average R² values for water and sediment discharge are 0·82 and 0·83, respectively. The sensitivity of the model to changes in the proportion of LMP was assessed for a single rain event by considering three scenarios of the Roujan catchment management with vegetative filters: 0% (Scenario 1), 18% (Scenario 2, real case) and 100% (Scenario 3). Compared to Scenario 2 (real case), soil losses decreased for Scenario 3 by 65% on the agricultural plot scale, 62% on the sub‐catchment scale and 45% at the outlet of the catchment and increased for Scenario 1 by 0% on the plot scale, 26% on the sub‐catchment scale and 18% at the outlet of the catchment. Copyright © 2010 John Wiley & Sons, Ltd.     

Sources and distribution of organic matter in northern Patagonia fjords,Chile (∼44–47°S): A multi-tracer approach for carbon cycling assessment

We investigated the provenance of organic matter in the inner fjord area of northern Patagonia, Chile (～44–47°S), by studying the elemental (organic carbon, total nitrogen), isotopic (δ¹³C, δ¹⁵N), and biomarker (n-alkanoic acids from vascular plant waxes) composition of surface sediments as well as local marine and terrestrial organic matter. Average end-member values of N/C, δ¹³C, and δ¹⁵N from organic matter were 0.127±0.010, ?19.8±0.3‰, and 9.9±0.5‰ for autochthonous (marine) sources and 0.040±0.018, ?29.3±2.1‰, and 0.2±3.0‰ for allochthonous (terrestrial) sources. Using a mixing equation based on these two end-members, we calculated the relative contribution of marine and terrestrial organic carbon from the open ocean to the heads of fjords close to river outlets. The input of marine-derived organic carbon varied widely and accounted for 13–96% (average 61%) of the organic carbon pool of surface sediments. Integrated regional calculations for the inner fjord system of northern Patagonia covered in this study, which encompasses an area of ～4280 km², suggest that carbon accumulation may account for between 2.3 and 7.8×10⁴ ton C yr^?1. This represents a storage capacity of marine-derived carbon between 1.8 and 6.2×10⁴ ton yr^?1, which corresponds to an assimilation rate of CO₂ by marine photosynthesis between 0.06 and 0.23×10⁶ ton yr^?1. This rate suggests that the entire fjord system of Patagonia, which covers an area of ～240,000 km², may represent a potentially important region for the global burial of marine organic matter and the sequestration of atmospheric CO₂.     

Monitoring in the Western Pacific region shows evidence of seagrass decline in line with global trends

Seagrass systems of the Western Pacific region are biodiverse habitats, providing vital services to ecosystems and humans over a vast geographic range. SeagrassNet is a worldwide monitoring program that collects data on seagrass habitats, including the ten locations across the Western Pacific reported here where change at various scales was rapidly detected. Three sites remote from human influence were stable. Seagrasses declined largely due to increased nutrient loading (4 sites) and increased sedimentation (3 sites), the two most common stressors of seagrass worldwide. Two sites experienced near-total loss from of excess sedimentation, followed by partial recovery once sedimentation was reduced. Species shifts were observed at every site with recovering sites colonized by pioneer species. Regulation of watersheds is essential if marine protected areas are to preserve seagrass meadows. Seagrasses in the Western Pacific experience stress due to human impacts despite the vastness of the ocean area and low development pressures.     

Mechanical stratigraphy as a factor controlling the development of a sandbox transfer zone: a three-dimensional analysis

In thrust belts, fold–fault terminations are common features of the structural architecture and can pose complicated problems to unravel, in particular when two or more terminations are in close proximity. Such terminations usually reflect pre-existing attributes. Amongst the many factors, lateral variations in the mechanical stratigraphy can control along-strike geometry and kinematics of fault-related folds.A displacement transfer zone was produced in a compressional sandbox model by means of two adjacent, mechanically different stratigraphic domains. The experiment allowed two discrete chains to develop in the different domains, so that a complex structural setting occurred in the connecting area. Periclinal folds, oblique thrust fronts and oblique ramps developed in the resulting transfer zone. The interaction between periclines in the transfer zone produced lateral culminations in the folded structures. The analysis of displacement across the structural domains revealed that a significant loss of slip along the faults occurred in the relay zone. In this area, imbricate faulting was partially replaced by layer-parallel shortening. A linear relationship appears to exist between the bed length of the thrust sheet and the related fault slip.     

Assessment of meteorological climate models as inputs for coastal studies

Modeling studies of future changes in coastal hydrodynamics, in terms of storm surges and wave climate, need appropriate wind and atmospheric forcings, a necessary requirement for the realistic reproduction of the statistics and the resolution of small scale features. This work compares meteorological results from different climate models in the Mediterranean area, with a focus on the Adriatic Sea, in order to assess their capability to reproduce coastal meteorological features and their possibility to be used as forcings for hydrodynamic simulations. Five meteorological datasets are considered. They are obtained from two regional climate models, implemented with different spatial resolutions and setups and are downscaled from two different global climate models. Wind and atmospheric pressure fields are compared with measurements at four stations along the Italian Adriatic coast. The analysis is carried out both on simulations of the control period 1960–1990 and on the A1B Intergovernmental Panel for Climate Change scenario projections (2070–2100), highlighting the ability of each model in reproducing the statistical coastal meteorological behavior and possible changes. The importance of simulated global- and regional-scale meteorological processes, in terms of correct spatial resolution of the phenomena, is also discussed. Within the Adriatic Sea, the meteorological climate is influenced by the local orography that controls the strengthening of north-eastern katabatic winds like Bora. Results show indeed that the increase in spatial resolution provides a more realistic wind forcing for the hydrodynamic simulations. Moreover, the chosen setup and the global climate models that drive the regional downscalings appear to play an important role in reproducing correct atmospheric pressure fields. The comparison between scenario and control simulations shows a small increase in the mean atmospheric pressure values, while a decrease in mean wind speed and in extreme wind events is observed, particularly for the datasets with higher spatial resolution. Finally, results suggest that an ensemble of downscaled climate models is likely to provide the most suitable climatic forcings (wind and atmospheric pressure fields) for coastal hydrodynamic modeling.