华东沿海高标准农田大气湍流特征的观测分析.pdf

6Advances in Meteorological Science and Technology 气象科技进展 13（3）-2023华东沿海高标准农田大气湍流特征的观测分析胡小旭1,2,5 王伟1,2,3 徐敏4 徐敬争5 陆琛6（1南京信息工程大学中国气象局生态系统碳源汇开放重点实验室，南京210044；2南京信息工程大学江苏省农业气象重点实验室，南京210044；3南京信息工程大学气象灾害预报预警与评估协同创新中心，南京210044；4江苏省气候中心，南京210019；5航天新气象科技有限公司，无锡214028；6江苏省海安市气象局，海安226600）摘要：为评估Monin-Obukhov相似理论和涡度相关系统在高标准农田的适用性，利用南通高标准农田水稻下垫面涡度相关观测数据，分析了大气稳定度、湍流方差相似性、湍流谱特征、湍流强度和湍流动能等。结果表明：该高标准稻田上方大气呈现白天不稳定、夜晚稳定的昼夜变化特征；Monin-Obukhov相似理论在该高标准农田适用，三维风速归一化标准差随大气稳定度的变化符合“1/3”次方规律，温度、比湿、CO2密度归一化标准差随大气稳定度的变化在大气不稳定时符合“1/3”次方规律；在惯性子区，三维风速的湍流能谱符合“2/3”次方斜率，垂直风速与温度、湿度和CO2密度的协谱符合“4/3”次方斜率，可见涡度相关系统能够观测该高标准农田的水热和CO2通量；三个方向上的湍流强度均随风速增大而减小，且水平方向湍流强度大于垂直方向；湍流动能随风速呈二次函数增大，呈现昼高夜低的日变化特征，大气中性时达到最大，且以动力湍流贡献为主。研究结果可为揭示华东地区高标准农田生态系统与大气之间的物质和能量交换特征以及机制提供参考。关键词：高标准农田，大气湍流特征，湍流强度，湍流谱，湍流动能DOI：10.3969/j.issn.2095-1973.2023.00.001Observation and Analysis on Atmospheric Turbulence Characteristics over a High Standard Farmland along the East China CoastHu Xiaoxu1,2,5,Wang Wei1,2,3,Xu Min4,Xu Jingzheng5,Lu Chen6(1 Key Laboratory of Ecosystem Carbon Source and Sink,China Meteorological Administration(ECSS-CMA),Nanjing University of Information Science&Technology,Nanjing 210044;2 Jiangsu Key Laboratory of Agricultural Meteorology,Nanjing University of Information Science&Technology,Nanjing 210044;3 Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters(CIC-FEMD),Nanjing University of Information Science&Technology,Nanjing 210044;4 Climate Center of Jiangsu Province,Nanjing 210019;5 Aerospace New Weather Technology Co.,Ltd,Wuxi 214028;6 Haian Meteorological Bureau of Jiangsu Province,Haian 226600)Abstract:In order to evaluate the applicability of Monin-Obukhov similarity theory and the eddy covariance system in high-standard farmland,this paper analyzed the atmospheric characteristics,including atmospheric stability,turbulent variance similarity,power spectra and cospectrum,turbulent intensity and turbulence kinetic energy,using the observation data of the eddy covariance system over the underlying surface of a high-standard farmland in Nantong,Jiangsu Province.The results show that the atmosphere over the farmland is unstable during the daytime and stable at night.The variation of normalized standard deviation of three-dimensional wind components with atmospheric stability obeys the“1/3”power law.The normalized standard deviation of temperature,specific humidity and CO2 density with atmospheric stability follows the“1/3”power law when atmosphere is unstable.Therefore,Monin-Obukhov similarity theory is applicable at such a high-standard farmland.In the inertia subregion,the power spectra of three-dimensional wind speeds follows the“2/3”slope relation,and the cospectrum of vertical wind speed and scalar quantities follows the“4/3”slope relation.So,the sensible heat,latent heat and CO2 fluxes over the high-standard farmland can be correctly observed by the eddy covariance system.The three-dimensional turbulence intensities decrease with the increase of wind speed,showing horizontal results stronger than those in vertical direction.The turbulent kinetic energy(TKE)increases with the increase of wind speed,following a quadratic function.TKE is high during the daytime and low at night,reaching its maximum under neutral atmosphere.Compared to thermal turbulence,mechanical turbulence contributes more to TKE.The research can advance our understanding of the energy and mass exchanges between high-standard farmland and atmosphere in East China.Keywords:high-standard farmland,atmospheric turbulence characteristics,turbulence intensity,turbulence energy spectra,turbulent kinetic energy收稿日期：2023年3月2日；修回日期：2023年4月12日第一作者：胡小旭（1987），Email：通信作者：王伟（1986），Email：资助信息：国家重点研发计划项目（2019YFA0607202）7Advances in Meteorological Science and Technology 气象科技进展 13（3）-2023Cover Story 封面报道0引言0引言农田占据了陆表约13%的面积1，它与大气之间的物质和能量交换影响着区域和全球气候，并在全球碳收支中占有重要地位。首先，卫星观测显示农田生态系统会加剧区域和全球变暖，在中国、印度和巴西等农业大国，其白天增暖效应超过5 oC2。其次，灌溉农田通过蒸散增湿和降低大气边界层高度显著增强湿热胁迫3，在中国华北、印度和巴基斯坦等地尤为明显。此外，农田生态系统是最活跃的碳库4，固碳量约占中国陆地生态系统的12%5。上述农田与大气之间的物质和能量交换均以湍流方式进行，观测分析大气湍流特征是准确量化农田物质和能量通量的基础。Monin-Obukhov相似理论是边界层湍流研究的理论基础，用于分析边界层中的外部参数对湍流扩散的影响。该理论大大推动了大气湍流观测试验的开展和观测仪器的研发4。测量湍流通量最直接的方式是涡度相关法6，该方法因其直接测定、理论假设少、观测信息全面等优点，被广泛用于观测农田上方的大气湍流特征。描述大气湍流特征的常见指标有湍流方差相似性、湍流强度、湍流能谱等4,7-8。余锦华等9 分析了常熟农业生态试验站的通量观测数据，发现在非中性条件下，水平风速标准差与大气稳定度的关系并不满足“1/3”次方定律，在平流输送明显的绿洲农田也观测到类似结果10，而“1/3”次方定律在云贵高原西部11和淮河流域12农田均成立。在湍流强度方面，杨智等11 观测发现云贵高原西部农田湍流强度比华北平原强，而李英等13 发现同为西部地区的成都平原的农田湍流强度却与华北平原观测结果接近，且农田湍流强度随风速的变化特征尚不明确14。郭建侠15 分析华北玉米下垫面的涡度相关资料发现，动量谱和感热谱在惯性子区均较为离散，不符合“4/3”次方规律。可见，由于农田下垫面类型多样且存在高度异质性，导致不同研究观测得到的农田大气湍流特征存在明显差异，使得大气边界层湍流研究基础理论之一Monin-Obukhov相似理论在农田的适用性存在挑战，亟须在开阔、平坦且均一的标准化农田开展大气湍流特征的观测分析。农业农村部经国务院批复同意于2021年9月6日印发的全国高标准农田建设规划（20212030年）指出，到2030年建成12亿亩（1亩666.7 m2）高标准农田，并改造提升现有高标准农田2.8亿亩，以此稳定保障1.2万亿斤（1斤0.5 kg）以上粮食产能。高标准农田是指田块平整、集中连片、设施完善、节水高效、农电配套、宜机作业、土壤肥沃、生态友好、抗灾能力强，与现代农业生产和经营方式相适应的旱涝保收、稳产高产的耕地16。高标准农田经过“田、水、路、林、村”综合整治，实现了“小田并大田”，下垫面变得均一、平坦，是通量观测的理想下垫面，理论上Monin-Obukhov相似理论更为适用，但这一推论尚未得到观测证实。南通市从2008年起建设高标准农田，目前拥有江苏省内单体最大的“万顷良田”建设示范区17。农田生态系统作为陆地生态系统的重要组成部分，明确高标准农田生态系统的碳收支特征对减少农田碳排放、增加土壤碳存储意义重大。本文基于南通高标准农田水稻下垫面2020年涡度相关观测数据，分析其大气稳定度分布、湍流方差相似性、湍流谱特征、湍流强度和湍流动能等，旨在论证Monin-Obukhov相似理论在高标准农田下垫面的适用性，评估基于湍流交换理论发展起来的涡度相关技术观测高标准农田碳通量的可靠性，为揭示华东地区典型农田生态系统与大气之间的物质和能量交换特征提供参考。1观测概况与研究方法1观测概况与研究方法1.1 观测概况1.1 观测概况观测站点位于江苏省南通市雅周镇“万亩良田”工程内，属北亚热带海洋性季风气候，四季分明。根据海安国家气象站的多年观测数据，19812010年，年平均气温为15.3，年平均日照时数为2000 h，年平均降水量为1015.1 mm（国家气象科学数据中心http:/ m半径范围内地势平坦，均为农田18。地下10100 cm的土壤类型为细沙土。种植方式为稻（610月）麦（11月翌年5月）轮作。整个观测系统由开路式涡度相关系统、小气候观测系统、数据采集器、通信设备和供电设备组成（图1）。涡度相关系统由开路式红外CO2/H2O分析仪（LI-7500DS，美国LI-COR公司）和三维超声风速仪（Windmaster Pro，英国Gill公司）组成，分别测量大气中的CO2/H2O浓度、三维风速和超声虚温，架设高度为2.5 m。小气候观测系统包括2路温湿度传感器（DHC2，航天新气象科技有限公司），用于测量大气温度和相对湿度，架设高度分别为2 m和3 m；翻斗式雨量传感器（SL3-1，上海气象仪器厂），架设高度为1.5 m；四分量净辐射传感器（FS-J1，航天新气象科技有限公司），用于测量向下短波、向上短波、向下长波和向上长波辐射通量密度，架设高度为2 m；5路光合有效辐射传感器（FS-PR，航天新气象科技8Advances in Meteorological Science and Technology 气象科技进展 13（3）-2023有限公司），1路架设高度为2 m，其他4路均架设在0.5 m高度；3路土壤热通量传感器（HFP01-05，荷兰Hukseflow公司），均埋设在土壤5 cm深度；5路土壤温湿度传感器，分别观测5、10、20、40和60 cm处的土壤温度和体积含水量。涡度相关系统的采样频率为10 Hz，小气候系统的采样频率为1 Hz。观测数据均由数据采集器（CR6000，美国Campbell公司）采集和存储。整个观测系统由4块100 W的太阳能板和4个120 Ah的蓄电池供电。1.2 研究方法1.2 研究方法1.2.1通量数据处理和质量控制1.2.1通量数据处理和质量控制选取2020年5月1日10月30日水稻生长期的观测数据，使用美国LI-COR公司的涡度相关数据处理软件EddyPro 6.1对10 Hz原始通量数据进行处理，得到30 min平均数据。数据后处理包括野点剔除、超声虚温订正、两次坐标旋转和密度效应订正等4。将30 min平均数据分为3个质量等级：0（最高）、1（中等）和2（最低）19，剔除1级和2级数据，仅保留0级数据分析大气湍流特征。通量贡献区描述的是来自通量传感器上风方向地表源对观测垂直通量的相对贡献6。为了明确通量观测信号的来源，使用美国LI-COR公司的涡度相关数据分析软件Tovi 2.8.1进行通量贡献源区分析20。经上述数据处理和质量控制后，共有4278条0级数据（占比53.5%）可用于湍流特征分析。1.2.2湍流特征量计算1.2.2湍流特征量计算（1）大气稳定度参数大气稳定度参数（）综合考虑了大气湍流的热力（浮力）和动力（风切变）生成机制，计算如下21：，（1）式中，z为观测高度；d为零平面位移；L为Obukhov长度。当时，大气为中性；当时，大气不稳定；当时，大气稳定22。，（2）式中，为垂直风速脉动（瞬时值与平均值之差）；为超声虚温脉动；为位温；k为von Karman常数（0.4）；g为重力加速度（9.8 m/s2）；和分别代表动力湍流和热力湍流的贡献。（2）湍流强度湍流强度（I）指风速标准差与平均风速的比值，用于描述湍流脉动量的相对波动程度。计算如下：，（3）式中，、分别为u（纬向风速）、v（经向风速）、w（垂直风速）方向上风速的标准差（m/s）；平均风速。（3）湍流动能湍流动能（）表示单位质量空气脉动速度所具有的动能，计算如下6：，（4）式中，、分别为、和方向上风速的方差。1.2.3湍流方差相似性1.2.3湍流方差相似性由Monin-Obukhov相似理论可知，在开阔、平坦且均一的下垫面上，当湍流充分发展时，近地层的三图1 研究区域位置（a，红色三角形为观测站点位置）和观测系统（b，农田生态系统小气候、通量观测系统的实景照片）Fig.1 Location of the research area(a,the red triangle representing the location of the observation station)and the observation system(b,the flux observation system for the microclimate of farmland ecosystem)(b)(a)117 E118 E119 E120 E121 E122 E31 N32 N33 N34 N35 N9Advances in Meteorological Science and Technology 气象科技进展 13（3）-2023Cover Story 封面报道维风速、温度、湿度和CO2密度等物理量的标准差经特征尺度参数无量纲化后，均可以表示为大气稳定度的函数23，即：，（5）式中，当x为三维风速u、v、w时，为摩擦风速；当x为温度（T）、湿度（q）和CO2密度（c）时，分别为、和。和分别为H2O和CO2的质量密度。对于风速，普适函数用下式拟合24：，（6）式中，和为常数；当时，前取负号，当时，前取正号。对于温度、湿度和CO2密度，普适函数用下式拟合25：，（7），（8）式中，、和为常数。1.2.4湍流能谱分析1.2.4湍流能谱分析湍流是由尺度大小不同的湍涡组成，按照频率（或波长）来研究各种尺度湍涡间的能量分布称为能谱分析。湍流能谱从低频到高频依次为含能涡区、惯性子区和耗散区。经Kolmogorov4证明，近地边界层内小尺度湍涡各向同性，在惯性子区内能量既不产生也不消耗，以“5/3”次方规律向更小尺度传递，可用于检验涡度相关系统的响应频率是否满足观测要求。，（9）式中，Fx为湍流能谱函数；x为三维风速u、v、w；为x方向上的无量纲Kolmogorov常数；E为湍流耗散率；波数，n为自然频率，为平均风速。当自然频率转换为归一化频率时，湍流能谱与归一化频率在惯性子区符合“2/3”次方规律，2个物理量之间的协谱与归一化频率在惯性子区符合“4/3”次方规律。2结果与分析2结果与分析2.1 通量贡献源区分析2.1 通量贡献源区分析图2是2020年水稻生长期的通量贡献源区图，等值线从内到外依次为10%80%的通量贡献源区范围，颜色表示空间上某一点对通量观测值的贡献率。由图2可知，水稻生长期内80%的通量贡献源区范围分布在距离涡度相关系统70 m以内的2块水稻田内，且主要通量贡献来源于涡度相关系统东侧的稻田，与研究时段内的主导风向一致。贡献率1.6 1.4 1.2 1.0 0.8 0.6 0.4 0.2 0103100 50 0 50 100y/mx/m50 m50050100图2 2020年水稻生长期涡度相关观测的通量贡献源区范围（底图来源于Google Earth，不同颜色表示每个点对涡度相关观测值的贡献率）Fig.2 Flux footprint area of eddy covariance system during rice growing period in 2020(The base image is from Google Earth.Different colors indicate the contribution of each point to eddy covariance observation.)2.2 风向、风速和大气稳定度特征2.2 风向、风速和大气稳定度特征如图3a所示，2020年水稻生长期内主导风向为东东南（ESE，13.1%），其次为东东北（ENE，12.3%）、东（E，12.1%）和东南（SE，10.1%）。观测高度2.5 m处平均风速在04 m/s，02 m/s和24 m/s的风速频率分别为53.4%和43.5%。由大气稳定度参数的概率密度分布图（图3b）可知，水稻生长期内，稻田上方大气以稳定状态为主，占研究时段的57.3%，大气不稳定频数占比为42.2%，大气中性频数占比为0.5%。稻田上方大气稳定度呈现显著的昼夜变化特征（图3c），白天（07：0016：00，北京时）不稳定，夜晚（17：00次日06：00，北京时，下同）稳定。2.3 湍流方差相似性2.3 湍流方差相似性2.3.1风速归一化标准差随大气稳定度参数的变化2.3.1风速归一化标准差随大气稳定度参数的变化由图4可知，三维风速归一化标准差随大气稳定度参数的变化均符合“1/3”次方规律，大气不稳定时的拟合效果更佳，尤其是在v方向上。大气不稳定时，v方向上的拟合值与观测值之间的相关系数R最大，为0.61；u方向上次之，R=0.60；w方向上最低，R=0.52。大气稳定时，三维风速的归一化标准差在15变化，变化幅度较大气不稳定时减小，“1/3”次方规律较大气不稳定时减弱。此时，拟合效果以w方向上最好，R=0.50；u和v方向上R分别为0.41和0.35。大气中性时，三维风速的归一化标准差都趋于常数，即，。10Advances in Meteorological Science and Technology 气象科技进展 13（3）-2023图3 2020年水稻农田风向、风速和大气稳定度（a）风向玫瑰图；（b）大气稳定度概率密度分布；（c）大气稳定度日变化（N、E、S、W分别表示北、东、南、西）Fig.3 Wind direction,wind speed and atmospheric stability over paddy field in 2020(a)Wind rose;(b)Probability density distribution of atmospheric stability parameter;(c)The diurnal variation of atmospheric stability parameter(N,E,S,W indicating wind direction of north,east,south,and west,respectively)图4 u、v、w风速分量归一化标准差（w、v、u）随大气稳定度参数（）的变化（ac）不稳定条件下；（df）稳定条件下（菱形为纵坐标在某个区间范围内的平均值，不稳定条件下区间为 102101100101102103104，稳定条件下区间为105104103102101100101）Fig.4 Normalized standard deviation of u,v,w wind components variation with atmospheric stability parameter (a-c)Unstable conditions;(d-f)Stable conditions(The diamond represents the bin average within a certain interval;interval is 102 to 101 to 100 to 101 to 102 to 103 to 104 under unstable conditions;interval is 105 to 104 to 103 to 102 to 101 to 100 to 101 under stable conditions)NNNENEENEESESESSESSWSWWSWWWNWNWNNWSE风向频率%8 68 46 24 02风速/(m/s)(a)14 12 10 8 6 4 2 0 2 4 6 8 10 12 14(b)概率分布/%大气稳定度参数0.80.60.40.200.20.40.60.81.01.21.450 45 40 35 30 25 20 15 10 5 0时间(c)000204060810121416182022240.25 0.20 0.15 0.10 0.05 0 0.05 0.10大气稳定度参数y2.55(14.01x)1/3,R0.60y2.623(11.048x)1/3,R0.41y2.115(11.779x)1/3,R0.35y1.163(11.318x)1/3,R0.50y2.40(18.91x)1/3,R0.61y1.1(12.55x)1/3,R0.52(a)(d)(b)(e)(c)(f)20 15 10 5 054321105032130 20 10 0864201021011001011021031041051041031021011001012.3.2标量归一化标准差随大气稳定度参数的变化2.3.2标量归一化标准差随大气稳定度参数的变化图5为温度、湿度和CO2密度归一化标准差随大气稳定度参数的变化。温度的归一化标准差随大气稳定度的变化在大气不稳定时符合“1/3”次方规律，但在大气稳定时符合“1”次方规律。比湿的归一化标准差随大气稳定度参数的变化在大气不稳定时符合“1/3”次方规律，大气稳定时趋于常数3.49。CO2密度归一化标准差与大气稳定度参数的关系在大气不稳定11Advances in Meteorological Science and Technology 气象科技进展 13（3）-2023Cover Story 封面报道图5 温度（T）、湿度（q）和CO2密度（c）的归一化标准差（t、h、c）随大气稳定度参数（）的变化（ac）不稳定条件下；（df）稳定条件下（菱形为纵坐标在某个区间范围内平均值，不稳定条件下区间为 101100101102103104，稳定条件下区间为104103102101100）Fig.5 Relationship between normalized standard deviation of T,q,c and atmospheric stability parameter (a-c)Unstable conditions;(d-f)Stable conditions(The diamond represents the ordinate average value within a certain interval;interval is 101 to 100 to 101 to 102 to 103 to 104 under unstable conditions;interval is 104 to 103 to 102 to 101 to 100 under stable conditions)和稳定条件虽然可以用“1/3”次方拟合，但结果较为离散，相关系数仅分别为0.13和0.05。2.4 湍流能谱分析2.4 湍流能谱分析图6为三维风速分量的归一化功率谱和垂直风速与温度、湿度、CO2的协谱。在惯性子区（010 Hz）内，u、v的功率谱均符合“2/3”次方规律，峰值分别出现在0.002 Hz和0.003 Hz附近；w的功率谱斜率略小于2/3，峰值出现在0.3 Hz附近。三维风速的功率谱在大于0.1 Hz的高频区间均上翘，说明观测系统受到高频噪音的影响。垂直风速与温度、湿度和CO2密度的协谱均与标准谱线26一致，在惯性子区内均符合“4/3”斜率特征，峰值都出现在0.1 Hz附近。可见，涡度相关系统能够有效地观测该高标准农田上方垂直风速与温度、湿度和CO2密度的协方差。2.5 湍流统计量分析2.5 湍流统计量分析2.5.1湍流强度2.5.1湍流强度由图7a可见，u和v方向上的湍流强度概率密度分布一致，分别有97.4%和96.1%的湍流强度集中分布在0.20.6，峰值均出现在0.45附近。与u、v方向相比，w方向上的湍流强度明显较低，概率密度峰值出现在0.2附近，0.050.4范围内的湍流强度占99.6%。三个方向上的湍流强度均随风速增大而减小（图7b）。风速小于2 m/s时，湍流强度随风速增加迅速降低；当风速超过3 m/s时，湍流强度趋于常数，即，。从区间均值来看，当风速超过1 m/s时，u方向上的湍流强度开始大于v方向的结果；当风速超过5 m/s时，。2.5.2湍流动能2.5.2湍流动能由图8a可知，大气中性时，湍流动能最大；当大气趋向稳定和不稳定时，湍流动能迅速下降。热力湍流对湍流动能的贡献（104）明显小于动力湍流。中性条件下，动力湍流对湍流动能的贡献最大。大气稳定时，逆温对应的负浮力抑制湍流发展，对湍流动能的贡献为负值。湍流动能随平均风速增大而增加（图8b），可用二次函数拟合两者之间的关系，相关系数R0.86。湍流动能呈现白天高、夜晚低的日变化特征（图8c），从06：00开始增大，13：00达到峰值0.96 m2/s2，随后下降，03：00降至最低值0.22 m2/s2。湍流动能日变化与风速一致，可见湍流动能主要受动力湍流控制。3讨论3讨论已有研究发现，长三角常熟地区水稻在非中性条件下水平风速标准差与大气稳定度之间的关系较分散，难以满足“1/3”次方定律9。华北玉米下垫面垂直风速与温度的协谱在惯性子区内不符合“4/3”次方斜率规律15。与上述研究相比，本研究选取的高标y3.416(x)1/3(x 0),R0.477y1.442(0.1121x)1/3,R0.23y0.4518x12.665,R0.75y3.492y2.913(0.002274x)1/3,R0.045y2.199(0.001866x)1/3,R0.128(a)(d)(b)(e)(c)(f)105100105104102100102100102102101100105100105104102100104103102101100101100101102103104c/|c*|h/|h*|t/|t*|12Advances in Meteorological Science and Technology 气象科技进展 13（3）-2023图6 三维风速分量的归一化功率谱（a，虚线为2/3次方斜率线）和垂直风速（w）与温度（T）、比湿（q）、CO2密度（c）的协谱（b，虚线为4/3次方斜率线）（黑色实线为Kaimal标准曲线）Fig.6 Normalized power spectra of three-dimensional wind components（a,dash line representing the 2/3 slope）and normalized cospectrum of vertical wind speed(w)with temperature(T),specific humidity(q)and CO2 density(c)（b,dash line representing the 4/3 slope)(solid line denoting the standard spectra of Kaimal)图7 三维风方向上湍流强度的概率密度分布（a）和湍流强度随风速的变化特征（b）（标记为0.2区间平均值）Fig.7 Probability density distribution of turbulence intensity of three-dimensional wind speeds(a)and variation of turbulence intensity with wind speed(b)(Marks are bin averages with width of 0.2)图8 湍流动能（）随大气稳定度（）、风速（）和时间的变化特征（a）湍流动能、动力湍流和热力湍流贡献随大气稳定度的变化；（b）湍流动能随风速的变化特征；（c）湍流动能和风速的日变化特征 Fig.8 Variations of turbulent kinetic energy()with atmospheric stability(),wind speed()and time(a)Turbulent kinetic energy,mechanical turbulence and thermal turbulence varying with atmospheric stability;(b)Turbulent kinetic energy varying with wind speed;(c)Diurnal variations of turbulent kinetic energy and wind speed(a)(b)2/34/3nSx/2xnCwx/cov(wx)fn(zd)/Ufn(zd)/UuvwwTwcwq1001011021001011021031041021001021041021001021009080706050403020100IuIvIw00.20.40.60.81.01.21.41.6(a)概率密度/%湍流强度IuIvIwmean(Iu)mean(Iv)mean(Iw)01234567(b)风速/(m/s)湍流强度101100101(b)3.0 2.5 2.0 1.5 1.0 0.5 0/(m/s)/(m2/s2)01234567y0.0817x20.0319x0.0536R0.86(c)2.5 2.4 2.3 2.2 2.1 2.0 1.9 1.8 1.7 1.60.9 0.8 0.7 0.6 0.5 0.4 0.3/(m/s)湍流动能 风速/(m2/s2)时间000408121620242.0 1.5 1.00.500.51.0(a)3.0 2.5 2.0 1.5 1.0 0.5 00.35 0.30 0.25 0.20 0.15 0.10 0.05/(m2/s2)u/(m/s)10414 12 10 8 6 4 2 0 213Advances in Meteorological Science and Technology 气象科技进展 13（3）-2023Cover Story 封面报道准农田中，三维风速归一化标准差随大气稳定度的变化均符合“1/3”次方规律，且三维风速的湍流谱在惯性子区中均符合“2/3”次方关系，垂直风速与标量的协谱在惯性子区中都符合“4/3”次方规律。可见，Monin-Obukhov相似理论在高标准农田更为适用，涡度相关系统能够准确地观测该高标准农田上的感热、水汽和CO2通量。大气近中性时，三维风速的归一化标准差都趋于常数。不同农田在中性条件下得到的 风 速 分 量 归 一 化 标 准 差 值 基 本 符 合的规律。本研究区域与常熟同属于长三角地区，本区域（、）与常熟稻田的结果（、）9接近，但 大 于 成 都 地 区 农 田 结 果（、）13，可能与周围地形和农作物差异有关。从湍流强度概率密度峰值来看，本研究中高标准农田的湍流强度峰值（、）约为华北平原观测结果（、）的2倍15，也明显大于成都平原的观测结果（、）13，与云贵高原大理地区结果（、）11相当。4结论4结论基于南通高标准农田水稻下垫面2020年510月涡度相关观测数据，分析了大气稳定度分布、湍流方差相似性、湍流谱特征、湍流强度和湍流动能等，得到以下结论。1）水稻生长期内，80%的通量贡献源区范围分布在距离涡度相关系统70 m