各向异性与非线性介质中弹性波场传播特征研究

The topic of this study is about the propagation features of elastic waves in the anisotropic and nonlinear media by numerical methods with high accuracy and stability. The main achievements of this paper are as followings: Firstly, basing on the third order elastic energy formula, principle of energy conservation and circumvolved matrix method, we firstly reported the equations of non-linear elastic waves with two dimensions and three components in VTI media. Secondly, several conclusions about some numerical methods have been obtained in this paper. Namely, the minimum suitable sample stepth in space is about 1/8-1/12 of the main wavelength in order to distinctly reduce the numerical dispersion resulted from the numerical mehtod, at the same time, the higher order conventional finite difference (CFD) schemes will give little contribution to avoid the numerical solutions error accumulating with time. To get the similar accuracy with the fourth order center finite difference method, the half truncation length of SFFT should be no less than 7. The FDFCT method can present with the numerical solutions without obvious dispersion when the paprameters of FCT is suitable (we think they should be in the scope from 0.0001 to 0.07). Fortunately, the NADM method not only can reported us with the higher order accuracy solutions (higher than that of the fourth order finite difference method and lower than that of the sixth order finite difference method), but also can distinctly reduce the numerical dispersion. Thirdly, basing on the numerial and theoretical analysis, we reported such nonlinear response accumulating with time as waveform aberration, harmonic generation and resonant peak shift shown by the propagation of one- and two-dimensional non-linear elasticwaves in this paper. And then, we drew the conclusion that these nonlinear responses are controlled by the product between nonlinear strength (SN) and the amplitude of the source. At last, the modified FDFCT numerical method presented by this paper is used to model the two-dimensional non-linear elastic waves propagating in VTI media. Subsequently, the wavelet analysis and polarization are adopted to investigate and understand the numerical results. And then, we found the following principles (attention: the nonlinear strength presented by this paper is weak, the thickness of the -nonlinear media is thin (200m), the initial energy of the source is weak and the anisotropy of the media is weak too): The non-linear response shown by the elastic waves in VTI media is anisotropic too; The instantaneous main frequency sections of seismic records resulted from the media with a non-linear layer have about 1/4 to 1/2 changes of the initial main frequency of source with that resulted from the media without non-linear layer; The responses shown by the elasic waves about the anisotropy and nonlinearity have obvious mutual reformation, namely, the non-linear response will be stronger in some directions because of the anisotropy and the anisotropic strength shown by the elastic waves will be stronger when the media is nonlinear.

Observational Evidence for Mode Coupling in the Chromospheric Network

Oscillations in network bright points (NBPs) are studied at a variety of chromospheric heights. In particular, the three-dimensional variation of NBP oscillations is studied using image segmentation and cross-correlation analysis between images taken in light of Ca II K3, Ha core, Mg I b2, and Mg I b1-0.4 Å. Wavelet analysis is used to isolate wave packets in time and to search for height-dependent time delays that result from upward- or downward-directed traveling waves. In each NBP studied, we find evidence for kink-mode waves (1.3, 1.9 mHz), traveling up through the chromosphere and coupling with sausage-mode waves (2.6, 3.8 mHz). This provides a means for depositing energy in the upper chromosphere. We also find evidence for other upward- and downward-propagating waves in the 1.3-4.6 mHz range. Some oscillations do not correspond to traveling waves, and we attribute these to waves generated in neighboring regions.

White-light oscillations during a flare on II Peg

We analyse the intensity oscillations observed in the gradual phase of a white-light flare on the RS CV n binary II Peg. Fast Fourier Transform power spectra and Wavelet analysis reveal a period of 220 s. The reliability of the oscillation is tested using several criteria. Oscillating coronal loop models are used to derive physical parameters such as temperature, electron density and magnetic field strength associated with the coronal loop. The derived parameters are consistent with the near-simultaneous X-ray observations of the flare. There is no evidence for oscillations in the quiescent state of the binary.

Propagating waves and magnetohydrodynamic mode coupling in the quiet-sun network

High-cadence multiwavelength optical observations were taken with the Dunn Solar Telescope at the National Solar Observatory, Sacramento Peak, accompanied by Advanced Stokes Polarimeter vector magnetograms. A total of 11 network bright points (NBPs) have been studied at different atmospheric heights using images taken in wave bands centered on Mg I b(1) - 0.4 Angstrom, Halpha, and Ca II K-3. Wavelet analysis was used to study wave packets and identify traveling magnetohydrodynamic waves. Wave speeds were estimated through the temporal cross-correlation of signals, in selected frequency bands of wavelet power, in each wavelength. Four mode-coupling cases were identified, one in each of four of the NBPs, and the variation of the associated Fourier power with height was studied. Three of the detected mode-coupling, transverse-mode frequencies were observed in the 1.2-1.6 mHz range (mean NBP apparent flux density magnitudes over 99-111 Mx cm(-2)), with the final case showing 2.0-2.2 mHz (with 142 Mx cm(-2)). Following this, longitudinal-mode frequencies were detected in the range 2.6-3.2 mHz for three of our cases, with 3.9-4.1 mHz for the remaining case. After mode coupling, two cases displayed a decrease in longitudinal-mode Fourier power in the higher chromosphere.

Eclipse observations of high-frequency oscillations in active region coronal loops

One of the mechanisms proposed for heating the corona above solar active regions is the damping of magnetohydrodynamic (MHD) waves. Continuing on previous work, we provide observational evidence for the existence of high-frequency MHD waves in coronal loops observed during the August 1999 total solar eclipse. A wavelet analysis is used to identify twenty 4 x 4 arcsec(2) areas showing intensity oscillations. All detections lie in the frequency range 0.15 - 0.25 Hz (7 - 4 s), last for at least 3 periods at a confidence level of more than 99% and arise just outside known coronal loops. This leads us to suggest that they occur in low emission-measure or different temperature loops associated with the active region.

Twisting flux tubes as a cause of micro-flaring activity

High-cadence optical observations of an H-alpha blue-wing bright point near solar AR NOAA 10794 are presented. The data were obtained with the Dunn Solar Telescope at the National Solar Observatory/Sacramento Peak using a newly developed camera system, the RAPID DUAL IMAGER. Wavelet analysis is undertaken to search for intensity-related oscillatory signatures, and periodicities ranging from 15 to 370 s are found with significance levels exceeding 95%. During two separate microflaring events, oscillation sites surrounding the bright point are observed to twist. We relate the twisting of the oscillation sites to the twisting of physical flux tubes, thus giving rise to reconnection phenomena. We derive an average twist velocity of 8.1 km s(-1) and detect a peak in the emitted flux between twist angles of 180 degrees and 230 degrees.

High-frequency oscillations in a solar active region observed with the RAPID DUAL IMAGER

High-cadence, synchronized, multiwavelength optical observations of a solar active region (NOAA 10794) are presented. The data were obtained with the Dunn Solar Telescope at the National Solar Observatory/Sacramento Peak using a newly developed camera system: the rapid dual imager. Wavelet analysis is undertaken to search for intensity related oscillatory signatures, and periodicities ranging from 20 to 370 s are found with significance levels exceeding 95%. Observations in the H-α blue wing show more penumbral oscillatory phenomena when compared to simultaneous G-band observations. The H-α oscillations are interpreted as the signatures of plasma motions with a mean velocity of 20 km s-1. The strong oscillatory power over H-α blue-wing and G-band penumbral bright grains is an indication of the Evershed flow with frequencies higher than previously reported.

Alfvén Waves in the Lower Solar Atmosphere

The flow of energy through the solar atmosphere and the heating of the Sun's outer regions are still not understood. Here, we report the detection of oscillatory phenomena associated with a large bright-point group that is 430,000 square kilometers in area and located near the solar disk center. Wavelet analysis reveals full-width half-maximum oscillations with periodicities ranging from 126 to 700 seconds originating above the bright point and significance levels exceeding 99%. These oscillations, 2.6 kilometers per second in amplitude, are coupled with chromospheric line-of-sight Doppler velocities with an average blue shift of 23 kilometers per second. A lack of cospatial intensity oscillations and transversal displacements rules out the presence of magneto-acoustic wave modes. The oscillations are a signature of Alfvén waves produced by a torsional twist of ±22 degrees. A phase shift of 180 degrees across the diameter of the bright point suggests that these torsional Alfvén oscillations are induced globally throughout the entire brightening. The energy flux associated with this wave mode is sufficient to heat the solar corona.

Chromospheric velocities of a C-class flare

Aims. We use high spatial and temporal resolution observations from the Swedish Solar Telescope to study the chromospheric velocities of a C-class flare originating from active region NOAA 10969.
Methods. A time-distance analysis is employed to estimate directional velocity components in Hα and Ca ii K image sequences. Also, imaging spectroscopy has allowed us to determine flare-induced line-of-sight velocities. A wavelet analysis is used to analyse the periodic nature of associated flare bursts.
Results. Time-distance analysis reveals velocities as high as 64 km s-1 along the flare ribbon and 15 km s-1 perpendicular to it. The velocities are very similar in both the Hα and Ca ii K time series. Line-of-sight Hα velocities are red-shifted with values up to 17 km s-1. The high spatial and temporal resolution of the observations have allowed us to detect velocities significantly higher than those found in earlier studies. Flare bursts with a periodicity of ≈60 s are also detected. These bursts are similar to the quasi-periodic oscillations observed at hard X-ray and radio wavelength data.
Conclusions. Some of the highest velocities detected in the solar atmosphere are presented. Line-of-sight velocity maps show considerable mixing of both the magnitude and direction of velocities along the flare path. A change in direction of the velocities at the flare kernel has also been detected which may be a signature of chromospheric evaporation.

Search for Rapid Changes in the Visible-Light Corona during the 21 June 2001 Total Solar Eclipse

Some 8000 images obtained with the Solar Eclipse Coronal Imaging System (SECIS) fast-frame CCD camera instrument located at Lusaka, Zambia, during the total eclipse of 21 June 2001 have been analysed to search for short-period oscillations in intensity that could be a signature of solar coronal heating mechanisms by MHD wave dissipation. Images were taken in white-light and Fe xiv green-line (5303 ) channels over 205 seconds (frame rate 39 s(-1)), approximately the length of eclipse totality at this location, with a pixel size of four arcseconds square. The data are of considerably better quality than those that we obtained during the 11 August 1999 total eclipse (Rudawy et al.: Astron. Astrophys. 416, 1179, 2004), in that the images are much better exposed and enhancements in the drive system of the heliostat used gave a much improved image stability. Classical Fourier and wavelet techniques have been used to analyse the emission at 29 518 locations, of which 10 714 had emission at reasonably high levels, searching for periodic fluctuations with periods in the range 0.1 -aEuro parts per thousand 17 seconds (frequencies 0.06 -aEuro parts per thousand 10 Hz). While a number of possible periodicities were apparent in the wavelet analysis, none of the spatially and time-limited periodicities in the local brightness curves was found to be physically important. This implies that the pervasive Alfv,n wave-like phenomena (Tomczyk et al.: Science 317, 1192, 2007) using polarimetric observations with the Coronal Multi-Channel Polarimeter (CoMP) instrument do not give rise to significant oscillatory intensity fluctuations.

Influence of the Pacific Decadal Oscillation, El Niño-Southern Oscillation and solar forcing on climate and primary productivity changes in the northeast Pacific

Evidence of 11-year Schwabe solar sunspot cycles, El Niño-Southern Oscillation (ENSO) and the Pacific Decadal Oscillation (PDO) were detected in an annual record of diatomaceous laminated sediments from anoxic Effingham Inlet, Vancouver Island, British Columbia. Radiometric dating and counting of annual varves dates the sediments from AD 1947-1993. Intact sediment slabs were X-rayed for sediment structure (lamina thickness and composition based on gray-scale), and subsamples were examined for diatom abundances and for grain size. Wavelet analysis reveals the presence of ~2-3, ~4.5, ~7 and ~9-12-year cycles in the diatom record and an w11e13 year record in the sedimentary varve thickness record. These cycle lengths suggest that both ENSO and the sunspot cycle had an influence on primary productivity and sedimentation patterns. Sediment grain size could not be correlated to the sunspot cycle although a peak in the grain size data centered around the mid-1970s may be related to the 1976-1977 Pacific climate shift, which occurred when the PDO index shifted from negative (cool conditions) to positive (warm conditions). Additional evidence of the PDO regime shift is found in wavelet and cross-wavelet results for Skeletonema costatum, a weakly silicified variant of S. costatum, annual precipitation and April to June precipitation. Higher spring (April/May) values of the North Pacific High pressure index during sunspot minima suggest that during this time, increased cloud cover and concomitant suppression of the Aleutian Low (AL) pressure system led to strengthened coastal upwelling and enhanced diatom production earlier in the year. These results suggest that the 11-year solar cycle, amplified by cloud cover and upwelling changes, as well as ENSO, exert significant influence on marine primary productivity in the northeast Pacific. The expression of these cyclic phenomena in the sedimentary record were in turn modulated by the phase of PDO, as indicated by the change in period of ENSO and suppression of the solar signal in the record after the 1976-1977 regime shift. © 2013 Elsevier Ltd and INQUA. All rights reserved.

Centennial-scale climate change in Ireland during the Holocene

We examine mid- to late Holocene centennial-scale climate variability in Ireland using proxy data from peatlands, lakes and a speleothem. A high degree of between-record variability is apparent in the proxy data and significant chronological uncertainties are present. However, tephra layers provide a robust tool for correlation and improve the chronological precision of the records. Although we can find no statistically significant coherence in the dataset as a whole, a selection of high-quality peatland water table reconstructions co-vary more than would be expected by chance alone. A locally weighted regression model with bootstrapping can be used to construct a ‘best-estimate’ palaeoclimatic reconstruction from these datasets. Visual comparison and cross-wavelet analysis of peatland water table compilations from Ireland and Northern Britain show that there are some periods of coherence between these records. Some terrestrial palaeoclimatic changes in Ireland appear to coincide with changes in the North Atlantic thermohaline circulation and solar activity. However, these relationships are inconsistent and may be obscured by chronological uncertainties. We conclude by suggesting an agenda for future Holocene climate research in Ireland. ©2013 Elsevier B.V. All rights reserved.

Ultrasound entropy may be a new non-invasive measure of pre-clinical vascular damage in young hypertensive patients

Background: The identification of pre-clinical microvascular damage in hypertension by non-invasive techniques has proved frustrating for clinicians. This proof of concept study investigated whether entropy, a novel summary measure for characterizing blood velocity waveforms, is altered in participants with hypertension and may therefore be useful in risk stratification.

Methods: Doppler ultrasound waveforms were obtained from the carotid and retrobulbar circulation in 42 participants with uncomplicated grade 1 hypertension (mean systolic/diastolic blood pressure (BP) 142/92 mmHg), and 26 healthy controls (mean systolic/diastolic BP 116/69 mmHg). Mean wavelet entropy was derived from flow-velocity data and compared with traditional haemodynamic measures of microvascular function, namely the resistive and pulsatility indices.

Results: Entropy, was significantly higher in control participants in the central retinal artery (CRA) (differential mean 0.11 (standard error 0.05 cms(-1)), CI 0.009 to 0.219, p 0.017) and ophthalmic artery (0.12 (0.05), CI 0.004 to 0.215, p 0.04). In comparison, the resistive index (0.12 (0.05), CI 0.005 to 0.226, p 0.029) and pulsatility index (0.96 (0.38), CI 0.19 to 1.72, p 0.015) showed significant differences between groups in the CRA alone. Regression analysis indicated that entropy was significantly influenced by age and systolic blood pressure (r values 0.4-0.6). None of the measures were significantly altered in the larger conduit vessel.

Conclusion: This is the first application of entropy to human blood velocity waveform analysis and shows that this new technique has the ability to discriminate health from early hypertensive disease, thereby promoting the early identification of cardiovascular disease in a young hypertensive population.

Essays on energy derivatives pricing and financial risk management

This thesis consists of an introductory chapter (essay I) and five more empirical essays on electricity markets and CO2 spot price behaviour, derivatives pricing analysis and hedging. Essay I presents the structure of the thesis and electricity markets functioning and characteristics, as well as the type of products traded, to be analyzed on the following essays. In the second essay we conduct an empirical study on co-movements in electricity markets resorting to wavelet analysis, discussing long-term dynamics and markets integration. Essay three is about hedging performance and multiscale relationships in the German electricity spot and futures markets, also using wavelet analysis. We concentrate the investigation on the relationship between coherence evolution and hedge ratio analysis, on a time-frequency-scale approach, between spot and futures which conditions the effectiveness of the hedging strategy. Essays four, five and six are interrelated between them and with the other two previous essays given the nature of the commodity analyzed, CO2 emission allowances, traded in electricity markets. Relationships between electricity prices, primary energy fuel prices and carbon dioxide permits are analyzed on essay four. The efficiency of the European market for allowances is examined taking into account markets heterogeneity. Essay five analyzes stylized statistical properties of the recent traded asset CO2 emission allowances, for spot and futures returns, examining also the relation linking convenience yield and risk premium, for the German European Energy Exchange (EEX) between October 2005 and October 2009. The study was conducted through empirical estimations of CO2 allowances risk premium, convenience yield, and their relation. Future prices from an ex-post perspective are examined to show evidence for significant negative risk premium, or else a positive forward premium. Finally, essay six analyzes emission allowances futures hedging effectiveness, providing evidence for utility gains increases with investor’s preference over risk. Deregulation of electricity markets has led to higher uncertainty in electricity prices and by presenting these essays we try to shed new lights about structuring, pricing and hedging in this type of markets.

Fractional dynamics in the trajectory control of redundant manipulators

Under the pseudoinverse control, robots with kinematical redundancy exhibit an undesirable chaotic joint motion which leads to an erratic behavior. This paper studies the complexity of fractional dynamics of the chaotic response. Fourier and wavelet analysis provides a deeper insight, helpful to know better the lack of repeatability problem of redundant manipulators. This perspective for the study of the chaotic phenomena will permit the development of superior trajectory control algorithms.