滑动控制理论在海洋机器人动态定位系统中的应用

本文论述了一种新的鲁棒性控制方法——滑动控制理论.并将该方法应用于海洋机器人的动态定位系统设计.仿真结果表明该方法是相当成功和有效的,对系统的参数变化有强鲁棒性.

电离层行进式扰动的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.

空间电流体系地磁反演及理论研究

The space currents definitely take effects on electromagnetic environment and also are scientific highlight in the space research. Space currents as a momentum and energy provider to Geospace Storm, disturb the varied part of geomagnetic field, distort magnetospheric configuration and furthermore take control of the coupling between magnetosphere and ionosphere. Due to both academic and commercial objectives above, we carry on geomagnetic inverse and theoretical studies about the space currents by using geomagnetic data from INTERMAGNET. At first, we apply a method of Natural Orthogonal Components (NOC) to decomposition the solar daily variation, especially for (solar quiet variation). NOC is just one of eign mode analysis, the most advantage of this method is that the basic functions (BFs) were not previously designated, but naturally came from the original data so that there are several BFs usually corresponding to the process really happened and have more physical meaning than the traditional spectrum analysis with the fixed BFs like Fourier trigonometric functions. The first two eign modes are corresponding to the and daily variation and their amplitudes both have the seasonal and day-to-day trend, that will be useful for evaluating geomagnetic activity indices. Because of the too strict constraints of orthogonality, we try to extend orthogonal contraints to the non-orthogonal ones in order to give more suitable and appropriate decomposition of the real processes when the most components did not satisfy orthogonality. We introduce a mapping matrix which can transform the real physical space to a new mathematical space, after that process, the modified components which associated with the physical processes have satisfied the orthogonality in the new mathematical space, furthermore, we can continue to use the NOC decomposition in the new mathematical space, and then all the components inversely transform back to original physical space, so that we would have finished the non-orthogonal decomposition which more generally in the real world. Secondly, geomagnetic inverse of the ring current’s topology is conducted. Configurational changes of the ring current in the magnetosphere lead to different patterns of disturbed ground field, so that the global configuration of ring current can be inferred from its geomagnetic perturbations. We took advantages of worldwide geomagnetic observatories network to investigate the disturbed geomagnetic field which produced by ring current. It was found that the ring current was not always centered at geomagnetic equator, and significantly deviated off the equator during several intense magnetic storms. The deviation owing to the tilting and latitudinal shifting of the ring current with respect to the earth’s dipole can be estimated from global geomagnetic survey. Furthermore those two configurational factors which gave a quantitative description of the ring current configuration, will be helpful to improve the Dst calibration and understand the dependence of ring current’s configuration on the plasma sheet location relative to the equator when magnetotail field warped. Thirdly, the energization and physical acceleration process of ring current during magnetic storm has been proposed. When IMF Bz component increase, the enhanced convection electric field drive the plasma injection into the inner magnetosphere. During the transport process, a dynamic heating is happened which make the particles more ‘hot’ when the injection is more deeply inward. The energy gradient along the injection path is equivalent to a kind of force, which resist the plasma more earthward injection, as a diamagnetic effect of the magnetosphere anti and repellent action to the exotically injected plasma. The acceleration efficiency has a power law form. We use analytical way to quantitatively describe the dynamical process by introducing a physical parameter: energization index, which will be useful to understand how the particle is heated. At the end, we give a scheme of how to get the from storm time geomagnetic data. During intense magnetic storms, the lognormal trend of geomagnetic Dst decreases depend on the heating dynamic of magnetosphere controlling ring current. The descending pattern of main phase is governed by the magnetospheric configuration, which can be describled by the energization index. The amplitude of Dst correlated with convection electric field or south component of the solar wind. Finally, the Dst index is predicted by upstream solar wind parameter. As we known space weather have posed many chanllenges and impacts on techinal system, the geomagnetic index for evaluating the activity space weather. We review the most popular Dst prediction method and repeat the Dst forecasting model works. A concise and convnient Key Points model of the polar region is also introduced to space weather. In summary, this paper contains some new quantitative and physical description of the space currents with special focus on the ring current. Whatever we do is just to gain a better understanding of the natural world, particularly the space environment around Earth through analytical deduction, algorithm designing and physical analysis, to quantitative interpretation. Applications of theoretical physics in conjunction with data analysis help us to understand the basic physical process govering the universe.

磁层-电离层耦合现象的分析研究

The magnetosphere-ionosphere coupling is mainly manifested by the trans- porting processes of energy into the ionosphere , the energy is carried by solar wind and firstly accumulate at the magnetosphere, and the coupling processes also significantly include the interaction between the magnetosphere and ionosphere for mass and energy. At the quiet condition, energy is delivered by the large-scale convection of the geomagnetic field; the huge energy from solar wind bulk will be injected into and consumed at the near magnetosphere and ionosphere by the geomagnetic storm and substorm activities. Aurorae and FACs (Field-aligned currents) are the important phenomena in the coupling processes. In the present work, firstly, we analyze the activity characteristics of auroral precipitating particle, secondly, we study the distribution characters of large-scale field aligned currents (LS FACs) at storm-time using the observations from different satellites at different altitudes. Finally, we investigate the evolution of the geomagnetic field configuration at the nightside sector on the onset of the expansion phase in a substorm event, the substorm event happened at 0430UT to 0630UT on 8th Nov. 2004. The main results as follows： At the first, the data of the estimated power input (EPI) of auroral particles from NOAA/POES (Polar orbiting environmental satellite) for some 30 years have been analyzed. The variation tendencies of the EPI generally coincide with aa, AE and Dst indices. The annual variation of EPI shows equinox peaks and an asymmetric-activity with a higher peak in the winter-hemisphere than in the summer-hemisphere. The diurnal UT variations are different from north and south hemisphere: for north hemisphere, the peak appears at 1200UT, and the relative deviation is 22% to the daily average of the north hemisphere. For south hemisphere, the maximal deviation is 22% at 2000UT. So the diurnal variation of EPI is more dominant than the annual variation which maximal deviation is 3% to 12% for different seasons. Studies on correlations of the hourly average of EPI, Pa, with AE and Dst indices show a correlation coefficient r=0.74 of Pa and AE, and r=-0.55 of Pa and Dst. The hourly EPIs for north and south polar regions, NPa and SPa, show a north-south asymmetry with a higher correlation of SPa and AE (or Dst). Time delays of EPI with respect to magnetic indices are examined, the maximum correlation coefficient of Pa with AE (r=0.78) occurs when the time delay =0, suggesting a synchronous activity of auroral electrojet and auroral precipitating particles, while =1-2h, the correlation coefficient of Pa with Dst is maximum (r=0.57), suggesting that the activity of auroral particle precipitating may influence the ring current on some extent. Sencondly, we use the high-resolution magnetic field vector data of the CHAMP satellite to investigate the distribution of large-scale FACs during the great magnetic storm on 7th to 8th Nov. 2004. The results show that, whether in the northern or southern hemisphere, the number and density of large-scale FACs during the main-phase are more and bigger than these during the recover-phase, and the number of large-scale FACs in morning sector obviously is more than that in afternoon sector. In terms of the magnetic indices, we find that large-scale FACs in morning sector significantly affected by the substorm activities, while in afternoon sector the large-scale FACs mainly indicate the fluctuations of the ring-current in storm time. Accordingly to the former studies, similarly, we find that in the morning sector, the scale of the large-scale FACs move to the high-latitude region, and in the afternoon sector, large-scale FACs distinctly expand to the low-latitude region. During the time periods that the NOAA/POES auroral precipitating particle power data temporally correspond to the large-scale FACs, the more the power of auroral particle is, the more and bigger the number and density of FACs are. At the same time, we use the magnetic field vector data of POLAR obtain a good form of region 1, region 2, and three pieces of cusp FACs during a single transit at 1930UT-2006UT on 07th. And the characteristics of simultaneous electric field and energy particles observations on Polar are coincide with the five FACs pieces. Finally, by means of the observation of Cluster 4 and Goes 10、 Goes 12, we analyze the evolution process of the change of the magnetic field configuration at night sector at the expansion phase of a substorm event which happened during 0430UT to 0630UT on 8th Nov. 2004, we find that the times of the beginning of the polarizations of magnetic field are observed from Goes 10 to Goes 12 then to Cluster 4. So, at the synchronous orbit ( 6.6 RE) to 10RE distance scale of the neutral sheet, the current disruption spread tailward. Simultaneously, the strengthen of the FACs deduced from these satellites’ magnetic field observations are almost consistent with the times of polarizations, as well as the high energy particles injection and the electric field dominant variation. The onset times determined by the magnetic field polarizations from these satellites are all ahead of the onset time that confirmed from the auroral electrojet indices. So, these characters of different observations can be used as the criterions to determine the onset time for the substorms of such type as we studied.

二氧化碳在孔隙咸水层中多相流动数值模拟研究

There has been a growing concern about the use of fossil fuels and its adverse effects on the atmospheric greenhouse and ecological environment. A reduction in the release rate of CO2 into the atmosphere poses a major challenge to the land ecology of China. The most promising way of achieving CO2 reduction is to dispose of CO2 in deep saline aquifers. Deep aquifers have a large potential for CO2 sequestration in geological medium in terms of volume and duration. Through the numerical simulation of multiphase flow in a porous media, the transformation and motion of CO2 in saline aquifers has been implemented under various temperature and hydrostatic pressure conditions, which plays an important role to the assessment of the reliability and safety of CO2 geological storage. As expected, the calculated results can provide meaningful and scientific information for management purposes. The key problem to the numerical simulation of multiphase flow in a porous media is to accurately capture the mass interface and to deal with the geological heterogeneity. In this study, the updated CE/SE (Space and time conservation element and solution element) method has been proposed, and the Hybrid Particle Level Set method (HPLS) has extended for multiphase flows in porous medium, which can accurately trace the transformation of the mass interface. The benchmark problems have been applied to evaluate and validate the proposed method. In this study, the reliability of CO2 storage in saline aquifers in Daqingzi oil field in Sunlong basin has been discussed. The simulation code developed in this study takes into account the state for CO2 covering the triple point temperature and pressure to the supercritical region. The geological heterogeneity has been implemented, using the well known geostatistical model (GSLIB) on the base of the hard data. The 2D and 3D model have been set up to simulate the CO2 multiphase flow in the porous saline aquifer, applying the CE/SE method and the HPLS method .The main contents and results are summarized as followings. (1) The 2D CE/SE method with first and second –order accuracy has been extended to simulate the multiphase flow in porous medium, which takes into account the contribution of source and sink in the momentum equation. The 3D CE/SE method with the first accuracy has been deduced. The accuracy and efficiency of the proposed CE/SE method have been investigated, using the benchmark problems. (2) The hybrid particle level set method has been made appropriate and extended for capturing the mass interface of multiphase flows in porous media, and the numerical method for level set function calculated has been formulated. (3) The closed equations for multiphase flow in porous medium has been developed, adept to both the Darcy flow and non-Darcy flow, getting over the limitation of Reynolds number to the calculation. It is found that Darcy number has a decisive influence on pressure as well as velocity given the Darcy number. (4) The new Euler scheme for numerical simulations of multiphase flows in porous medium has been proposed, which is efficient and can accurately capture the mass interface. The artificial compressibility method has been used to couple the velocities and pressure. It is found that the Darcy number has determinant effects on the numerical convergence and stability. In terms of the different Darcy numbers, the coefficient of artificial compressibility and the time step have been obtained. (5) The time scale of the critical instability for critical CO2 in the saline aquifer has been found, which is comparable with that of completely CO2 dissolved saline aquifer. (6) The concept model for CO2 multiphase flows in the saline aquifer has been configured, based on the temperature, pressure, porosity as well as permeability of the field site .Numerical simulation of CO2 hydrodynamic trapping in saline aquifers has been performed, applying the proposed CE/SE method. The state for CO2 has been employed to take into account realistic reservoir conditions for CO2 geological sequestration. The geological heterogeneity has been sufficiently treated , using the geostatistical model. (7) It is found that the Rayleigh-Taylor instability phenomenon, which is associated with the penetration of saline fluid into CO2 fluid in the direction of gravity, has been observed in CO2 multiphase flows in the saline aquifer. Development of a mushroom-type spike is a strong indication of the formation of Kelvin-Helmholtz instability due to the developed short wavelength perturbations present along the interface and parallel to the bulk flow. Additional key findings: the geological heterogeneity can distort the flow convection. The ascending of CO2 can induce the persistent flow cycling effects. The results show that boundary conditions of the field site have determinant effects on the transformation and motion of CO2 in saline aquifers. It is confirmed that the proposed method and numerical model has the reliability to simulate the process of the hydrodynamic trapping, which is the controlling mechanism for the initial period of CO2 storage at time scale of 100 years.

波动方程偏移速度分析

On the subject of oil and gas exploration, migration is an efficacious technique for imagining structures underground. Wave-equation migration (WEM) dominates over other migration methods in accuracy, despite of higher computational cost. However, the advantages of WEM will emerge as the progress of computer technology. WEM is sensitive to velocity model more than others. Small velocity perturbations result in grate divergence in the image pad. Currently, Kirrchhoff method is still very popular in the exploration industry for the reason of difficult to provide precise velocity model. It is very urgent to figure out a way to migration velocity modeling. This dissertation is mainly devoted to migration velocity analysis method for WEM: 1. In this dissertation, we cataloged wave equation prestack depth migration. The concept of migration is introduced. Then, the analysis is applied to different kinds of extrapolate operator to demonstrate their accuracy and applicability. We derived the DSR and SSR migration method and apply both to 2D model. 2. The output of prestack WEM is in form of common image gathers (CIGs). Angle domain common image gathers (ADCIGs) gained by wave equation are proved to be free of artifacts. They are also the most potential candidates for migration velocity analysis. We discussed how to get ADCIGs by DSR and SSR, and obtained ADCIGs before and after imaging separately. The quality of post stack image is affected by CIGs, only the focused or flattened CIGs generate the correct image. Based on wave equation migration, image could be enhanced by special measures. In this dissertation we use both prestack depth residual migration and time shift imaging condition to improve the image quality. 3. Inaccurate velocities lead to errors of imaging depth and curvature of coherent events in CIGs. The ultimate goal of migration velocity analysis (MVA) is to focus scattered event to correct depth and flatten curving event by updating velocities. The kinematic figures are implicitly presented by focus depth aberration and kinetic figure by amplitude. The initial model of Wave-equation migration velocity analysis (WEMVA) is the output of RMO velocity analysis. For integrity of MVA, we review RMO method in this dissertation. The dissertation discusses the general ideal of RMO velocity analysis for flat and dipping events and the corresponding velocity update formula. Migration velocity analysis is a very time consuming work. Respect to computational convenience, we discus how RMO works for synthetic source record migration. In some extremely situation, RMO method fails. Especially in the areas of poorly illuminated or steep structure, it is very difficult to obtain enough angle information for RMO. WEMVA based on wave extrapolate theory, which successfully overcome the drawback of ray based methods. WEMVA inverses residual velocities with residual images. Based on migration regression, we studied the linearized scattering operator and linearized residual image. The key to WEMVA is the linearized residual image. Residual image obtained by Prestack residual migration, which based on DSR is very inefficient. In this dissertation, we proposed obtaining residual migration by time shift image condition, so that, WEMVA could be implemented by SSR. It evidently reduce the computational cost for this method.

浑善达克沙地全新世气候变化与风沙活动—夏日淖湖泊沉积记录

Atmospheric dust as an important factor for climate feedbacks is mainly derived from and drove by the aeolian activities of the semi-arid regions in Asia. Therefore, reconstructing the history of aeolian activities in this region has a great significance in understanding the dynamics of past and future climate changes. This paper made a systematical sedimentology, geochemisty and high-resolution chronology analysis on the sediments from the Xiarinao Lake, which located in Hunshandake Desert, Inner Mongolia, and compared with the meteorological records for identifying the indicator of aeolian activities and climate changes. Based on above work, the evolution history of the aeolian activities and climate changes since 11,000 a BP was reconstructed and the dynamics of the aeolian activities and climate changes in the Hunshandake Desert was discussed. The indicator of aeolian activities was established through the comparison of the clastic and chemical records with the meteorological registers in recent 50 years: 1. The sediments of Xiarinao Lake are mainly derived from eaolian clastic materials and composed of two major components: one is the sand fraction (<4 ) similar to the composition of sand dunes, the other is the silt fraction (>4 ) similar to the atmospheric dust. The sand content and the median grain size, particularly the sand content, show a close association with the intensity of wind, indicating that the particle-size composition of the sediments reflect the variations of the aeolian activities. 2. The proportion of soluble salts relative to the total carbonates in the sediments was correlated to the arid degree registered in meteorological records, suggesting that arid degree might be associated with changes of wind speed and the relative proportion of soluble salts to the total carbonate could be used as the indicator of aridity. 3. The δ18O of authigenic calcite in the Xiarinao Lake show a well correlation with both the atmosphere precipitation and variation of the moisture indicated by the concentrations of soluble salts, suggesting that δ18O of anthigenic calcite is a reliable indicator of humidity. Based on above studies, the evolution history of the aeolian activities and climate changes since 11,000 a BP was reconstructed. At the same time, their association and underlying dynamics was discussed: 1. The evolution history of the eaolian activity was reconstructed by the sand content and median grain size. The result showed that aeolian activities had experienced several different climate periods in the Holocene: the aeolian activities fluctuated with a 500-year cycles during the interval from 10,900 to 8200 a BP; no significant aeolian activities had been developed from 8200 to 6300 a BP; during the interval of 6300 to 2600 a BP, sand content and media grain size had increased gradually, the sand content increased from about 5% to 25%, indicating the intensity of aeolian activity increased in the Hunshandake Desert; since 2600 a BP, intensity of aeolian activity has become stronger and activated the sand dunes. 2. The relative arid degree indicated by the proportion of soluble salts relative to the total carbonate show that both the wind speed and dry degree had being increased since 11,000 a BP. During the interval of 10900 to 6300 a BP, no detectable soluble salt was found in the sediments, indicating the humidity condition; the proportion increased to 20% from 6300 to 4200 a BP, showing the climate became drier; the little change of proportion (20%) suggest that the climate was relative stable during the interval of 4200 to 2600 a BP; after 2600 a BP, the proportion increased, showing the climate became further dry. 3. The δ18O evolution of the lake water was established through analysis of authigenetic calcite δ18O for revealing the variations of relative effective humidity in the Hunshandake Desert. δ18O maintained around -6 ‰ during the interval of 10900 to 8200 a BP, showing the climate was a little dry relative to that of the whole Holocnene; in the period of 8200 to 6300 a BP, δ18O had the most negative values, indicating that it was the most humid interval in Holocene; from 6300 to 4200 a BP, δ18O increased from about -7.5 ‰ to about -3 ‰ gradually, suggesting increased aridity; since 4200 a BP, the climate has become stable and dry. On the whole, the records in this study show that the variation of the relative humidity and aridity is consistent with that of summer monsoon in the Hunshandake Desert. The relationships among the grain size, soluble salts relative content, and the authigenic calcite δ18O indicate aeolian activities is controlled by the aridity in the studied area, and the increased aeolian activities is closely associated with the temperature changes under the relative dry condition. The high wind speed under the low temperature leads to the increased aeolian activity, and the temperature change associated with the aeolian activities might be related with the perturbations of the Siberian High Pressure.

中新世秦安黄土一古土壤序列记录的古气候变化周期

Analysis of periodic oscillations of climate is very important in understanding the behavior of the climate system. Milankovitch hypothesis, which holds that the glacial-interglacial climatic cycles during the Quaternary were primarily driven by variations in orbital parameters, has been supported by substantial geological evidence. Continuous long-term and high-resolution records are crucial to detect how variations of Earth's orbital parameters affected climate before the Quaternary when the boundary conditions were significantly different. Qinan loess formed in the Miocene is nearly continuous aeolian deposit in northern China. Previous study has established a constrained chronology, which provides a basis to examine long-term climatic variations. One of important issues to untangle the mechanisms behind major climate changes is the investigation of climate cycles recorded in Qinan loess. In this paper, two climatic proxies, magnetic susceptibility and redness, are analyzed for QA-I section to evaluate climate cycles using Maximum entropy spectral analysis and Blackman-Tuckey method. Main conclusions are presented as following: Results exhibit significant peaks at periods of 100 ka, 64 ka, 41 ka, 30 ka and 23 ka, but also 1000 ka, 600 ka and 400 ka. These peaks correspond to the dominant periods of the Earth's orbit parameters, which indicates that the formation of the aeolian sediment in northern China might be primarily driven by variations in orbital parameters. Fluctuations with different cycles respectively dominated in different periods. Major shifts in the dominant cycles occurred at 20.3, 19.0, 17.9, 15.2, 12.5 and 11.3 Myr ago. The transition that happened at 17.9 Myr ago was synchronous with the uplift of the Tibetan Plateau, while others at 15.2, 12.5 and 11.3 Myr ago were in good agreement with the timing of the development of Antarctic ice sheet. Therefore we inferred that these shifts might be related to changes in global ice volume and/or the Tibetan uplift. 3. The strong period of 100 ka is observed between 17.9 and 15.2, and 12.5 and 11.3 Myr ago. Ice sheet-climate models that have been used to explain the cause of the 100 ka period since the middle Pleistocene couldn't be responsible for driving the 100 ka climate cycle in the Miocene in Northern China because of the different boundary of climatic conditions between the Quaternary and Miocene. Further investigation is needed to understand how this cycle became dominant in Qinan loess records during these two time segments.

Formation of the V center on H-type zeolites with thermal treatment

The effect of thermal treatment on H-MCM-22 and H-ZSM-5 zeolites was investigated using the electron spin resonance technique. A six-line signal (denoted as A, g = 2.048, A = 22. 15 G) was detected on H-MCM-22 after He purging at high temperatures, whose intensities increased with the treating temperature. The same signal was also found on H-ZSM-5 zeolites with different crystal sizes. The paramagnetic center was identified as a V center, namely, a hole of an electron trapped on an oxygen atom bonding to a nearby aluminum atom. These signals appeared only on a dealuminated sample or a sample concomitantly with dealumination. The formation of the hole might involve an electron transferring from the lattice oxygen to a nonframework aluminum species, and the hyperfine splitting is caused by the interaction between the electron hole locating on the p orbit of oxygen and the framework aluminum bonding with the oxygen. The signal disappeared after the sample was exposed to air or oxygen at room temperature. However, the process was reversible. A new set of signals (denoted as B, g(1) = 2.008, g(2) = 2.003, g(3) = 1.9985) was observed after oxygen adsorption on the H-MCM-22 pretreated with He at 973 K or He purging at 973 K on the H-MCM-22 pretreated with oxygen at 813 K, which was attributed to the O- species.

Photodissociation of 1-bromo-2-chloroethane at 266 nm

The photodissociation of CH2BrCH2Cl at 266 nm has been investigated on the universal crossed molecular beam machine. The primary dissociation step leads exclusively to the formation of CH2CH2Cl radicals and Br atoms in the electronic ground state as well as in the spin-orbit excited state, with a branching ratio 2 +/- 1:8 +/- 1. Photofragment total c.m. translational energy distribution P(E-t) has been obtained and about 64% of the available energy is partitioned into translational energy for Br channel and about 28.5% of the available energy is partitioned into translational energy for Br* channel. The anisotropy parameters are determined to be beta(Br*) = 0.8 +/- 0.2 and beta(Br) = -0.6 +/- 0.2, respectively. Some CH2CH2Cl radicals with large internal excitation (corresponding to formation of ground state Br channel) may undergo secondary dissociation to form CH2CH2 +/- Cl. The experimental results are discussed in terms of a model that involves the initial excitation of two repulsive electronic states: one from an parallel transition to the (3)Q(0) state, and the other from a perpendicular transition to the (3)Q(1), (1)Q states. (C) 1999 Elsevier Science B.V. All rights reserved.

Object Recognition By Alignment Using Invariant Projections of Planar Surfaces

In order to recognize an object in an image, we must determine the best transformation from object model to the image. In this paper, we show that for features from coplanar surfaces which undergo linear transformations in space, there exist projections invariant to the surface motions up to rotations in the image field. To use this property, we propose a new alignment approach to object recognition based on centroid alignment of corresponding feature groups. This method uses only a single pair of 2D model and data. Experimental results show the robustness of the proposed method against perturbations of feature positions.

Theories of Comparative Analysis

Comparative analysis is the problem of predicting how a system will react to perturbations in its parameters, and why. For example, comparative analysis could be asked to explain why the period of an oscillating spring/block system would increase if the mass of the block were larger. This thesis formalizes the task of comparative analysis and presents two solution techniques: differential qualitative (DQ) analysis and exaggeration. Both techniques solve many comparative analysis problems, providing explanations suitable for use by design systems, automated diagnosis, intelligent tutoring systems, and explanation based generalization. This thesis explains the theoretical basis for each technique, describes how they are implemented, and discusses the difference between the two. DQ analysis is sound; it never generates an incorrect answer to a comparative analysis question. Although exaggeration does occasionally produce misleading answers, it solves a larger class of problems than DQ analysis and frequently results in simpler explanations.

Electrons in feldspar II: a consideration of the influence of conduction band-tail states on luminescence processes

Poolton, Nigel; Ozanyan, K.B.; Wallinga, J.; Murray, A.S., (2002) 'Electrons in feldspar II: a consideration of the influence of conduction band-tail states on luminescence processes', Physics and Chemistry of Minerals 29(3) pp.217-225 RAE2008

Predicting interpretability of metabolome models based on behavior, putative identity, and biological relevance of explanatory signals

Enot, D. P., Beckmann, M., Overy, D., Draper, J. (2006). Predicting interpretability of metabolome models based on behavior, putative identity, and biological relevance of explanatory signals. Proceedings of the National Academy of Sciences of the USA, 103(40), 14865-14870. Sponsorship: BBSRC RAE2008

Neural Modeling and Imaging of the Cortical Interactions Underlying Syllable Production

This paper describes a neural model of speech acquisition and production that accounts for a wide range of acoustic, kinematic, and neuroimaging data concerning the control of speech movements. The model is a neural network whose components correspond to regions of the cerebral cortex and cerebellum, including premotor, motor, auditory, and somatosensory cortical areas. Computer simulations of the model verify its ability to account for compensation to lip and jaw perturbations during speech. Specific anatomical locations of the model's components are estimated, and these estimates are used to simulate fMRI experiments of simple syllable production with and without jaw perturbations.