Peri-Gondwanan Ordovician crustal fragments in the high-grade basement of the Eastern Rhodope Massif,Bulgaria: evidence from U-Pb LA-ICP-MS zircon geochronology and geochemistry

Field, geochemical, and geochronologic data of high-grade basement metamafic and evolved rocks are used to identify the nature and timing of pre-Alpine crustal growth of the Rhodope Massif. These rocks occur intrusive into clastic-carbonate metasedimentary succession. Petrography and mineral chemistry show compositions consistent with Alpine amphibolite-facies metamorphism that obliterated the original igneous textures of the protoliths. Bulk-rock geochemistry identifies low-Ti tholeiitic to calc-alkaline gabbroic-basaltic and plagiogranite precursors, with MORB-IAT supra-subduction zone signature and trace elements comparable to modern back-arc basalts. The U-Pb zircon dating revealed a mean age of 455 Ma for the magmatic crystallization of the protoliths that contain inherited Cambrian (528–534 Ma) zircons. Carboniferous, Jurassic, and Eocene metamorphic events overprinted the Ordovician protoliths. The radiometric results of the metamorphic rocks demonstrate that Ordovician oceanic crust was involved in the build-up of the Rhodope high-grade basement. Dating of Eocene-Oligocene volcanic rocks overlying or cross-cutting the metamorphic rocks supplied Neoproterozoic, Ordovician and Permo-Carboniferous xenocrystic zircons that were sampled en route to the surface from the basement. The volcanic rocks thus confirm sub-regionally present Neoproterozoic and Paleozoic igneous and metamorphic basement. We interpret the origin of the Middle-Late Ordovician oceanic magmatism in a back-arc rift-spreading center propagating along peri-Gondwanan Cadomian basement terrane related to the Rheic Ocean widening. The results highlight the presence of elements of Cadomian northern Gondwana margin in the high-grade basement and record of Rheic Ocean evolution. The eastern Rhodope Massif high-grade basement compared to adjacent terranes with Neoproterozoic and Cambro-Ordovician evolution shares analogous tectono-magmatic record providing a linkage among basement terranes incorporated in the Alpine belt of the north Aegean region.     

Sea Ice Dynamics Induced by External Stochastic Fluctuations

The influence of stochastic fluctuations in the atmosphere and in the ocean caused by different occasional phenomena (noises) on dynamic processes of sea ice growth with a mushy layer is studied. It is shown that atmospheric temperature variances substantially increase the sea ice thickness, whereas dispersion variations of turbulent flows in the ocean to a great extent decrease the ice content produced by false bottom evolution.     

The tidal and wind induced hydrodynamics of the composite system Adriatic Sea/Lagoon of Venice

A shallow water hydrostatic 2D hydrodynamic numerical model, based on the boundary conforming coordinate system, was used to simulate aspects of both general and small scale oceanic features occurring in the composite system constituted by the Adriatic Sea and the Lagoon of Venice (Italy), under the influence of tide and realistic atmospheric forcing. Due to a specific technique for the treatment of movable lateral boundaries, the model is able to simulate efficiently dry up and flooding processes within the lagoon. Firstly, a model calibration was performed by comparing the results of the model, forced using tides and ECMWF atmospheric pressure and wind fields, with observations collected for a set of 33 mareographic stations uniformly distributed in the Adriatic Sea and in the Lagoon of Venice. A second numerical experiment was then carried out by considering only the tidal forcing. Through a comparison between the results obtained in the two experiments it was possible to assess the reliability of the estimated parameter through the composite forcing. Model results were then verified by comparing simulated amplitude and phase of each tidal constituent as well as tidal velocities simulated at the inlets of the lagoon and in the Northern Adriatic Sea with the corresponding observed values. The model accurately reproduces the observed harmonics: mean amplitude differences rarely exceed 1 cm, while phase errors are commonly confined below 15°. Semidiurnal and diurnal currents were correctly reproduced in the northern basin and a good agreement was obtained with measurements carried out at the lagoon inlets. On this basis, the outcomes of the hydrodynamic model were analyzed in order to investigate: (i) small-scale coastal circulation features observed at the interface between the adjoining basins, which consist often of vortical dipoles connected with the tidal flow of Adriatic water entering and leaving the Lagoon of Venice and with along-shore current fields connected with specific wind patterns; (ii) residual oscillations, which are often connected to meteorological forcing over the basin. In particular, it emerges that small-scale vortical features generated near the lagoon inlet can be efficiently transported toward the open sea, thus contributing to the water exchange between the two marine regions, and a realistic representation of observed residual oscillations in the area would require a very detailed knowledge of atmospheric as well as remote oceanic forcing.     

Geochemistry of late Cenozoic lavas on Kunashir Island, Kurile Arc

Middle Miocene to Quaternary lavas on Kunashir Island in the southern zone of the Kurile Arc were examined for major, trace, and Sr–Nd–Pb isotope compositions. The lavas range from basalt through to rhyolite and the mafic lavas show typical oceanic island arc signatures without significant crustal or sub-continental lithosphere contamination. The lavas exhibit across-arc variation, with increasingly greater fluid-immobile incompatible element contents from the volcanic front to the rear-arc; this pattern, however, does not apply to some other incompatible elements such as B, Sb, and halogens. All Sr–Nd–Pb isotope compositions reflect a depleted source with Indian Ocean mantle domain characteristics. The Nd and Pb isotope ratios are radiogenic in the volcanic front, whereas Sr isotope ratios are less radiogenic. These Nd isotope ratios covary with incompatible element ratios such as Th/Nd and Nb/Zr, indicating involvement of a slab-derived sediment component by addition of melt or supercritical fluid capable of mobilizing these high field-strength elements and rare earth elements from the slab. Fluid mobile elements, such as Ba, are also elevated in all basalt suites, suggesting involvement of slab fluid derived from altered oceanic crust. The Kurile Arc lavas are thus affected both by slab sediment and altered basaltic crust components. This magma plumbing system has been continuously active from the Middle Miocene to the present.     

A comparison of continuous mass‐lumped finite elements with finite differences for 3‐D wave propagation

The finite‐difference method on rectangular meshes is widely used for time‐domain modelling of the wave equation. It is relatively easy to implement high‐order spatial discretization schemes and parallelization. Also, the method is computationally efficient. However, the use of finite elements on tetrahedral unstructured meshes is more accurate in complex geometries near sharp interfaces. We compared the standard eighth‐order finite‐difference method to fourth‐order continuous mass‐lumped finite elements in terms of accuracy and computational cost. The results show that, for simple models like a cube with constant density and velocity, the finite‐difference method outperforms the finite‐element method by at least an order of magnitude. Outside the application area of rectangular meshes, i.e., for a model with interior complexity and topography well described by tetrahedra, however, finite‐element methods are about two orders of magnitude faster than finite‐difference methods, for a given accuracy.     

Global warming mitigation by sulphur loading in the stratosphere: dependence of required emissions on allowable residual warming rate

An approach to mitigate global warming via sulphur loading in the stratosphere (geoengineering) is studied, employing a large ensemble of numerical experiments with the climate model of intermediate complexity IAP RAS CM. The model is forced by the historical+SRES A1B anthropogenic greenhouse gases+tropospheric sulphates scenario for 1860–2100 with additional sulphur emissions in the stratosphere in the twenty-first century. Different ensemble members are constructed by varying values of the parameters governing mass, horizontal distribution and radiative forcing of the stratospheric sulphates. It is obtained that, given a global loading of the sulphates in the stratosphere, among those studied in this paper latitudinal distributions of geoengineering aerosols, the most efficient one at the global basis is that peaked between 50°N and 70°N and with a somewhat smaller burden in the tropics. Uniform latitudinal distribution of stratospheric sulphates is a little less efficient. Sulphur emissions in the stratosphere required to stop the global temperature at the level corresponding to the mean value for 2000–2010 amount to more than 10 TgS/year in the year 2100. These emissions may be reduced if some warming is allowed to occur in the twenty-first century. For instance, if the global temperature trend S _g in every decade of this century is limited not to exceed 0.10 K/decade (0.15 K/decade), geoengineering emissions of 4–14 TgS/year (2–7 TgS/year) would be sufficient. Even if the global warming is stopped, temperature changes in different regions still occur with a magnitude up to 1 K. Their horizontal pattern depends on implied latitudinal distribution of stratospheric sulphates. In addition, for the stabilised global mean surface air temperature, global precipitation decreases by about 10%. If geoengineering emissions are stopped after several decades of implementation, their climatic effect is removed within a few decades. In this period, surface air temperature may grow with a rate of several Kelvins per decade. The results obtained with the IAP RAS CM are further interpreted employing a globally averaged energy–balance climate model. With the latter model, an analytical estimate for sulphate aerosol emissions in the stratosphere required climate mitigation is obtained. It is shown that effective vertical localisation of the imposed radiative forcing is important for geoengineering efficiency.     

Harmonic analysis of tides from standard observations of the sea level

Peculiarities of harmonic analysis of tides by using standard observation data of the sea level with a 6-hourly time step are considered. Masking frequencies causing difficulties when interpreting the results of the data analysis are presented. A method of control and editing of standard observations by using a non-recursive filter is proposed. The data interpolation approach based on direct and inverse Fourier transformation is considered.