淮河流域农田近地层湍流通量特征研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
通过分析2007年9月~2008年8月年在淮河流农田下垫面采用涡动相关法(Eddy Covariance,简称EC)测得的地气间湍流通量以及辐射平衡分量和廓线资料,研究淮河流域农田地气相互作用日、季变化特征,分析了农田不同下垫面以及作物不同生长期地表反照率、水热通量、CO2通量、总体输送系数以及能量平衡比率等的变化规律,并且将涡度法与梯度法测得的感热、潜热和动量通量进行对比分析,以期能够全面深入了解农田生态系统的湍流输送特征。得到的主要结果有:
     (1)不同高度的风速、温度和湿度日、季变化均很明显。风速夜间小于白天,风速随高度的增加而增大;温度廓线夜间呈逆温状态,白天温度随高度的增加而减小;湿度基本随高度的增加呈递减趋势,但也出现了逆湿现象。10 cm以上土壤温度有明显的日、季变化特征,而10 cm以下土壤温度只存在明显月变化。农田土壤湿度与降水及下垫面状况有关。
     (2)农田下垫面近地层辐射平衡各分量和能量通量随着季节、天气和下垫面的变化,其变化的幅度、形态存在着明显的差异。淮河流域热量资源较丰富,太阳辐射年总量约为4878.8 MJ·m-2;净辐射在各能量间的分配和地表反照率主要受下垫面状况的影响,感热通量在下垫面为裸地时最大,占净辐射比例可达30%~50%;潜热一直是农田能量消耗的主要形式,占净辐射比例为50%~80%,土壤热通量占净辐射的比例一般均小于10%;反照率变化表现为:裸地<麦田<稻田<雪地。农田下垫面的总体输送过程中动力因子的作用大于热力作用;CO2通量受下垫面的影响较大,在不同下垫面以及作物不同的生长期,CO2通量均有明显的变化,CO2净吸收在作物抽穗期最大可达2.86 g/(m2·h),研究表明近中性条件是最有利于水热和CO2通量垂直传输的大气状态。
     (3)地表能量平衡分析表明:农田能量平衡处于不闭合状态,全天能量闭合率平均为0.89,各下垫面的能量平衡闭合率表现为:裸地>小麦>水稻>雪地,裸地闭合率最高可达0.91,而雪地只有0.46,白天能量平衡闭合程度明显大于夜间,白天能量闭合率平均为0.78,夜间平均仅有0.24;农田能量不闭合与夜间能量不平衡、下垫面部分能量吸收项的忽略以及湍流强度对湍流通量输送的影响有关,整个观测期涡度相关仪器直接观测的湍流能量(LE+H)始终小于有效能( Rn ? G),涡度相关仪器直接观测的湍流通量可能有被低估的趋势。
     (4)农田中梯度观测和涡度法观测的感热通量有较好的一致性,两方法所观测通量的复相关系数均超过了0.6,且通过信度检验,可认为梯度法与涡度法总体具有较好的相关性。
     (5)对采用MDV方法插补前后的涡度数据在显著性α=0.05水平上进行Z检验,除CO2通量正值区插补前后数据的分布差异显著外,其它数据插补前后差异均不显著。
The features of surface albedo, water and heat flux, CO2 flux, bulk transfer coefficients and energy balance in different underlying surface and growing period have been investigated based on the data of turbulent fluxes observed by Eddy Covariance method, radiation balance components and profile data over cropland in Huaihe River Basin from Sept.2007 to Aug.2008 in order to complete insight into the features of turbulent transfers over the cropland ecosystem. The seasonal and diurnal variation of the interaction between surface and atmosphere were also elucidated from these data. Furthermore, the sensible heat, latent heat and momentum fluxes are calculated by the eddy correlation and vertical gradient method are compared. The main results can be listed as follows:
     (1)It is found that diurnal and seasonal variation features of wind speed, temperature and humidity in different height is obvious, Wind velocity at daytime is larger than night which increase with elevate observation height. The air temperature decrease as elevate observation height at daytime and temperature inversion state appear at night, Humidity decrease frequently as elevate observation height, but humidity inversion state is also observed. The diurnal and seasonal variation features of soil temperature from 0cm layer to 10cm layer is obvious, however, only monthly variation regulation exist surpass that depth. The soil humidity of cropland is dependent on precipitate and underlying surface state.
     (2)The analysis of radiation balance components and energy flux in Huaihe River basin shows that: the range and pattern variation of radiation balance components and energy flux are influenced by the change of season, weather and underlying surface state. Annual total solar radiation amount in this area reaches around 4878.8 MJ·m-2, indicated that the Huaihe River basin has an abundance heat resource. We also found that allocation of net radiation among various energy and surface albedo were strongly influenced by underlying surface state. The heat fluxes of sensible and latent are mainly budget items of net radiation. The largest ratio of the sensible heat flux to the net radiation is about 30% to 50% on bare soil. The Latent heat flux is a prevailing item of the net radiation over the cropland, which is 50% to 80% of net radiation, whereas soil heat flux account for less than 10%. The results from the surface albedo investigation suggest that: bare field< wheat field< paddy field < snow field. The bulk transfer coefficients of momentum is larger than the bulk transfer coefficients of sensible heat at the same observation height. To some extent, CO2 flux is determined by underlying surface state, and the variation of CO2 flux with different underlying surface and crop growing stage is obvious. The net CO2 assimilation reached maximum value of 2.86 g/(m2·h) in jointing stage during the crop growing period, and study further shows that water, heat and CO2 flux vertical transportation was in favor of near neutral conditions.
     (3)The analysis of the surface energy balance shows that the phenomenon of energy imbalance exists over cropland. The average energy closure rate was 0.91 in bald field, whereas only 0.46 in snow field and the study shows that: bare field> wheat field> paddy field> snow field during the observed period. The average energy closure rate was 0.89 during the observed period, which 0.78 at daytime and 0.24 at night. The energy in-closure of cropland was influenced by energy imbalance at night, the ignorance of some energy item and turbulent intensity impact on turbulent transfers. During the observation period, the turbulent fluxes(LE+H)that observed by eddy covariance method always less then available energy( Rn ? G),which indicated that turbulent fluxes might be underestimated.
     (4)Sensible heat, latent heat, momentum flux observed by vertical gradient method was consistent with eddy covariance method. The multiple correlation coefficients between data acquired by these two methods were greater than 0.6 and reliability test indicated that there was a fine correlation between two methods.
     (5)Using Z-test at 0.05 level to analysis the difference of EC data before and after interpolate by MDV method, it shows that the difference of CO2 positive value region was obvious.
引文
[1]于贵瑞,李正泉,孙晓敏等《.陆地生态系统通量观测的原理与方法》[M].高等教育出版社, 2006, 13:391-392
    [2]杨晓光,沈彦俊等.麦田水热传输的控制效应及非线性特征分析[J].中国农业大学学报, 2000, 01: 018
    [3] Harazono, Y.J. Kim, A. Miyata, T. Choi, J.I. Yun and J.W. Kim, Measurement of energybudget components during the International Rice Experiment (IREX) in Japan[J].HydrologicalProcesses, 12: 2018-2092
    [4]王文杰,石福臣等.陆地生态系统二氧化碳通量网的建设和发展[J].东北林业大学学报, 2002, 30: 57-61
    [5]周秀骥.长江三角洲低层大气与生态系统相互作用研究.北京气象出版社[J].2004: 157-172
    [6] Winbank, baldoc-chi, Swinbank W C. Measurement of vertical transfer of heat and water vapor by eddies in the lower atmosphere .Journal of Meteorology, 1951, 8: 135-145
    [7] Onue, Ionue I.An aerodynamic measurement of photosynthesis over a paddy field[J].In Proceeding of the 7th Japan National Congress of applied echaics, 1958, 211-214
    [8] Kaimal J C, Wyngaard J C. The Kansas and Minnesota experiments Boundary-Layer[J]. Meteorology, 1990, 50: 31-47
    [9] Dennison M S, Berry J F. Wetlands. Guide to science, law, and technology [M]. Park Ridge, NJ, USA. Noyes publications, 1989, 439
    [10]张桂华,王艳秋等.气候变暖对黑龙江作物生长影响及对策[J].自然灾害学报, 2004, 13 (3): 95-100
    [11] Wood S, Sebastian K. Pilot analysis of Global Ecosystem Agroecosystems Technical Report. Washington D C Word Resources Institute and International Food Policy Research Institute. 2000
    [12] Lal R, Kimbel J M, Follet R F, etal. The potential of U S copland to sequestration carbon and mitigation to the greenhouse effect chelsea M I[J].Sleeping Bear press, 1998, 55-87
    [13] Cole C V, Cerri C, Minami K, etal. Agricultural options for mitigation of greenhouse gas emission, Climate Change 1995: Impactism Adaptations and Mitigation.of Climate Change. 1996, 23: 745-771
    [14]涂钢.半干旱区不同下垫面地气相互作用研究.中国科学院研究生院博士学位论文, 2007年1月: 9
    [15]牛国跃,洪钟祥,孙菽芬.陆面过程研究的现状与发展趋势[J].地球科学进展, 1997, 12(1): 20 -25
    [16] Joon Kim, Jinil Yun.Measuring fluxes energy, water and C02 over forests, paddies and coastal areasin Korea[J].International Workshop for Advanced Flux Network and Flux Evaluation Proceedings. Hokkaido University, Sapporo, Japan, 2001, 73-78
    [17]胡雯,王东勇,张苏等.淮河流域旱涝易发区农田的陆气相互作用―CO2通量、能量交换和水汽输送的季节变化特征[J].自然灾害学报, 2006, 15(6): 92-100
    [18] Schmugge T J, Andre J C. Land surface evaporation measurement and parameterization [M].Springer Verlag, New York, 1991: 424
    [19] Sellers P J, Hall F G, Asrar C, etal. The first ISLSCP field experiment (FIFE) [J]. Bull.Amer. Meteor. Soc, 1988, 69 (1): 22-27
    [20] Goutorbe J P, Lebel T, Tinga A, etal. A large scale study of land-atmosphere interaction in the semi-arid tropics (HAPEX-Sahel) [J]. Ann. Geophys, 1994, 12: 53-64
    [21]陈海山,孙照渤.陆气相互作用及陆面模式的研究进展[J].南京气象学院学报, 2002, 25(2): 277-288
    [22] Bolle H J, Andre J C, Arrue J Letal. The European field experiment in a desertification threatened area [J]. Ann. Geophys, 1993, (2): 173-189
    [23] Lundin L C, Halldin S. NOPEX Experiment Plan: Part I and 11. NOPEX cenUal Office Uppsala, Sweden, 1994
    [24] Hall F GS Sellers P J, Apps M. BOREAS: Boreal Ecosystem Atmosphere Study [J]. IEEE.Geosc. Remote. Sens. Soc. Newsletter, Mar, 1993, 9-17
    [25] Leuning R, Raupach M R, Coppin P A, et al. Spatial and temporal variations in fluxes of energy, water vapour and carbon dioxide during OASIS 1994 and 1995[J].Boundary-Layer Meteorology, 2004, 110: 3-38
    [26] Lundin L C, Halldin S. NOPEX Experiment Plan: Part I and 11. NOPEX cenUal Office Uppsala, Sweden, 1994
    [27] Law ford R G. Amidterm report on the GEWEX Continental-scale International Project (GCIP) [J]. J. Geophys. Res., 1999, 104(D16): 19279-19292
    [28] JSC. Report of the 15th Session of the Joint Scientific Committee of WCRP [C]. WMO/I'DNo. 632. WMO, Geneva, 1994: 97
    [29] Sellers P J, C A Nobre, D J Fitjarrald, etal.LAMBADA-BATERISTA: A preliminary science plan for a large scale biosphere-atmosphere field experiment in the Amazon Basin[J].ISLSCP IGPO Washington DC, 1993: 49
    [30] Charney J G, W J Quirk, S M Chew, etal.A comparative study of the effects of albedo change on drought in semi-arid region[J].J AtmosSci, 1977, 34(9): 1366-1385
    [31] Walker J, P R Rowntree.The effect of soil moisture on circulation and rainfall in a tropical model[J].Q.J R Meteorol Soc, 1977(103): 29-46
    [32] Shukal J, Y Mintz.The influence of land-surface evapotranspiration on Earth’s Climate [J]. Science, 1982 (215): 1498-1501.
    [33] Huang J, H M Van den Dool,K P Georgakakos.Analysis of model-calculated soil moisture over the United States(1931–1993) and applications to long-range temperature fore-casts[J].J Climate, 1996(9): 1350-1362.
    [34] Durre I, J M Wallace, D P Lettenmaier. Dependence of extreme daily maximum temperatures on antecedent soil moisture in the contiguous United States during summer [J].Climate, 2000(13): 2641-2651
    [35] Randal D, Koster.Regions of Strong Coupling Between Soil Moisture and Precipitation [J].Science, 2004, 305: 1138-1140
    [36] Kim J, Verma S B. Comoponents of surface energy balance in a temperate grassland ecosystem [J].Boundary-Layer Meteorology, 1990, 51(4): 401-417
    [37] George L Vourlitis, Nicolau Priante Filho. Seasonal variations in the evapotranspiration of a transitional tropical forest of Mato Grosso [J]. Brazil. Water Resources Research, 2002, 38(10): 10-29
    [38] Sellers, P J, Hall F, Margolis H, etal. The boreal ecosystem-atmosphere study (BORERS): an overview and early results from the 1994 field year [J]. Bull. Amer. Met. Soc .1995, 76: 1549-1577
    [39] Shannon Spencer, George Huczek, Bradley Muir. Developing a student-scientist partnership: Boreal Forest Watch [J]. J. Sci. Educ. Technol., 1998, 7 (1): 31-43
    [40] Halldin S, Gryning S-E, Gottschallc L Energy water and carbon exchange in a boreal forest landscape-NOPEX experience [J]. Agr. Forest. Meteorol., 1999, 98-99: 5-29
    [41] IPCC Climate Change 1994.In: Houghton J T, L G Meira Filho, J.Bruce etal.Radiative forcing of climate change and an evolution of the IPCC 1922 Emission Scenarios[M].U K: Cambridge University Press, 1995
    [42]王介民.陆面过程实验和地气相互作用研究―从HEIFE到IMGRASS和GAME-Tibet/ TIPEX[J].高原气象, 1999, 18(3): 280-294
    [43]胡隐樵,高由禧.黑河试验(HEIFE)—对干旱地区陆面过程的一些新认识[J].气象学报, 1994, 52 (3) : 285- 296
    [44]胡隐樵,高由禧,王介民等.黑河实验((HEIFE)的一些研究成果[J].高原气象, 1994, 13 (3): 225-236.
    [45]王介民,邱华盛.中日合作亚洲季风实验一青藏高原实验(GAME-Tibet) [J].中国科学院院刊, 2000: 386-388
    [46]张雁.淮河流域能量与水分循环试验和研究(HUBEX)项目进展[J].气象科技, 2000: 11-15
    [47] Lu Daren. Inner Mongolia semi-arid grassland soil-vegitation-atmosphere interaction (IMGRASS)[J].Global Change News Letter, 1997, 31: 4-5
    [48]张强,黄荣辉,王胜等.西北干旱区陆一气相互作用试验(NWC-ALIEX)及其研究进展[J].地球科学进展, 2005, 20 (4 ): 427-441
    [49]胡国权,丁一汇.HUBEX试验期间不同地面的能量收支研究[J].气候与环境研究, 2001, 6(2): 228-233
    [50] Halldin S, Gryning S-E, Gottschallc L Energy water and carbon exchange in a boreal forest landscape-NOPEX experience [J]. Agr. Forest. Meteorol., 1999, 98-99: 5-29
    [51]吕达仁,陈佐忠,王庚辰等.内蒙古半干旱草原气候一生态相互作用问题—IMGRASS计划初步结果[J].地学前缘, 2002, 9 (2): 307-320
    [52]马耀明,戴有学,马伟强等.干旱半干旱区非均匀地表区域能量通量的卫星遥感参数化[J].高原气象, 2004. 23 (2 ): 139-146
    [53]胡隐樵.黑河实验(HEIFE)能量平衡和水汽输送研究进展[ J ].地球科学进展, 1994, 9 (4): 30-34
    [54]黄宝霞.半干旱地区地表能量平衡特征研究.研究生学位论文, 2003: 12
    [55]李家伦,洪钟祥,罗卫东.青藏高原改则地区近地层通量观测研究[ J ].大气科学, 1999, 23 (2): 142-151
    [56]卞林根,陆龙骅,程彦杰,等.青藏高原东南部昌都地区近地层湍流输送地观测研究[J].应用气象学报, 2001, 12 (1): 1-13
    [57] Gao Z, Wang J, Ma Y, Kim J, Choi T, Lee H, Assanuma J, Su Z. Calculation ofNear2surface Layer Turbulent Transport and Analysis of Surface Thermal Equilibrium Features in Nagqu of Tibet[ J ]. Phys. Chem. Earth (B), 2000, 25 (2) : 135-139
    [58] Gao Z, Bian L, Zhou X. Measurements of Turbulent Transfer in the Near surface Layer over a Rice Paddy in China [ J ]. J. Geophys. Res. 2003, 108 (D13): 4387
    [59]王澄海,黄宝霞,杨兴国.陇中黄土高原植被覆盖和裸露下垫面地表通量和总体输送系数研究[J].高原气象, 2007年2月Vol. 26: 30-38
    [60]李国平,段廷扬,巩远发.青藏高原近地层通量特征的合成分析[J ] .气象学报, 2002, 60 (4): 453-460
    [61]李国平,段廷扬,吴贵芬.青藏高原西部的地面热源强度及地面热量平衡[J ].地理科学, 2003, 23 (1) : 13-18
    [62]余锦华,刘晶淼,丁裕国.青藏高原西部地表通量的年、日变化特征[J ].高原气象, 2004, 23 (3): 353-359
    [63]梁玲,吕世华,柳媛普.黄土高原植被变化对环境影响的数值模拟[J].高原气象, 2006, 25 (4): 575-582
    [64] Randal D K, Suarez M J. Soil moisture memory in climate models [J].Journal of Hydrometeorology, 2001, 2: 558 - 570
    [65] Li Jikang, Islam S. Estimation of root zone soil moisture and surface fluxes partitioning using near surface soil moisture measurements [J]. Journal of Hydrometeorology, 2002, 225: 1-14
    [66] Budyko M I. The Heat Balance of the Earth’s Surface [M].Leningrad, Russia: Gidrometeoizdat, 1956
    [67] Manabe S. Climate and the ocean circulation.1.The atmospheric circulation and the hydrology of the earth’s surface [J]. Monthly Weather Review, 1969, 97: 739-805
    [68] Deardorff J W. Efficient prediction of ground surface temperature and moisture, with inclusion of a layer of vegetation [J].Journal of Geophysical Research, 1978, 83 (C4): 1889-1903
    [69] Dickinson R E, Henderson S A, Kennedy P J, et al. Biosphere Atmosphere Transfer Scheme (BATS) for NCAR Community Climate Model [R]. NCAR/TN-275+STR, 1986
    [70] Dickinson R E, Henderson S A, Kennedy P J. Biosphere Atmosphere Transfer Scheme (BATS) version E as coupled to the NCAR Community Climate Model. NCAR, Tech Note[R]. NCAR/TN-387+STR, 1993
    [71] Sellers P J, Mintz Y, Sud Y C. A simple biosphere model (SiB) for use within general circulation models [J]. Journal of Atmospheric Science, 1986, 43(6): 505-531
    [72] Chen T H, Henderson S A, Milly P, et al. Cabauw experimental results from the project for intercomparison of landsurface parameterization scheme [ J]. Journal of Climate, 1997, 10: 1194-1215
    [73] Shao Y, Henderson S A. Validation of soil moisture simulation in land surface parameterization schemes with HAPEX data [J]. Global and Planetary Change, 1996, 13: 11-46
    [74] Lohmann D, Lettenmaier D P, Liang Xu, etal. The project for intercomparison of land-surface parameterization schemes (PILPS) phase(c) Red-Arkansas River basin experiment: Spatial and temporal analysis of water fluxes [J]. Global andPlanetary Change, 1998, 19: 161-179
    [75]雍斌,张万昌等.水文模型与陆面模式耦合研究进展[J].冰川冻土. 2006, 28(6): 961-970
    [76]牛国跃,洪钟祥,孙淑芬.陆面过程研究现状与发展趋势[J].地理科学进展, 1997, 12(1): 20-25
    [77]薛根元,周锁铨,孙照渤等.陆面过程研究的新进展[J].科技通报, 2005, 21(4): 378-385
    [78]朱治林,孙晓敏,张仁华.用微气象方法估算淮河流域能量平衡(HUBEX/IOP1998/99)的统计分析和比较研究[J].大气科学, 2003, 20(2): 285-291
    [79] Swinbank W C. The measurement of vertical transfers of heat and water vapor by eddies in the lower atmosphere [J]. J. Meteor, 1951, 8: 135-145
    [80]陈发祖.微气象学研究的若干进展一兼评涡度相关技术的作用[M].中国农业小气候研究进展,北京:气象出版社, 1993: 18-25
    [81] Dyer A 1, Hicks B B. Flux gradient transport of heat and water in an unstable atmosphere [J].Quart.J, Roy. Met. Soc., 1970, 9: 501-508
    [82] Hicks B B. Wind profile relationships from the Wangara experiment [J]. Q J R Meteor Soc, 1976, 102: 535-551
    [83]刘树华,黄子深,刘立超等.空气动力学法在湍流通量计算中的误差分析[J].气象, 21(3): 3-7
    [84] Eva Falge, Dennis Baldocchi, Richard Olson, etal. Gap filling strategies for long term energy flux date sets [J]. Agricultural and Forest Meteorology, 2001(107): 71-77
    [85] Webb E.K., Pearman G.I.and Leuning R.Correction of flux measurements for density effects due to heat and water vapor transfer [J].Quart J. Roy. Meteorology, 1980, 106: 85-100
    [86]高志球,卞林根,陆龙骅等.水稻不同生长期稻田能量收支、CO2通量模拟研究[J].应用气象学报, 2004, 15(2): 129-141
    [87]赵逸舟,马耀明,马伟强等.藏北高原土壤温湿变化特征分析[J] .冰川冻土, 2007, 29(4): 578-583
    [88]李国平,赵邦杰,卢敬华.青藏高原总体输送系数的特征[J].气象学报, 2002, 60(1): 60-67.
    [89]张述文,邱崇践,张卫东.用变分方法估算淮河流域总体输送系数[J].高原气象, 2004, 23(4): 506-511
    [90] Monteith J.I., Unsworth M. H. Principles of environmental physics. 2nded [J].London: Arnold 1990
    [91]卞林根,高志球,陆龙骅等.长江下游农业生态区CO2通量的观测研究[J].应用气象学报, 2005, 16(6): 828-834
    [92]郭建侠,卞林根,戴永久等.在华北玉米生育期观测的16 m高度CO2浓度及通量特征[J].大气科学, 2007, 31(4) 695-707
    [93]郭晓峰,康凌,蔡旭晖等.华南农田下垫面地气交换和能量收支的观测研究[J].大气科学, 2006, 30(3)453-463
    [94]刘树华,麻益民.农田近地面层CO2和湍流通量特征研究[J].1997, 55(2): 187-199
    [95]李正泉,于贵瑞,温学发等.中国通量观测网络(ChinaFLUX)能量平衡闭合状况的评价[J].地球科学2004, 34: 46-56
    [96] Schmid H P. Experimental design for flux measurements: matching scales of observations and fluxes [J]. Agric For Meteorol, 1997, 87: 179-200
    [97] Mayocchi C L, Bristow K L. Soil surface heat flux: some general questions and comments on measurements[J].Agric For Meteorol, 1995, 75: 43-50
    [98] Verhoef A, Hurk B J M, Adrie F G, et al. Thermal properties forvineyard (EFEDA-I) and savanna (HAPEX-Sahel) sites[J].Agric For Meteorol, 1996, 78: 1
    [99] Stannard D I, Blanford J H, Kustas W P. Interpretation of surface flux measurements in heterogeneous terrain during the Monsoon experiment [J].Water Resour Res, 1994, 30(5):1227-1239.
    [100] Mc Caughey J H. Energy balance storage terms in a mature mixed forest at Petawawa Ontario- a case study [J]. Bound-Lay Meteorol, 1985, 31: 89-101
    [101] Moore C J. Frequency response corrections for eddy correlation systems [J].Bound-Lay Meteorol, 1986, 37: 17-35
    [102]陆龙弊,卞林根,程彦杰等.冬季北京城市近地层微气象特征[J].应用气象学报, 2002, 13: 34-41
    [103]朱治林,孙晓敏,温学发等.中国通量网(ChinaFlux)夜间CO2涡度相关通量数据处理方法[J].中国科学D辑地球科学, 2006, 36: 34-44
NGLC 2004-2010.National Geological Library of China All Rights Reserved.
Add:29 Xueyuan Rd,Haidian District,Beijing,PRC. Mail Add: 8324 mailbox 100083
For exchange or info please contact us via email.