Models for the study of neuronal activity during cognitive tasks

Assessment of brain connectivity among different brain areas during cognitive or motor tasks is a crucial problem in neuroscience today. Aim of this research study is to use neural mass models to assess the effect of various connectivity patterns in cortical EEG power spectral density (PSD), and investigate the possibility to derive connectivity circuits from EEG data. To this end, two different models have been built. In the first model an individual region of interest (ROI) has been built as the parallel arrangement of three populations, each one exhibiting a unimodal spectrum, at low, medium or high frequency. Connectivity among ROIs includes three parameters, which specify the strength of connection in the different frequency bands. Subsequent studies demonstrated that a single population can exhibit many different simultaneous rhythms, provided that some of these come from external sources (for instance, from remote regions). For this reason in the second model an individual ROI is simulated only with a single population. Both models have been validated by comparing the simulated power spectral density with that computed in some cortical regions during cognitive and motor tasks. Another research study is focused on multisensory integration of tactile and visual stimuli in the representation of the near space around the body (peripersonal space). This work describes an original neural network to simulate representation of the peripersonal space around the hands, in basal conditions and after training with a tool used to reach the far space. The model is composed of three areas for each hand, two unimodal areas (visual and tactile) connected to a third bimodal area (visual-tactile), which is activated only when a stimulus falls within the peripersonal space. Results show that the peripersonal space, which includes just a small visual space around the hand in normal conditions, becomes elongated in the direction of the tool after training, thanks to a reinforcement of synapses.

Shallow and deep deformation in northern Apennines region using seismological data

For its particular position and the complex geological history, the Northern Apennines has been considered as a natural laboratory to apply several kinds of investigations. By the way, it is complicated to joint all the knowledge about the Northern Apennines in a unique picture that explains the structural and geological emplacement that produced it. The main goal of this thesis is to put together all information on the deformation - in the crust and at depth - of this region and to describe a geodynamical model that takes account of it. To do so, we have analyzed the pattern of deformation in the crust and in the mantle. In both cases the deformation has been studied using always information recovered from earthquakes, although using different techniques. In particular the shallower deformation has been studied using seismic moment tensors information. For our purpose we used the methods described in Arvidsson and Ekstrom (1998) that allowing the use in the inversion of surface waves [and not only of the body waves as the Centroid Moment Tensor (Dziewonski et al., 1981) one] allow to determine seismic source parameters for earthquakes with magnitude as small as 4.0. We applied this tool in the Northern Apennines and through this activity we have built up the Italian CMT dataset (Pondrelli et al., 2006) and the pattern of seismic deformation using the Kostrov (1974) method on a regular grid of 0.25 degree cells. We obtained a map of lateral variations of the pattern of seismic deformation on different layers of depth, taking into account the fact that shallow earthquakes (within 15 km of depth) in the region occur everywhere while most of events with a deeper hypocenter (15-40 km) occur only in the outer part of the belt, on the Adriatic side. For the analysis of the deep deformation, i.e. that occurred in the mantle, we used the anisotropy information characterizing the structure below the Northern Apennines. The anisotropy is an earth properties that in the crust is due to the presence of aligned fluid filled cracks or alternating isotropic layers with different elastic properties while in the mantle the most important cause of seismic anisotropy is the lattice preferred orientation (LPO) of the mantle minerals as the olivine. This last is a highly anisotropic mineral and tends to align its fast crystallographic axes (a-axis) parallel to the astenospheric flow as a response to finite strain induced by geodynamic processes. The seismic anisotropy pattern of a region is measured utilizing the shear wave splitting phenomenon (that is the seismological analogue to optical birefringence). Here, to do so, we apply on teleseismic earthquakes recorded on stations located in the study region, the Sileny and Plomerova (1996) approach. The results are analyzed on the basis of their lateral and vertical variations to better define the earth structure beneath Northern Apennines. We find different anisotropic domains, a Tuscany and an Adria one, with a pattern of seismic anisotropy which laterally varies in a similar way respect to the seismic deformation. Moreover, beneath the Adriatic region the distribution of the splitting parameters is so complex to request an appropriate analysis. Therefore we applied on our data the code of Menke and Levin (2003) which allows to look for different models of structures with multilayer anisotropy. We obtained that the structure beneath the Po Plain is probably even more complicated than expected. On the basis of the results obtained for this thesis, added with those from previous works, we suggest that slab roll-back, which created the Apennines and opened the Tyrrhenian Sea, evolved in the north boundary of Northern Apennines in a different way from its southern part. In particular, the trench retreat developed primarily south of our study region, with an eastward roll-back. In the northern portion of the orogen, after a first stage during which the retreat was perpendicular to the trench, it became oblique with respect to the structure.

Substructuring approache in state space models for dynamic system parameters identification

In the recent decade, the request for structural health monitoring expertise increased exponentially in the United States. The aging issues that most of the transportation structures are experiencing can put in serious jeopardy the economic system of a region as well as of a country. At the same time, the monitoring of structures is a central topic of discussion in Europe, where the preservation of historical buildings has been addressed over the last four centuries. More recently, various concerns arose about security performance of civil structures after tragic events such the 9/11 or the 2011 Japan earthquake: engineers looks for a design able to resist exceptional loadings due to earthquakes, hurricanes and terrorist attacks. After events of such a kind, the assessment of the remaining life of the structure is at least as important as the initial performance design. Consequently, it appears very clear that the introduction of reliable and accessible damage assessment techniques is crucial for the localization of issues and for a correct and immediate rehabilitation. The System Identification is a branch of the more general Control Theory. In Civil Engineering, this field addresses the techniques needed to find mechanical characteristics as the stiffness or the mass starting from the signals captured by sensors. The objective of the Dynamic Structural Identification (DSI) is to define, starting from experimental measurements, the modal fundamental parameters of a generic structure in order to characterize, via a mathematical model, the dynamic behavior. The knowledge of these parameters is helpful in the Model Updating procedure, that permits to define corrected theoretical models through experimental validation. The main aim of this technique is to minimize the differences between the theoretical model results and in situ measurements of dynamic data. Therefore, the new model becomes a very effective control practice when it comes to rehabilitation of structures or damage assessment. The instrumentation of a whole structure is an unfeasible procedure sometimes because of the high cost involved or, sometimes, because it’s not possible to physically reach each point of the structure. Therefore, numerous scholars have been trying to address this problem. In general two are the main involved methods. Since the limited number of sensors, in a first case, it’s possible to gather time histories only for some locations, then to move the instruments to another location and replay the procedure. Otherwise, if the number of sensors is enough and the structure does not present a complicate geometry, it’s usually sufficient to detect only the principal first modes. This two problems are well presented in the works of Balsamo [1] for the application to a simple system and Jun [2] for the analysis of system with a limited number of sensors. Once the system identification has been carried, it is possible to access the actual system characteristics. A frequent practice is to create an updated FEM model and assess whether the structure fulfills or not the requested functions. Once again the objective of this work is to present a general methodology to analyze big structure using a limited number of instrumentation and at the same time, obtaining the most information about an identified structure without recalling methodologies of difficult interpretation. A general framework of the state space identification procedure via OKID/ERA algorithm is developed and implemented in Matlab. Then, some simple examples are proposed to highlight the principal characteristics and advantage of this methodology. A new algebraic manipulation for a prolific use of substructuring results is developed and implemented.

Entwicklung von computergestützten Hangstabilitätsmodellen zur Erstellung einer Gefahrenhinweiskarte für die Region Rheinhessen

Das Ziel der Arbeit war die Entwicklung computergestützter Methoden zur Erstellung einer Gefahrenhinweiskarte für die Region Rheinhessen, zur Minimierung der Hangrutschungsgefährdung. Dazu wurde mit Hilfe zweier statistischer Verfahren (Diskriminanzanalyse, Lo­gistische Regression) und einer Methode aus dem Bereich der Künstlichen Intelligenz (Fuzzy Logik) versucht, die potentielle Gefährdung auch solcher Hänge zu klassifizieren, die bis heute noch nicht durch Massenbewegungen aufgefallen sind. Da ingenieurgeologische und geotechnische Hanguntersuchungen aus Zeit und Kostengründen im regionalen Maßstab nicht möglich sind, wurde auf punktuell vorhandene Datenbestände zu einzelnen Rut­schungen des Winters 1981/82, die in einer Rutschungsdatenbank zu­sammengefaßt sind, zurückgegriffen, wobei die daraus gewonnenen Erkenntnisse über Prozeßmechanismen und auslösende Faktoren genutzt und in das jeweilige Modell integriert wurden. Flächenhafte Daten (Lithologie, Hangneigung, Landnutzung, etc.), die für die Berechnung der Hangstabilität notwendig sind, wurden durch Fernerkundungsmethoden, dem Digitalisieren von Karten und der Auswertung von Digitalen Geländemodellen (Reliefanalyse) gewonnen. Für eine weiterführende Untersuchung von einzelnen, als rutschgefährdet klassifi­zierten Bereichen der Gefahrenhinweiskarte, wurde am Beispiel eines Testgebietes, eine auf dem infinite-slope-stability Modell aufbauende Me­thode untersucht, die im Maßstabsbereich von Grundkarten (1:5000) auch geotechnische und hydrogeologische Parameter berücksichtigt und damit eine genauere, der jeweiligen klimatischen Situation angepaßte, Gefahrenabschätzung ermöglicht.

Conditional and transgenic mouse models as tools to study the role of TGFbeta, BMP signaling in skeletal development

Die TGFbeta/BMP Signaltransduktionskaskade ist wichtig für viele Entwicklungsprozesse fast aller embryonaler sowie extraembryonaler Gewebe und sie ist ebenso essentiell bei der Aufrechterhaltung der Homöostase im adulten Organismus. In vielen Mausmodellen und Zellkulturversuchen wurde gezeigt, dass Liganden dieses Signalweges in verschiedene Stadien der Knorpel- und Knochenentwicklung involviert sind. BMPs sind beispielsweise maßgeblich an der frühen Kondensation und Bildung des Knorpels und später an Proliferation und Hypertrophie der Chondrozyten beteiligt. BMPs können ektopisch Knochenbildung auslösen und das Expressionsmuster der Liganden und spezifischen Rezeptoren in der Wachstumsfuge lässt auf eine wichtige Rolle der BMPs in der Wachstumsfuge schließen. Der gezielte knock out der BMP-Rezeptoren Bmpr1a und Bmpr1b in proliferierenden Chondrozyten führt zur Ausbildung einer generellen Chondrodysplasie. Smad1, Smad5 und Smad8 sind die Mediatoren der BMP-Signalkaskade. Im Rahmen der vorliegenden Arbeit sollte die Rolle und Funktion der Smad1- und Smad5-Proteine in der Wachstumsfuge untersucht werden. Hierzu wurden konditionale Smad1-knock out-Mäuse mit einer transgenen Mauslinie gekreuzt, die die Cre-Rekombinase spezifisch in proliferierenden Chondrozyten exprimiert. Diese Mäuse wurden mit und ohne heterozygotem Smad5-Hintergrund charakterisiert. Bei einem knock out von Smad1 allein konnte ein leichte Verkürzung der Wachstumsfuge beobachtet werden, wobei prähypertrophe und hypertrophe Zone gleichermaßen betroffen waren. Dieser Phänotyp war verstärkt in Mäusen mit zusätzlichem heterozygotem Smad5-Hintergrund. Eine Verringerung der Proliferationsrate konnte zusammen mit einer verminderten Ihh-Expression nachgewiesen werden. Zusätzlich konnte anhand von Röntgenaufnahmen eine Dysorganisation der nasalen Region und ein fehlendes nasales Septum beobachtet werden. Produktion und Mineralisation der extrazellulären Matrix waren nicht beeinträchtigt. Um die Rolle der BMP- und TGFbeta-Signalkaskaden während der endochondralen Ossifikation zu vergleichen, wurden transgene Mäuse generiert, in denen die TGFbeta-Signalkaskade spezifisch in proliferierenden Chondrozyten gestört war. Zwei Mauslinien, die ähnliche Phänotypen zeigten, wurden untersucht. Esl1 ist ein TGFbeta-bindendes Protein, von dem man annimmt, dass es die TGFbeta-Signalkaskade inhibieren kann. Esl1-knock out-Mäuse sind kleiner als Wildtypmäuse und die Überexpression von Esl1 in proliferierenden Chondrozyten führt zu einer Verlängerung der Wachstumsfuge und einer verstärkten Proliferationsrate. Knorpelmarker, wie Col2a1 und Sox9 sind in diesen Mäusen herunterreguliert, während Col10a1 und Ihh als Marker für die hypertrophe und prähypertrophe Zone herunterreguliert waren. Dies führt zu der Annahme, dass mehr Zellen in die terminale Differenzierung eintreten. Bei transgenen Mäusen, in denen ein dominant-negativer (dn) TGFbeta-Rezeptor in proliferierenden Chondrozyten überexprimiert wurde, konnte eine verlängerte prähypertrophe Zone, eine erhöhte Ihh-Expression, sowie eine verstärkte Proliferationsrate beobachtet werden. Zusätzlich konnte in homozygoten Tieren ein craniofacialer Phänotyp beschrieben werden, der zu Problemen bei der Nahrungsaufnahme und damit zu einer starken Wachstumsbeeinträchtigung führte. Die BMP- und TGFbeta-Signalkaskaden haben möglicherweise antagonistische Effekte in der Wachstumsfuge. Während der Ausfall von BMP in proliferierenden Chondrozyten aufgrund einer gesunkenen Proliferationsrate zu einer Verkürzung der Wachstumsfuge führte, kann man in Mäusen mit einer Störung der TGFbeta-Signalkaskade eine verstärkte Proliferation in einer daher verlängerten Wachstumsfuge beobachten. Ein weiteres Ziel dieser Arbeit war die Generation einer transgenen Mauslinie, die die Cre-Rekombinase spezifisch in hypertrophen Chondrozyten exprimiert. Promoterstudien mit transgenen Mäusen weisen darauf hin, dass ein putatives AP1-Element, etwa 4 kb vor dem ersten Exon des Col10a1 gelegen, wichtig für die spezifische Expression in hypertrophen Chondrozyten ist. Ein Konstrukt, dass vier Kopien dieses Elements und den basalen Promoter enthält, wurde benutzt, um die Cre-Rekombinase spezifisch zu exprimieren. Diese Mauslinie befindet sich in der Testphase und erste Daten deuten auf eine spezifische Expression der Cre-Rekombinase in hypertrophen Chondrozyten hin.

Lattice preferred orientation in omphacite: Examples from the Tauern Window, Austria and the Western Gneiss Region, Norway

It is lively debated how eclogites find their way from deep to mid-crustal levels during exhumation. Different exhumation models for high-pressure and ultrahigh-pressure rocks were suggested in previous studies, based mainly on field observations and less on microstructural studies on the exhumed rocks. The development and improvement of electron microscopy techniques allows it, to focus interest on direct investigations of microstructures and crystallographic properties in eclogites. In this case, it is of importance to study the applicability of crystallographic measurements on eclogites for exhumation processes and to unravel which processes affect eclogite textures. Previous studies suggested a strong relationship between deformation and lattice preferred orientation (LPO) in omphacite but it is still unclear if the deformation is related to the exhumation of eclogites. This study is focused on the questions which processes affect omphacite LPO and if textural investigations of omphacite are applicable for studying eclogite exhumation. Therefore, eclogites from two examples in the Alps and in the Caledonides were collected systematically and investigated with respect to omphacite LPO by using the electron backscattered diffraction (EBSD) technique. Omphacite textures of the Tauern Window (Austria) and the Western Gneiss Region (Norway) were studied to compare lattice preferred orientation with field observations and suggested exhumation models from previous studies. The interpretation of omphacite textures, regarding the deformation regime is mainly based on numerical simulations in previous studies. Omphacite LPO patterns of the Eclogite Zone are clearly independent from any kind of exhumation process. The textures were generated during omphacite growth on the prograde path of eclogite development until metamorphic peak conditions. Field observations in the Eclogite Zone show that kinematics in garnet mica schist, surrounding the eclogites, strongly indicate an extrusion wedge geometry. Stretching lineations show top-N thrusting at the base and a top-S normal faulting with a sinistral shear component at the top of the Eclogite Zone. The different shear sense on both sides of the unit does not affect the omphacite textures in any way. The omphacite lattice preferred orientation patterns of the Western Gneiss Region can not be connected with any exhumation model. The textures were probably generated during the metamorphic peak and reflect the change from subduction to exhumation. Eclogite Zone and Western Gneiss Region differ significantly in size and especially in metamorphic conditions. While the Eclogite Zone is characterized by constant P-T conditions (600-650°C, 20-25 kbar), the Western Gneiss Region contains a wide P-T range from high- to ultrahigh pressure conditions (400-800°C, 20-35 kbar). In contrast to this, the omphacite textures of both units are very similar. This means that omphacite LPO is independent from P-T conditions and therefore from burial depth. Further, in both units, omphacite LPO is independent from grain and subgrain size as well as from any shape preferred orientation (SPO) on grain and subgrain scale. Overall, omphacite lattice preferred orientation are generated on the prograde part of omphacite development. Therefore, textural investigations on omphacite LPO are not applicable to study eclogite exhumation.

Determinants of cesarean delivery, a population-based study in the Emilia Romagna Region

Cesarean Delivery (CD) rates are rising in many parts of the world. In order to define strategies to reduce them, it is important to explore the role of clinical and organizational factors. This thesis has the objective to describe the contemporary CD practice and study clinical and organizational variables as determinants of CD in all women who gave birth between 2005 and June 2010 in the Emilia Romagna region (Italy). All hospital discharge abstracts of women who delivered between 2005 and mid 2010 in the region were selected and linked with birth certificates. In addition to descriptive statistics, in order to study the role of clinical and organizational variables (teaching or non-teaching hospital, birth volumes, time and day of delivery) multilevel Poisson regression models and a classification tree were used. A substantial inter-hospital variability in CD rate was found, and this was only partially explained by the considered variables. The most important risk factors of CD were: previous CD (RR 4,95; 95%CI: 4,85-5,05), cord prolapse (RR 3,51; 95% CI:2,96-4,16), and malposition/malpresentation (RR 2,72; 95%CI: 2,66-2,77). Delivery between 7 pm and 7 am and during non working days protect against CD in all subgroups including those with a small number of elective CDs while delivery at a teaching hospital and birth volumes were not statistically significant risk factors. The classification tree shows that previous CD and malposition/malpresentation are the most important variables discriminating between high and low risk of CD. These results indicate that other not considered factors might explain CD variability and do not provide clear evidence that small hospitals have a poor performance in terms of CD rate. Some strategies to reduce CD could be found by focusing on the differences in delivery practice between day and night and between working and no-working day deliveries.

Kinematic models of interseismic deformation from inversion of GPS and InSAR measurements to estimate fault parameters and coupling degree

We have used kinematic models in two Italian regions to reproduce surface interseismic velocities obtained from InSAR and GPS measurements. We have considered a Block modeling, BM, approach to evaluate which fault system is actively accommodating the occurring deformation in both considered areas. We have performed a study for the Umbria-Marche Apennines, obtaining that the tectonic extension observed by GPS measurements is explained by the active contribution of at least two fault systems, one of which is the Alto Tiberina fault, ATF. We have estimated also the interseismic coupling distribution for the ATF using a 3D surface and the result shows an interesting correlation between the microseismicity and the uncoupled fault portions. The second area analyzed concerns the Gargano promontory for which we have used jointly the available InSAR and GPS velocities. Firstly we have attached the two datasets to the same terrestrial reference frame and then using a simple dislocation approach, we have estimated the best fault parameters reproducing the available data, providing a solution corresponding to the Mattinata fault. Subsequently we have considered within a BM analysis both GPS and InSAR datasets in order to evaluate if the Mattinata fault may accommodate the deformation occurring in the central Adriatic due to the relative motion between the North-Adriatic and South-Adriatic plates. We obtain that the deformation occurring in that region should be accommodated by more that one fault system, that is however difficult to detect since the poor coverage of geodetic measurement offshore of the Gargano promontory. Finally we have performed also the estimate of the interseismic coupling distribution for the Mattinata fault, obtaining a shallow coupling pattern. Both of coupling distributions found using the BM approach have been tested by means of resolution checkerboard tests and they demonstrate that the coupling patterns depend on the geodetic data positions.

Numerical and semi-analytical models of sliding masses: application to the 1783 Scilla tsunamigenic landslide

The Scilla rock avalanche occurred on 6 February 1783 along the coast of the Calabria region (southern Italy), close to the Messina Strait. It was triggered by a mainshock of the Terremoto delle Calabrie seismic sequence, and it induced a tsunami wave responsible for more than 1500 casualties along the neighboring Marina Grande beach. The main goal of this work is the application of semi-analtycal and numerical models to simulate this event. The first one is a MATLAB code expressly created for this work that solves the equations of motion for sliding particles on a two-dimensional surface through a fourth-order Runge-Kutta method. The second one is a code developed by the Tsunami Research Team of the Department of Physics and Astronomy (DIFA) of the Bologna University that describes a slide as a chain of blocks able to interact while sliding down over a slope and adopts a Lagrangian point of view. A wide description of landslide phenomena and in particular of landslides induced by earthquakes and with tsunamigenic potential is proposed in the first part of the work. Subsequently, the physical and mathematical background is presented; in particular, a detailed study on derivatives discratization is provided. Later on, a description of the dynamics of a point-mass sliding on a surface is proposed together with several applications of numerical and analytical models over ideal topographies. In the last part, the dynamics of points sliding on a surface and interacting with each other is proposed. Similarly, different application on an ideal topography are shown. Finally, the applications on the 1783 Scilla event are shown and discussed.

PM(10) remote sensing from geostationary SEVIRI and polar-orbiting MODIS sensors over the complex terrain of the European Alpine region

The subject of this study is to investigate the capability of spaceborne remote sensing data to predict ground concentrations of PM10 over the European Alpine region using satellite derived Aerosol Optical Depth (AOD) from the geostationary Spinning Enhanced Visible and InfraRed Imager (SEVIRI) and the polar-orbiting MODerate resolution Imaging Spectroradiometer (MODIS). The spatial and temporal resolutions of these aerosol products (10 km and 2 measurements per day for MODIS, ∼ 25 km and observation intervals of 15 min for SEVIRI) permit an evaluation of PM estimation from space at different spatial and temporal scales. Different empirical linear relationships between coincident AOD and PM10 observations are evaluated at 13 ground-based PM measurement sites, with the assumption that aerosols are vertically homogeneously distributed below the planetary Boundary Layer Height (BLH). The BLH and Relative Humidity (RH) variability are assessed, as well as their impact on the parameterization. The BLH has a strong influence on the correlation of daily and hourly time series, whilst RH effects are less clear and smaller in magnitude. Despite its lower spatial resolution and AOD accuracy, SEVIRI shows higher correlations than MODIS (rSEV∼ 0.7, rMOD∼ 0.6) with regard to daily averaged PM10. Advantages from MODIS arise only at hourly time scales in mountainous locations but lower correlations were found for both sensors at this time scale (r∼ 0.45). Moreover, the fraction of days in 2008 with at least one satellite observation was 27% for SEVIRI and 17% for MODIS. These results suggest that the frequency of observations plays an important role in PM monitoring, while higher spatial resolution does not generally improve the PM estimation. Ground-based Sun Photometer (SP) measurements are used to validate the satellite-based AOD in the study region and to discuss the impact of aerosols' micro-physical properties in the empirical models. A lower error limit of 30 to 60% in the PM10 assessment from space is estimated in the study area as a result of AOD uncertainties, variability of aerosols properties and the heterogeneity of ground measurement sites. It is concluded that SEVIRI has a similar capacity to map PM as sensors on board polar-orbiting platforms, with the advantage of a higher number of observations. However, the accuracy represents a serious limitation to the applicability of satellites for ground PM mapping, especially in mountainous areas.

Improved water vapour spectroscopy in the 4174-4300 cm-1 region and its impact on SCIAMACHY HDO/H2O measurements

The relative abundance of the heavy water isotopologue HDO provides a deeper insight into the atmospheric hydrological cycle. The SCanning Imaging Absorption spectroMeter for Atmospheric CartograpHY (SCIAMACHY) allows for global retrievals of the ratio HDO/H2O in the 2.3 micron wavelength range. However, the spectroscopy of water lines in this region remains a large source of uncertainty for these retrievals. We therefore evaluate and improve the water spectroscopy in the range 4174–4300 cm−1 and test if this reduces systematic uncertainties in the SCIAMACHY retrievals of HDO/H2O. We use a laboratory spectrum of water vapour to fit line intensity, air broadening and wavelength shift parameters. The improved spectroscopy is tested on a series of ground-based high resolution FTS spectra as well as on SCIAMACHY retrievals of H2O and the ratio HDO/H2O. We find that the improved spectroscopy leads to lower residuals in the FTS spectra compared to HITRAN 2008 and Jenouvrier et al. (2007) spectroscopy, and the retrievals become more robust against changes in the retrieval window. For both the FTS and SCIAMACHY measurements, the retrieved total H2O columns decrease by 2–4% and we find a negative shift of the HDO/H2O ratio, which for SCIAMACHY is partly compensated by changes in the retrieval setup and calibration software. The updated SCIAMACHY HDO/H2O product shows somewhat steeper latitudinal and temporal gradients and a steeper Rayleigh distillation curve, strengthening previous conclusions that current isotope-enabled general circulation models underestimate the variability in the near-surface HDO/H2O ratio.

The Optimal Confidence Region for a Random Parameter

Under a two-level hierarchical model, suppose that the distribution of the random parameter is known or can be estimated well. Data are generated via a fixed, but unobservable realization of this parameter. In this paper, we derive the smallest confidence region of the random parameter under a joint Bayesian/frequentist paradigm. On average this optimal region can be much smaller than the corresponding Bayesian highest posterior density region. The new estimation procedure is appealing when one deals with data generated under a highly parallel structure, for example, data from a trial with a large number of clinical centers involved or genome-wide gene-expession data for estimating individual gene- or center-specific parameters simultaneously. The new proposal is illustrated with a typical microarray data set and its performance is examined via a small simulation study.

Shear band evolution in Zr/Hf-based bulk metallic glasses under static and dynamic indentations

Bulk metallic glasses (BMGs) exhibit superior mechanical properties as compared with other conventional materials and have been proposed for numerous engineering and technological applications. Zr/Hf-based BMGs or tungsten reinforced BMG composites are considered as a potential replacement for depleted uranium armor-piercing projectiles because of their ability to form localized shear bands during impact, which has been known to be the dominant plastic deformation mechanism in BMGs. However, in conventional tensile, compressive and bending tests, limited ductility has been observed because of fracture initiation immediately following the shear band formation. To fully investigate shear band characteristics, indentation tests that can confine the deformation in a limited region have been pursued. In this thesis, a detailed investigation of thermal stability and mechanical deformation behavior of Zr/Hf-based BMGs is conducted. First, systematic studies had been implemented to understand the influence of relative compositions of Zr and Hf on thermal stability and mechanical property evolution. Second, shear band evolution under indentations were investigated experimentally and theoretically. Three kinds of indentation studies were conducted on BMGs in the current study. (a) Nano-indentation to determine the mechanical properties as a function of Hf/Zr content. (b) Static Vickers indentation on bonded split specimens to investigate the shear band evolution characteristics beneath the indention. (c) Dynamic Vickers indentation on bonded split specimens to investigate the influence of strain rate. It was found in the present work that gradually replacing Zr by Hf remarkably increases the density and improves the mechanical properties. However, a slight decrease in glass forming ability with increasing Hf content has also been identified through thermodynamic analysis although all the materials in the current study were still found to be amorphous. Many indentation studies have revealed only a few shear bands surrounding the indent on the top surface of the specimen. This small number of shear bands cannot account for the large plastic deformation beneath the indentations. Therefore, a bonded interface technique has been used to observe the slip-steps due to shear band evolution. Vickers indentations were performed along the interface of the bonded split specimen at increasing loads. At small indentation loads, the plastic deformation was primarily accommodated by semi-circular primary shear bands surrounding the indentation. At higher loads, secondary and tertiary shear bands were formed inside this plastic zone. A modified expanding cavity model was then used to predict the plastic zone size characterized by the shear bands and to identify the stress components responsible for the evolution of the various types of shear bands. The applicability of various hardness—yield-strength ( H −σγ ) relationships currently available in the literature for bulk metallic glasses (BMGs) is also investigated. Experimental data generated on ZrHf-based BMGs in the current study and those available elsewhere on other BMG compositions were used to validate the models. A modified expanding-cavity model, employed in earlier work, was extended to propose a new H −σγ relationship. Unlike previous models, the proposed model takes into account not only the indenter geometry and the material properties, but also the pressure sensitivity index of the BMGs. The influence of various model parameters is systematically analyzed. It is shown that there is a good correlation between the model predictions and the experimental data for a wide range of BMG compositions. Under dynamic Vickers indentation, a decrease in indentation hardness at high loading rate was observed compared to static indentation hardness. It was observed that at equivalent loads, dynamic indentations produced more severe deformation features on the loading surface than static indentations. Different from static indentation, two sets of widely spaced semi-circular shear bands with two different curvatures were observed. The observed shear band pattern and the strain rate softening in indentation hardness were rationalized based on the variations in the normal stress on the slip plane, the strain rate of shear and the temperature rise associated with the indentation deformation. Finally, a coupled thermo-mechanical model is proposed that utilizes a momentum diffusion mechanism for the growth and evolution of the final spacing of shear bands. The influence of strain rate, confinement pressure and critical shear displacement on the shear band spacing, temperature rise within the shear band, and the associated variation in flow stress have been captured and analyzed. Consistent with the known pressure sensitive behavior of BMGs, the current model clearly captures the influence of the normal stress in the formation of shear bands. The normal stress not only reduces the time to reach critical shear displacement but also causes a significant temperature rise during the shear band formation. Based on this observation, the variation of shear band spacing in a typical dynamic indentation test has been rationalized. The temperature rise within a shear band can be in excess of 2000K at high strain rate and high confinement pressure conditions. The associated drop in viscosity and flow stress may explain the observed decrease in fracture strength and indentation hardness. The above investigations provide valuable insight into the deformation behavior of BMGs under static and dynamic loading conditions. The shear band patterns observed in the above indentation studies can be helpful to understand and model the deformation features under complex loading scenarios such as the interaction of a penetrator with armor. Future work encompasses (1) extending and modifying the coupled thermo-mechanical model to account for the temperature rise in quasistatic deformation; and (2) expanding this model to account for the microstructural variation-crystallization and free volume migration associated with the deformation.

Critical evaluation of diameter increment modelling in the Great Lakes region

Simulations of forest stand dynamics in a modelling framework including Forest Vegetation Simulator (FVS) are diameter driven, thus the diameter or basal area increment model needs a special attention. This dissertation critically evaluates diameter or basal area increment models and modelling approaches in the context of the Great Lakes region of the United States and Canada. A set of related studies are presented that critically evaluate the sub-model for change in individual tree basal diameter used in the Forest Vegetation Simulator (FVS), a dominant forestry model in the Great Lakes region. Various historical implementations of the STEMS (Stand and Tree Evaluation and Modeling System) family of diameter increment models, including the current public release of the Lake States variant of FVS (LS-FVS), were tested for the 30 most common tree species using data from the Michigan Forest Inventory and Analysis (FIA) program. The results showed that current public release of the LS-FVS diameter increment model over-predicts 10-year diameter increment by 17% on average. Also the study affirms that a simple adjustment factor as a function of a single predictor, dbh (diameter at breast height) used in the past versions, provides an inadequate correction of model prediction bias. In order to re-engineer the basal diameter increment model, the historical, conceptual and philosophical differences among the individual tree increment model families and their modelling approaches were analyzed and discussed. Two underlying conceptual approaches toward diameter or basal area increment modelling have been often used: the potential-modifier (POTMOD) and composite (COMP) approaches, which are exemplified by the STEMS/TWIGS and Prognosis models, respectively. It is argued that both approaches essentially use a similar base function and neither is conceptually different from a biological perspective, even though they look different in their model forms. No matter what modelling approach is used, the base function is the foundation of an increment model. Two base functions – gamma and Box-Lucas – were identified as candidate base functions for forestry applications. The results of a comparative analysis of empirical fits showed that quality of fit is essentially similar, and both are sufficiently detailed and flexible for forestry applications. The choice of either base function in order to model diameter or basal area increment is dependent upon personal preference; however, the gamma base function may be preferred over the Box-Lucas, as it fits the periodic increment data in both a linear and nonlinear composite model form. Finally, the utility of site index as a predictor variable has been criticized, as it has been widely used in models for complex, mixed species forest stands though not well suited for this purpose. An alternative to site index in an increment model was explored, using site index and a combination of climate variables and Forest Ecosystem Classification (FEC) ecosites and data from the Province of Ontario, Canada. The results showed that a combination of climate and FEC ecosites variables can replace site index in the diameter increment model.