Seismic sequences analysis for estimation of earthquake source parameters: corner frequency, stress drop, and seismic moment observations

The present study has been carried out with the following objectives: i) To investigate the attributes of source parameters of local and regional earthquakes; ii) To estimate, as accurately as possible, M0, fc, Δσ and their standard errors to infer their relationship with source size; iii) To quantify high-frequency earthquake ground motion and to study the source scaling. This work is based on observational data of micro, small and moderate -earthquakes for three selected seismic sequences, namely Parkfield (CA, USA), Maule (Chile) and Ferrara (Italy). For the Parkfield seismic sequence (CA), a data set of 757 (42 clusters) repeating micro-earthquakes (0 ≤ MW ≤ 2), collected using borehole High Resolution Seismic Network (HRSN), have been analyzed and interpreted. We used the coda methodology to compute spectral ratios to obtain accurate values of fc , Δσ, and M0 for three target clusters (San Francisco, Los Angeles, and Hawaii) of our data. We also performed a general regression on peak ground velocities to obtain reliable seismic spectra of all earthquakes. For the Maule seismic sequence, a data set of 172 aftershocks of the 2010 MW 8.8 earthquake (3.7 ≤ MW ≤ 6.2), recorded by more than 100 temporary broadband stations, have been analyzed and interpreted to quantify high-frequency earthquake ground motion in this subduction zone. We completely calibrated the excitation and attenuation of the ground motion in Central Chile. For the Ferrara sequence, we calculated moment tensor solutions for 20 events from MW 5.63 (the largest main event occurred on May 20 2012), down to MW 3.2 by a 1-D velocity model for the crust beneath the Pianura Padana, using all the geophysical and geological information available for the area. The PADANIA model allowed a numerical study on the characteristics of the ground motion in the thick sediments of the flood plain.

Toward a Magnitude Scale for Low Frequency Seismicity The LP energy for Campi Flegrei caldera

The energy released during a seismic crisis in volcanic areas is strictly related to the physical processes in the volcanic structure. In particular Long Period seismicity, that seems to be related to the oscillation of a fluid-filled crack (Chouet , 1996, Chouet, 2003, McNutt, 2005), can precedes or accompanies an eruption. The present doctoral thesis is focused on the study of the LP seismicity recorded in the Campi Flegrei volcano (Campania, Italy) during the October 2006 crisis. Campi Flegrei Caldera is an active caldera; the combination of an active magmatic system and a dense populated area make the Campi Flegrei a critical volcano. The source dynamic of LP seismicity is thought to be very different from the other kind of seismicity ( Tectonic or Volcano Tectonic): it’s characterized by a time sustained source and a low content in frequency. This features implies that the duration–magnitude, that is commonly used for VT events and sometimes for LPs as well, is unadapted for LP magnitude evaluation. The main goal of this doctoral work was to develop a method for the determination of the magnitude for the LP seismicity; it’s based on the comparison of the energy of VT event and LP event, linking the energy to the VT moment magnitude. So the magnitude of the LP event would be the moment magnitude of a VT event with the same energy of the LP. We applied this method to the LP data-set recorded at Campi Flegrei caldera in 2006, to an LP data-set of Colima volcano recorded in 2005 – 2006 and for an event recorded at Etna volcano. Experimenting this method to lots of waveforms recorded at different volcanoes we tested its easy applicability and consequently its usefulness in the routinely and in the quasi-real time work of a volcanological observatory.

Palaeo wind system reconstruction of the last glacial period over Europe, using high resolution proxy data and model-data-comparison

The atmosphere is a global influence on the movement of heat and humidity between the continents, and thus significantly affects climate variability. Information about atmospheric circulation are of major importance for the understanding of different climatic conditions. Dust deposits from maar lakes and dry maars from the Eifel Volcanic Field (Germany) are therefore used as proxy data for the reconstruction of past aeolian dynamics.rnrnIn this thesis past two sediment cores from the Eifel region are examined: the core SM3 from Lake Schalkenmehren and the core DE3 from the Dehner dry maar. Both cores contain the tephra of the Laacher See eruption, which is dated to 12,900 before present. Taken together the cores cover the last 60,000 years: SM3 the Holocene and DE3 the marine isotope stages MIS-3 and MIS-2, respectively. The frequencies of glacial dust storm events and their paleo wind direction are detected by high resolution grain size and provenance analysis of the lake sediments. Therefore two different methods are applied: geochemical measurements of the sediment using µXRF-scanning and the particle analysis method RADIUS (rapid particle analysis of digital images by ultra-high-resolution scanning of thin sections).rnIt is shown that single dust layers in the lake sediment are characterized by an increased content of aeolian transported carbonate particles. The limestone-bearing Eifel-North-South zone is the most likely source for the carbonate rich aeolian dust in the lake sediments of the Dehner dry maar. The dry maar is located on the western side of the Eifel-North-South zone. Thus, carbonate rich aeolian sediment is most likely to be transported towards the Dehner dry maar within easterly winds. A methodology is developed which limits the detection to the aeolian transported carbonate particles in the sediment, the RADIUS-carbonate module.rnrnIn summary, during the marine isotope stage MIS-3 the storm frequency and the east wind frequency are both increased in comparison to MIS-2. These results leads to the suggestion that atmospheric circulation was affected by more turbulent conditions during MIS-3 in comparison to the more stable atmospheric circulation during the full glacial conditions of MIS-2.rnThe results of the investigations of the dust records are finally evaluated in relation a study of atmospheric general circulation models for a comprehensive interpretation. Here, AGCM experiments (ECHAM3 and ECHAM4) with different prescribed SST patterns are used to develop a synoptic interpretation of long-persisting east wind conditions and of east wind storm events, which are suggested to lead to an enhanced accumulation of sediment being transported by easterly winds to the proxy site of the Dehner dry maar.rnrnThe basic observations made on the proxy record are also illustrated in the 10 m-wind vectors in the different model experiments under glacial conditions with different prescribed sea surface temperature patterns. Furthermore, the analysis of long-persisting east wind conditions in the AGCM data shows a stronger seasonality under glacial conditions: all the different experiments are characterized by an increase of the relative importance of the LEWIC during spring and summer. The different glacial experiments consistently show a shift from a long-lasting high over the Baltic Sea towards the NW, directly above the Scandinavian Ice Sheet, together with contemporary enhanced westerly circulation over the North Atlantic.rnrnThis thesis is a comprehensive analysis of atmospheric circulation patterns during the last glacial period. It has been possible to reconstruct important elements of the glacial paleo climate in Central Europe. While the proxy data from sediment cores lead to a binary signal of the wind direction changes (east versus west wind), a synoptic interpretation using atmospheric circulation models is successful. This shows a possible distribution of high and low pressure areas and thus the direction and strength of wind fields which have the capacity to transport dust. In conclusion, the combination of numerical models, to enhance understanding of processes in the climate system, with proxy data from the environmental record is the key to a comprehensive approach to paleo climatic reconstruction.rn

Towards the 3D attenuation imaging of active volcanoes: methods and tests on real and simulated data

The purpose of my PhD thesis has been to face the issue of retrieving a three dimensional attenuation model in volcanic areas. To this purpose, I first elaborated a robust strategy for the analysis of seismic data. This was done by performing several synthetic tests to assess the applicability of spectral ratio method to our purposes. The results of the tests allowed us to conclude that: 1) spectral ratio method gives reliable differential attenuation (dt*) measurements in smooth velocity models; 2) short signal time window has to be chosen to perform spectral analysis; 3) the frequency range over which to compute spectral ratios greatly affects dt* measurements. Furthermore, a refined approach for the application of spectral ratio method has been developed and tested. Through this procedure, the effects caused by heterogeneities of propagation medium on the seismic signals may be removed. The tested data analysis technique was applied to the real active seismic SERAPIS database. It provided a dataset of dt* measurements which was used to obtain a three dimensional attenuation model of the shallowest part of Campi Flegrei caldera. Then, a linearized, iterative, damped attenuation tomography technique has been tested and applied to the selected dataset. The tomography, with a resolution of 0.5 km in the horizontal directions and 0.25 km in the vertical direction, allowed to image important features in the off-shore part of Campi Flegrei caldera. High QP bodies are immersed in a high attenuation body (Qp=30). The latter is well correlated with low Vp and high Vp/Vs values and it is interpreted as a saturated marine and volcanic sediments layer. High Qp anomalies, instead, are interpreted as the effects either of cooled lava bodies or of a CO2 reservoir. A pseudo-circular high Qp anomaly was detected and interpreted as the buried rim of NYT caldera.

Non-thermal emission in High Frequency Peaked blazars towards the Square Kilometer Array era

In this thesis work we will explore and discuss the properties of the gamma-ray sources included in the first Fermi-LAT catalog of sources above 10 GeV (1FHL), by considering both blazars and the non negligible fraction of still unassociated gamma-ray sources (UGS, 13%). We perform a statistical analysis of a complete sample of hard gamma-ray sources, included in the 1FHL catalog, mostly composed of HSP blazars, and we present new VLBI observations of the faintest members of the sample. The new VLBI data, complemented by an extensive search of the archives for brighter sources, are essential to gather a sample as large as possible for the assessment of the significance of the correlation between radio and very high energy (E>100 GeV) emission bands. After the characterization of the statistical properties of HSP blazars and UGS, we use a complementary approach, by focusing on an intensive multi-frequency observing VLBI and gamma-ray campaign carried out for one of the most remarkable and closest HSP blazar Markarian 421.

The effect of arts speech therapy on cerebral oxygenation and low frequency hemodynamic oscillations measured using functional near-infrared spectroscopy

Introduction: As a previous study revealed, arts speech therapy (AST) affects cardiorespiratory interaction [1]. The aim of the present study was to investigate whether AST also has effects on brain oxygenation and hemodynamics measured non-invasively using near-infrared spectroscopy (NIRS). Material and methods: NIRS measurements were performed on 17 subjects (8 men and 9 women, mean age: 35.6 ± 12.7 y) during AST. Each measurement lasted 35 min, comprising 8 min pre-baseline, 10 min recitation and 20 min post-baseline. For each subject, measurements were performed for three different AST recitation tasks (recitation of alliterative, hexameter and prose verse). Relative concentration changes of oxyhemoglobin (Δ[O2Hb]) and deoxyhemoglobin (Δ[HHb]) as well as the tissue oxygenation index (TOI) were measured using a Hamamatsu NIRO300 NIRS device and a sensor placed on the subjects forehead. Movement artifacts were removed using a novel method [2]. Statistical analysis (Wilcoxon test) was applied to the data to investigate (i) if the recitation causes changes in the median values and/or in the Mayer wave power spectral density (MW-PSD, range: 0.07–0.13 Hz) of Δ[O2Hb], Δ[HHb] or TOI, and (ii) if these changes vary between the 3 recitation forms. Results: For all three recitation styles a significant (p < 0.05) decrease in Δ[O2Hb] and TOI was found, indicating a decrease in blood flow. These decreases did not vary significantly between the three styles. MW-PSD increased significantly for Δ[O2Hb] when reciting the hexameter and prose verse, and for Δ[HHb] and TOI when reciting alliterations and hexameter, representing an increase in Mayer waves. The MW-PSD increase for Δ[O2Hb] was significantly larger for the hexameter verse compared to alliterative and prose verse Conclusion: The study showed that AST affects brain hemodynamics (oxygenation, blood flow and Mayer waves). Recitation caused a significant decrease in cerebral blood flow for all recitation styles as well as an increase in Mayer waves, particularly for the hexameter, which may indicate a sympathetic activation. References 1. D. Cysarz, D. von Bonin, H. Lackner, P. Heusser, M. Moser, H. Bettermann. Am J Physiol Heart Circ Physiol, 287 (2) (2004), pp. H579–H587 2. F. Scholkmann, S. Spichtig, T. Muehlemann, M. Wolf. Physiol Meas, 31 (5) (2010), pp. 649–662

Frequency of and predictors for patient-reported medical and medication errors in Switzerland

OBJECTIVES: To analyse the frequency of and identify risk factors for patient-reported medical errors in Switzerland. The joint effect of risk factors on error-reporting probability was modelled for hypothetical patients. METHODS: A representative population sample of Swiss citizens (n = 1306) was surveyed as part of the Commonwealth Fund’s 2010 lnternational Survey of the General Public’s Views of their Health Care System’s Performance in Eleven Countries. Data on personal background, utilisation of health care, coordination of care problems and reported errors were assessed. Logistic regression analysis was conducted to identify risk factors for patients’ reports of medical mistakes and medication errors. RESULTS: 11.4% of participants reported at least one error in their care in the previous two years (8% medical errors, 5.3% medication errors). Poor coordination of care experiences was frequent. 7.8% experienced that test results or medical records were not available, 17.2% received conflicting information from care providers and 11.5% reported that tests were ordered although they had been done before. Age (OR = 0.98, p = 0.014), poor health (OR = 2.95, p = 0.007), utilisation of emergency care (OR = 2.45, p = 0.003), inpatient-stay (OR = 2.31, p = 0.010) and poor care coordination (OR = 5.43, p <0.001) are important predictors for reporting error. For high utilisers of care that unify multiple risk factors the probability that errors are reported rises up to p = 0.8. CONCLUSIONS: Patient safety remains a major challenge for the Swiss health care system. Despite the health related and economic burden associated with it, the widespread experience of medical error in some subpopulations also has the potential to erode trust in the health care system as a whole.

Subconductance states of Cx30 gap junction channels: data from transfected HeLa cells versus data from a mathematical model

Human HeLa cells expressing mouse connexin30 were used to study the electrical properties of gap junction channel substates. Experiments were performed on cell pairs using a dual voltage-clamp method. Single-channel currents revealed discrete levels attributable to a main state, a residual state, and five substates interposed, suggesting the operation of six subgates provided by the six connexins of a gap junction hemichannel. Substate conductances, gamma(j,substate), were unevenly distributed between the main-state and the residual-state conductance (gamma(j,main state) = 141 pS, gamma(j,residual state) = 21 pS). Activation of the first subgate reduced the channel conductance by approximately 30%, and activation of subsequent subgates resulted in conductance decrements of 10-15% each. Current transitions between the states were fast (<2 ms). Substate events were usually demarcated by transitions from and back to the main state; transitions among substates were rare. Hence, subgates are recruited simultaneously rather than sequentially. The incidence of substate events was larger at larger gradients of V(j). Frequency and duration of substate events increased with increasing number of synchronously activated subgates. Our mathematical model, which describes the operation of gap junction channels, was expanded to include channel substates. Based on the established V(j)-sensitivity of gamma(j,main state) and gamma(j,residual state), the simulation yielded unique functions gamma(j,substate) = f(V(j)) for each substate. Hence, the spacing of subconductance levels between the channel main state and residual state were uneven and characteristic for each V(j).

Higher frequency of chromosomal aberrations in ovarian endometriosis compared to extragonadal endometriosis: A possible link to endometrioid adenocarcinoma

Endometriosis may progress to invasive endometrioid adenocarcinoma, particularly in the ovary. Up to now, little is known of the molecular mechanisms possibly involved in the malignant transformation of endometriosis. Therefore, in this study, extragonadal endometriosis (n = 10), ovarian endometriosis without malignancy (n = 10), ovarian endometriosis with direct transition into endometrioid adenocarcinoma (n = 8), and normal endometrium (n = 12) were investigated for numerical chromosomal aberrations by fluorescence in situ hybridization using centromere enumeration probes. The proportions of cells with aneusomies were semiquantitatively assessed. Trisomies 1 and 7, and monosomies 9 and 17 were found in endometriosis, ovarian endometrioid adenocarcinoma, and normal endometrium. The proportions of aneusomic cells were significantly higher in ovarian endometrioid carcinoma compared with ovarian endometriosis (P < 0.001), and in ovarian endometriosis compared with extragonadal endometriosis and normal endometrium (P < 0.001). The data provide new evidence of a common lineage of endometriosis and ovarian endometrioid carcinoma. The higher frequency of chromosomal aberrations in endometrioid carcinoma than in endometriosis may reflect an expansion of aberrant cell clones already present in endometriosis during the progression to cancer. The higher frequency of chromosomal aberrations in ovarian endometriosis than in extragonadal endometriosis suggests a role of the ovarian stromal milieu in the induction of genetic changes, which may eventually lead to invasive cancer.

Alternative protocol to initiate high-frequency oscillatory ventilation: an experimental study

INTRODUCTION: The objective was to study the effects of a novel lung volume optimization procedure (LVOP) using high-frequency oscillatory ventilation (HFOV) upon gas exchange, the transpulmonary pressure (TPP), and hemodynamics in a porcine model of surfactant depletion. METHODS: With institutional review board approval, the hemodynamics, blood gas analysis, TPP, and pulmonary shunt fraction were obtained in six anesthetized pigs before and after saline lung lavage. Measurements were acquired during pressure-controlled ventilation (PCV) prior to and after lung damage, and during a LVOP with HFOV. The LVOP comprised a recruitment maneuver with a continuous distending pressure (CDP) of 45 mbar for 2.5 minutes, and a stepwise decrease of the CDP (5 mbar every 5 minute) from 45 to 20 mbar. The TPP level was identified during the decrease in CDP, which assured a change of the PaO2/FIO2 ratio < 25% compared with maximum lung recruitment at CDP of 45 mbar (CDP45). Data are presented as the median (25th-75th percentile); differences between measurements are determined by Friedman repeated-measures analysis on ranks and multiple comparisons (Tukey's test). The level of significance was set at P < 0.05. RESULTS: The PaO2/FiO2 ratio increased from 99.1 (56.2-128) Torr at PCV post-lavage to 621 (619.4-660.3) Torr at CDP45 (CDP45) (P < 0.031). The pulmonary shunt fraction decreased from 51.8% (49-55%) at PCV post-lavage to 1.03% (0.4-3%) at CDP45 (P < 0.05). The cardiac output and stroke volume decreased at CDP45 (P < 0.05) compared with PCV, whereas the heart rate, mean arterial pressure, and intrathoracic blood volume remained unchanged. A TPP of 25.5 (17-32) mbar was required to preserve a difference in PaO2/FIO2 ratio < 25% related to CDP45; this TPP was achieved at a CDP of 35 (25-40) mbar. CONCLUSION: This HFOV protocol is easy to perform, and allows a fast determination of an adequate TPP level that preserves oxygenation. Systemic hemodynamics, as a measure of safety, showed no relevant deterioration throughout the procedure.

Childhood extraordinary daytime urinary frequency-a case series and a systematic literature review

Childhood extraordinary daytime urinary frequency is likely a common but underreported condition characterized by daytime frequent voiding and typically not linked with complaints of burning, urinary incontinence, altered urinary stream, changes in the nighttime voiding pattern, excessive fluid intake and excessive urinary volume. To determine the features and outcome of extraordinary daytime urinary frequency, we report our experience with 14 children and the results of a formal systematic analysis of peer-reviewed English-language literature on this topic. Nineteen case series were found (together with 16 mostly pertinent comments), with each case series providing details on from one to 119 children. On the basis of our experience and the findings of our systematic analysis, we conclude that, in general practice, extraordinary daytime urinary frequency is a common cause of urinary frequency, that the age of such patients is, on average, 6 years and that the micturation abnormalities persist for an average of 6 months. The results of this review must be viewed with an understanding of the limitations of the analysis process, which incorporated data exclusively from case series.

Frequency of spontaneous ventricular tachycardia in a pediatric population

Few data exist on the incidence of spontaneously occurring ventricular tachycardia (VT) in an unselected pediatric population. The aim of this study was to define the incidence and outcomes of VT in a general pediatric population. A retrospective analysis was performed of all documented episodes of VT in children referred to a single center during a 10-year study period ending in December 2005. The study center drains a stable referral area with 252,000 children aged <16 years, with no other pediatric cardiologic or pediatric intensive care services available. Twenty-seven patients with spontaneously occurring episodes of VT were observed, accounting for a VT incidence of 1.1 episodes/100,000 childhood years. Thirteen patients had VT in the absence of structural heart disease, and 14 had VT in the presence of a wide range of underlying cardiac disease. Overall mortality was 5 of 27 patients (19%), but mortality was seen exclusively in patients with underlying heart disease; for this subgroup of patients, mortality was 36%. Idiopathic VT in children with structurally normal hearts carried a good prognosis, and treatment was required in a minority (20%) of these patients. In conclusion, this study highlights that VT in childhood is rare, and outcomes are highly dependent on the underlying pathologic substrate.

A comparison of ground cover and frequency estimation methods for post-harvest soil monitoring

The amount and type of ground cover is an important characteristic to measure when collecting soil disturbance monitoring data after a timber harvest. Estimates of ground cover and bare soil can be used for tracking changes in invasive species, plant growth and regeneration, woody debris loadings, and the risk of surface water runoff and soil erosion. A new method of assessing ground cover and soil disturbance was recently published by the U.S. Forest Service, the Forest Soil Disturbance Monitoring Protocol (FSDMP). This protocol uses the frequency of cover types in small circular (15cm) plots to compare ground surface in pre- and post-harvest condition. While both frequency and percent cover are common methods of describing vegetation, frequency has rarely been used to measure ground surface cover. In this study, three methods for assessing ground cover percent (step-point, 15cm dia. circular and 1x5m visual plot estimates) were compared to the FSDMP frequency method. Results show that the FSDMP method provides significantly higher estimates of ground surface condition for most soil cover types, except coarse wood. The three cover methods had similar estimates for most cover values. The FSDMP method also produced the highest value when bare soil estimates were used to model erosion risk. In a person-hour analysis, estimating ground cover percent in 15cm dia. plots required the least sampling time, and provided standard errors similar to the other cover estimates even at low sampling intensities (n=18). If ground cover estimates are desired in soil monitoring, then a small plot size (15cm dia. circle), or a step-point method can provide a more accurate estimate in less time than the current FSDMP method.

Soft-decision equalization techniques for frequency selective MIMO channels

Multi-input multi-output (MIMO) technology is an emerging solution for high data rate wireless communications. We develop soft-decision based equalization techniques for frequency selective MIMO channels in the quest for low-complexity equalizers with BER performance competitive to that of ML sequence detection. We first propose soft decision equalization (SDE), and demonstrate that decision feedback equalization (DFE) based on soft-decisions, expressed via the posterior probabilities associated with feedback symbols, is able to outperform hard-decision DFE, with a low computational cost that is polynomial in the number of symbols to be recovered, and linear in the signal constellation size. Building upon the probabilistic data association (PDA) multiuser detector, we present two new MIMO equalization solutions to handle the distinctive channel memory. With their low complexity, simple implementations, and impressive near-optimum performance offered by iterative soft-decision processing, the proposed SDE methods are attractive candidates to deliver efficient reception solutions to practical high-capacity MIMO systems. Motivated by the need for low-complexity receiver processing, we further present an alternative low-complexity soft-decision equalization approach for frequency selective MIMO communication systems. With the help of iterative processing, two detection and estimation schemes based on second-order statistics are harmoniously put together to yield a two-part receiver structure: local multiuser detection (MUD) using soft-decision Probabilistic Data Association (PDA) detection, and dynamic noise-interference tracking using Kalman filtering. The proposed Kalman-PDA detector performs local MUD within a sub-block of the received data instead of over the entire data set, to reduce the computational load. At the same time, all the inter-ference affecting the local sub-block, including both multiple access and inter-symbol interference, is properly modeled as the state vector of a linear system, and dynamically tracked by Kalman filtering. Two types of Kalman filters are designed, both of which are able to track an finite impulse response (FIR) MIMO channel of any memory length. The overall algorithms enjoy low complexity that is only polynomial in the number of information-bearing bits to be detected, regardless of the data block size. Furthermore, we introduce two optional performance-enhancing techniques: cross- layer automatic repeat request (ARQ) for uncoded systems and code-aided method for coded systems. We take Kalman-PDA as an example, and show via simulations that both techniques can render error performance that is better than Kalman-PDA alone and competitive to sphere decoding. At last, we consider the case that channel state information (CSI) is not perfectly known to the receiver, and present an iterative channel estimation algorithm. Simulations show that the performance of SDE with channel estimation approaches that of SDE with perfect CSI.

Cardiovascular mortality and exposure to extremely low frequency magnetic fields: a cohort study of Swiss railway workers

BACKGROUND: Exposure to intermittent magnetic fields of 16 Hz has been shown to reduce heart rate variability, and decreased heart rate variability predicts cardiovascular mortality. We examined mortality from cardiovascular causes in railway workers exposed to varying degrees to intermittent 16.7 Hz magnetic fields. METHODS: We studied a cohort of 20,141 Swiss railway employees between 1972 and 2002, including highly exposed train drivers (median lifetime exposure 120.5 muT-years), and less or little exposed shunting yard engineers (42.1 muT-years), train attendants (13.3 muT-years) and station masters (5.7 muT-years). During 464,129 person-years of follow up, 5,413 deaths were recorded and 3,594 deaths were attributed to cardio-vascular diseases. We analyzed data using Cox proportional hazards models. RESULTS: For all cardiovascular mortality the hazard ratio compared to station masters was 0.99 (95%CI: 0.91, 1.08) in train drivers, 1.13 (95%CI: 0.98, 1.30) in shunting yard engineers, and 1.09 (95%CI: 1.00, 1.19) in train attendants. Corresponding hazard ratios for arrhythmia related deaths were 1.04 (95%CI: 0.68, 1.59), 0.58 (95%CI: 0.24, 1.37) and 10 (95%CI: 0.87, 1.93) and for acute myocardial infarction 1.00 (95%CI: 0.73, 1.36), 1.56 (95%CI: 1.04, 2.32), and 1.14 (95%CI: 0.85, 1.53). The hazard ratio for arrhythmia related deaths per 100 muT-years of cumulative exposure was 0.94 (95%CI: 0.71, 1.24) and 0.91 (95%CI: 0.75, 1.11) for acute myocardial infarction. CONCLUSION: This study provides evidence against an association between long-term occupational exposure to intermittent 16.7 Hz magnetic fields and cardiovascular mortality.