The variation of upper ocean structure and paleoproductivity in the Kuroshio source region during the last 200 kyr

A core from the source region of the Kuroshio warm current (east of the Luzon Island) was analyzed using several proxies in order to study the variability of the Western Pacific Warm Pool (WPWP) during the last two glacial-interglacial cycles. Primary productivity (PP) variations were deduced from variations in the coccolith flora. Primary productivity was higher during glacial periods (the end of Marine Isotope Stage [MIS] 3, some periods in MIS 2 and 6), and decreased during interglacial periods (MIS 7, MIS Se and probably MIS 5c-5d), with the lowest PP in MIS 5e. variations in the delta C-13 difference in benthic and bulk carbonate, thus in the vertical gradient of delta C-13 in dissolved inorganic carbon (Delta delta C-13(c). (wuellerstorfi-N. dutertrei) and Delta delta C-13(c.) (wuellerstorfi-coccolith)) Coincided With the PP Changes, showing that export productivity was low during interglacial periods (MIS 7, MIS 5e and Holocene) and high during glacial periods (MIS 6, probably MIS 5c-5d, late MIS 4 and late MIS 3). Comparison of foraminiferal carbonate dissolution indicators and PP changes reveals that nannofossil assemblage in core Ph05-5 is not sensitive to carbonate dissolution intensity. The depth of the thermocline (DOT) was estimated from planktonic forminiferal assemblages, and was relatively greater during interglacial periods (MIS 7, MIS 5e, probably MIS 5c and Holocene) than during glacials (middle MIS 6, probably MIS 5b and 5d, some periods in MIS 4, MIS 3 and MIS 2). Good coherence between the paleoproductivity records and the DOT suggests that the DOT changes could be the primary control factor in changes of paleoproductivity, and the glacial high productivity in the Kuroshio source region could be associated with a global increase of nutrient concentration in the intermediate waters that upwelled into the photic zone. The low CO2 values derived for intervals of high productivity and a relatively shallow DOT suggest that the changes in biological productivity and DOT in the equatorial Pacific could have modified atmospheric CO2 concentrations. High Sea Surface Temperatures (SSTs) during the warm MIS 5e in combination with intensified monsoonal rain fall could have resulted in a more intense stratification of the upper waters, resulting in low nutrient supply to the surface waters and a resulting decrease in productivity. (C) 2010 Elsevier B.V. All rights reserved.

Effect of low-frequency Rossby wave on thermal structure of the upper southwestern tropical Indian Ocean

We investigate the influence of low-frequency Rossby waves on the thermal structure of the upper southwestern tropical Indian Ocean (SWTIO) using Argo profiles, satellite altimetric data, sea surface temperature, wind field data and the theory of linear vertical normal mode decomposition. Our results show that the SWTIO is generally dominated by the first baroclinic mode motion. As strong downwelling Rossby waves reach the SWTIO, the contribution of the second baroclinic mode motion in this region can be increased mainly because of the reduction in the vertical stratification of the upper layer above thermocline, and the enhancement in the vertical stratification of the lower layer under thermocline also contributes to it. The vertical displacement of each isothermal is enlarged and the thermal structure of the upper level is modulated, which is indicative of strong vertical mixing. However, the cold Rossby waves increase the vertical stratification of the upper level, restricting the variability related to the second baroclinic mode. On the other hand, during decaying phase of warm Rossby waves, Ekman upwelling and advection processes associated with the surface cyclonic wind circulation can restrain the downwelling processes, carrying the relatively colder water to the near-surface, which results in an out-of-phase phenomenon between sea surface temperature anomaly (SSTA) and sea surface height anomaly (SSHA) in the SWTIO.

El Nino phenomenon in simple ocean data assimilation data

To study how the air and sea interact with each other during El Nino/La Nina onsets, extended associate pattern analysis (EAPA) is adopted with the simple ocean data assimilation (SODA) data. The results show that as El Nino/La Nina's parents their behaviors are quite different, there does not exist a relatively independent tropical atmosphere but does exist a relatively independent tropical Pacific Ocean because the air is heated from the bottom surface instead of the top surface and of much stronger baroclinic instability than the sea and has a very large inter-tropical convergence zone covering the most tropical Pacific Ocean. The idea that it is the wester burst and wind convergence, coming from middle latitudes directly that produce the seawater eastward movement and meridional convergence in the upper levels and result in the typical El Nino sea surface temperature warm signal is confirmed again.

A hybrid coupled ocean-atmosphere model and ENSO prediction study

A hybrid coupled ocean-atmosphere model is designed, which consists of a global AGCM and a simple anomaly ocean model in the tropical Pacific. Retroactive experimental predictions initiated in each season from 1979 to 1994 are performed. Analyses indicate that: (1) The overall predictive capability of this model for SSTA over the central-eastern tropical Pacific can reach one year, and the error is not larger than 0.8 degrees C. (2) The prediction skill depends greatly on the season when forecasts start. However, the phenomenon of SPB (spring prediction barrier) is not found in the model. (3) The ensemble forecast method can effectively improve prediction results. A new initialization scheme is discussed.

Effect of Stokes drift on upper ocean mixing

Stokes drift is the main source of vertical vorticity in the ocean mixed layer. In the ways of Coriolis - Stokes forcing and Langmuir circulations, Stokes drift can substantially affect the whole mixed layer. A modified Mellor-Yamada 2.5 level turbulence closure model is used to parameterize its effect on upper ocean mixing conventionally. Results show that comparing surface heating with wave breaking, Stokes drift plays the most important role in the entire ocean mixed layer, especially in the subsurface layer. As expected, Stokes drift elevates both the dissipation rate and the turbulence energy in the upper ocean mixing. Also, influence of the surface heating, wave breaking and wind speed on Stokes drift is investigated respectively. Research shows that it is significant and important to assessing the Stokes drift into ocean mixed layer studying. The laboratory observations are supporting numerical experiments quantitatively.

An approach to prediction of the South China Sea summer monsoon onset

In the present paper, correlation between the South China Sea summer monsoon (SCSSM) onset and heat content in the upper layer of the warm pool in the western Pacific Ocean is examined using the Scripps Institution of Oceanography dataset for the period of 1955-1998 and an approach to prediction the SCSSM onset is proposed. Correlation showes that there exists interdecadal variability of the SCSSM onset demarcated by 1970 with the largest correlation coefficient in the area west of the center of the warm pool rather than near its centers, implying certain effect from other factors involved besides ENSO. As the correlation is poor for the period before 1970, the heat content anomaly of the warm pool after 1970 is used to indicate early or late onset of the SCSSM beforehand. An ideal representative area (1A degrees x1A degrees) for the warm pool heat content was determined with its center at 3A degrees N/138A degrees E. The nearest TAO (TAO-Tropical Atmosphere Ocean-array) mooring to the center is at 2A degrees N/137A degrees E, and chosen to calculate the heat content for prediction. It is suggested that the TAO mooring at 2A degrees N/137A degrees E could be used to predict the SCSSM onset with the heat content in the upper layer, if the correlation between the SCSSM onset and the heat content of the warm pool runs like that of after 1970. On the other hand, if the situation does like the one before 1970, the representative station is determined at 13A degrees S/74A degrees E with relatively poor correlation, meaning that the warm pool in the western Pacific Ocean plays more important role in the SCSSM onset than the Indian Ocean.

Mixed-layer water oscillations in tropical Pacific for ENSO cycle

The main modes of interannal variabilities of thermocline and sea surface wind stress in the tropical Pacific and their interactions are investigated, which show the following results. (1) The thermocline anomalies in the tropical Pacific have a zonal dipole pattern with 160 W as its axis and a meridional seesaw pattern with 6-8 degrees N as its transverse axis. The meridional oscillation has a phase lag of about 90 to the zonal oscillation, both oscillations get together to form the El Nino/La Nina cycle, which behaves as a mixed layer water oscillates anticlockwise within the tropical Pacific basin between equator and 12 degrees N. (2) There are two main patterns of wind stress anomalies in the tropical Pacific, of which the first component caused by trade wind anomaly is characterized by the zonal wind stress anomalies and its corresponding divergences field in the equatorial Pacific, and the abnormal cross- equatorial flow wind stress and its corresponding divergence field, which has a sign opposite to that of the equatorial region, in the off-equator of the tropical North Pacific, and the second component represents the wind stress anomalies and corresponding divergences caused by the ITCZ anomaly. (3) The trade winds anomaly plays a decisive role in the strength and phase transition of the ENSO cycle, which results in the sea level tilting, provides an initial potential energy to the mixed layer water oscillation, and causes the opposite thermocline displacement between the west side and east side of the equator and also between the equator and 12 degrees N of the North Pacific basin, therefore determines the amplitude and route for ENSO cycle. The ITCZ anomaly has some effects on the phase transition. (4) The thermal anomaly of the tropical western Pacific causes the wind stress anomaly and extends eastward along the equator accompanied with the mixed layer water oscillation in the equatorial Pacific, which causes the trade winds anomaly and produces the anomalous wind stress and the corresponding divergence in favor to conduce the oscillation, which in turn intensifies the oscillation. The coupled system of ocean-atmosphere interactions and the inertia gravity of the mixed layer water oscillation provide together a phase-switching mechanism and interannual memory for the ENSO cycle. In conclusion, the ENSO cycle essentially is an inertial oscillation of the mixed layer water induced by both the trade winds anomaly and the coupled ocean-atmosphere interaction in the tropical Pacific basin between the equator and 12 degrees N. When the force produced by the coupled ocean-atmosphere interaction is larger than or equal to the resistance caused by the mixed layer water oscillation, the oscillation will be stronger or maintain as it is, while when the force is less than the resistance, the oscillation will be weaker, even break.

Carbon Dioxide Exchange Between the Atmosphere and an Alpine Shrubland Meadow During the Growing Season on the Qinghai-Tibetan Plateau

In the present study, we used the eddy covariance method to measure CO2 exchange between the atmosphere and an alpine shrubland meadow ecosystem (37°36'N, 101°18'E; 3 250 m a.s.l.) on the Qinghai-Tibetan Plateau, China, during the growing season in 2003, from 20 April to 30 September. This meadow is dominated by formations of Potentilla fruticosa L. The soil is Mol-Cryic Cambisols. During the study period, the meadow was not grazed. The maximum rates of CO2 uptake and release derived from the diurnal course of CO2 flux were -9.38 and 5.02 μmol•m-2•s-1, respectively. The largest daily CO2 uptake was 1.7 g C•m-2•d-1 on 14 July, which is less than half that of an alpine Kobresia meadow ecosystem at similar latitudes. Daily CO2 uptake during the measurement period indicated that the alpine shrubland meadow ecosystem may behave as a sink of atmospheric CO2 during the growing season. The daytime CO2 uptake was correlated exponentially or linearly with the daily photo synthetic photon flux density each month. The daytime average water use efficiency of the ecosystem was 6.47 mg CO2/g H2O. The efficiency of the ecosystem increased with a decrease in vapor pressure deficit.

GEO-ECOLOGICAL TRANSECT STUDIES IN NORTHEAST TIBET (QINGHAI, CHINA) REVEAL HUMAN-MADE MID-HOLOCENE ENVIRONMENTAL CHANGES IN THE UPPER YELLOW RIVER CATCHMENT CHANGING FOREST TO GRASSLAND

Geo-ecological transect studies in the pastures of the upper catchment of the HuangHe (99 degrees 30'-100 degrees 00'E/35 degrees 30'-35 degrees 40'N'; 3,000-4,000 in a.s.l., Qinghai province, China) revealed evidence that pastures replace forests. Plot-based vegetation records and fenced grazing exclosure experiments enabled the identification of grazing indicator plants for the first time. The mapping of vegetation patterns of pastures with isolated juniper and Spruce forests raise questions as to the origin of the grasslands, which arc widely classified as "natural" at present. Soil investigations and charcoal fragments of Juniperus (8,153 +/- 63 uncal BP) and Picea (6,665 +/- 59 uncal BP) provide evidence of the wider presence of forests. As temperatures and rainfall records undoubtedly represent a forest climate, it is assumed that the present pastures have replaced forests. Circumstantial evidence arising from investigations into the environmental history of the Holocene effectively substantiates this theory.

Effects of temperature on CO2 exchange between the atmosphere and an alpine meadow

The alpine meadow ecosystem on the Qinghai-Tibetan Plateau is characterized by low temperatures because of its high elevation. The low-temperature environment may limit both ecosystem photosynthetic CO2 uptake and ecosystem respiration, and thus affect the net ecosystem CO2 exchange (NEE). We clarified the low-temperature constraint on photosynthesis and respiration in an alpine meadow ecosystem on the northern edge of the plateau using flux measurements obtained by the eddy covariance technique in two growing seasons. When we compared NEE during the two periods, during which the leaf area index and other environmental parameters were similar but the mean temperature differed, we found that NEE from 9 August to 10 September 2001, when the average temperature was low, was greater than that during the same period in 2002, when the average temperature was high, but the ecosystem gross primary production was similar during the two periods. Further analysis showed that ecosystem respiration was significantly higher in 2002 than in 2001 during the study period, as estimated from the relationship between temperature and nighttime ecosystem respiration. The results suggest that low temperature controlled the NEE mainly through its influence on ecosystem respiration. The annual NEE, estimated from 15 January 2002 to 14 January 2003, was about 290 g CO2 m(-2) year(-1). The optimum temperature for ecosystem NEE under light-saturated conditions was estimated to be around 15 degrees C.

Energy exchange between the atmosphere and a meadow ecosystem on the Qinghai-Tibetan Plateau

To reveal the potential contribution of grassland ecosystems to climate change, we examined the energy exchange over an alpine Kobresia meadow on the northeastern Qinghai-Tibetan Plateau. The annual pattern of energy exchange showed a clear distinction between periods of frozen soil with the daily mean soil temperature at 5 cm (T-s5 ≤ 0 ° C) and non-frozen soil (T-s5 > 0 ° C). More than 80% of net radiation was converted to sensible heat (H) during the frozen soil period, but H varied considerably with the change in vegetation during the non-frozen soil period. Three different sub-periods were further distinguished for the later period: (1) the pre-growth period with Bowen ratio (β) > 1 was characterized by a high β of 3.0 in average and the rapid increase of net radiation associated with the increases of H, latent heat (LE) and soil heat; (2) during the Growth period when β ≤ 1, the LE was high but H fluxes was low with β changing between 0.3 and 0.4; (3) the post-growth period with average β of 3.6 when H increased again and reached a second maximum around early October. The seasonal pattern suggests that the phenology of the vegetation and the soil water content were the major factors affecting the energy partitioning in the alpine meadow ecosystem. © 2005 Elsevier B.V. All rights reserved.

Carbon dioxide exchange between the atmosphere and an alpine meadow ecosystem on the Qinghai-Tibetan Plateau, China

We used the eddy covariance method to measure the M exchange between the atmosphere and an alpine meadow ecosystem (37degrees29-45'N, 101degrees12-23'E, 3250m a.s.l.) on the Qinghai-Tibetan Plateau, China in the 2001 and 2002 growing seasons. The maximum rates Of CO2 uptake and release derived from the diurnal course Of CO2 flux (FCO2) were -10.8 and 4.4 mumol m(-2) s(-1), respectively, indicating a relatively high net carbon sequestration potential as compared to subalpine coniferous forest at similar elevation and latitude. The largest daily CO2 uptake was 3.9 g cm(-2) per day on 7 July 2002, which is less than half of those reported for lowland grassland and forest at similar latitudes. The daily CO2 uptake during the measurement period indicated that the alpine ecosystem might behave as a sink of atmospheric M during the growing season if the carbon lost due to grazing is not significant. The daytime CO2 uptake was linearly correlated with the daily photosynthetic photon flux density each month. The nighttime averaged F-CO2 showed a positive exponential correlation with the soil temperature, but apparently negative correlation with the soil water content. (C) 2004 Elsevier B.V. All rights reserved.

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

中高层大气对低层大气声重波响应的分析

Acoustic Gravity waves (AGW) play an important role in balancing the atmospheric energy and momentum budget. Propagation of gravity wave in the atmosphere is one of the important factors of changing middle and upper atmosphere and ionosphere. The purpose of this dissertation is to study the propagation of gravity wave in a compression atmosphere whit means of numerical simulation and to analyze the response of middle and upper atmosphere to pulse disturbance from lower atmosphere. This work begins with the establishment of 2-D fully nonlinear compressible atmospheric dynamic model in polar coordinate, which is used ton numerically study gravity wave propagation. Then the propagation characteristics of acoustic gravity wave packets are investigated and discussed. We also simulate the response of middle and upper atmosphere to pulse disturbance of lower atmosphere in background winds or without background winds by using this model and analyze the data we obtained by using Fourier Transform (FT), Short-time Fourier Transform (STFT) and Empirical Mode Decomposition (EMD) method which is an important part of Hilbert-Huang Transform (HHT). The research content is summarized in the following: 1. By using a two-dimensional full-implicit-continuous-Eulerian (FICE) scheme and taking the atmospheric basic motion equations as the governing equations, a numerical model for nonlinear propagation of acoustic gravity wave disturbance in two-dimensional polar coordinates is solved. 2. Then the propagation characteristics of acoustic gravity wave packets are investigated and discussed. Results of numerical simulation show that the acoustic gravity wave packets propagate steadily upward and keep its shape well after several periods. 3. We simulate the response of middle and upper atmosphere to pulse disturbance of lower atmosphere in background winds or without background winds by using this model, and obtain the distribution of a certain physical quantity in time and space from earth’s surface to 300km above. The results reveal that the response of ionosphere occurs at a large horizontal distance from the source and the disturbance becomes greater with increasing of height. The situation when the direction of the background wind is opposite to or the same as the direction of disturbed velocity of gravity-wave is studied. The results show that gravity wave propagating against winds is easier than those propagating along winds and the background wind can accelerate gravity wave propagation. Just upon the source, an acoustic wave component with period of 6 min can be found. These images of simulation are similar to observations of the total electron content (TEC) disturbances caused by the great Sumatra-Andaman earthquake on December 26 in 2004. 4. Using the EMD method the disturbed velocity data of a certain physical quantity in time and space can be decomposed into a series of intrinsic mode function (IMF) and a trend mode respectively. The results of EMD reveal impact of the gravity wave frequency under the background winds.

太阳耀斑和日食期间电离层行为的统计和模拟研究

The ionosphere is the ionized component of the Earth's upper atmosphere. Solar EUV radiation is the source of ionospheric ionization. Thus the ionosphere is affected strongly by the variations in solar radiation. Solar flares and solar eclipses can induce remarkable short time changes in solar radiation: the solar radiation would increase suddenly during solar flares and decrease significantly during solar eclipses. Solar flare and eclipse events not only affect directly the photochemical processes, but also affect the dynamic processes, and even affect the neutral atmosphere, which is strongly coupled with the ionosphere. The study on the ionospheric response to solar flares and eclipses can advance our knowledge on the ionosphere and its photochemical and dynamic processes and help us to evaluate the ionospheric parameters (such as ion loss coefficients). In addition, the study on the ionospheric responses to solar flares and eclipses is an important part of the ionospheric space weather, which can provide guides for space weather monitoring. This thesis devotes to the study on the ionospheric responses to solar flares and solar eclipses. I have developed two models to simulate the variations of solar EUV radiation during solar flares and solar eclipses, and involved in developing a 2D mid- and low-latitude ionospheric model. On the basis of some observed data and the ionospheric model, I study the temporal and spatial variations of the ionosphere during solar flares and eclipses, and investigate the influences of solar activity, solar zenith angle, neutral gas density, and magnetic dip angle on the ionospheric responses to solar flares and solar eclipses. The main points of my works and results are summarized as follows. 1. The ionospheric response to the X17.2 solar flare on October 28, 2003 was modeled via using a one-dimension theoretical ionospheric model. The simulated variation of TEC is in accordance with the observations, though there are some differences in the amplitude of the variation. Then I carried out a series of simulations to explore the local time and seasonal dependences of the ionospheric responses to solar flares. These calculations show that the ionospheric responses are largely related with the solar zenith angle (SZA). During the daytime (small SZA), most of the increases in electron density occur at altitudes below 300 km with a peak at around 115 km; whereas around sunrise and sunset (SZA>90°), the strongest ionospheric responses occur at much higher altitudes. The TEC increases slower at sunrise than at sunset, which is caused by the difference in the evolution of SZA at sunrise and sunset: SZA decreases with time at sunrise and increase with time at sunset. The ionospheric response is largest in summer and smallest in winter, which is also related to the seasonal difference of SZA. 2. Based on the observations from the ionosondes in Europe and the ionospheric model, I investigated the differences of the ionosphere responses to solar eclipses between the E-layer and F1-layer. Both the observation and simulation show that the decrease in foF1 due to the solar eclipses is larger than that in foE. This effect is due to that the F1 region locates at the transition height between the atomic ion layer and the molecular ion layer. With the revised model of solar radiation during solar flares, our model calculates the radiations from both the inside and outside of photosphere. Large discrepancy can be found between the observations and the calculations with an unrevised model, while the calculations with the revised model consist with the observations. 3. I also explore the effects of the F2-layer height, local time, solar cycle, and magnetic dip angle on the ionospheric responses to solar eclipses via using an ionospheric model and study on the solar zenith angle and the dip dependences by analyzing the data derived from 23 ionosonde stations during seven eclipse events. Both the measured and simulated results show that these factors have significant effect on the ionospheric response. The larger F2-layer height causes the smaller decrease in foF2, which is because that the electron density response decreases with height. The larger dip results in the smaller eclipse effect on the F2 layer, because the larger dip would cause the more diffusion from the top ionosphere which can make up for the plasma loss. The foF2 response is largest at midday and decreases with the increasing SZA. The foF2 response is larger at high solar activity than at low solar activity. The simulated results show that the local time and solar activity discrepancy of the eclipse effect mainly attribute to the difference of the background neutral gas density. 4. I carried out a statistical study on the latitudinal dependence of the ionospheric response to solar eclipses and modeled this latitudinal dependence by the ionospheric model. Both the observations and simulations show that the foF2 and TEC responses have the same latitudinal dependence: the eclipse effects on foF2 and TEC are smaller at low latitudes than at middle latitudes; at the middle latitudes (>40°), the eclipse effect decreases with increasing latitude. In addition, the simulated results show the change in electron temperature at the heights of above 300 km of low latitudes is much smaller than that at the same heights of middle latitudes. This is due to the smaller decrease in photoelectron production rate at its conjugate low heights. 5. By analyzing the observed data during the October 3, 2005 solar eclipse, I find some significant disturbances in the conjugate region of the eclipse region, including a decrease in Te, an increase in foF2 and TEC, and an uprising in hmF2. I also simulated the ionosphere behavior during this eclipse using a mid-low latitude ionospheric model. The simulations reproduce the measured ionospheric disturbances mentioned above in the conjugated hemisphere. The simulations show that the great loss of arriving photoelectron heat from the eclipse region is the principal driving source for the disturbances in the conjugate hemisphere.