Thermochronologic and geodynamic evolution of the Pontides: Sakarya terrane (Karakaya Complex) and Istanbul terrane, Turkey.

An integrated array of analytical methods -including clay mineralogy, vitrinite reflectance, Raman spectroscopy on carbonaceous material, and apatite fission-track analysis- was employed to constrain the thermal and thermochronological evolution of selected portions of the Pontides of northern Turkey. (1) A multimethod investigation was applied for the first time to characterise the thermal history of the Karakaya Complex, a Permo-Triassic subduction-accretion complex cropping out throughout the Sakarya Zone. The results indicate two different thermal regimes: the Lower Karakaya Complex (Nilüfer Unit) -mostly made of metabasite and marble- suffered peak temperatures of 300-500°C (greenschist facies); the Upper Karakaya Complex (Hodul and the Orhanlar Units) –mostly made of greywacke and arkose- yielded heterogeneous peak temperatures (125-376°C), possibly the result of different degree of involvement of the units in the complex dynamic processes of the accretionary wedge. Contrary to common belief, the results of this study indicate that the entire Karakaya Complex suffered metamorphic conditions. Moreover, a good degree of correlation among the results of these methods demonstrate that Raman spectroscopy on carbonaceous material can be applied successfully to temperature ranges of 200-330°C, thus extending the application of this method from higher grade metamorphic contexts to lower grade metamorphic conditions. (2) Apatite fission-track analysis was applied to the Sakarya and the İstanbul Zones in order to constrain the exhumation history and timing of amalgamation of these two exotic terranes. AFT ages from the İstanbul and Sakarya terranes recorded three distinct episodes of exhumation related to the complex tectonic evolution of the Pontides. (i) Paleocene - early Eocene ages (62.3-50.3 Ma) reflect the closure of the İzmir-Ankara ocean and the ensuing collision between the Sakarya terrane and the Anatolide-Tauride Block. (ii) Late Eocene - earliest Oligocene (43.5-32.3 Ma) ages reflect renewed tectonic activity along the İzmir-Ankara. (iii) Late Oligocene- Early Miocene ages reflect the onset and development of the northern Aegean extension. The consistency of AFT ages, both north and south of the tectonic contact between the İstanbul and Sakarya terranes, suggest that such terranes were amalgamated in pre-Cenozoic times. (3) Fission-track analysis was also applied to rock samples from the Marmara region, in an attempt to constrain the inception and development of the North Anatolian Fault system in the region. The results agree with those from the central Pontides. The youngest AFT ages (Late Oligocene - early Miocene) were recorded in the western portion of the Marmara Sea region and reflect the onset and development of northern Aegean extension. Fission-track data from the eastern Marmara Sea region indicate rapid Early Eocene exhumation induced by the development of the İzmir-Ankara orogenic wedge. Thermochronological data along the trace of the Ganos Fault –a segment of the North Anatolian Fault system- indicate the presence of a tectonic discontinuity active by Late Oligocene time, i.e. well before the arrival of the North Anatolian Fault system in the area. The integration of thermochronologic data with preexisting structural data point to the existence of a system of major E-W-trending structural discontinuities active at least from the Late Oligocene. In the Early Pliocene, inception of the present-day North Anatolian Fault system in the Marmara region occurred by reactivation of these older tectonic structures. 


Low-temperature thermochronological evolution of the Menderes and Alanya massifs (Turkey)

The application of two low-temperature thermochronometers [fission-track analysis and (U-Th)/He analyses, both on apatite] to various tectonostratigraphic units of the Menderes and Alanya Massifs of Turkey has provided significant new constraints to the understanding of their structural evolution. The Menderes Massif of western Anatolia is one of the largest metamorphic core complexes on Earth. The integration of the geochronometric dataset presented in this dissertation with preexisting ones from the literature delineates three groups of samples within the Menderes Massif. In the northern and southern region the massif experienced a Late Oligocene-Early Miocene tectonic denudation and surface uplift; whereas data from the central region are younger, with most ages ranging between the Middle-Late Miocene. The results of this study are consistent with the interpretation for a symmetric exhumation of the Menderes Massif. The Alanya Massif of SW Anatolia presents a typical nappe pile consisting of thrust sheets with contrasting metamorphic histories. Petrological and geochronological data clearly indicate that the tectonometamorphic evolution Alanya started from Late Cretaceous with the northward subduction of an ‘Alanya ocean’ under the Tauride plate. As an effect of the closure of the İzmir–Ankara–Erzincan ocean, northward backthrusting during the Paleocene-Early Eocene created the present stacking order. Apatite fission-track ages from this study range from 31.8 to 26.8 Ma (Late Rupelian-Early Chattian) and point to a previously unrecognized mid-Oligocene cooling/exhumation episode. (U-Th)/He analysis on zircon crystals obtained from the island of Cyprus evidentiate that the Late Cretaceous trondhjemites of the Troodos Massif not recorded a significant cooling event. Instead results for the Late Triassic turbiditic sandstones of the Vlambouros Formation show that the Mamonia mélange was never buried enough to reach the closure temperature of the ZHe radiometric system (ca. 200°C), thus retaining the Paleozoic signature of a previous sedimentary cycle.

Rapporti tra tettonica e clima durante l'esumazione nelle Alpi Centrali. Evidenze dalla termocronologia e dall'analisi strutturale lungo il traforo ferroviario del Sempione

(U-Th)/He and fission-track analyses of apatite along deep-seated tunnels crossing high-relief mountain ranges offer the opportunity to investigate climate and tectonic forcing on the topographic evolution. In this study, the thermochronologic analysis of a large set of samples collected in the Simplon railway tunnel (western-central Alps; Italy and Switzerland) and along its surface trace, coupled with kinematic and structural analysis of major fault zones intersecting the tunnel, constrains the phenomena controlling the topographic and structural evolution, during the latest stage of exhumation of the Simplon Massif, and the timing in which they operated. The study area is located at the western margin of the Lepontine metamorphic dome where a complex nappe-stack pertaining to the Penninic and Ultrahelvetic domains experienced a fast exhumation from the latest Oligocene onward. The exhumation was mainly accommodated by a west-dipping low-angle detachment (the Simplon Fault Zone) which is located just 8 km to the west of the tunnel. However, along the section itself several faults related to two principal phases both with important dip-slip kinematics have been detected. Cooling rates derived from our thermocronological data vary from about 10 °C/Ma at about 10 Ma to about 35 °C/Ma in the last 5 Ma. Such increase in the cooling rate corresponds to the most important climatic change recorded in the northern hemisphere in the last 10 Ma, i.e. the shift to wetter conditions at the end of the Messinian salinity crisis and the inception of glacial cycles in the northern hemisphere. In addition, (U-Th)/He and fission-track age patterns lack of important correlation with the topography suggesting that the present-day relief morphology is the result of recent erosional dynamics. More in details, the (U-Th)/He tunnel ages show an impressive uniformity at 2 Ma, whereas cooling rates calculated at 1 Ma increase towards the two major valleys. This indicates a focusing of erosive processes in the valleys which led to the shaping of present-day topography. Structural analysis documents the presence of two phases of brittle deformation postdating the metamorphic phases in the area. The first one is directly related to the last phase of activity along the Simplon Fault Zone and is characterized by extension towards SO and vertical shortening. The young one is characterized by extension towards NO and horizontal shortening in a along the NE-SO direction. Structures related to the first phase of brittle deformation generate important variations in the older ages' dataset, until 3 Ma, suggesting that tectonics controlled rocks exhumation up to that age. Structures related to the second phase generate some variations also in the younger age dataset, highlighting the activity of faults bordering the massif and suggesting a continuous activity also after 2 Ma. However, most of (U-Th)/He tunnel ages, varying slightly around 2 Ma, document that the Simplon area has experienced primarily erosional exhumation in this time span. In conclusion, all our data suggest that in the central Italian Alps the climatic signal gradually overrode the tectonic effects after about 5 Ma, as a consequence of the climatic instability started at end of Messinian salinity crisis and improved by the onset of glaciations in the northern hemisphere.

Erosione e sollevamento nell'arco calabro-peloritano

The Calabrian-Peloritani arc represents key site to unravel evolution of surface processes on top of subducting lithosphere. During the Pleistocene, in fact the arc uplifted at rate of the order of about 1mm/yr, forming high-standing low-relief upland (figure 2). Our study is focused on the relationship between tectonic and land evolution in the Sila Massif, Messina strait and Peloritani Mts. Landforms reflect a competition between tectonic, climatic, and surficial processes. Many landscape evolution models that explore feedbacks between these competing processes, given steady forcing, predict a state of erosional equilibrium, where the rates of river incision and hillslope erosion balance rock uplift. It has been suggested that this may be the final constructive stage of orogenic systems. Assumptions of steady erosion and incision are used in the interpretation of exhumation and uplift rates from different geologic data, and in the formulation of fluvial incision and hillslope evolution models. In the Sila massif we carried out cosmogenic isotopes analysis on 24 samples of modern fluvial sediments to constrain long-term (~103 yr) erosion rate averaged on the catchment area. 35 longitudinal rivers profiles have been analyzed to study the tectonic signal on the landscape evolution. The rivers analyzed exhibit a wide variety of profile forms, diverging from equilibrium state form. Generally the river profiles show at least 2 and often 3 distinct concave-up knickpoint-bounded segments, characterized by different value of concavity and steepness indices. River profiles suggest three main stages of incision. The values of ks and θ in the lower segments evidence a decrease in river incision, due probably to increasing uplift rate. The cosmogenic erosion rates pointed out that old landscape upland is eroding slowly at ~0.1 mm/yr. In the contrary, the flanks of the massif is eroding faster with value from 0.4 to 0.5 mm/yr due to river incision and hillslope processes. Cosmogenic erosion rates mach linearly with steepness indices and with average hillslope gradient. In the Messina area the long term erosion rate from low-T thermochronometry are of the same order than millennium scale cosmogenic erosion rate (1-2 mm/yr). In this part of the chain the fast erosion is active since several million years, probably controlled by extensional tectonic regime. In the Peloritani Mts apatite fission-track and (U-Th)/He thermochronometry are applied to constraint the thermal history of the basement rock. Apatite fission-track ages range between 29.0±5.5 and 5.5±0.9 Ma while apatite (U-Th)/He ages vary from 19.4 to 1.0 Ma. Most of the AFT ages are younger than the overlying terrigenous sequence that in turn postdates the main orogenic phase. Through the coupling of the thermal modelling with the stratigraphic record, a Middle Miocene thermal event due to tectonic burial is unravel. This event affected a inner-intermediate portion of the Peloritani belt confined by young AFT data (<15 Ma) distribution. We interpret this thermal event as due to an out-of–sequence thrusting occurring in the inner portion of the belt. Young (U-Th)/He ages (c. 5 Ma) record a final exhumation stage with increasing rates of denudation since the Pliocene times due to postorogenic extensional tectonics and regional uplift. In the final chapter we change the spatial scale to insert digital topography analysis and field data within a geodynamic model that can explain surface evidence produced by subduction process.

Analysis of optimal control problems for the incompressible MHD equations and implementation in a finite element multiphysics code

This thesis deals with the study of optimal control problems for the incompressible Magnetohydrodynamics (MHD) equations. Particular attention to these problems arises from several applications in science and engineering, such as fission nuclear reactors with liquid metal coolant and aluminum casting in metallurgy. In such applications it is of great interest to achieve the control on the fluid state variables through the action of the magnetic Lorentz force. In this thesis we investigate a class of boundary optimal control problems, in which the flow is controlled through the boundary conditions of the magnetic field. Due to their complexity, these problems present various challenges in the definition of an adequate solution approach, both from a theoretical and from a computational point of view. In this thesis we propose a new boundary control approach, based on lifting functions of the boundary conditions, which yields both theoretical and numerical advantages. With the introduction of lifting functions, boundary control problems can be formulated as extended distributed problems. We consider a systematic mathematical formulation of these problems in terms of the minimization of a cost functional constrained by the MHD equations. The existence of a solution to the flow equations and to the optimal control problem are shown. The Lagrange multiplier technique is used to derive an optimality system from which candidate solutions for the control problem can be obtained. In order to achieve the numerical solution of this system, a finite element approximation is considered for the discretization together with an appropriate gradient-type algorithm. A finite element object-oriented library has been developed to obtain a parallel and multigrid computational implementation of the optimality system based on a multiphysics approach. Numerical results of two- and three-dimensional computations show that a possible minimum for the control problem can be computed in a robust and accurate manner.

Gait analysis using a single wearable inertial measurement unit

Procedures for quantitative walking analysis include the assessment of body segment movements within defined gait cycles. Recently, methods to track human body motion using inertial measurement units have been suggested. It is not known if these techniques can be readily transferred to clinical measurement situations. This work investigates the aspects necessary for one inertial measurement unit mounted on the lower back to track orientation, and determine spatio-temporal features of gait outside the confines of a conventional gait laboratory. Apparent limitations of different inertial sensors can be overcome by fusing data using methods such as a Kalman filter. The benefits of optimizing such a filter for the type of motion are unknown. 3D accelerations and 3D angular velocities were collected for 18 healthy subjects while treadmill walking. Optimization of Kalman filter parameters improved pitch and roll angle estimates when compared to angles derived using stereophotogrammetry. A Weighted Fourier Linear Combiner method for estimating 3D orientation angles by constructing an analytical representation of angular velocities and allowing drift free integration is also presented. When tested this method provided accurate estimates of 3D orientation when compared to stereophotogrammetry. Methods to determine spatio-temporal features from lower trunk accelerations generally require knowledge of sensor alignment. A method was developed to estimate the instants of initial and final ground contact from accelerations measured by a waist mounted inertial device without rigorous alignment. A continuous wavelet transform method was used to filter and differentiate the signal and derive estimates of initial and final contact times. The technique was tested with data recorded for both healthy and pathologic (hemiplegia and Parkinson’s disease) subjects and validated using an instrumented mat. The results show that a single inertial measurement unit can assist whole body gait assessment however further investigation is required to understand altered gait timing in some pathological subjects.

Power Transient Analysis of Experimental Devices for Jules Horowitz Material Testing Reactor (JHR)

The objective of this thesis is the power transient analysis concerning experimental devices placed within the reflector of Jules Horowitz Reactor (JHR). Since JHR material testing facility is designed to achieve 100 MW core thermal power, a large reflector hosts fissile material samples that are irradiated up to total relevant power of 3 MW. MADISON devices are expected to attain 130 kW, conversely ADELINE nominal power is of some 60 kW. In addition, MOLFI test samples are envisaged to reach 360 kW for what concerns LEU configuration and up to 650 kW according to HEU frame. Safety issues concern shutdown transients and need particular verifications about thermal power decreasing of these fissile samples with respect to core kinetics, as far as single device reactivity determination is concerned. Calculation model is conceived and applied in order to properly account for different nuclear heating processes and relative time-dependent features of device transients. An innovative methodology is carried out since flux shape modification during control rod insertions is investigated regarding the impact on device power through core-reflector coupling coefficients. In fact, previous methods considering only nominal core-reflector parameters are then improved. Moreover, delayed emissions effect is evaluated about spatial impact on devices of a diffuse in-core delayed neutron source. Delayed gammas transport related to fission products concentration is taken into account through evolution calculations of different fuel compositions in equilibrium cycle. Provided accurate device reactivity control, power transients are then computed for every sample according to envisaged shutdown procedures. Results obtained in this study are aimed at design feedback and reactor management optimization by JHR project team. Moreover, Safety Report is intended to utilize present analysis for improved device characterization.

Dynamic analysis of the motorcycle chattering behaviour by means of symbolic multibody modelling

Aim of this research is the development and validation of a comprehensive multibody motorcycle model featuring rigid-ring tires, taking into account both slope and roughness of road surfaces. A novel parametrization for the general kinematics of the motorcycle is proposed, using a mixed reference-point and relative-coordinates approach. The resulting description, developed in terms of dependent coordinates, makes it possible to efficiently include rigid-ring kinematics as well as road elevation and slope. The equations of motion for the multibody system are derived symbolically and the constraint equations arising from the dependent-coordinate formulation are handled using a projection technique. Therefore the resulting system of equations can be integrated in time domain using a standard ODE algorithm. The model is validated with respect to maneuvers experimentally measured on the race track, showing consistent results and excellent computational efficiency. More in detail, it is also capable of reproducing the chatter vibration of racing motorcycles. The chatter phenomenon, appearing during high speed cornering maneuvers, consists of a self-excited vertical oscillation of both the front and rear unsprung masses in the range of frequency between 17 and 22 Hz. A critical maneuver is numerically simulated, and a self-excited vibration appears, consistent with the experimentally measured chatter vibration. Finally, the driving mechanism for the self-excitation is highlighted and a physical interpretation is proposed.