Palaeoenvironmental changes from pollen record in deep sea core PC-1 from northern Okinawa Trough, East China Sea during the past 24 ka

A pollen record of core PC-1 from the northern Okinawa Trough, East China Sea (ECS), provides information on vegetation and climate changes since 24 cal. kaBP. A total of 103 samples were palynologically analyzed at 8 cm intervals with a time resolution of 230 a. Four pollen zones are recognized: zone I (812-715 cm, 24.2-21.1 cal. kaBP), zone II (715-451 cm, 21.1-15.2 cal. kaBP), zone III (451-251 cm, 15.2-10.8 cal. kaBP), zone IV (251-0 cm, 10.8-0.3 cal. kaBP), corresponding to Late MIS 3, Last Glacial Maximum (LGM), deglaciation and Holocene, respectively. The LGM is characterized by the dominance of herbs, mainly Artemisia, and high pollen influx, implying an open vegetation on the exposed continental shelf and a cool and dry climate. The deglaciation is a climate warming stage with Pinus percentage increased and Artemisia percentage decreased and a rapid sea-level rise. The Holocene is characterized by predominance of tree pollen with rapid increase in Castanea-Castanopsis indicating the development of mixed evergreen and deciduous broad-leaved forest and a warm, humid climate. Low pollen influx during the Holocene probably implies submergence of the continental shelf and retreat of the pollen source area. The vegetation indicated by pollen assemblage found in this upper zone is consistent with the present vegetation found in Kyushu, Japan. Originating from the humid mountain area of North Luzon of the Philippines, Tasmania and New Zealand, Phyllocladus with sporadic occurrence throughout PC-1 core probably suggests the influence of Palaeo-Kuroshio Current or intense summer monsoon. The observed changes in Pinus and Herbs percentage indicate fluctuations of the sea level, and high Pinus percentage corresponds to high sea level. Spectrum analysis of the pollen percentage record reveals many millennial-scale periodicities, such as periodicities of 6.8, 3.85 2.2, 1.6 ka.

History and variability of Asian interior aridity recorded by eolian flux in the Chinese Loess Plateau during the past 7 Ma

Eolian flux in the Chinese Loess Plateau was reconstructed by measuring the dry bulk density and CaCO3 content of the late Cenozoic loess-paleosol-red clay sequences in the Lingtai profile. Comparison of eolian flux variation between the Lingtai profile and the ODP sites 885/886 in the North Pacific shows a significant wet-dry variability in addition to a gradual drying trend in the dust source regions in interior Asia. Especially, the increase of eolian fluxes from both continental and pelagic eolian sediments indicates a sharp drying of the dust source regions between 3.6 and 2.6 MaBP, which might be attributed to the tectonic uplift of the Tibetan Plateau, which cut down the moisture input to the interior Asia. The average value and variability of eolian flux are higher after 2.6 MaBP than before, which may be related to the Quaternary climatic fluctuations on the glacial-interglacial timescale after the commencement of major Northern Hemisphere Glaciations. The eolian fluxes of the Lingtai profile and Core V21-146 in northwest Pacific show a synchronous variation on the 10(4)-10(5) a timescale, indicating that the flux variations from both continental and marine records are closely correlated to the Quaternary climatic fluctuation forced by the ice volume changes on a global scale.

Evolution and variability of the East Asian summer monsoon during the Pliocene: Evidence from clay mineral records of the South China Sea

The Late Pliocene is thought to be characterized by the simultaneous intensification of both the East Asian winter monsoon (EAWM) and East Asian summer monsoon (EASM). However, the evolution of the EASM during the Pliocene remains still controversial and only little is known about the dynamics of the EASM during the Pliocene on orbital time scales. Here we use clay mineral assemblages in sediments from Ocean Drilling Program (ODP) Site 1143 in the southern South China Sea (SCS) to obtain proxy records of past changes in the EASM climate during the Pliocene. Provenance analysis suggests that illite, chlorite and kaolinite originated mainly from the Mekong River drainage area. Smectite was derived mainly from the Indonesian islands. The kaolinite/illite ratio and the chemical index of alteration (CIA) of siliciclastic sediments allowed us to reconstruct the history of chemical weathering and physical erosion of the Mekong River drainage area and thus, the evolution of,the EASM during the Pliocene. Our clay minerals proxy data suggests a stronger EASM during the Early Pliocene than during the Late Pliocene. We propose that the long-term evolution of the EASM has been driven by global cooling rather than the uplift of the Tibetan Plateau. Spectral analysis of kaolinite/ illite ratio displays a set of strong periodicities at 100 ka, 30 ka, 28 ka, 25 ka, and 22 ka. with no clear obliquityrelated signal. Our study suggests that the Pliocene EASM intensity on orbital time scales is not only controlled by the Northern Hemisphere summer insolation, but also strongly influenced by equatorial Pacific ENSO-like ocean atmosphere dynamics. (C) 2010 Elsevier B.V. All rights reserved.

Coherence between solar activity and the East Asian winter monsoon variability in the past 8000 years from Yangtze River-derived mud in the East China Sea

AMS(14)C dating and grain-size analysis for Core PC-6, located in the middle of a mud area on the inner shelf of the East China Sea (ECS), were used to rebuild the Holocene history of the East Asian winter monsoon (EAWM). The 7.5-m core recorded the history of environmental changes during the postglacial transgression. The core's mud section (the upper 450 cm) has been formed mainly by suspended sediment delivered from the Yangtze River mouth by the ECS Winter Coastal Current (ECSWCC) since 7.6 kyr BP. Using a mathematical method called "grain size vs. standard deviatioW', we can divide the Core PC-6's grain-size distribution into two populations at about 28 mu m. The fine population (< 28 mu m) is considered to be transported by the ECSWCC as suspended loads. Content of the fine population changes little and represents a stable sedimentary environment in accord with the present situation. Thus, variation of mean grain-size from the fine population would reflect the strength of ECSWCC, which is mainly controlled by the East Asian winter monsoon. Abrupt increasing mean grain size in the mud section is inferred to be transported by sudden strengthened ECSWCC, which was caused by the strengthened EAWM. Thus, the high resolution mean grain-size variation might serve as a proxy for reconstruction of the EAWM. A good correlation between sunspot change and the mean grain-size of suspended fine population suggests that one of the primary controls on centennial- to decadal-scale changes of the EAWM in the past 8 ka is the variations of sun irradiance, i.e., the EAWM will increase in intensity when the number of sunspots decreases. Spectral analyses of the mean grain-size time series of Core PC-6 show statistically significant periodicities centering on 2463, 1368, 128, 106, 100, 88-91, 7678, and 70-72 years. The EAWM and the East Asian summer monsoon (EASM) agree with each other well on these cycles, and the East Asian Monsoon (EAM) and the Indian Monsoon also share in concurrent cycles in Holocene, which are in accord with the changes of the sun irradiance. (c) 2005 Elsevier B.V. All rights reserved.

Dynamics of intraseasonal sea level and thermocline variability in the equatorial Atlantic during 2002-03

Satellite and in situ observations in the equatorial Atlantic Ocean during 2002-03 show dominant spectral peaks at 40-60 days and secondary peaks at 10-40 days in sea level and thermocline within the intraseasonal period band (10-80 days). A detailed investigation of the dynamics of the intraseasonal variations is carried out using an ocean general circulation model, namely, the Hybrid Coordinate Ocean Model (HYCOM). Two parallel experiments are performed in the tropical Atlantic Ocean basin for the period 2000-03: one is forced by daily scatterometer winds from the Quick Scatterometer (QuikSCAT) satellite together with other forcing fields, and the other is forced by the low-passed 80-day version of the above fields. To help in understanding the role played by the wind-driven equatorial waves, a linear continuously stratified ocean model is also used. Within 3 degrees S-3 degrees N of the equatorial region, the strong 40-60-day sea surface height anomaly (SSHA) and thermocline variability result mainly from the first and second baroclinic modes equatorial Kelvin waves that are forced by intraseasonal zonal winds, with the second baroclinic mode playing a more important role. Sharp 40-50-day peaks of zonal and meridional winds appear in both the QuikSCAT and Pilot Research Moored Array in the Tropical Atlantic (PIRATA) data for the period 2002-03, and they are especially strong in 2002. Zonal wind anomaly in the central-western equatorial basin for the period 2000-06 is significantly correlated with SSHA across the equatorial basin, with simultaneous/ lag correlation ranging from-0.62 to 0.74 above 95% significance. Away from the equator (3 degrees-5 degrees N), however, sea level and thermocline variations in the 40-60-day band are caused largely by tropical instability waves (TIWs). On 10-40-day time scales and west of 10 degrees W, the spectral power of sea level and thermocline appears to be dominated by TIWs within 5 degrees S-5 degrees N of the equatorial region. The wind-driven circulation, however, also provides a significant contribution. Interestingly, east of 10 W, SSHA and thermocline variations at 10 40- day periods result almost entirely from wind-driven equatorial waves. During the boreal spring of 2002 when TIWs are weak, Kelvin waves dominate the SSHA across the equatorial basin (2 degrees S-2 degrees N). The observed quasi-biweekly Yanai waves are excited mainly by the quasi-biweekly meridional winds, and they contribute significantly to the SSHA and thermocline variations in 1 degrees-5 degrees N and 1 degrees-5 degrees S regions.

The variability of eddy kinetic energy in the South China Sea deduced from satellite altimeter data

We used fifteen years (1993-2007) of altimetric data, combined from different missions (ERS-1/2, TOPEX/Poseidon, Jason-1, and Envisat), to analyze the variability of the eddy kinetic energy (EKE) in the South China Sea (SCS). We found that the EKE ranged from 64 cm(2)/s(2) to 1 390 cm(2)/s(2) with a mean value of 314 cm(2)/s(2). The highest EKE center was observed to the east of Vietnam (with a mean value of 509 cm(2)/s(2)) and the second highest EKE region was located to the southwest of Taiwan Island (with a mean value of 319 cm(2)/s(2)). We also found that the EKE structure is the consequence of the superposition of different variability components. First, interannual variability is important in the SCS. Spectral analysis of the EKE interannual signal (IA-EKE) shows that the main periodicities of the IA-EKE to the east of Vietnam, to the southwest of Taiwan Island, and in the SCS are 3.75, 1.87, and 3.75 years, respectively. It is to the south of Taiwan Island that the IA-EKE signal has the most obvious impact on EKE variability. In addition, the IA-EKE exhibit different trends in different regions. An obvious positive trend is observed along the east coast of Vietnam, while a negative trend is found to the southwest of Taiwan Island and in the east basin of Vietnam. Correlation analysis shows that the IA-EKE has an obvious negative correlation with the SSTA in Nio3 (5A degrees S-5A degrees N, 90A degrees W-150A degrees W). El Nio-Southern Oscillation (ENSO) affects the IA-EKE variability in the SCS through an atmospheric bridge-wind stress curl over the SCS. Second, the seasonal cycle is the most obvious timescale affecting EKE variability. The locations of the most remarkable EKE seasonal variabilities in the SCS are to the east of Vietnam, to the southwest of Taiwan, and to the west of Philippines. To the east of Vietnam, the seasonal cycle is the dominant mechanism controlling EKE variability, which is attributed primarily to the annual cycle there of wind stress curl. In this area, the maximum EKE is observed in autumn. To the southwest of Taiwan Island, the EKE is enlarged by the stronger SCS circulation, which is caused by the intrusion branch from the Kuroshio in winter. Finally, intra-annual and mesoscale variability, although less important than the former, cannot be neglected. The most obvious intra-annual and mesoscale variability, which may be the result of baroclinic instability of the background flow, are observed to the southwest of Taiwan Island. Sporadic events can have an important effect on EKE variability.

Analysis of the sea-level variability along the Chinese coast and estimation of the impact of a CO2-perturbed atmospheric circulation

The relationship between monthly sea-level data measured at stations located along the Chinese coast and concurrent large-scale atmospheric forcing in the period 1960-1990 is examined. It is found that sea-level varies quite coherently along the whole coast, despite the geographical extension of the station set. A canonical correlation analysis between sea-level and sea-level pressure (SLP) indicates that a great part of the sea-level variability can be explained by the action of the wind stress on the ocean surface. The relationship between sea-level and sea-level pressure is analyzed separately for the summer and winter half-years. In winter, one factor affecting sea-level variability at all stations is the SLP contrast between the continent and the Pacific Ocean, hence the intensity of the winter Monsoon circulation. Another factor that affects coherently all stations is the intensity of the zonal circulation at mid-latitudes. In the summer half year, on the other hand, the influence of SLP on sea-level is spatially less coherent: the stations in the Yellow Sea are affected by a more localized circulation anomaly pattern, whereas the rest of the stations is more directly connected to the intensity of the zonal circulation. Based on this analysis, statistical models (different for summer and winter) to hindcast coastal sealevel anomalies from the large-scale SLP field are formulated. These models have been tested by fitting their internal parameters in a test period and reproducing reasonably the sea-level evolution in an independent period. These statistical models are also used to estimate the contribution of the changes of the atmospheric circulation on sea-level along the Chinese coast in an altered climate. For this purpose the ouput of 150 year-long experiment with the coupled ocean-atmosphere model ECHAM1-LSG has been analyzed, in which the atmospheric concentration of greenhouse gases was continuously increased from 1940 until 2090, according to the Scenario A projection of the Intergovermental Panel on Climate Change. In this experiment the meridional (zonal) circulation relevant for sea-level tends to become weaker (stronger) in the winter half year and stronger (weaker) in summer. The estimated contribution of this atmospheric circulation changes to coastal sea-level is of the order of a few centimeters at the end of the integration, being in winter negative in the Yellow Sea and positive in the China Sea with opposite signs in the summer half-year.

Intraseasonal variability of Indian Ocean sea surface temperature during boreal winter: Madden-Julian Oscillation versus submonthly forcing and processes

[ 1] Intraseasonal variability of Indian Ocean sea surface temperature (SST) during boreal winter is investigated by analyzing available data and a suite of solutions to an ocean general circulation model for 1998 - 2004. This period covers the QuikSCAT and Tropical Rainfall Measuring Mission (TRMM) observations. Impacts of the 30 - 90 day and 10 - 30 day atmospheric intraseasonal oscillations (ISOs) are examined separately, with the former dominated by the Madden-Julian Oscillation (MJO) and the latter dominated by convectively coupled Rossby and Kelvin waves. The maximum variation of intraseasonal SST occurs at 10 degrees S - 2 degrees S in the wintertime Intertropical Convergence Zone (ITCZ), where the mixed layer is thin and intraseasonal wind speed reaches its maximum. The observed maximum warming ( cooling) averaged over ( 60 degrees E - 85 degrees E, 10 degrees S - 3 degrees S) is 1.13 degrees C ( - 0.97 degrees C) for the period of interest, with a standard deviation of 0.39 degrees C in winter. This SST change is forced predominantly by the MJO. While the MJO causes a basin-wide cooling ( warming) in the ITCZ region, submonthly ISOs cause a more complex SST structure that propagates southwestward in the western-central basin and southeastward in the eastern ocean. On both the MJO and submonthly timescales, winds are the deterministic factor for the SST variability. Short-wave radiation generally plays a secondary role, and effects of precipitation are negligible. The dominant role of winds results roughly equally from wind speed and stress forcing. Wind speed affects SST by altering turbulent heat fluxes and entrainment cooling. Wind stress affects SST via several local and remote oceanic processes.

IMPACTS OF THE TROPICAL PACIFIC COUPLED PROCESS ON THE INTERANNUAL VARIABILITY IN THE INDIAN OCEAN

The basic features of climatology and interannual variations of tropical Pacific and Indian Oceans were analyzed using a coupled general circulation model (CGCM), which was constituted with an intermediate 2.5-layer ocean model and atmosphere model ECHAM4. The CGCM well captures the spatial and temporal structure of the Pacific El Nino-Southern Oscillation (ENSO) and the variability features in the tropical Indian Ocean. The influence of Pacific air-sea coupled process on the Indian Ocean variability was investigated carefully by conducting numerical experiments. Results show that the occurrence frequency of positive/negative Indian Ocean Dipole (IOD) event will decrease/increase with the presence/absence of the coupled process in the Pacific Ocean. Further analysis demonstrated that the air-sea coupled process in the Pacific Ocean affects the IOD variability mainly by influencing the zonal gradient of thermocline via modulating the background sea surface wind.

Variability of surface circulation in the South China Sea from satellite altimeter data

11-year satellite altimeter sea surface height (SSH) anomaly data from January 1993 to December 2003 are used to present the dominant spatial patterns and temporal variations of the South China Sea (SCS) surface circulation through Empirical Orthogonal Function (EOF) analysis. The first three EOF modes show the obvious seasonal variations of SSH in the SCS. EOF mode one is generally characterized by a basin-wide circulation. Mode two describes the double-cell basin scale circulation structure. The two cells were located off west of the Luzon Island and southeast of Vietnam, respectively. EOF mode three presents the mesoscale eddy structure in the western SCS, which develops into a strong cyclonic eddy rapidly from July to September. EOF mode one and mode three are also embedded with interannual signals, indicating that the SCS surface circulation variation is influenced by El Nino events prominently. The strong El Nino of 1997/98 obviously changed the SCS circulation structure. This study also shows that there existed a series of mesoscale eddies in the western SCS, and their temporal variation indicates intra-seasonal and interannual signals.

Seasonal and interannual variability of primary and export production in the South China Sea: a three-dimensional physical-biogeochemical model study

To investigate the seasonal and interannual variations in biological productivity in the South China Sea (SCS), a Pacific basin-wide physical - biogeochemical model has been developed and used to estimate the biological productivity and export flux in the SCS. The Pacific circulation model, based on the Regional Ocean Model Systems (ROMS), is forced with daily air-sea fluxes derived from the NCEP (National Centers for Environmental Prediction) reanalysis between 1990 and 2004. The biogeochemical processes are simulated with a carbon, Si(OH)(4), and nitrogen ecosystem (CoSiNE) model consisting of silicate, nitrate, ammonium, two phytoplankton groups (small phytoplankton and large phytoplankton), two zooplankton grazers (small micrograzers and large mesozooplankton), and two detritus pools. The ROMS-CoSiNE model favourably reproduces many of the observed features, such as ChI a, nutrients, and primary production (PP) in the SCS. The modelled depth-integrated PP over the euphotic zone (0-125 m) varies seasonally, with the highest value of 386 mg C m (-2) d (-1) during winter and the lowest value of 156 mg C m (-2) d (-1) during early summer. The annual mean value is 196 mg C m (-2) d (-1). The model-integrated annual mean new production (uptake of nitrate), in carbon units, is 64.4 mg C m (-2) d (-1) which yields an f-ratio of 0.33 for the entire SCS. The modelled export ratio (e-ratio: the ratio of export to PP) is 0.24 for the basin-wide SCS. The year-to-year variation of biological productivity in the SCS is weaker than the seasonal variation. The large phytoplankton group tends to dominate over the smaller phytoplankton group, and likely plays an important role in determining the interannual variability of primary and new production.

Phytoplankton biomass size structure and its regulation in the Southern Yellow Sea (China): Seasonal variability

Phytoplankton size structure plays a significant role in controlling the carbon flux of marine pelagic ecosystems. The mesoscale distribution and seasonal variation of total and size-fractionated phytoplankton biomass in surface waters. as measured by chlorophyll a (Chl a), was studied in the Southern Yellow Sea using data from four cruises during 2006-2007. The distribution of Chl a showed a high degree of spatial and temporal variation in the study area. Chl a concentrations were relatively high in the summer and autumn, with a mean of 142 and 1.27 mg m(-3), respectively. Conversely, in the winter and spring. the average Chl a levels were only 098 and 0.99 mg m(-3) Total Chl a showed a clear decreasing gradient from coastal areas to the open sea in the summer, autumn and winter cruises. Patches of high Chl a were observed in the central part of the Southern Yellow Sea in the spring due to the onset of the phytoplankton bloom. The eutrophic coastal waters contributed at least 68% of the total phytoplankton biomass in the surface layer. Picophytoplankton showed a consistent and absolute dominance in the central region of the Southern Yellow Sea (>40%) in all of the cruises, while the proportion of microphytoplankton was the highest in coastal waters The relative proportions of pico- and nanophytoplankton decreased with total biomass, whereas the proportion of the micro-fraction increased with total biomass. Relationships between phytoplankton biomass and environmental factors were also analysed. The results showed that the onset of the spring bloom was highly dependent on water column stability. Phytoplankton growth was limited by nutrient availability in the summer due to the strong thermocline. The combined effects of P-limitation and vertical mixing in the autumn restrained the further increase of phytoplankton biomass in the Surface layer. The low phytoplankton biomass in winter was caused by vertical dispersion due to intense mixing. Compared with the availability of nutrients. temperature did not seem to cause direct effects on phytoplankton biomass and its size structure. Although interactions of many different environmental factors affected phytoplankton distributions. hydrodynamic conditions seemed to be the dominant factor. Phytoplankton size structure was determined mainly by the size-differential capacity in acquiring resource. Short time scale events, such as the spring bloom and the extension of Yangtze River plume, can have substantial influences, both on the total Chl a concentration and on the size structure of the phytoplankton. (C) 2009 Elsevier Ltd. All rights reserved.

Temporal and spatial variability in nitrogen uptake kinetics during harmful dinoflagellate blooms in the East China Sea

Temporal and spatial variability in the kinetic parameters of uptake of nitrate (NO3-), ammonium (NH4+), urea, and glycine was measured during dinoflagellate blooms in Changjiang River estuary and East China Sea coast, 2005. Karenia mikimotoi was the dominant species in the early stage of the blooms and was succeeded by Prorocentrum donghaiense. The uptake of nitrogen (N) was determined using N-15 tracer techniques. The results of comparison kinetic parameters with ambient nutrients confirmed that different N forms were preferentially taken up during different stages of the bloom. NO3- (V-max 0.044 h(-1); K-s 60.8 mu M-N) was an important N source before it was depleted. NH4+ (V-max 0.049 h(-1); K-s 2.15 mu M-N) was generally the preferred N. Between the 2 organic N sources, urea was more preferred when K. mikimotoi dominated the bloom (V-max 0.020 h(-1); K-s 1.35 mu M-N) and glycine, considered as a dominant amino acid, was more preferred when P. donghaiense dominated the bloom (V-max 0.025 h(-1); K-s 1.76 mu M-N). The change of N uptake preference by the bloom-forming algae was also related to the variation in ambient N concentrations. Published by Elsevier B.V.

Variability in radiocarbon ages of biochemical compound classes of high molecular weight dissolved organic matter in estuaries

High molecular weight dissolved organic matter (HMW-DOM, > 1000 Da) represents a major fraction (> 30%) of dissolved organic carbon (DOC) in the ocean and thus plays an important role in the global biogeochemical cycling of carbon and many other elements. Its organic sources and formation mechanisms, however, are still not well understood especially in estuarine and coastal regions where multiple natural and anthropogenic sources contribute to total HMW-DOM. In this paper we report our measurements of natural radiocarbon (C-14) abundances and stable carbon isotope (C-13) compositions of the major biochemical compound classes: amino acids, carbohydrates and lipids separated from eight HMW-DOM samples collected from five US estuaries as part of our on-going study of sources, distribution and transport of chromophoric dissolved organic matter (CDOM) in estuarine and coastal waters. Distinct differences in both C-14 and C-13 values were found among the bulk HMW-DOM samples as well as the individual compound classes. Radiocarbon ages of the major compound classes varied by as much as 27,000 years in a single sample. The calculated average radiocarbon ages of the compound fractions of HMW-DOM indicate that the total lipid fraction is very "old", while the acid-insoluble fraction is slightly younger. Total amino acid and carbohydrate fractions, however, have relatively modern apparent C-14 ages. The significant variability in C-14 ages among the compound classes indicates not only multiple organic carbon sources but also different formation and turnover pathways controlling the cycling of different biochemical components of HMW-DOM in estuarine and coastal waters. (c) 2006 Elsevier Ltd. All rights reserved.

Spatial variability in biomass and production of heterotrophic bacteria in the East China Sea and the Yellow Sea

The distributions of heterotrophic bacterial abundance and production were investigated in the East China Sea and the Yellow Sea during the autumn of 2000 and spring of 2001. Bacterial abundance varied in the range 3.2-15.7 (averaging 5.7) x 10(5) and 2.3-13.6 (averaging 6.2) x 10(5) cells cm(-3) in the spring and autumn, respectively. During autumn, bacterial production (BP) (0.27-7.77 mg C m(-3) day(-1)) was on average 3 fold that in spring (0.001-2.04 mg C m(-3) day(-1)). Bacterial average turnover rate (ratio of bacterial production:bacterial biomass, mu=0.21 day(-1)) in autumn was 3 times as high as in spring (0.07 day(-1)). The ratio of integrated bacterial biomass to integrated phytoplankton biomass in the euphotic zone ranged from 4 to 101% (averaging 35%) in spring and 24 to 556% (averaging 121%) in autumn. The results indicate that the distributions of heterotrophic bacteria were controlled generally by temperature in spring and additionally by substrate supply in autumn. (C) 2010 Elsevier Ltd. All rights reserved.