电离层行进式扰动的GPS台网监测研究

In this dissertation, we investigated two types of traveling ionospheric disturbances (TIDs)/gravity waves (GWs) triggered separately by auroral energy input during super geomagnetic storms and solar terminator (ST) under quiet geomagnetic conditions (kp<3+) using TEC measurements from the global network of GPS receivers. Research into the generation and propagation of TIDs/GWs during storms greatly enhance our understandings on the evolution processes of energy transportation from the high-latitude’s magnetosphere to the low-latitude ionosphere and the conjugated effect of TIDs propagation between the northern and southern hemispheres. Our results revealed that the conjugacy of propagation direction between the northern and southern hemispheres was subject to the influence of Coriolis force. We also figure out the evolution processes of ionospheric disturbances at the global scale. These are important topics that had not been well addressed previously. In addition, we also obtained thee wave structures of medium scale TIDs excited by the solar terminator (ST) moving over the northern America and physical mechanisms involved. Our observations confirm that the ST is a stable and repetitive source of ionospheric wave disturbances and the evidence of solar terminator generated disturbances has been demonstrated experimentally via the GPS TEC measurement. The main researches and results of this dissertation are as follows. First, the global traveling ionospheric disturbances (TIDs) during the drastic magnetic storms of October 29–31, 2003 were analyzed using the Global Position System (GPS) total electron content (TEC) data observed in the Asian-Australian, European and North American sectors. We collected the most comprehensive set of the TEC data from more than 900 GPS stations on the International GNSS Services (IGS) website and introduce here a strategy that combines polynomial fitting and multi-channel maximum entropy spectral analysis to obtain TID parameters. Moreover, in collaboration with my thesis advisor, I have developed an imaging technique of 2-dimensional map of TIDs structures to obtain spatial and temporal maps of large scale traveling ionospheric disturbances (LSTIDs). The clear structures of TEC perturbations map during the passage of TIDs were displayed. The results of our study are summarized as follows: (1) Large-scale TIDs (LSTIDs) and medium-scale TIDs (MSTIDs) were detected in all three sectors after the sudden commencement (SC) of the magnetic storm, and their features showed longitudinal and latitudinal dependences. The duration of TIDs was longer at higher latitudes than at middle latitudes, with a maximum of about 16 h. The TEC variation amplitude of LSTIDs was larger in the North American sector than in the two other sectors. At the lower latitudes, the ionospheric perturbations were more complicated, and their duration and amplitude were relatively longer and larger. (2) The periods and phase speeds of TIDs were different in these three sectors. In Europe, the TIDs propagated southward; in North America and Asia, the TIDs propagated southwestward; in the near-equator region, the disturbances propagated with the azimuth (the angle of the propagation direction of the LSTIDs measured clockwise from due north with 0°) of 210° showing the influence of Coriolis force; in the Southern Hemisphere, the LSTIDs propagated conjugatedly northwestward. Both the southwestward and northeastward propagating LSTIDs are found in the equatorial region. These results mean that the Coriolis effect cannot be ignored for the wave propagation of LSTIDs and that the propagation direction is correlated with the polar magnetic activity. (3) The day (day of year: 301) before the SC (sudden commencement) of magnetic storm, we observed a sudden TEC skip disturbances (±10 TECU). It should be a response for the high flux of proton during the solar flare event, but not the magnetic storms. Next, the most comprehensive and dense GPS network’s data from North-America region were used in this paper to analyze the medium scale traveling ionospheric disturbances (MSTIDs) which were generated by the moving solar terminator during the quiet days in 2005. We applied the multi-channel maximum entropy spectral analysis to calculated TID parameters, and found that the occurrence of ST-MSTIDs depends on the seasonal variations. The results of our study are summarized as follows: (1) MSTIDs stimulated by the moving ST (ST-MSTIDs) are detected at mid-latitudes after the passage of the solar terminator with the life time of 2～3 hours and the variation amplitude of 0.2～0.8 TECU. Spectral analysis indicated that the horizontal wavelength, average period, horizontal phase velocity of the MSTIDs are around 300±150 km，150±80 m/s and 25±15 min, respectively. In addition, ST-MSTIDs have wave fronts elongating the moving ST direction and almost parallel to ST. (2) The statistical results demonstrate that the dusk MSTIDs stimulated by ST is more obvious than the dawn MSTIDs in summer. On the contrary, the more-pronounced dawn MSTIDs occurs in winter. (3) Further analysis indicates that the seasonal variations of ST-MSTIDs occurrence frequency are most probably related to the seasonal differences of the variations of EUV flux in the ionosphere region and recombination process during sunrise and sunset period at mid-latitudes. Statistical study of occurrence characteristics of TIDs using the GPS network in North-American and European during solar maximum, In conclusion, statistical studies of the propagation characteristics of TIDs, which excited by the two common origins including geomagnetic storms and moving solar terminator, were involved with global GPS TEC databasein this thesis. We employed the multichannel maximum entropy spectral analysis method to diagnose the characteristics of propagation and evolvement of ionospheric disturbances, also, the characteristics of their regional distribution and climatological variations were revealed by the statistic analysis. The results of these studies can improve our knowledge about the energy transfer in the solar-terrestrial system and the coupling process between upper and lower atmosphere (thermosphere-ionosphere-mesosphere). On the other hand, our results of the investigation on TIDs generated by particular linear origin such as ST are important for developing ionospheric irregularity physics and modeling the transionosphere radio wave propagation. Besides, the GPS TEC representation of the ST-generated ionospheric structure suggests a better possibility for investigating this phenomenon. Subsequently, there are scientific meaning of the result of this dissertation to deeply discuss the energy transfer and coupling in the ionosphere, as well as realistic value to space weather forecast in the ionosphere region.

POLAR NETWORK INDEX AS A MAGNETIC PROXY FOR THE SOLAR CYCLE STUDIES

The Sun has a polar magnetic field which oscillates with the 11 yr sunspot cycle. This polar magnetic field is an important component of the dynamo process which operates in the solar convection zone and produces the sunspot cycle. We have direct systematic measurements of the Sun's polar magnetic field only from about the mid-1970s. There are, however, indirect proxies which give us information about this field at earlier times. The Ca-K spectroheliograms taken at the Kodaikanal Solar Observatory during 1904-2007 have now been digitized with 4k x 4k CCD and have higher resolution (similar to 0.86 arcsec) than the other available historical data sets. From these Ca-K spectroheliograms, we have developed a completely new proxy (polar network index, hereafter PNI) for the Sun's polar magnetic field. We calculate PNI from the digitized images using an automated algorithm and calibrate our measured PNI against the polar field as measured by the Wilcox Solar Observatory for the period 1976-1990. This calibration allows us to estimate the polar fields for the earlier period up to 1904. The dynamo calculations performed with this proxy as input data reproduce reasonably well the Sun's magnetic behavior for the past century.

Electric and magnetic losses and gains in determining the sign of refractive index

Within the framework of classic electromagnetic theories, we have studied the sign of refractive index of optical medias with the emphases on the roles of the electric and magnetic losses and gains. Starting from the Maxwell equations for an isotropic and homogeneous media, we have derived the general form of the complex refractive index and its relation with the complex electric permittivity and magnetic permeability, i.e. n = root epsilon mu, in which the intrinsic electric and magnetic losses and gains are included as the imaginary parts of the complex permittivity and permeability, respectively, as epsilon = epsilon(r) + i(epsilon i) and mu = mu(r) + i mu(i). The electric and magnetic losses are present in all passive materials, which correspond, respectively, to the positive imaginary permittivity and permeability epsilon(i) > 0 and mu(i) > 0. The electric and magnetic gains are present in materials where external pumping sources enable the light to be amplified instead of attenuated, which correspond, respectively, to the negative imaginary permittivity and permeability epsilon(i) < 0 and mu(i) < 0. We have analyzed and determined uniquely the sign of the refractive index, for all possible combinations of the four parameters epsilon(r), mu(r), epsilon(i), and mu(i), in light of the relativistic causality. A causal solution requires that the wave impedance be positive Re {Z} > 0. We illustrate the results for all cases in tables of the sign of refractive index. One of the most important messages from the sign tables is that, apart from the well-known case where simultaneously epsilon < 0 and mu < 0, there are other possibilities for the refractive index to be negative n < 0, for example, for epsilon(r) < 0, mu(r) > 0, epsilon(i) > 0, and mu(i) > 0, the refractive index is negative n < 0 provided mu(i)/epsilon(i) > mu(r)/vertical bar epsilon(r)vertical bar. (c) 2006 Elsevier B.V. All rights reserved.

Random walk of magnetic field lines for different values of the energy range spectral index

An analytical nonlinear description of field-line wandering in partially statistically magnetic systems was proposed recently. In this article the influence of the wave spectrum in the energy range onto field-line random walk is investigated by applying this formulation. It is demonstrated that in all considered cases we clearly obtain a superdiffusive behavior of the field-lines. If the energy range spectral index exceeds unity a free-streaming behavior of the field-lines can be found for all relevant length-scales of turbulence. Since the superdiffusive results obtained for the slab model are exact, it seems that superdiffusion is the normal behavior of field-line wandering.

X-ray fluorescence scanning of discrete samples on the Schwalbenberg II loess-paleosol sequence, raw data, magnetic susceptibility, organic carbon and U-ratio as grainsize index

The Schwalbenberg II loess-paleosol sequence (LPS) denotes a key site for Marine Isotope Stage (MIS 3) in Western Europe owing to eight succeeding cambisols, which primarily constitute the Ahrgau Subformation. Therefore, this LPS qualifies as a test candidate for the potential of temporal high-resolution geochemical data obtained X-ray fluorescence (XRF) scanning of discrete samplesproviding a fast and non-destructive tool for determining the element composition. The geochemical data is first contextualized to existing proxy data such as magnetic susceptibility (MS) and organic carbon (Corg) and then aggregated to element log ratios characteristic for weathering intensity [LOG (Ca/Sr), LOG (Rb/Sr), LOG (Ba/Sr), LOG (Rb/K)] and dust provenance [LOG (Ti/Zr), LOG (Ti/Al), LOG (Si/Al)]. Generally, an interpretation of rock magnetic particles is challenged in western Europe, where not only magnetic enhancement but also depletion plays a role. Our data indicates leaching and top-soil erosion induced MS depletion at the Schwalbenberg II LPS. Besides weathering, LOG (Ca/Sr) is susceptible for secondary calcification. Thus, also LOG (Rb/Sr) and LOG (Ba/Sr) are shown to be influenced by calcification dynamics. Consequently, LOG (Rb/K) seems to be the most suitable weathering index identifying the Sinzig Soils S1 and S2 as the most pronounced paleosols for this site. Sinzig Soil S3 is enclosed by gelic gleysols and in contrast to S1 and S2 only initially weathered pointing to colder climate conditions. Also the Remagen Soils are characterized by subtle to moderate positive excursions in the weathering indices. Comparing the Schwalbenberg II LPS with the nearby Eifel Lake Sediment Archive (ELSA) and other more distant German, Austrian and Czech LPS while discussing time and climate as limiting factors for pedogenesis, we suggest that the lithologically determined paleosols are in-situ soil formations. The provenance indices document a Zr-enrichment at the transition from the Ahrgau to the Hesbaye Subformation. This is explained by a conceptual model incorporating multiple sediment recycling and sorting effects in eolian and fluvial domains.

Magnetic resonance imaging analysis of the upper cervical spine extensor musculature in an asymptomatic cohort: an index of fat within muscle

AIM: To establish a simple method to quantify muscle/fat constituents in cervical muscles of asymptomatic women using magnetic resonance imaging (MRI), and to determine whether there is an age effect within a defined age range. MATERIALS AND METHODS: MRI of the upper cervical spine was performed for 42 asymptomatic women aged 18-45 years. The muscle and fat signal intensities on axial spin echo T1-weighted images were quantitatively classified by taking a ratio of the pixel intensity profiles of muscle against those of intermuscular fat for the rectus capitis posterior major and minor and inferior obliquus capitis muscles bilaterally. Inter- and intra-examiner agreement was scrutinized. RESULTS: The average relative values of fat within the upper cervical musculature compared with intermuscular fat indicated that there were only slight variations in indices between the three sets of muscles. There was no significant correlation between age and fat indices. There were significant differences for the relative fat within the muscle compared with intermuscular fat and body mass index for the right rectus capitis posterior major and right and left inferior obliquus capitis muscles (p = 0.032). Intraclass correlation coefficients for intraobserver agreement ranged from 0.94 to 0.98. Inter-rater agreement of the measurements ranged from 0.75 to 0.97. CONCLUSION: A quantitative measure of muscle/fat constituents has been developed, and results of this study indicate that relative fatty infiltration is not a feature of age in the upper cervical extensor muscles of women aged 18-45 years. (C) 2005 The Royal College of Radiologists. Published by Elsevier Ltd. All rights reserved.

Multifractal analysis of geomagnetic storm and solar flare indices and their class dependence

The multifractal properties of two indices of geomagnetic activity, D st (representative of low latitudes) and a p (representative of the global geomagnetic activity), with the solar X-ray brightness, X l , during the period from 1 March 1995 to 17 June 2003 are examined using multifractal detrended fluctuation analysis (MF-DFA). The h(q) curves of D st and a p in the MF-DFA are similar to each other, but they are different from that of X l , indicating that the scaling properties of X l are different from those of D st and a p . Hence, one should not predict the magnitude of magnetic storms directly from solar X-ray observations. However, a strong relationship exists between the classes of the solar X-ray irradiance (the classes being chosen to separate solar flares of class X-M, class C, and class B or less, including no flares) in hourly measurements and the geomagnetic disturbances (large to moderate, small, or quiet) seen in D st and a p during the active period. Each time series was converted into a symbolic sequence using three classes. The frequency, yielding the measure representations, of the substrings in the symbolic sequences then characterizes the pattern of space weather events. Using the MF-DFA method and traditional multifractal analysis, we calculate the h(q), D(q), and τ (q) curves of the measure representations. The τ (q) curves indicate that the measure representations of these three indices are multifractal. On the basis of this three-class clustering, we find that the h(q), D(q), and τ (q) curves of the measure representations of these three indices are similar to each other for positive values of q. Hence, a positive flare storm class dependence is reflected in the scaling exponents h(q) in the MF-DFA and the multifractal exponents D(q) and τ (q). This finding indicates that the use of the solar flare classes could improve the prediction of the D st classes.

Magnetic resonance spectroscopy and neurocognitive dysfunction in obstructive sleep apnea before and after CPAP treatment

STUDY OBJECTIVES: To determine whether cerebral metabolite changes may underlie abnormalities of neurocognitive function and respiratory control in OSA. DESIGN: Observational, before and after CPAP treatment. SETTING: Two tertiary hospital research institutes. PARTICIPANTS: 30 untreated severe OSA patients, and 25 age-matched healthy controls, all males free of comorbidities, and all having had detailed structural brain analysis using voxel-based morphometry (VBM). MEASUREMENTS AND RESULTS: Single voxel bilateral hippocampal and brainstem, and multivoxel frontal metabolite concentrations were measured using magnetic resonance spectroscopy (MRS) in a high resolution (3T) scanner. Subjects also completed a battery of neurocognitive tests. Patients had repeat testing after 6 months of CPAP. There were significant differences at baseline in frontal N-acetylaspartate/choline (NAA/Cho) ratios (patients [mean (SD)] 4.56 [0.41], controls 4.92 [0.44], P = 0.001), and in hippocampal choline/creatine (Cho/Cr) ratios (0.38 [0.04] vs 0.41 [0.04], P = 0.006), (both ANCOVA, with age and premorbid IQ as covariates). No longitudinal changes were seen with treatment (n = 27, paired t tests), however the hippocampal differences were no longer significant at 6 months, and frontal NAA/Cr ratios were now also significantly different (patients 1.55 [0.13] vs control 1.65 [0.18] P = 0.01). No significant correlations were found between spectroscopy results and neurocognitive test results, but significant negative correlations were seen between arousal index and frontal NAA/Cho (r = -0.39, corrected P = 0.033) and between % total sleep time at SpO(2) < 90% and hippocampal Cho/Cr (r = -0.40, corrected P = 0.01). CONCLUSIONS: OSA patients have brain metabolite changes detected by MRS, suggestive of decreased frontal lobe neuronal viability and integrity, and decreased hippocampal membrane turnover. These regions have previously been shown to have no gross structural lesions using VBM. Little change was seen with treatment with CPAP for 6 months. No correlation of metabolite concentrations was seen with results on neurocognitive tests, but there were significant negative correlations with OSA severity as measured by severity of nocturnal hypoxemia. CITATION: O'Donoghue FJ; Wellard RM; Rochford PD; Dawson A; Barnes M; Ruehland WR; Jackson ML; Howard ME; Pierce RJ; Jackson GD. Magnetic resonance spectroscopy and neurocognitive dysfunction in obstructive sleep apnea before and after CPAP treatment.

Ventilation distribution in rats : part 2 – a comparison of electrical impedance tomography and hyperpolarised helium magnetic resonance imaging

Background: Hyperpolarised helium MRI (He3 MRI) is a new technique that enables imaging of the air distribution within the lungs. This allows accurate determination of the ventilation distribution in vivo. The technique has the disadvantages of requiring an expensive helium isotope, complex apparatus and moving the patient to a compatible MRI scanner. Electrical impedance tomography (EIT) a non-invasive bedside technique that allows constant monitoring of lung impedance, which is dependent on changes in air space capacity in the lung. We have used He3MRI measurements of ventilation distribution as the gold standard for assessment of EIT. Methods: Seven rats were ventilated in supine, prone, left and right lateral position with 70% helium/30% oxygen for EIT measurements and pure helium for He3 MRI. The same ventilator and settings were used for both measurements. Image dimensions, geometric centre and global in homogeneity index were calculated. Results: EIT images were smaller and of lower resolution and contained less anatomical detail than those from He3 MRI. However, both methods could measure positional induced changes in lung ventilation, as assessed by the geometric centre. The global in homogeneity index were comparable between the techniques. Conclusion: EIT is a suitable technique for monitoring ventilation distribution and inhomgeneity as assessed by comparison with He3 MRI.

Spinal inflammation in the absence of sacroiliac joint inflammation on magnetic resonance imaging in patients with active nonradiographic axial spondyloarthritis

Objective: To evaluate the presence of spinal inflammation with and without sacroiliac (SI) joint inflammation on magnetic resonance imaging (MRI) in patients with active nonradiographic axial spondyloarthritis (SpA), and to compare the disease characteristics of these subgroups. Methods: ABILITY-1 is a multicenter, randomized, controlled trial of adalimumab versus placebo in patients with nonradiographic axial SpA classified using the Assessment of SpondyloArthritis international Society axial SpA criteria. Baseline MRIs were centrally scored independently by 2 readers using the Spondyloarthritis Research Consortium of Canada (SPARCC) method for the SI joints and the SPARCC 6-discovertebral unit method for the spine. Positive evidence of inflammation on MRI was defined as a SPARCC score of >2 for either the SI joints or the spine. Results: Among patients with baseline SPARCC scores, 40% had an SI joint score of >2 and 52% had a spine score of >2. Forty-nine percent of patients with baseline SI joint scores of <2, and 58% of those with baseline SI joint scores of >2, had a spine score of >2. Comparison of baseline disease characteristics by baseline SI joint and spine scores showed that a greater proportion of patients in the subgroup with a baseline SPARCC score of >2 for both SI joints and spine were male, and patients with spine and SI joint scores of <2 were younger and had shorter symptom duration. SPARCC spine scores correlated with baseline symptom duration, and SI joint scores correlated negatively with the baseline Bath Ankylosing Spondylitis Disease Activity Index, but neither correlated with the baseline Ankylosing Spondylitis Disease Activity Score, total back pain, the patient's global assessment of disease activity, the Bath Ankylosing Spondylitis Functional Index, morning stiffness, nocturnal pain, or C-reactive protein level. Conclusion: Assessment by experienced readers showed that spinal inflammation on MRI might be observed in half of patients with nonradiographic axial SpA without SI joint inflammation.

Normalized wall index specific and MRI-based stress analysis of atherosclerotic carotid plaques: A study comparing acutely symptomatic and asymptomatic patients

Background: Biomechanical stresses play an important role in determining plaque stability. Quantification of these simulated stresses can be potentially used to assess plaque vulnerability and differentiate different patient groups. Methods and Results: 54 asymptomatic and 45 acutely symptomatic patients underwent in vivo multicontrast magnetic resonance imaging (MRI) of the carotid arteries. Plaque geometry used for finite element analysis was derived from in vivo MRI at the sites of maximum and minimum plaque burden. In total, 198 slices were used for the computational simulations. A pre-shrink technique was used to refine the simulation. Maximum principle stress at the vulnerable plaque sites (ie, critical stress) was extracted for the selected slices and a comparison was performed between the 2 groups. Critical stress in the slice with maximum plaque burden is significantly higher in acutely symptomatic patients as compared to asymptomatic patients (median, inter quartile range: 198.0 kPa (119.8-359.0 kPa) vs 138.4 kPa (83.8-242.6 kPa), P=0.04). No significant difference was found in the slice with minimum plaque burden between the 2 groups (196.7 kPa (133.3-282.7 kPa) vs 182.4 kPa (117.2-310.6 kPa), P=0.82). Conclusions: Acutely symptomatic carotid plaques have significantly high biomechanical stresses than asymptomatic plaques. This might be potentially useful for establishing a biomechanical risk stratification criteria based on plaque burden in future studies.

Association between biomechanical structural stresses of atherosclerotic carotid plaques and subsequent ischaemic cerebrovascular events - A longitudinal in vivo magnetic resonance imaging-based finite element study

Background: High-resolution magnetic resonance (MR) imaging has been used for MR imaging-based structural stress analysis of atherosclerotic plaques. The biomechanical stress profile of stable plaques has been observed to differ from that of unstable plaques; however, the role that structural stresses play in determining plaque vulnerability remains speculative. Methods: A total of 61 patients with previous history of symptomatic carotid artery disease underwent carotid plaque MR imaging. Plaque components of the index artery such as fibrous tissue, lipid content and plaque haemorrhage (PH) were delineated and used for finite element analysis-based maximum structural stress (M-C Stress) quantification. These patients were followed up for 2 years. The clinical end point was occurrence of an ischaemic cerebrovascular event. The association of the time to the clinical end point with plaque morphology and M-C Stress was analysed. Results: During a median follow-up duration of 514 days, 20% of patients (n=12) experienced an ischaemic event in the territory of the index carotid artery. Cox regression analysis indicated that M-C Stress (hazard ratio (HR): 12.98 (95% confidence interval (CI): 1.32-26.67, pZ0.02), fibrous cap (FC) disruption (HR: 7.39 (95% CI: 1.61e33.82), p Z 0.009) and PH (HR: 5.85 (95% CI: 1.27e26.77), p Z 0.02) are associated with the development of subsequent cerebrovascular events. Plaques associated with future events had higher M-C Stress than those which had remained asymptomatic (median (interquartile range, IQR): 330 kPa (229e494) vs. 254 kPa (166-290), p Z0.04). Conclusions: High biomechanical structural stresses, in addition to FC rupture and PH, are associated with subsequent cerebrovascular events.

Lens shape and refractive index distribution in type 1 diabetes

Purpose: To compare lens dimensions and refractive index distributions in type 1 diabetes and age-matched control groups. Methods: There were 17 participants with type 1 diabetes, consisting of two subgroups (7 young [23 ± 4 years] and 10 older [54 ± 4 years] participants), with 23 controls (13 young, 24 ± 4 years; 10 older, 55 ± 4 years). For each participant, one eye was tested with relaxed accommodation. A 3T clinical magnetic resonance imaging scanner was used to image the eye, employing a multiple spin echo (MSE) sequence to determine lens dimensions and refractive index profiles along the equatorial and axial directions. Results: The diabetes group had significantly smaller lens equatorial diameters and larger lens axial thicknesses than the control group (diameter mean ± 95% confidence interval [CI]: diabetes group 8.65 ± 0.26 mm, control group 9.42 ± 0.18 mm; axial thickness: diabetes group 4.33 ± 0.30 mm, control group 3.80 ± 0.14 mm). These differences were also significant within each age group. The older group had significantly greater axial thickness than the young group (older group 4.35 ± 0.26 mm, young group 3.70 ± 0.25 mm). Center refractive indices of diabetes and control groups were not significantly different. There were some statistically significant differences between the refractive index fitting parameters of young and older groups, but not between diabetes and control groups of the same age. Conclusions: Smaller lens diameters occurred in the diabetes groups than in the age-matched control groups. Differences in refractive index distribution between persons with and without diabetes are too small to have important effects on instruments measuring axial thickness.