SIMPLE SIMULATION OF THE ANNUAL VARIATION OF THE SPECIFIC PHOTOSYNTHESIS RATE IN JIAOZHOU BAY

A simple diagnostic simulation of the annual cycling of the surface specific photosynthesis rate (SPR) in Jiaazhou Bay is described in this paper. Light intensity, temperature and nutrients(nitrate ammonia, phosphate) were considered as main factors controlling photosynthesis of phytoplankton and were introduced into the model by different function equations. The simulated variation ofspecific photosynthesis rate coincided with the measured data. Analysis of the effect of every factor onphotosynthesis indicated that the variation of photosynthesis rate was controlled by all these three factors,while temperature showed good correlation with SPR as measurement showed. This diagnostic simulationyielded the values of some parameter relating with the photosynthesis in Jiaozhou Bay.     

EXAMINATION OF SILICATE LIMITATION OF PRIMARY PRODUCTION IN JIAOZHOU BAY, CHINA——I. SILICATE BEING A LIMITING FACTOR OF PHYTOPLANKTON PRIMARY PRODUCTION

Jiaozhou Bay data collected from May 1991 to February 1994, in 12 seasonal investigations, and provided the authors by the Ecological Station of Jiaozhou Bay, were analyzed to determine the spatiotemporal variations in temperature, light, nutrients (NO3^--N, NO2^--N, NH4^ -N, SIO3^2--Si, PO4^3--P), phytoplankton, and primary production in Jiaozhou Bay. The results indicated that only silicate correlated well in time and space with, and had important effects on, the characteristics, dynamic cycles and trends of, primary production in Jiaozhou Bay. The authors developed a corresponding dynamic model of primary production and silicate and water temperature. Eq. ( 1 ) of the model shows that the primary production variation is controlled by the nutrient Si and affected by water temperature; that the main factor controlling the primary production is Si; that water temperature affects the composition of the structure of phytoplankton assemblage; that the different populations of the phytoplankton assemblage occupy different ecological niches for C, the apparent ratio of conversion of silicate in seawater into phytoplankton biomas and D, the coefficient of water temperature‘s effect on phytoplankton biomass. The authors researched the silicon source of Jiaozhou Bay, the biogeochemical sediment process of the silicon, the phytoplankton predominant species and the phytoplankton structure. The authors considered silicate a limiting factor of primary production in Jiaozhou Bay, whose decreasing concentration of silicate from terrestrial source is supposedly due to dilution by current and uptake by phytoplankton; quantified the silicate assimilated by phytoplankton, the intrinsic ratio of conversion of silicon into phytoplankton biomass, the proportion of silicate uptaken by phytoplankton and diluted by current; and found that the primary production of the phytoplankton is determined by the quantity of the silicate assimilated by them. The phenomenon of apparently high plant-nutrient concentTations but low phytoplankton biomass in some waters is reasonably explained in this paper.     

Influence of Seawater Temperature on Phytoplankton Growth in Jiaozhou Bay, China

The phytoplankton reproduction capacity (PRC), as a new concept regarding chlorophyll-a and primary production (PP) is described. PRC is different from PP, carbon assimilation number (CAN) or photosynthetic rate ( P^B ) . PRC quantifies phytoplankton growth with a special consideration of the effect of seawater temperature. Observation data in Jiaozhou Bay, Qingdao, China, collected from May 1991 to February 1994 were used to analyze the horizontal distribution and seasonal variation of the PRC in Jiaozhou Bay in order to determine the characteristics, dynamic cycles and trends of phytoplankton growth in Jiaozhou Bay; and to develop a corresponding dynamic model of seawater temperature vs. PRC. Simulation curves showed that seawater temperature has a dual function of limiting and enhancing PRC. PRC‘s periodicity and fluctuation are similar to those of the seawater temperature. Nutrient silicon in Jiaozhou Bay satisfies phytoplankton growth from June 7 to November 3. When nutrients N, P and Si satisfy the phytoplankton growth and solar irradiation is sufficient, the PRC would reflect the influence of seawater temperature on phytoplankton growth. Moreover, the result quantitatively explains the scenario of one-peak or two-peak phytoplankton reproduction in Jiaozhou Bay, and also quantitatively elucidates the internal mechanism of the one- or two-peak phytoplankton reproduction in the global marine areas.     

COUPLED PHYSICAL-ECOLOGICAL MODELLING IN THE CENTRAL PART OF JIAOZHOU BAY Ⅱ.COUPLED WITH AN ECOLOGICAL MODEL

Sharpies‘ 1-D physical rrozlel maploying tide-wind driven turbulence closure and surface heating-cooling physics, was coupled with an eculogical rnodet with 9-biochemical components: phytoplankton, zooplankton, shellfish, autotmphic and heterotrophic bacterioplankton, dissolved organic carbon (DOC), suspended detritus and sinking particles to simulate the armual evolution of ecosystem in thecentral part of Jiaozhou Bay. The coupled modeling results showed that the phytoplankton shading effectcould reduce seawater temperamre by 2℃, so that photosynthesis efficiency should be less than 8% ; that the loss of phytoplankton by zooplankton grazing in winter tended to be compensated by phytoplankton advection and diffusion from the otrtside of the Bay; that the incidem irradiance intensity could be the mostimportant factor for phytoplankton grcr, wth rate; and that it was the bacterial secondary prnduction that maintained the maximum zooplankton biomass in winter usually observed in the 1990s, indicating that themicrobial food loop was extremely important for ecosystem study of Jiaozhou Bay.     

Silicon limitation on primary production and its destiny in Jiaozhou Bay,China Ⅷ: The variation of atmospheric carbon caused by both phytoplankton and human

Statistical analysis on data collected in the Jiaozhou Bay (Shandong, China) from May 1991 to February 1994 and those collected in Hawaii from March 1958 to December 2007 shows dynamic and cyclic changes in atmospheric carbon in the Northern Pacific Ocean (NPO), as well as the variation in space-time distribution of phytoplankton primary production and atmospheric carbon in the study regions. The study indicates that the human beings have imposed an important impact on the changing trends of the atmospheric carbon. Primary production in the Jiaozhou Bay presents a good example in this regard. In this paper, dynamic models of the atmospheric carbon in the NPO, the cyclic variations in the atmospheric carbon, and primary production in the Jiaozhou Bay are studied with simulation curves presented. A set of equations were established that able to calculate the rate and acceleration of increasing carbon discharged anthropologically into the atmosphere and the conversion rate of phytoplankton to atmospheric carbon. Our calculation shows that the amount of atmospheric carbon absorbed by one unit of primary production in the Jiaozhou Bay is (3.21–9.74)×10-9/(mgC·m-2d-1), and the amount of primary production consumed by a unit of atmospheric carbon is 102.66–311.52 (mgC·m-2d-1/10-6). Therefore, we consider that the variation of atmospheric carbon is a dynamic process controlled by the increase of carbon compound and its cyclic variation, and those from anthropologic discharge, and phytoplankton growth.     

Seasonal change in body length of important small copepods and relationship with environmental factors in Jiaozhou Bay,China

Differences among species in prosome length and in species’ response to environmental factors do exist. Therefore, it is useful to examine prosome length for different copepod species in variable environments. Seasonal variations in prosome length of four small copepods and their copepodite stages in the Jiaozhou Bay were compared and the relative influence of temperature, salinity, and chlorophyll concentration were examined. Two peaks were found in the mean prosome length of Paracalanus parvus (in early winter and May). For Acartia bifilosa, the maximum values of all copepodites occurred mainly from February to April, and decreased to the bottom in July. Prosome length of Acartia pacifica peaked when it first appeared in June, then reached to the minimum in July. Parvocalanus crassirostris only appeared from late summer to autumn and the mean prosome length showed no clear changes. Correlations of adult prosome length with environmental factors were evaluated. For the four species, temperature was negatively correlated to prosome length except for P. crassirostris. But the different species varied markedly in their responds to temperature. A. bifilosa showed a more definite trend of size variation with temperature than P. parvus and A. pacifica. Correlations of prosome length with salinity were significantly positive for almost all the small copepods. The relationship between chlorophyll concentration and prosome length was complicated for these copepods, but for P. parvus, chlorophyll concentration was also an important affecting factor. Furthermore, investigation needs to be done on food quality for some copepod. These results are essential to estimate the biomass and the production, and to understand these small copepods’ population dynamics in this human-affected bay.     

Pelagic Nitrogen Cycling in Jiaozhou Bay, a Model Study I： the Conceptual Model

A zero-dimensional box model (PNCMjzb) with six state variables (ammonium, nitrate, dissolved organic nitrogen, phytoplankton, zooplankton and detritus) was developed to study nitrogen cycling in the Jiaozhou Bay pelagic ecosystem. The dominant processes within these compartments are considered with nitrogen as flow currency. Phytoplankton and zooplankton are treated as separate state variables, assuming that the species composition was dominated by two or three species the dynamic constants of which are similar and that they represent the entire plankton community. The microbial loop has not been integrated explicitly in the model. The turnover of bacteria is included implicitly in processes such as detritus decomposition, DON remineralization, pelagic nitrification and denitrification. The model is driven by two forcing variables, viz. water temperature and light intensity. Historical data from the1980s and 1990s were compiled and used for model calibration. In this paper (part Ⅰ), the consideration of every main compartment in the model is interpreted in detail. And the applied equations and parameters are presented. The main results from the simulations together with discussion about phytoplankton dynamics and primary production in Jiaozhou Bay are presented in the next paper (part Ⅱ).     

Steady increase in water clarity in Jiaozhou Bay in the Yellow Sea from 2000 to 2018: Observations from MODIS

The Moderate Resolution Imaging Spectroradiometer(MODIS) surface reflectance data were used to analyze the temporal and spatial distribution characteristics of water clarity( Z_(sd)) in the Jiaozhou Bay, Qingdao, China, in the Yellow Sea from 2000 to 2018. Z_(sd)retrieval models were regionally optimized using in-situ data with coincident MODIS images, and then were used to retrieve the Z_(sd) products in Jiaozhou Bay from 2000–2018. The analysis of the Z_(sd) results suggests that the spatial distribution of relative Z_(sd) spatial characteristics in Jiaozhou Bay was stable, being higher Z_(sd) in the southeast and a lower Z_(sd) in the northwest. The annual mean Z_(sd) in Jiaozhou Bay showed a significant upward trend, with an annual increase of approximately 0.02 m. Water depth and wind speed were important factors affecting the spatial distribution and annual variation of Z_(sd) in Jiaozhou Bay, respectively.     

Sensitivity analysis of thermal equilibrium parameters of MIKE 11 model: A case study of Wuxikou Reservoir in Jiangxi Province of China

Sensitivity analysis of thermal equilibrium parameters in the reservoir module of MIKE 11 model was conducted for the Wuxikou Reservoir in Jiangxi Province of China in order to apply the module to the environmental impact assessment to accurately predict water temperature of reservoirs.Results showed that radiation parameter A and evaporation-first parameter were much more sensitive than other parameters.The values of the radiation parameter A ranged from 0.10 to 0.34.The values of evaporation-first parameter varied from 0 to 10.The sensitivity of solar absorption parameters was less than that of evaporation parameter,of which light attenuation values ranged from 0.5 to 0.7,and this parameter would not impact model results if it was more than 2.Constants in Beer’s law ranged from 0.2 to 0.7.Radiation parameter B was not more sensitive than evaporation parameter and its reasonable range was higher than 0.48.The fitting curves showed consistent changing tendency for these parameters within the reasonable ranges.Additionally,all the thermal equilibrium parameters had much more important effects on surface water temperature than deep water temperature.Moreover,if no observed data could be obtained,the local empirical value would be used to input to the MIKE 11 model to simulate the changes in the discharged outflow-water temperature qualitatively.     

Seasonal variations of group-specific phytoplankton cell death in Xiamen Bay,China

The importance of phytoplankton cell death is being increasingly recognized,however,there are still no published reports on this in Xiamen Bay.In this study,the proportion of dead phytoplankton cells associated with environmental factors was assessed at a station in Xiamen Bay from December2012 to December 2013,using a cell digestion assay,which is an effective method to analyze dead/living cells in complex natural phytoplankton communities.The percentages of dead cells(%DC) in the total phytoplankton in summer(16%±6%) were lower than those in winter(27%±16%).Six groups of phytoplankton(G1-G6) were categorized by flow cytometry.These phytoplankton communities with diverse seasonal variations in%DC had different responses to environmental constraints.The main factors affecting mortality were temperature and salinity,while nutrient concentration showed little influence on phytoplankton death.Additionally,our results provide evidence that chlorophyll a concentrations had an inverse relationship with total phytoplankton%DC and viable cell abundance was more meaningful than total cells to explain variations in environmental parameters(such as Chl a).Moreover,the lowest mean%DC in total phytoplankton was 16%±6%at our sample site,which is in a subtropical area with high water temperatures,full solar radiation,and rich nutrients.This indicates that phytoplankton cell death is a process that cannot be ignored.In summary,phytoplankton cell death is important in understanding the dynamics of phytoplankton communities and the functioning of subtropical ecosystems.     

Annual Occurrence of Calanus sinicus Carcasses in the Jiaozhou Bay and the Adjacent Waters

The dead proportion of a calanoid copepod Calanus sinicus population was investigated monthly with neutral red staining method at 12 stations in the Jiaozhou Bay and the adjacent waters from December 2008 to November 2009. C. sinicus could be observed through the whole year, with an evident numerical peak from February to April and an extremely low abundance in summer. Annual variation of the dead proportion differed greatly between adult females and copepodites. Dead females were observed only in the southern part and outside the bay sporadically in five months, when the egg production was the most active. Dead proportion of females was usually less than 5%. Dead copepodites could be observed in all but three months in summer in the whole study area, and the proportion varied between 3.8% – 18.2%. Death percentage of copepodites increased significantly in the northern part in January and June, when water temperature exceeded it's favorable limits; however, the fluctuated salinity had no significant effect on the survival of both females and copepodites. Thus the dead proportion of C. sinicus was relatively low in the Jiaozhou Bay, while difference could be observed between females and copepodites. Dead copepodites were common in all samples, but dead females were presented mainly as ‘reproduction cost' during generation alternation.     

Spatial and seasonal variations in the trophic spectrum of demersal fish assemblages in Jiaozhou Bay,China

Trophic structure of fish communities is fundamental for ecosystem-based fisheries management, and trophic spectrum classifies fishes by their positions in food web, which provides a simple summary on the trophic structure and ecosystem function. In this study, both fish biomass and abundance trophic spectra were constructed to study the spatial and seasonal variations in the trophic structure of demersal fish assemblages in Jiaozhou Bay, China. Data were collected from four seasonal bottom trawl surveys in Jiaozhou Bay from February to November in 2011. Trophic levels(TLs) of fishes were determined by nitrogen stable isotope analysis. This study indicated that most of these trophic spectra had a single peak at trophic level(TL) of 3.4–3.7, suggesting that demersal fish assemblages of Jiaozhou Bay were dominated by secondary consumers(eg. Pholis fangi and Amblychaeturichthys hexanema). The spatial and seasonal variations of trophic spectra of Jiaozhou Bay reflected the influence of fish reproduction, fishing pressure and migration of fishes. Two-way analysis of variance(ANOVA) showed that seasonal variations in trophic spectra in Jiaozhou Bay were significant(P 0.05), but variations among different areas were not significant( P 0.05). The trophic spectrum has been proved to be a useful tool to monitor the trophic structure of fish assemblages. This study highlighted the comprehensive application of fish biomass and abundance trophic spectra in the study on trophic structure of fish assemblages.     

Numerical Simulation of Nutrient and Phytoplankton Dynamics in Guangxi Coastal Bays,China

The increasing riverine pollutants have resulted in nutrient enrichment and deterioration of water quality in the coastal water of Guangxi Province, China. However, the quantitative relationship between nutrient loads and water quality responses, which is crucial for developing eutrophication control strategies, is not well studied. In this study, the riverine fluxes of nutrients were quan- tified and integrated with nutrient cycling and phytoplankton dynamics by using box models for Guangxi coastal bays. The model concepts and biogeochemical equations were the same; while most model parameters were specific for each bay. The parameters were calibrated with seasonal observations during 2006--2007, and validated with yearly averaged measurements in 2009. The gen-eral features of nutrient and phytoplankton dynamics were reproduced, and the models were proved feasible under a wide range of bay conditions. Dissolved inorganic nitrogen was depleted during the spring algal bloom in Zhenzhu Bay and Fangcheng Bay with relatively less nutrient inputs. Phosphorus concentration was high in spring, which decreased then due to continuous phytoplankton consumption. Chlorophyll-a concentration reached its annual maximum in summer, but was the minimum in winter. Eutrophication was characterized by both an increase in nutrient concentrations and phytoplankton biomass in Lianzhou Bay. Either about 80% re-duction of nitrogen or 70% reduction of phosphorus was required to control the algal bloom in Lianzhou Bay. Defects of the models were discussed and suggestions to the environmental protection of Guangxi coastal bays were proposed.     

Numerical simulation of nutrient and phytoplankton dynamics in Guangxi coastal bays,China

The increasing riverine pollutants have resulted in nutrient enrichment and deterioration of water quality in the coastal water of Guangxi Province, China. However, the quantitative relationship between nutrient loads and water quality responses, which is crucial for developing eutrophication control strategies, is not well studied. In this study, the riverine fluxes of nutrients were quantified and integrated with nutrient cycling and phytoplankton dynamics by using box models for Guangxi coastal bays. The model concepts and biogeochemical equations were the same; while most model parameters were specific for each bay. The parameters were calibrated with seasonal observations during 2006–2007, and validated with yearly averaged measurements in 2009. The general features of nutrient and phytoplankton dynamics were reproduced, and the models were proved feasible under a wide range of bay conditions. Dissolved inorganic nitrogen was depleted during the spring algal bloom in Zhenzhu Bay and Fangcheng Bay with relatively less nutrient inputs. Phosphorus concentration was high in spring, which decreased then due to continuous phytoplankton consumption. Chlorophyll-a concentration reached its annual maximum in summer, but was the minimum in winter. Eutrophication was characterized by both an increase in nutrient concentrations and phytoplankton biomass in Lianzhou Bay. Either about 80% reduction of nitrogen or 70% reduction of phosphorus was required to control the algal bloom in Lianzhou Bay. Defects of the models were discussed and suggestions to the environmental protection of Guangxi coastal bays were proposed.     

REPRODUCTION ECOLOGY OF FISHES IN JIAOZHOU BAY WATERS

Observations on gonad development of 84 species of marine fishes in Jiaozhou Bay waters showed the breeding periods of more than 50 of these species are mostly from April to August, with peaks between May to June. A great many of the species spawn in the warm season, only a few species in the cold season. The number of spawning species increases with the rise and decreases with the fall in water temperature. The absolute fecundity of 50 teleost species, which usually increases with body length and weight, was estimated to range from about 300 to 9,000,000 eggs. The reproductive patterns of fishes in Jiaozhou Bay waters are discussed.     

Autumn photoproduction of carbon monoxide in Jiaozhou Bay,China

Carbon monoxide(CO) plays a significant role in global warming and atmospheric chemistry. Global oceans are net natural sources of atmospheric CO. CO at surface ocean is primarily produced from the photochemical degradation of chromophoric dissolved organic matter(CDOM). In this study, the effects of photobleaching, temperature and the origin(terrestrial or marine) of CDOM on the apparent quantum yields(AQY) of CO were studied for seawater samples collected from Jiaozhou Bay. Our results demonstrat that photobleaching, temperature and the origin of CDOM strongly affected the efficiency of CO photoproduction. The concentration, absorbance and fluorescence of CDOM exponentially decreased with increasing light dose. Terrestrial riverine organic matter could be more prone to photodegradation than the marine algae-derived one. The relationships between CO AQY and the dissolved organic carbon-specific absorption coefficient at 254 nm for the photobleaching study were nonlinear, whereas those of the original samples were strongly linear. This suggests that: 1) terrestrial riverine CDOM was more efficient than marine algae-derived CDOM for CO photoproduction; 2) aromatic and olefinic moieties of the CDOM pool were affected more strongly by degradation processes than by aliphatic ones. Water temperature and the origin of CDOM strongly affected the efficiency of CO photoproduction. The photoproduction rate of CO in autumn was estimated to be 31.98 μmol m-2 d-1 and the total DOC photomineralization was equivalent to 3.25%- 6.35% of primary production in Jiaozhou Bay. Our results indicate that CO photochemistry in coastal areas is important for oceanic carbon cycle.     

Effects of spatio-temporal and environmental factors on distribution and abundance of wintering anchovy Engraulis japonicus in central and southern Yellow Sea

We investigated the spatio-temporal and environmental factors that affect the distribution and abundance of wintering anchovy and quantifi es the infl uences of these factors. Generalized additive models(GAMs) were developed to examine the variation in species distribution and abundance with a set of spatiotemporal and oceanographic factors, using data collected by bottom trawl surveys and remote sensing in the central and southern Yellow Sea during 2000–2011. The fi nal model accounted for 28.21% and 41.03% of the variance in anchovy distribution and abundance, respectively. The results of a two-step GAM showed that hour, longitude, latitude, temperature gradient(TGR), and chlorophyll a(Chl- a) concentration best explained the anchovy distribution(presence/absence) and that a model including year, longitude, latitude, depth, sea surface temperature(SST), and TGR best described anchovy abundance(given presence). Longitude and latitude were the most important factors affecting both distribution and abundance, but the area of high abundance tended to be east and south of the area where anchovy were most likely to be present. Hour had a signifi cant effect on distribution, but year was more important for anchovy abundance, indicating that the anchovy catch ratio varied across the day but abundance had an apparent interannual variation. With respect to environmental factors, TGR and Chl- a concentration had effects on distribution, while depth, SST, and TGR affected abundance. Changes in SST between two successive years or between any year and the 2000–2011 mean were not associated with changes in anchovy distribution or abundance. This fi nding indicated that short- and long-term water temperature changes during 2000–2011 were not of suffi cient magnitude to give rise to variation in wintering anchovy distribution or abundance in the study area. The results of this study have important implications for fi sheries management.     

Composition of stable isotope in precipitation and its influences by different vapor sources in the eastern Qilian Mountains

To better understand the process of precipitation and water cycle, the composition of stable isotope in precipitation and its influences by different vapor sources in the eastern of Qilian Mountains were conducted from June 2013 to May 2014. The total of 100 precipitation samples were collected in Wushaoling national meteorological station located in the eastern of Qilian Mountains. The analysis indicates that the slope of Local Meteoric Water Line is lower than that of Global Meteoric Water Line. The average values of δ18 O and δD in precipitation are higher in summer but lower in winter. Except for negative correlation with relative humidity, the stable isotope values in precipitation are positive correlations with temperature, precipitation and water vapor pressure. Influenced by water vapor source, the values of d-excess are lower for the Westerly wind and the South Asia Monsoon onJuly and the Westerly wind and the East Asia Monsoon on August, but they are higher for the Westerly wind on other months, that they are also influenced by the weather conditions in rainfall process. The variation of stable isotope in precipitation exhibited significant temperature effect, and there is also some precipitation amount effect in spring and summer.     

Carbon biomass,carbon-to-chlorophyll a ratio and the growth rate of phytoplankton in Jiaozhou Bay,China

Carbon biomass,carbon-to-chlorophyll a ratio(C:Chl a),and the growth rate of phytoplankton cells were studied during four seasonal cruises in 2017 and 2018 in Jiaozhou B ay,China.Water samples were collected from 12 stations,and phytoplankton carbon biomass(phyto-C) was estimated from microscopemeasured cell volumes.The phyto-C ranged from 5.05 to 78.52 μg C/L in the bay,and it constituted a mean of 38.16% of the total particulate organic carbon in the bay.High phyto-C values appeared mostly in the northern or northeastern bay.Diatom carbon was predominant during all four cruises.Dinoflagellate carbon contributed much less(30%) to the total phyto-C,and high values appeared often in the outer bay.The C:Chl a of phytoplankton cells varied from 11.50 to 61.45(mean 3 1.66),and high values appeared in the outer bay during all four seasons.The phyto-C was also used to calculate the intrinsic growth rates of phytoplankton cells in the bay,and phytoplankton growth rates ranged from 0.56 to 1.96/d;the rate was highest in summer(mean 1.79/d),followed by that in fall(mean 1.24/d) and spring(mean 1.17/d),and the rate was lowest in winter(mean 0.77/d).Temperature and silicate concentration were found to be the determining factors of phytoplankton growth rates in the bay.To our knowledge,this study is the first report on phytoplankton carbon biomass and C:Chl a based on water samples in Jiaozhou B ay,and it will provide useful information for studies on carbon-based food web calculations and carbon-based ecosystem models in the bay.     

Pb distribution and translocation in Jiaozhou Bay

The trends of distribution, translocation and seasonal change of heavy metal Pb were studied based on the surface and bottom water sampling in Jiaozhou Bay in 1979, and compared with those in 1990＇s. The results showed that the source of Pb in the bay was from wastewater and sewage in the east of Jiaozhou Bay from ocean vessels. Pb concentration was higher in spring and lower in summer and autumn, and remained stable through sedimentation in the bottom layer. The overall water quality was good in 1970＇s. Compared with the environmental monitoring data of 1995-1999, Pb pollution had become serious. Therefore, more efforts should be made to protect the bay from Pb pollution.