Localized geodetic deformation of an approximately arctangent shape around locked strike-slip faults is widely reported, but there are also important exceptions showing distributed interseismic deformation. Understanding the controlling mechanism is important to the interpretation of geodetic observations for hazard assessment and geodynamic analysis. In this thesis, I use simple finite element models to separately study the two major contributors to the deformation: far-field loading and previous earthquakes. The models feature a vertical strike-slip fault in an elastic layer overlying a viscoelastic substrate of Maxwell or Burgers rheology, with or without weaknesses representing extensions of the fault either along strike or to greater depth. If the locked fault is loaded only from the far field without the effects of previous earthquakes, localized deformation occurs only if local mechanical weaknesses below the fault and/or somewhere along strike are introduced. I first show that the effects of far-field loading are rather limited even in the presence of extreme weaknesses. Then I use idealized earthquake cycle models to investigate the effects of past seismic events in a viscoelastic Earth. I demonstrate that, after a phase of fast postseismic deformation just after the earthquake, the localization of interseismic deformation is controlled mainly by the recurrence interval of past earthquakes. Given viscosity, shorter recurrence leads to greater interseismic localization, regardless of the rheological model used. The presence of a low-viscosity deep fault zone does not change this conclusion, although it tends to lessen localization by promoting faster postseismic stress relaxation. Distributed interseismic deformation, although less reported in the literature, is a natural consequence of very long recurrence and in theory should be as common as localized deformation. The apparent propensity of the latter is likely associated with the much greater quantity and better quality of geodetic observations from higher-rate and shorter-recurrence faults. Using viscoelastic earthquake-cycle models, I also explore the role of nearby earthquakes and creeping segments along the same fault. For faults of relatively short recurrence, frequent ruptures of nearby segments, modelled using a migrating rupture sequence with or without temporal clustering, further enhance localization. For faults of very long recurrence, faster near-fault deformation induced by a recent earthquake may give a false impression of localized interseismic deformation. / Graduate

With the recent increase in the threat of biological terrorism as well as the continual risk of other diseases, the research in public health surveillance and disease monitoring has grown tremendously. There is an abundance of data available in all sorts of forms. Hospitals, federal and local governments, and industries are all collecting data and developing new methods to be used in the detection of anomalies. Many of these methods are developed, applied to a real data set, and incorporated into software. This research, however, takes a different view of the evaluation of these methods. We feel that there needs to be solid statistical evaluation of proposed methods no matter the intended area of application. Using proof-by-example does not seem reasonable as the sole evaluation criteria especially concerning methods that have the potential to have a great impact in our lives. For this reason, this research focuses on determining the properties of some of the most common anomaly detection methods. A distinction is made between metrics used for retrospective historical monitoring and those used for prospective on-going monitoring with the focus on the latter situation. Metrics such as the recurrence interval and time-to-signal measures are therefore the most applicable. These metrics, in conjunction with control charts such as exponentially weighted moving average (EWMA) charts and cumulative sum (CUSUM) charts, are examined. Two new time-to-signal measures, the average time-between-signal events and the average signal event length, are introduced to better compare the recurrence interval with the time-to-signal properties of surveillance schemes. The relationship commonly thought to exist between the recurrence interval and the average time to signal is shown to not exist once autocorrelation is present in the statistics used for monitoring. This means that closer consideration needs to be paid to the selection of which of these metrics to report. The properties of a commonly applied scan method are also studied carefully in the strictly temporal setting. The counts of incidences are assumed to occur independently over time and follow a Poisson distribution. Simulations are used to evaluate the method under changes in various parameters. In addition, there are two methods proposed in the literature for the calculation of the p-value, an adjustment based on the tests for previous time periods and the use of the recurrence interval with no adjustment for previous tests. The difference in these two methods is also considered. The quickness of the scan method in detecting an increase in the incidence rate as well as the number of false alarm events that occur and how long the method signals after the increase threat has passed are all of interest. These estimates from the scan method are compared to other attribute monitoring methods, mainly the Poisson CUSUM chart. It is shown that the Poisson CUSUM chart is typically faster in the detection of the increased incidence rate. / Ph. D.

We examine some of the methodologies implemented by the Centers for Disease Control and Prevention's (CDC) BioSense program. The program uses data from hospitals and public health departments to detect outbreaks using the Early Aberration Reporting System (EARS). The EARS method W2 allows one to monitor syndrome counts (W2count) from each source and the proportion of counts of a particular syndrome relative to the total number of visits (W2rate). We investigate the performance of the W2r method designed using an empiric recurrence interval (RI) in this dissertation research. An adaptive threshold monitoring method is introduced based on fitting sample data to the underlying distributions, then converting the current value to a Z-score through a p-value. We compare the upper thresholds on the Z-scores required to obtain given values of the recurrence interval for different sets of parameter values. We then simulate one-week outbreaks in our data and calculate the proportion of times these methods correctly signal an outbreak using Shewhart and exponentially weighted moving average (EWMA) charts. Our results indicate the adaptive threshold method gives more consistent statistical performance across different parameter sets and amounts of baseline historical data used for computing the statistics. For the power analysis, the EWMA chart is superior to its Shewhart counterpart in nearly all cases, and the adaptive threshold method tends to outperform the W2 rate method. Two modified W2r methods proposed in the dissertation also tend to outperform the W2r method in terms of the RI threshold functions and in the power analysis. / Ph. D.

Cette thèse examine les impacts sur la morphologie des tributaires du fleuve Saint-Laurent des changements dans leur débit et leur niveau de base engendrés par les changements climatiques prévus pour la période 2010–2099. Les tributaires sélectionnés (rivières Batiscan, Richelieu, Saint-Maurice, Saint-François et Yamachiche) ont été choisis en raison de leurs différences de taille, de débit et de contexte morphologique. Non seulement ces tributaires subissent-ils un régime hydrologique modifié en raison des changements climatiques, mais leur niveau de base (niveau d’eau du fleuve Saint-Laurent) sera aussi affecté. Le modèle morphodynamique en une dimension (1D) SEDROUT, à l’origine développé pour des rivières graveleuses en mode d’aggradation, a été adapté pour le contexte spécifique des tributaires des basses-terres du Saint-Laurent afin de simuler des rivières sablonneuses avec un débit quotidien variable et des fluctuations du niveau d’eau à l’aval. Un module pour simuler le partage des sédiments autour d’îles a aussi été ajouté au modèle. Le modèle ainsi amélioré (SEDROUT4-M), qui a été testé à l’aide de simulations à petite échelle et avec les conditions actuelles d’écoulement et de transport de sédiments dans quatre tributaires du fleuve Saint-Laurent, peut maintenant simuler une gamme de problèmes morphodynamiques de rivières. Les changements d’élévation du lit et d’apport en sédiments au fleuve Saint-Laurent pour la période 2010–2099 ont été simulés avec SEDROUT4-M pour les rivières Batiscan, Richelieu et Saint-François pour toutes les combinaisons de sept régimes hydrologiques (conditions actuelles et celles prédites par trois modèles de climat globaux (MCG) et deux scénarios de gaz à effet de serre) et de trois scénarios de changements du niveau de base du fleuve Saint-Laurent (aucun changement, baisse graduelle, baisse abrupte). Les impacts sur l’apport de sédiments et l’élévation du lit diffèrent entre les MCG et semblent reliés au statut des cours d’eau (selon qu’ils soient en état d’aggradation, de dégradation ou d’équilibre), ce qui illustre l’importance d’examiner plusieurs rivières avec différents modèles climatiques afin d’établir des tendances dans les effets des changements climatiques. Malgré le fait que le débit journalier moyen et le débit annuel moyen demeurent près de leur valeur actuelle dans les trois scénarios de MCG, des changements importants dans les taux de transport de sédiments simulés pour chaque tributaire sont observés. Ceci est dû à l’impact important de fortes crues plus fréquentes dans un climat futur de même qu’à l’arrivée plus hâtive de la crue printanière, ce qui résulte en une variabilité accrue dans les taux de transport en charge de fond. Certaines complications avec l’approche de modélisation en 1D pour représenter la géométrie complexe des rivières Saint-Maurice et Saint-François suggèrent qu’une approche bi-dimensionnelle (2D) devrait être sérieusement considérée afin de simuler de façon plus exacte la répartition des débits aux bifurcations autour des îles. La rivière Saint-François est utilisée comme étude de cas pour le modèle 2D H2D2, qui performe bien d’un point de vue hydraulique, mais qui requiert des ajustements pour être en mesure de pleinement simuler les ajustements morphologiques des cours d’eau. / This thesis investigates the impacts of climate-induced changes in discharge and base level on the morphology of Saint-Lawrence River tributaries for the period 2010–2099. The selected tributaries (Batiscan, Richelieu, Saint-Maurice, Saint-François and Yamachiche rivers) were chosen because of their differences in size, flow regime and morphological setting. Not only will these tributaries experience an altered hydrological regime as a consequence of climate change, but their base level (Saint-Lawrence River water level) will also change. A one-dimensional (1D) morphodynamic model (SEDROUT), originally developed for aggrading gravel-bed rivers, was adapted for the specific context of the Saint-Lawrence lowland tributaries in order to simulate sand-bed rivers with variable daily discharge and downstream water level fluctuations. A module to deal with sediment routing in channels with islands was also added to the model. The enhanced model (SEDROUT4-M), which was tested with small-scale simulations and present-day conditions in four tributaries of the Saint-Lawrence River, can now simulate a very wide range of river morphodynamic problems. Changes in bed elevation and bed-material delivery to the Saint-Lawrence River over the 2010–2099 period were simulated with SEDROUT4-M for the Batiscan, Richelieu and Saint-François rivers for all combinations of seven tributary hydrological regimes (present-day and those predicted using three global climate models (GCM) and two greenhouse gas emission scenarios) and three scenarios of how the base level provided by the Saint-Lawrence River will alter (no change, gradual decrease, step decrease). The effects on mean annual sediment delivery and bed elevation differ between GCM and seem to be related to whether the river is currently aggrading, degrading or in equilibrium, which highlights the importance of investigating several rivers using several climate models in order to determine trends in climate change impacts. Despite the fact that mean daily discharge and mean annual maximum discharge remain close to their current values in the three GCM scenarios for daily discharge, marked changes occur in the mean annual sediment transport rates in each simulated tributary. This is due to the important effect of more frequent large individual flood events under future climate as well as a shift of peak annual discharge from the spring towards the winter, which results in increased variability of bed-material transport rates. Some complications with the 1D modelling approach to capture the complex geometry of the Saint-Maurice and Saint-François rivers suggest that the use of a two-dimensional (2D) approach should be seriously considered to accurately simulate the discharge distribution at bifurcations around islands. The Saint-François River is used as a test case for the 2D model H2D2, which performs well from a hydraulics point of view but which needs to be adapted to fully simulate morphological adjustments in the channel.

No description available.