A Cancer Recognition Method Based on DNA Microarray

The accurate cancer classification is of great importance in clinical treatment. Recently, the DNA microarray technology provides a promising approach to the diagnosis and prognosis of cancer types. However, it has no perfect method for the multiclass classification problem. The difficulty lies in the fact that the data are of high dimensionality with small sample size. This paper proposed an automatic classification method of multiclass cancers based on Biomimetic pattern recognition (BPR). To the public GCM data set, the average correct classification rate reaches 80% under the condition that the correct rejection rate is 81%.

隆肛蛙属Feirana种群的系统学与动物地理学研究

对隆肛蛙属的物种构成进行了订正，建立新属肛刺蛙属Yerana gen. nov.；订正后的隆肛蛙属现仅隶2种， 即隆肛蛙F. quadrana和太行隆肛蛙F. taihangnicus。运用形态学分析探讨了隆肛蛙属物种及种群的形态差异和分类关系，通过分子系统学研究探讨了隆肛蛙属物种及种群的分类和系统发育关系，运用动物地理学方法结合系统发育关系探讨了隆肛蛙属种群的地理分布格局成因与历史过程。主要结果和推论如下： 1．隆肛蛙属物种构成的订正及一新属建立 建立新属肛刺蛙属，将隆肛蛙属中的原叶氏隆肛蛙F. yei归隶新属肛刺蛙属并更名为叶氏肛刺蛙Y. yei，，新属建立的主要依据为：（1）雄性肛部隆起，肛孔下方有两个布满黑刺的大的白色球形隆起，具单咽下内声囊， 第一指具婚刺；（2）形态量度分析表明叶氏肛刺蛙与隆肛蛙和太行隆肛蛙的形态差异远大于后两者之间的差异；（3）叶氏肛刺蛙的分布区与隆肛蛙和太行隆肛蛙的分布区距离较远且呈隔离状态；（4）分子系统学研究资料（Jiang et al.,2005）证明叶氏肛刺蛙与隆肛蛙和太行隆肛蛙非单系发生；叶氏肛刺蛙在第二支中位于基部。因此，隆肛蛙属现仅隶2种，即隆肛蛙和太行隆肛蛙。 2．隆肛蛙属种群形态学研究 对隆肛蛙属中隆肛蛙和太行隆肛蛙的15个地理种群565只标本的28项形态性状进行了测量，运用典型判别分析法对其分析的结果表明：（1）太行隆肛蛙与隆肛蛙形态差异明显，支持其为不同的物种；（2）原隆肛蛙河南伏牛山种群和山西中条山种群应为太行隆肛蛙的地理种群；（3）隆肛蛙不同地理种群之间形态差异明显，其中四川安县种群、陕西周至种群和湖北利川种群与模式产地重庆巫山种群的差异可能达到了亚种或亚种以上分化水平。对隆肛蛙属量度分析的15个种群进行定性形态分析表明其分为三种形态型，对应隆肛蛙、过渡型和太行隆肛蛙，其变异特征主要为内跗褶、雄性肛部隆起及疣粒分布、第五趾外侧缘膜等，这与量度分析结果相似。 3．隆肛蛙属种群分子系统学研究 测定隆肛蛙属Feirana的2种19种群的线粒体12S rRNA和16S rRNA基因片段、ND2基因的DNA序列，比对后共计1953bps。（1）遗传多样性与距离分析：结果表明，隆肛蛙属种群具很高的遗传多样性，19个种群样品表现出19种单倍型（遗传多样性指数Hd＝1.0）； ND2基因的进化信息含量远高于12SrRNA和16SrRNA。隆肛蛙属2种群组内的种群间的遗传距离远小于两种群组间的距离，种群在不同基因上的遗传距离表现的关系与对应的系统树一致。（2）系统发育关系分析：结果表明，不同基因片断基于不同方法构建的隆肛蛙属种群系统发育树结构基本一致，基本表明隆肛蛙属种群为单系发生；它们在系统树中分为两大支，分别对应于隆肛蛙和太行隆肛蛙；支持中条山种群（沁水、历山和济源种群）和伏牛山种群（栾川和内乡种群）为太行隆肛蛙的地理种群，而原隆肛蛙秦岭中东段的部分种群（柞水、宁陕、长安大坝沟种群）也应为太行隆肛蛙的地理种群。（3）亚种分化分析：根据遗传距离分析和系统发育关系分析结果，并考虑形态上的差异情况以及地理分布信息，隆肛蛙所隶种群组可分为2亚种，即隆肛蛙指名亚种F. quadrana quadrana包括四川盆地东缘大巴山东段-巫山-武陵山北麓种群和秦岭中段（周至板房子和长安广货街）种群，他们在系统关系树上聚为一支；安县亚种F. quadrana anxianensis包括四川盆地西缘岷山东麓-龙门山-大巴山和秦岭西段的种群（安县、青川、文县、南江和凤县种群），他们在系统关系树上聚为一支。太行隆肛蛙所隶种群组也可分为2亚种，即太行隆肛蛙指名亚种F. taihangnicus taihangnicus包括中条山的种群（沁水、历山和济源种群）和中东秦岭的部分种群（柞水、长安大坝沟和宁陕种群），他们在系统关系树上聚为一支；太行隆肛蛙伏牛亚种F. taihangnicus funiuensis，为伏牛山地区的种群（栾川和内乡种群），他们在系统关系树上聚为一支。 4．隆肛蛙属种群动物地理学研究 隆肛蛙属19种群的分歧年代分析： 以长江巫山段和黄河三门峡段的形成历史时期为参考点，根据已测隆肛蛙属19种群及其外群包括N. pleski、P. yunnanesis、P. robertingeri、F. limnocharis的1953bps DNA序列构建分子钟，获得各支系的分歧年代。结果表明：①棘蛙族在70Ma左右开始其独立演化历程，这与Roelants et al.（2004）的分析结果～60±15Ma左右开始分化基本一致，后者印证了本文的分子钟。②隆肛蛙属的起始分化年代较早，隆肛蛙和太行隆肛蛙两种群组的最近祖先种群大概在46Ma～50Ma左右；隆肛蛙和太行隆肛蛙种群组内的种群分化年代相对两种群组间晚得多， 隆肛蛙种群组内两亚种分化起始年代约为10Ma左右，而太行隆肛蛙种群组内两亚种分化起始年代约为6Ma。 隆肛蛙属种群分布格局形成过程分析： ①隆肛蛙属的系统关系与地理分布格局密切相关，大部分系统分支分级与地理距离成正比；②隆肛蛙属最近祖先种群的分化中心可能位于秦岭中部地区， 隆肛蛙属的种群分布格局的形成表现为隔离分化与扩散相结合的机制，由隔离分化产生的隆肛蛙祖先种群主要从秦岭中部向西南方向扩散，后隔离分化为两亚种；太行隆肛蛙祖先种群向东北方向扩散也分化为两亚种。 隆肛蛙属种群分布区域地质历史的探讨：本文所建分子钟和种群分化方式印证了该区域的几次主要地质事件，包括岷山-龙门山-西秦岭等地区的快速差异隆起、第四纪冰期等。 The specific composition of the genus Feirana should be revised. A new genus Yerana gen. nov.（Ranidae：Dicroglossinae）was established based on morphological data-set and molecular phylogeny, as a result, only two species F. quadrana and F. taihangnicus are classified into Feirana now. Morphological differences and taxonomy of populations of Feirana were investigated based on morphological and morphometric data; phylogenetic relationships and taxonomy of populations of Feirana were elucidated using molecular data, and then the proceeding of the distribution pattern of populations of Feirana were discussed. The main results and conclusions and proposals were presented as following: 1. Revising of the specific composition of the genus Feirana and establishment of a new genus The new genus Yerana, only containing the type species Y. yei, was established based on the following evidences: (1) In adult male, distinct up-heaved circular vesicle presents around the anal, and under anal there are two white balls on which black spines exist, black horny spines scatter on the upper side of first finger, and internal single subgular vocal sac presents; (2) there is obvious morphometric differences between Yerana and Feirana; (3) Yerana is distributed far from Feirana; (4) evidences of molecular phylogeny（Jiang et al.,2005）suggested that Yerana take a special phylogenetic clade which is different from other genus included in the tribe Paini. As a result, there are only two species in Feirana, i.e., F. quadrana and F. taihangnicus. 2. Morphological research of populations of Feirana Twenty-eight characters of 565 individuals of 15 populations of the genus Feirana were measured, the results of Canonical Discriminant analysis of the morphometric data-set indicated that: (1) there are very prominent differences between the two species F. quadrana and F. taihangnicus. The validity of species F. taihangnicus was approved here; (2) Mt. Funiu population and Mt. Zhongtiao population should belong to the species F. taihangnicus; (3) Obvious differences exist among 12 populations of F. quadrana, the differentiation among Zhouzhi population, Anxian population, Lichuan population, and Wushan population together with the others probably reach sub-specific or specific level. Result of morphological comparison between 15 different populations show that 3 morphological types are recogenized in according with F. quadrana, F. taihangnicus and intergradation, this result conform to the result of morphometric analysis. 3. Molecular phylogenetic study on populaions of Feirana Fragment of 12SrRNA and 16SrRNA genes, and ND2 gene of 19 populations of two species of Feirana were sequenced and aligned, from which 1953 bps were received. (1) analyses of genetic distance and hereditary diversity indicated that: genetic distance between populations in each group were less than distance between two groups of Feirana, 19 haplotypes were recognized from 19 samples of 19 populations, so the hereditary diversity of populations of Feirana was very high (Hd=1.0), phylogenetic information in ND2 gene is more than fragment sequence of 12SrRNA and 16SrRNA genes. (2) Result of molecular phylogeny indicate that the phylogenetic trees constructed using different methods based on different sequence data sets showed the revised genus Feirana is monophyletic since the 19 populations of Feirana were firstly clustered together as one large clade, which was further clustered into two major clades, corresponding to F. quadrana(GroupⅠ) and F. taihangnicus(GroupⅡ), respectively. So populations of Qinshui and Lishan in Mt. Zhongtiao, populations of Luanchuan and Neixiang in Mt. Funiu, and populations of Zhashui, Dabagou of Chang’an and Ningshan in eastern Mt. Qinling should belong to the species F. taihangnicus; (3) Subspecific differentiation. on the basis of genetic distance, phylogenetic trees and geographical distribution, F. quadrana should have two subspecies, i.e., F. quadrana qudadrana, consisting of the populations Guanghuojie of Chang’an and Zhouzhi in Mid-Mt. Qinling, populations in Wushan area and northern Mt. Wuling (Lichuan), and F. qudadrana anxianensis, consisting of the populations in eastern Mt. Ming shan-Mt. Longmen-western Mt. Daba-western Mt. Qinling (Anxian, Qingchuan, Wenxian, Nanjiang and Fengxian); F. taihangnicus should also has two subspecies, i.e., F. taihangnicus taihangnicus, consisting of the populations in Mt. Zhongtiao and eastern Mt. Qinling, and F. taihangnicus funiuensis, consisting of the populations in Mt. Funiu. 4. Zoogeography of populaions of Feirana Analysis for divergent time of 19 populations of Feirana: Using the dates of run-through of Wushan segment of Changjiang River as the time when the population of Lichuan started differentiated from the populations of Wushan and Shennongjia, and the dates of Sanmenxia segment of Yellow River as the time when the populations in Mt. Zhongtiao started differentiated from the population of Dabagou in Chang’an, molecular clock was established using sequences with 1953 bps of 19 populations of Feirana and outgroup including N. pleski, P. yunnanesis, P. robertingeri, F. limnocharis in order to estimate divergent time of all clades. Result of that indicated that: ① the tribe Paini started to evolve independently at about 70Ma when is in consistent with that estimated by Roelants et al.（2004）with result of about ～60±15Ma, they were corroborated by each other, this confirms the validity of this molecular clock; ② divergent time for speciation of Feriana is early, ancestral populations of F. quadrana and F. taihangnicus were found about 46Ma～50Ma; differentiation of populations within species is greatly late to the divergence of the two species, divergent time for F. quadrana is 10Ma and divergent time for F. taihangnicus is 6Ma. Proceeding of distribution pattern of Feirana. Phylogenetic relationships of populations of Feirana matched quite with distribution pattern of them, the relationships among clades showed in phylogenetic trees is direct ratio to geographical distance of them; the estimated date of speciation between two species of Feirana was as early as speciation of Paa yunnanesis and Nanara pleski; middle part of Mt. Qinling is the center of speciation of Feirana, combination of mult-events of dispersal and vicariance are probably the mechanism of speciation of Feirana, F. quadrana colonized the mid-Mt. Qinling and then differentiated into two subspecies in southwest direction, ancestral population of F. taihangnicus colonized the mid-Mt. Qinling and then differentiated into two subspecies in northeast direction. On geological history of the distribution of Feirana. According to molecular clock and speciation model of populations of Feirana, some geological events are confirmed, including special rise of Mt. Minshan- Mt. Longmen-western Mt. Qinling, glacial age.

区域水土流失研究中降雨量空间插值方法选取

在区域水土流失模型研究中,空间插值可提供每个计算栅格的气象要素资料。考虑到研究区域降雨与高程相关性很弱,不宜采用梯度距离反比法(GIDS),故采用距离反比法(IDW)和普通克里格法(Kriging),对延安示范区及其周围共50个站点2000—2003年的5—10月逐月降雨量进行插值。交叉验证结果表明:对2种插值方法,二者经过对数变换后平均相对误差(MRE)为8.30%和7.67%,分别比原始数据插值后的MRE下降了23.17%和23.50%,说明插值精度得到了提升,对研究区域某一年逐月降水的插值Kriging方法比IDW方法更加精确。

Impacts of LUCC processes on potential land productivity in China in the 1990s

Using meteorological data and RS dynamic land-use observation data set, the potential land productivity that is limited by solar radiation and temperature is estimated and the impacts of recent LUCC processes on it are analyzed in this paper. The results show that the influence of LUCC processes on potential land productivity change has extensive and unbalanced characteristics. It generally reduces the productivity in South China and increases it in North China, and the overall effect is increasing the total productivity by 26.22 million tons. The farmland reclamation and original farmlands losses are the primary causes that led potential land productivity to change. The reclamation mostly distributed in arable-pasture and arable-forest transitional zones and oasises in northwestern China has made total productivity increase by 83.35 million tons, accounting for 3.50% of the overall output. The losses of original farmlands driven by built-up areas invading and occupying arable land are mostly distributed in the regions which have rapid economic development, e.g. Huang-Huai-Hai plain, Yangtze River delta, Zhujiang delta, central part of Gansu, southeast coastal region, southeast of Sichuan Basin and Urumqi-Shihezi. It has led the total productivity to decrease 57.13 million tons, which is 2.40% of the overall output.

Environmental Kuznets Curve for SO_2 Emissions in China

This study examines the link between the economic growth and the environmental quality. Based on a panel data set, a N-shaped Environmental Kuzents Curve has been found for the sample period: a cubic relationship between per capita GDP and emissions of sulphur dioxide (SO2). We also find that energy consumption is an important determinant of environmental degradation. The empirical results suggest that we should promote environmental protection as soon as possible.

Efficient HOG human detection

While Histograms of Oriented Gradients (HOG) plus Support Vector Machine (SVM) (HOG+SVM) is the most successful human detection algorithm, it is time-consuming. This paper proposes two ways to deal with this problem. One way is to reuse the features in blocks to construct the HOG features for intersecting detection windows. Another way is to utilize sub-cell based interpolation to efficiently compute the HOG features for each block. The combination of the two ways results in significant increase in detecting humans-more than five times better. To evaluate the proposed method, we have established a top-view human database. Experimental results on the top-view database and the well-known INRIA data set have demonstrated the effectiveness and efficiency of the proposed method. (C) 2010 Elsevier B.V. All rights reserved.

Supervised Gaussian Process Latent Variable Model for Dimensionality Reduction

The Gaussian process latent variable model (GP-LVM) has been identified to be an effective probabilistic approach for dimensionality reduction because it can obtain a low-dimensional manifold of a data set in an unsupervised fashion. Consequently, the GP-LVM is insufficient for supervised learning tasks (e. g., classification and regression) because it ignores the class label information for dimensionality reduction. In this paper, a supervised GP-LVM is developed for supervised learning tasks, and the maximum a posteriori algorithm is introduced to estimate positions of all samples in the latent variable space. We present experimental evidences suggesting that the supervised GP-LVM is able to use the class label information effectively, and thus, it outperforms the GP-LVM and the discriminative extension of the GP-LVM consistently. The comparison with some supervised classification methods, such as Gaussian process classification and support vector machines, is also given to illustrate the advantage of the proposed method.

Comparison of different methods for variable selection

In chemistry for chemical analysis of a multi-component sample or quantitative structure-activity/property relationship (QSAR/QSPR) studies, variable selection is a key step. In this study, comparisons between different methods were performed. These methods include three classical methods such as forward selection, backward elimination and stepwise regression; orthogonal descriptors; leaps-and-bounds regression and genetic algorithm. Thirty-five nitrobenzenes were taken as the data set. From these structures quantum chemical parameters, topological indices and indicator variable were extracted as the descriptors for the comparisons of variable selections. The interesting results have been obtained. (C) 2001 Elsevier Science B.V. All rights reserved.

Selection for variables in quantitative structure-activity/property relationship studies - Comparison between methods of orthogonal descriptors and leaps-and-bounds regression analysis

In this paper, the comparison of orthogonal descriptors and Leaps-and-Bounds regression analysis is performed. The results obtained by using orthogonal descriptors are better than that obtained by using Leaps-and-Bounds regression for the data set of nitrobenzenes used in this study. Leaps-and-Bounds regression can be used effectively for selection of variables in quantitative structure-activity/property relationship(QSAR/QSPR) studies. Consequently, orthogonalisation of descriptors is also a good method for variable selection for studies on QSAR/QSPR.

Molecular phylogeny and species identification of pufferfish of the genus Takifugu (Tetraodontiformes, Tetraodontidae)

The phylogenetic relationships and species identification of pufferfishes of the genus Takifugu were examined by use of randomly amplified polymorphic DNA (RAPD) and sequencing of the amplified partial mitochondrial 16S ribosomal RNA genes. Amplifications with 200 ten-base primers under predetermined optimal reaction conditions yielded 1962 reproducible amplified fragments ranging from 200 to 3000 bp. Genetic distances between 5 species of Takifugu and Lagocephalus spadiceus as the outgroup were calculated from the presence or absence of the amplified fragments. Approximately 572 bp of the 16S ribosonial RNA gene was amplified, using universal primers, and used to determine the genetic distance values. Topological phylogenic trees for the 5 species of Takifugu and outgroup were generated from neighbor-joining analysis based on the data set of RAPD analysis and sequences of mitochondrial 16S rDNA. The genetic distance between Takifugu rubripes and Takifugu pseudommus was almost the same as that between individuals within cacti species, but much smaller than that between T. rubripes, T. pseudommus, and the other species. The molecular data gathered from both analysis of mitochondria and nuclear DNA strongly indicated that T. rubripes and T. pseudommus should be regarded as the same species. A fragment of approximately 900 bp was amplified from the genome of all 26 T. pseudommus individuals examined and 4 individuals of intermediate varieties between T. rubripes and T. pseudommus. Of the 32 T. rubripes individuals, only 3 had the amplified fragment. These results suggest that this fragment may be useful in distinguishing between T. rubripes and T. pseudommus.

Circulation in the western tropical Pacific Ocean and its seasonal variation

An assimilation data set based on the GFDL MOM3 model and the NODC XBT data set is used to examine the circulation in the western tropical Pacific and its seasonal variations. The assimilated and observed velocities and transports of the mean circulation agree well. Transports of the North Equatorial Current (NEC), Mindanao Current (MC), North Equatorial Countercurrent (NECC) west of 140degreesE and Kuroshio origin estimated with the assimilation data display the seasonal cycles, roughly strong in boreal spring and weak in autumn, with a little phase difference. The NECC transport also has a semi-annual fluctuation resulting from the phase lag between seasonal cycles of two tropical gyres' recirculations. Strong in summer during the southeast monsoon period, the seasonal cycle of the Indonesian throughflow (ITF) is somewhat different from those of its upstreams, the MC and New Guinea Coastal Current (NGCC), implying the monsoon's impact on it.

On the western boundary currents in the Philippine Sea

In this study we describe the velocity structure and transport of the North Equatorial Current (NEC), the Kuroshio, and the Mindanao Current (MC) using repeated hydrographic sections near the Philippine coast. A most striking feature of the current system in the region is the undercurrent structure below the surface flow. Both the Luzon Undercurrent and the Mindanao Undercurrent appear to be permanent phenomena. The present data set also provides an estimate of the mean circulation diagram (relative to 1500 dbar) that involves a NEC transport of 41 Sverdrups (Sv), a Kuroshio transport of 14 Sv, and a MC transport of 27 Sv, inducing a mass balance better than 1 Sv within the region enclosed by stations. The circulation diagram is insensitive to vertical displacements of the reference level within the depth range between 1500 and 2500 dbar. Transport fluctuations are, in general, consistent with earlier observations; that is, the NEC and the Kuroshio vary in the same phase with a seasonal signal superimposed with interannual variations, and the transport of the MC is dominated by a quasi-biennial oscillation. Dynamic height distributions are also examined to explore the dynamics of the current system.

COASTAL SEA-LEVEL AND THE LARGE-SCALE CLIMATE STATE - A DOWNSCALING EXERCISE FOR THE JAPANESE ISLANDS

A major problem which is envisaged in the course of man-made climate change is sea-level rise. The global aspect of the thermal expansion of the sea water likely is reasonably well simulated by present day climate models; the variation of sea level, due to variations of the regional atmospheric forcing and of the large-scale oceanic circulation, is not adequately simulated by a global climate model because of insufficient spatial resolution. A method to infer the coastal aspects of sea level change is to use a statistical ''downscaling'' strategy: a linear statistical model is built upon a multi-year data set of local sea level data and of large-scale oceanic and/or atmospheric data such as sea-surface temperature or sea-level air-pressure. We apply this idea to sea level along the Japanese coast. The sea level is related to regional and North Pacific sea-surface temperature and sea-level air pressure. Two relevant processes are identified. One process is the local wind set-up of water due to regional low-frequency wind anomalies; the other is a planetary scale atmosphere-ocean interaction which takes place in the eastern North Pacific.

Evolution of nutrient structure and phytoplankton composition in the Jiaozhou Bay ecosystem

The inventories of nutrients in the surface water and large phytoplankton( > 69 pm) were analyzed from the data set of JERS ecological database about a typical coastal waters, the Jiaozhou Bay, China, from 1960s for N, P and from 1980s; for Si. By examining long-term changes of nutrient concentration, calculating stoichiometric balance, and comparing diatom composition, Si limitation of diatom production was found to be more possible. The possibility of Si limitation was from 37% in 1980s to 50% in 1990s. Jiaozhou Bay ecosystem is becoming serious eutrophication, with notable increase of NO2-N, NO3-N and NH4-N from 0.1417 mumol/L, 0.5414 mumol/L, 1.7222 mumol/L in 1960s to 0.9551 mumol/L, 3.001 mumol/L, 8.0359 mumol/L in late 1990s respectively and prominent decrease of Si from 4.2614 mumol/L in 1980s to 1.5861 mumol/L in late 1990s; the nutrient structure is controlled by nitrogen; the main limiting nutrient is probably silicon; because of the Si limitation the phytoplankton community structure has changed drastically.

Role of the Jiaozhon Bay as a source/sink of CO2 over a seasonal cycle

The seasonal evolution of dissolved inorganic carbon (DIC) and CO2 air-sea fluxes in the Jiaozhou Bay was investigated by means of a data set from four cruises covering a seasonal cycle during 2003 and 2004. The results revealed that DIC had no obvious seasonal variation, with an average concentration of 2035 mu mol kg(-1) C in surface water. However, the sea surface partial pressure of CO2 changed with the season. pCO(2) was 695 mu atm in July and 317 mu atm in February. Using the gas exchange coefficient calculated with Wanninkhof's model, it was concluded that the Jiaozhou Bay was a source of atmospheric CO, in spring, summer, and autumn, whereas it was a sink in winter. The Jiaozhou Bay released 2.60 x 10(11) mmol C to the atmosphere in spring, 6.18 x 10(11) mmol C in summer, and 3.01 x 10(11) mmol C in autumn, whereas it absorbed 5.32 x 10(10) mmol C from the atmosphere in winter. A total of 1.13 x 10(11) mmol C was released to the atmosphere over one year. The behaviour as a carbon source/sink obviously varied in the different regions of the Jiaozhou Bay. In February, the inner bay was a carbon sink, while the bay mouth and the Outer bay were carbon sources. In June and July, the inner and Outer bay were carbon sources, but the strength was different, increasing from the inner to the outer bay. In November, the inner bay was a carbon source, but the bay Mouth was a carbon sink. The outer bay was a weaker CO2 Source. These changes are controlled by many factors, the most important being temperature and phytoplankton. Water temperature in particular was the main factor controlling the carbon dioxide system and the behaviour of the Jiaozhou Bay as a carbon source/sink. The Jiaozhou Bay is a carbon dioxide source when the water temperature is higher than 6.6 degrees C. Otherwise, it is a carbon sink. Phytoplankton is another controlling factor that may play an important role in behaviour as a carbon source or sink in regions where the source or sink nature is weaker.