柴达木盆地地震勘探关键技术研究

The rugged surface topography determined the seismic data acquisition construction conditions and the seismic wave explosive and receiver quality in Qaidam Basin. This dissertation systematically researched the seismic acquisition, imaging process and the attribute analysis techniques of complicated oil and gas reservoir. The main research achievements and cognitions are as follows: 1. Through the stimulation effects research and analysis from the aspect of lithologic water-containing differences, it’s specific that stable hydrous sand layer can effectively enhance the stimulation effects combined with the corresponding field tests. The seismic data S/N ratio has been improved due to the combination explosive stimulation. Through the fold number and maximum offset analyses of target horizon, the complicated geometry has been optimized and the S/N ratio of seismic data has been improved, which made an important basis for improvement of 3D seismic data. 2. It has been proved that the first arrival refraction static correction method under the model constraint of fine surface survey is suitable to the Qaidam Basin of western areas by the real seismic data processing. Although the refraction horizon of near surface has some changes in a certain extent, it’s steady basically. The refraction horizon can be continuously traced in sections, so it’s qualified for the refraction static correction method on the whole. 3. The research is based on the curved-ray pre-stack time migration techniques of rough topography, and improved the imaging precision of complex areas. This techniques adopted the constant and variable velocity scanning mode and enhanced the velocity analysis precision. The 3D pre-stack time migration techniques reasonably solved the imaging and velocity multiple solutions problems of steep-dip faults and the intersections of horizontal layers. What’s more, fine velocity analysis and mute are very important to enhance the imaging precision of the seismic data in complicated Wunan areas. 4. The 3D seismic data edge-preserving processing methods have been realized due to the image process techniques. Because this method uses the large range filter, it can attenuate the noise maximally. The faults, break points, lithologic pinchout points and lithologic body of small scale such as river will not be influenced by blur because of the edge-preserving characterization of the method which is really an effective assistant technique of low S/N ratio seismic data attribute analysis. 5. The use of spectral decomposition technique can effectively identify the reservoirs. The special geology body which will not be identified (or without obvious characters) in the seismic profile may be found through the details changes of different frequencies in the amplitude profiles.

中低纬度变化磁场的暴时特征研究

Magnetic storm is a kind of severe disturbances in the whole solar-earth electromagnetic space. It has significant effects on communication, electric power, oil transport pipe and human activities in space. Therefore, magnetic storms are worth for applications systems, not only being a favorable issue for scientists. In this paper, the spatial and temporal distributions of the magnetic fields produced by the magnetosphere-ionosphere current systems during storms are studied. Four parts are included in this paper decomposion of different disturbances with different origins, topological structure of the ring current, the asymmetric characteristics of the ring current, and the statistic peculiarities of the day-to-day variability (DTD) of Sq. 1 The decomposition of magnetic disturbances at mid-low latitudes and its evolutions during storms Transient variations in the geomagnetic field recorded at mid-low latitudes mainly include the storm-time variation (Dst), solar quiet daily variation (Sq) and disturbance daily variation (SD). With the data of the geomagnetic meridian chain observatories in China, 25 storms during the period of 1997 to 1999 have been analyzed. According to the features of different variations, a method of “three-steps decomposition” is developed by using the method of Natural Orthogonal Components (NOC), Correlation Analysis and Fourier Analysis to separate those three components in turn. The results show that, the first eigenmode by the MNOC clearly describing the special distribution and temporal evolution of storm-time variation, in addition, Correlation Analysis and Fourier Analysis offer a useful method to extract the Sq and SD variations. The latitudinal shift of the Sq current focus seems to be the principal reason of the day-to-day variaitons in the daily range of Sq. The magnitude of SD reaches a maximum during the main phase, and then gradually decreases. 2 The topology structure of the ring current during storms Both the mechanism of the ring current and the geomagnetic data suggest that the central plane of the ring current is declining to the geomagnetic equator plane with a tilt angle δ. Using the H and Z component data at two stations in a meridian chain, we deduce a new parameter describing the invariable peculiarity of different storms. Then the δ angle is calculated by using the data from a meridian chain and tested with the ERC model. Finally the deduced tilt angles are used to modify Dst index. 3 The asymmetric characteristics of the ring current during storms The variations of the geomagnetic field at mid-low latitudes show a significant dawn-dusk asymmetry, resulting from the superposition of the fields from the symmetric ring current and the partial ring current. On the basis of the data from the 20°E, 30°E meridian chains and 30°N latitudinal chain, the dawn-dusk asymmetry is investigated by using three methods, namely, statistic analysis, ring current model calculation and typical event analysis. This characteristic implies the asymmetry of the spatial distribution of the ring current. In addition, during the main phase after the sudden commencement (SC), H field increases and reaches maximum around noontime, implying the effect of the Chapman-Ferraro current. 4 The statistic characteristics of the day-to-day variability and its mechanism The day-to-day variability of the geomagnetic Sq field is studied by using the magnetic data from a meridian chain of magnetometers along 120° E longitude. The method of NOC is applied to separate the Sq variation from complicated disturbances. The first eigenmode with the largest eigenvalue represents fairly well the Sq variation with a conspicuous day-to-day variability in the daily range. For the stations on the same north- or south-side of the Sq current system focus, the day-to-day variations show a positive correlation. In contrast, for the stations on the different sides of the Sq focus, they show a negative correlation, suggesting an important role of latitudinal shift of the Sq current system focus to the day-to-day variability of the Sq daily range. The Sq daily range is correlated with the magnetic indices Ap and Dst in a peculiar way: on some severe disturbed days, noticeably enhancements of the Sq are observed, implying increases of the ionospheric conductivities and/or tidal wind velocities; on other severe disturbed days, however, dramatically reduced Sq variations occur, suggesting dominant effects of the ‘disturbance dynamo’ process.

基于匹配追踪方法的地震信号时频分析研究

Seismic signal is a typical non-stationary signal, whose frequency is continuously changing with time and is determined by the bandwidth of seismic source and the absorption characteristic of the media underground. The most interesting target of seismic signal’s processing and explaining is to know about the local frequency’s abrupt changing with the time, since this kind of abrupt changing is indicating the changing of the physical attributes of the media underground. As to the seismic signal’s instantaneous attributes taken from time-frequency domain, the key target is to search a effective, non-negative and fast algorithm time-frequency distribution, and transform the seismic signal into this time-frequency domain to get its instantaneous power spectrum density, and then use the process of weighted adding and average etc. to get the instantaneous attributes of seismic signal. Time-frequency analysis as a powerful tool to deal with time variant non-stationary signal is becoming a hot researching spot of modern signal processing, and also is an important method to make seismic signal’s attributes analysis. This kind of method provides joint distribution message about time domain and frequency domain, and it clearly plots the correlation of signal’s frequency changing with the time. The spectrum decomposition technique makes seismic signal’s resolving rate reach its theoretical level, and by the method of all frequency scanning and imaging the three dimensional seismic data in frequency domain, it improves and promotes the resolving abilities of seismic signal vs. geological abnormal objects. Matching pursuits method is an important way to realize signal’s self-adaptive decomposition. Its main thought is that any signal can be expressed by a series of time-frequency atoms’ linear composition. By decomposition the signal within an over completed library, the time-frequency atoms which stand for the signal itself are selected neatly and self-adaptively according to the signal’s characteristics. This method has excellent sparse decomposition characteristics, and is widely used in signal de-noising, signal coding and pattern recognizing processing and is also adaptive to seismic signal’s decomposition and attributes analysis. This paper takes matching pursuits method as the key research object. As introducing the principle and implementation techniques of matching pursuits method systematically, it researches deeply the pivotal problems of atom type’s selection, the atom dictionary’s discrete, and the most matching atom’s searching algorithm, and at the same time, applying this matching pursuits method into seismic signal’s processing by picking-up correlative instantaneous messages from time-frequency analysis and spectrum decomposition to the seismic signal. Based on the research of the theory and its correlative model examination of the adaptively signal decomposition with matching pursuit method, this paper proposes a fast optimal matching time-frequency atom’s searching algorithm aimed at seismic signal’s decomposition by frequency-dominated pursuit method and this makes the MP method pertinence to seismic signal’s processing. Upon the research of optimal Gabor atom’s fast searching and matching algorithm, this paper proposes global optimal searching method using Simulated Annealing Algorithm, Genetic Algorithm and composed Simulated Annealing and Genetic Algorithm, so as to provide another way to implement fast matching pursuit method. At the same time, aimed at the characteristics of seismic signal, this paper proposes a fast matching atom’s searching algorithm by means of designating the max energy points of complex seismic signal, searching for the most optimal atom in the neighbor area of these points according to its instantaneous frequency and instantaneous phase, and this promotes the calculating efficiency of seismic signal’s matching pursuit algorithm. According to these methods proposed above, this paper implements them by programmed calculation, compares them with some open algorithm and proves this paper’s conclusions. It also testifies the active results of various methods by the processing of actual signals. The problems need to be solved further and the aftertime researching targets are as follows: continuously seeking for more efficient fast matching pursuit algorithm and expanding its application range, and also study the actual usage of matching pursuit method.

太阳辐射对热层和电离层变化性的影响

Solar ultraviolet (UV) radiation at wavelengths less than 400 nm is an important source of energy for aeronomic processes throughout the solar system. Solar UV photons are absorbed in planetary atmospheres, as well as throughout the heliosphere, via photodissociation of molecules, photoionization of molecules and atoms, and photoexcitation toexcitation including resonance scattering. In this paper, the solar irradiances data measured by TIMED SEE, as well as the solar proxies such as F10.7 and Mg II, thermosphere neutral density of CHAMP measurements and topside ionospheric plasmas densities from DMSP, are used to analyze solar irradiance effects on the variabilities of the thermosphere and the ionosphere. First, thermosphere densities near 410 km altitude are analyzed for solar irradiance variability effects during the period 2002-2004. Correlations between the densities and the solar irradiances for different spectral lines and wavelength ranges reveal significantly different characteristics. The density correlates remarkably well with all the selected solar irradiances except the lower chromospheric O I (130.4 nm) emission. Among the chosen solar proxies, the Mg II core-to-wing ratio index, EUV (30-120 nm) and F10.7 show the highest correlations with the density for short-term (< ~27 days) variations. For both long- (> ~27 days) and short-term variations, linear correlation coefficients exhibit a decreasing trend from low latitudes towards high latitudes. The density variability can be effectively modeled (capturing 71% of the variance) using multiple solar irradiance indices, including F10.7, SEUV (the EUV 30-120 nm index), and SFUV (the FUV 120-193 nm index), in which a lag time of 1 day was used for both F10.7 and SEUV, and 5 days for SFUV. In our regression formulation SEUV has the largest contribution to the density variation (40%), with the F10.7 having the next largest contribution (32%) and SFUV accounting for the rest (28%). Furthermore, a pronounced period of about 27.2 days (mean period of the Sun's rotation) is present in both density and solar irradiance data of 2003 and 2004, and a pronounced period of about 54.4 days (doubled period of the solar rotation) is also revealed in 2004. However, soft X-ray and FUV irradiances did not present a pronounced 54.4 day period in 2004, in spite of their high correlation with the densities. The Ap index also shows 54-day periodicities in 2004, and magnetic activity, together with solar irradiance, affects the 54-day variation in density significantly. In addition, NRLMSISE00, DTM-2000 and JB2006 model predictions are compared with density measurements from CHAMP to assess their accuracy, and the results show that these models underestimate the response of the thermosphere to variations induced by solar rotation. Next, the equatorial topside ionospheric plasmas densities Ni are analyzed for solar irradiance variability effects during the period 2002-2005. Linear correlations between Ni and the solar irradiances for different wavelength ranges reveal significantly different characteristics. XUV (0-35 nm) and EUV (115-130 nm) show higher correlation with Ni for the long-term variations, whereas EUV (35-115 nm) show higher correlation for the short-term variations. Moreover, partial correlation analysis shows that the long-term variations of Ni are affected by both XUV (0-35 nm) and EUV (35-115 nm), whereas XUV (0-35 nm) play a more important role; the short-term variations of Ni are mostly affected by EUV (35-115 nm). Furthermore, a pronounced period of about 27 days is present in both Ni and solar irradiance data of 2003 and 2004, and a pronounced period of about 54 days is also revealed in 2004. Finally, prompted by previous studies that have suggested solar EUV radiation as a means of driving the semiannual variation, we investigate the intra-annual variation in thermosphere neutral density near 400 km during 2002-2005. The intra-annual variation, commonly referred to as the ‘semiannual variation’, is characterized by significant latitude structure, hemispheric asymmetries, and inter-annual variability. The magnitude of the maximum yearly difference, from the yearly minimum to the yearly maximum, varies by as much as 60% from year to year, and the phases of the minima and maxima also change by 20-40 days from year to year. Each annual harmonic of the intra-annual variation, namely, annual, semiannual, ter-annual and quatra-annual, exhibits a decreasing trend from 2002 through 2005 that is correlated with the decline in solar activity. In addition, some variations in these harmonics are correlated with geomagnetic activity, as represented by the daily mean value of Kp. Recent empirical models of the thermosphere are found to be deficient in capturing most of the latitude dependencies discovered in our data. In addition, the solar flux and geomagnetic activity proxies that we have employed do not capture some latitude and inter-annual variations detected in our data. It is possible that these variations are partly due to other effects, such as seasonal-latitudinal variations in turbopause altitude (and hence O/N2 composition) and ionosphere coupling processes that remain to be discovered in the context of influencing the intra-annual variations depicted here. Our results provide a new dataset to challenge and validate thermosphere-ionosphere general circulation models that seek to delineate the thermosphere intra-annual variation and to understand the various competing mechanisms that may contribute to its existence and variability. We furthermore suggest that the term “intra-annual” variation be adopted to describe the variability in thermosphere and ionosphere parameters that is well-captured through a superposition of annual, semiannual, ter-annual, and quatra-annual harmonic terms, and that “semiannual’ be used strictly in reference to a pure 6-monthly sinusoidal variation. Moreover, we propose the term “intra-seasonal” to refer to those shorter-term variations that arise as residuals from the above Fourier representation.

一种TEC地图构建方法与全球电子总含量的气候学研究

Global positioning system (GPS) can not only provide precise service for navigation and timing, but also be used to investigate the ionospheric variation. From the GPS observations, we can obtain total electron content (TEC), so-called GPS TEC, which is used to characterize the ionospheric structure. This thesis mainly concerns about GPS TEC data processing and ionospheric climatological analysis as follows. Firstly, develop an algorithm for high-resolution global ionospheric TEC mapping. According to current algorithms in global TEC mapping, we propose a practical way to calibrate the original GPS TEC with the existing GIM results. We also finish global/local TEC mapping by model fitting with the processed GPS TEC data; in practice, we apply it into the local TEC mapping in Southeast of China and obtain some initial results. Next, suggest a new method to calculate equivalent ionospheric global electron content (GEC). We calculate such an equivalent GEC with the TEC data along the geographic longitude 120°E. With the climatological analysis, we can see that GEC climatological variation is mainly composed of three factors: solar cycle, annual and semiannual variations. Solar cycle variation is dominant among them, which indicates the most prominent influence; both annual and semiannual variations play a secondary role and are modulated by solar activity. We construct an empirical GEC model driven by solar activity and seasonal factors on the basis of partial correlation analysis. Generally speaking, our researches not only show that GPS is advantageous in now-casting ionospheric TEC as an important observation, but also show that GEC may become a new index to describe the solar influence on the global ionosphere since the great correlation between GEC and solar activity factor indicates the close relationship between the ionosphere and solar activity.

Numerical solution of the incompressible Navier-Stokes equations with an upwind compact difference scheme

A new finite difference method for the discretization of the incompressible Navier-Stokes equations is presented. The scheme is constructed on a staggered-mesh grid system. The convection terms are discretized with a fifth-order-accurate upwind compact difference approximation, the viscous terms are discretized with a sixth-order symmetrical compact difference approximation, the continuity equation and the pressure gradient in the momentum equations are discretized with a fourth-order difference approximation on a cell-centered mesh. Time advancement uses a three-stage Runge-Kutta method. The Poisson equation for computing the pressure is solved with preconditioning. Accuracy analysis shows that the new method has high resolving efficiency. Validation of the method by computation of Taylor's vortex array is presented.

Experimental Study On Liquid Free Surface In Buoyant-Thermocapillary Convection

We investigate the surface deformations of buoyant-thermocapillary convection in a rectangular cavity clue to gravity and temperature gradient between the two sidewalls. The cavity is 52mm x 42mm in horizontal cross section, the thickness of liquid layer h is changed from 2.5mm to 6.5mm. Surface deformations of h = 3.5mm and 6.0mm are discussed and compared. Temperature difference is increased gradually, and the flow in the liquid layer will change from stable convection to unstable convection. Two kinds of optical diagnostic system with image processor are developed for study of the kinetics of buoyant-thermocapillary convection, they give out the information of liquid free surface. The quantitative results are calculated by Fourier transform and correlation analysis, respectively. With the increasing temperature gradient, surface deformations calculated are more declining. It is interesting phenomenon that the inclining directions of the convections in thin and thick liquid layers are different. For a thin layer, the convection is mainly controlled by thermocapillary effect. However, for a thick layer, the convection is mainly controlled by buoyancy effect. The surface deformation theoretically analysed is consistent with our experimental results. The present experiment proves that surface deformation is related to temperature gradient and thickness of the liquid layer. In other words, surface deformation lies on capillary convection and buoyancy convection.

The influence of buoyancy on the oscillatory thermocapillary convection with small bond number

This work is an experimental study of unidirectional bamboo-epoxy laminates of varying laminae number, in which tensile, compressive, flexural and interlaminar shear properties are evaluated. Further, the disposition of bamboo fibre, the parenchymatous tissue, and the resin matrix under different loading conditions are examined. Our results show that the specific strength and specific modulus of bamboo-epoxy laminates are adequate, the former being 3 to 4 times that of mild steel. Its mechanical properties are generally comparable to those of ordinary glass-fibre composites. The fracture behaviour of bamboo-epoxy under different loading conditions were observed using both acoustic emission techniques and scanning electron microscopy. The fracture mode varied with load, the fracture mechanism being similar to glass and carbon reinforced composites. Microstructural analyses revealed that natural bamboo is eligibly a fibre composite in itself; its inclusion in a plastic matrix will help solve the problems of cracking due to desiccation and bioerosion caused by insect pests. Furthermore, the thickness and shape of the composite can be tailored during fabrication to meet specific requirements, thereby enabling a wide spectrum of applications.

大口径八边形钕玻璃片支撑系统的光机集成分析

控制由机械装夹方式所引入的波前畸变以提高惯性约束聚变(ICF)输出光束的质量，是在大口径钕玻璃片主放大器结构设计中必须考虑的。提出了一种新的有限元变形结果与光学元件面形畸变之间的数据处理方式，并与传统方式进行了对比。基于新的数据处理接口，利用光机集成分析方法对大口径八边形钕玻璃片的支撑系统结构设计参数进行优化。优化的结果保证了由支撑系统引起的透过波前畸变小于十分之一波长，同时波前畸变与设计参数变动的相关性最小。

Effects of deposition rates on laser damage threshold of TiO2/SiO2 high reflectors

TiO2 single layers and TiO2/SiO2 high reflectors (HR) are prepared by electron beam evaporation at different TiO2 deposition rates. It is found that the changes of properties of TiO2 films with the increase of rate, such as the increase of refractive index and extinction coefficient and the decrease of physical thickness, lead to the spectrum shift and reflectivity bandwidth broadening of HR together with the increase of absorption and decrease of laser-induced damage threshold. The damages are found of different morphologies: a shallow pit to a seriously delaminated and deep crater, and the different amorphous-to-anatase-to-rutile phase transition processes detected by Raman study. The frequency shift of Raman vibration mode correlates with the strain in. film. Energy dispersive X-ray analysis reveals that impurities and non-stoichiometric defects are two absorption initiations resulting to the laser-induced transformation. (C) 2008 Elsevier B. V. All rights reserved.

Influence of leaf tissue structure on host feeding selection by pea leaf miner Liriomyza huidobrensis (Diptera : Agromyzidae)

Host feeding selection by the female pea leafminer, Liriomyza huidobrensis, on 47 species of plants was studied. The leaves were sectioned by microtome, and 15 characteristics of the leaf tissue structure were measured under a microscope. Correlation analysis between host feeding selection and leaf tissue structure indicated that the preference of host feeding selection was positively correlated with the percentage of moisture content of leaves and negatively with thickness of the epidermis wall, and densities of the palisade and spongy tissues of leaves. Leaf tissue structure was influential in feeding and probing behavior of female L. huidobrensis. So, thickness of epidermis wall, densities of the palisade and spongy tissues can act as a physical barrier to female oviposition. Furthermore, higher densities of palisade and spongy tissues can be considered a resistant trait which affects mining of leaf miner larvae as well. As a result, plants with lower leaf moisture content may not be suitable for the development of L. huidobrensis.

Plankton Community Succession in Artificial Systems Subjected to Cyanobacterial Blooms Removal using Chitosan-Modified Soils

Using artificial systems to simulate natural lake environments with cyanobacterial blooms, we investigated plankton community succession by polymerase chain reaction-denaturing gradient gel electrophoresis (PCR-DGGE) fingerprinting and morphological method. With this approach, we explored potential ecological effects of a newly developed cyanobacterial blooms removal method using chitosan-modified soils. Results of PCR-DGGE and morphological identification showed that plankton communities in the four test systems were nearly identical at the beginning of the experiment. After applying the newly developed and standard removal methods, there was a shift in community composition, but neither chemical conditions nor plankton succession were significantly affected by the cyanobacteria removal process. The planted Vallisneria natans successfully recovered after cyanobacteria removal, whereas that in the box without removal process did not. Additionally, canonical correspondence analysis indicated that other than for zooplankton abundance, total phosphorus was the most important environmental predictor of planktonic composition. The present study and others suggest that dealing with cyanobacteria removal using chitosan-modified soils can play an important role in controlling cyanobacterial blooms in eutrophicated freshwater systems.

Spatial distribution of benthic algae in the Gangqu River, Shangrila, China

This study consisted of sampling benthic algae at 32 sites in the Gangqu River, an important upstream tributary of the Yangtze River. Our aims were to characterize the benthic algae communities and relationships with environmental variables. Among the 162 taxa observed, Achnanthes linearis and Achnanthes lanceolata var. elliptica were the dominant species (17.10% and 14.30% of the total relative abundance, respectively). Major gradients and principal patterns of variation within the environmental variables were detected by principal component analysis (PCA). Then non-metric multidimensional scaling (NMS) divided all the sites into three groups, which were validated by multi-response permutation procedures (MRPP). Canonical correspondence analysis (CCA) indicated that three environmental variables (TN, TDS, and TP) significantly affected the distribution of benthic algae. Weighted averaging regression and cross-calibration produced strong models for predicting TN and TDS concentration, which enabled selection of algae taxa as potentially sensitive indicators of certain TN and TDS levels: for TN, Achnanthes lanceolata, Achnanthes lanceolata var. elliptica, and Cymbella ventricosa var. semicircularis; for TDS, Cocconeis placentula, Cymbella alpina var. minuta, and Fragilaria virescens. The present study represents an early step in establishing baseline conditions. Further monitoring is suggested to gain a better understanding of this region.

Response of macrozoobenthos communities to ecological engineering remediation in a hypertrophic urban lake

The relationship of macrozoobenthos communities with some environmental variables, and their response to the ongoing restoration measures were studied in a small hypertrophic urban lake near the Yangzte River, China. Twenty taxa including 9 oligochaetes, 7 insects, 2 mollusks and two other animals were found during March 2005 to May 2006. The reappearance of some indigenous macrozoobenthos species showed that the ecological engineering remediation carried out was helpful for the recovery of the macrozoobenthos communities. Through canonical correspondence analysis (CCA), it was detected that temperature (T), conductivity (COND), Secchi depth/deep (SD/Deep), TN, total suspended solids (SS) and chemical oxygen demand (CODcr) were significant environmental factors that influenced the pattern of macrozoobenthos. Limnodrilus hoffmeisteri, Tanypus sp. and Alocinma longicornis could be used as potential bio-indicators in monitoring the development of ecological restoration.

Controlling factors of spring-summer phytoplankton succession in Lake Taihu (Meiliang Bay, China)

The spring-summer successions of phytoplankton and crustacean zooplankton were examined weekly in Meiliang Bay of the subtropical Lake Taihu in 2004 and 2005. During the study period, the ecosystem of Meiliang Bay was characterized by (i) clearly declined nitrogen compounds (nitrate, TN, and ammonium) and slowly increased phosphorus compounds (TP and SRP), (ii) increased total phytoplankton density and rapid replacement of chlorophyta (mainly Ulothrix) by cyanobacteria (mainly Microcystis), and (iii) rapid replacement of large-sized crustaceans (Daphnia and Moina) by small-sized ones (Bosmina, Limnoithona, and Ceriodaphnia). Results from the CCA and correlation analysis indicate that the spring-summer phytoplankton succession was primarily controlled by abiotic factors. Cyanobacteria were mainly promoted by increased temperature and decreased concentrations of nitrogen compounds. The pure contribution of crustacean was low for the variation of phytoplankton suggesting a weak top-down control by crustacean zooplankton in the subtropical Lake Taihu.