黄土高原半干旱区农田生态系统蒸散与作物系数特征

蒸散是地表能量平衡和水分平衡的重要组分,与水循环过程密切相关.采用涡度相关技术,对黄土高原半干旱区农田生态系统2010年生长季(4-9月)蒸散特征进行观测,分析了农田作物系数与环境因子的关系.结果表明： 在观测期间,研究区各月潜热通量(LE)日变化呈近似“单峰型”曲线特征,最大峰值出现在8月(151.4 W·m^-2);日间能量分配方式存在明显季节差异,4-6月的日间能量分配表现为LE/R_n＜H/R_n(R_n为净辐射,H为感热通量),7-9月的日间能量分配方式(LE/R_n＞H/R_n)与4-6月相反.黄土高原半干旱区农田日蒸散率存在显著季节变化特征,最大日蒸散率为4.69 mm·d^-1.风速(W_s)、空气相对湿度(RH)、土壤含水量(θ)和饱和水汽压差(D)是作物系数(K_c)的主要影响因子.K_c随W_s增加呈指数降低趋势,随RH、θ增加呈指数增长趋势,随D增加呈线性降低趋势.     

黄土高原春小麦农田蒸散及其影响因素

蒸散与水循环、能量平衡密切相关,是黄土高原雨养农田生态系统最重要的水通量之一。准确测定半干旱区农田生态系统蒸散,对增强陆气相互作用的理解以及科学应对气候变化有重要意义。采用涡度相关技术对黄土高原春小麦农田生态系统蒸散进行了观测,利用气象梯度系统进行环境因子观测;分析了春小麦农田生态系统蒸散日、季动态及其环境影响因子。结果表明,黄土高原半干旱区春小麦农田生态系统蒸散呈早晚低、中午高的"单峰型"日变化特征;最大日峰值出现在8月(0.22mm/h)。生长季蒸散日峰值高于非生长季。春小麦农田最大日蒸散率值相对较低,这可能与该地区干旱少雨的气候特征有关。农田蒸散且具有明显的季节动态,与降水季节分布密切相关。7、8月份降水较多,月蒸散量较高。全年蒸散量(318.0 mm)略低于年降水量(332.3 mm);蒸散量与降水量比值为95.7%。非生长季蒸散量显著低于生长季(4—9月);二者之比为0.26。农田蒸散随土壤含水量和空气温度(低于26℃)增大呈指数增长趋势;随空气相对湿度、太阳辐射、风速增大呈先增大后降低的二次曲线变化趋势。净辐射是黄土高原半干旱区农田生态系统蒸散主要环境控制因子,土壤含水量次之。     

玉米农田蒸散过程及其对气候变化的响应模拟

应用基于生理生态学过程的EALCO模型,对玉米农田生态系统的蒸散(ET)过程进行了模拟,在模型检验基础上,使用该模型模拟了玉米农田生态系统ET过程对未来气候变化的响应。结果表明,EALCO模型中能量与水过程的动态耦合机制使模型能够较好地模拟农田蒸散过程,基于涡度相关法的观测值与模型模拟值在小时、日尺度上均吻合较好,模型可以解释67%的日蒸散的变化特征。对土壤蒸发与冠层蒸腾的分别模拟显示,生长季土壤蒸发约占ET的36%。温度的升高会引起ET与冠层蒸腾的增加,同时土壤蒸发减少;ET对降水减少的响应较为敏感,主要表现在土壤蒸发的下降。大气CO2浓度升高对冠层蒸腾影响显著,该情景下冠层蒸腾下降幅度最大。研究所假设的2100年气候情景下,该农田生态系统生长季蒸散将减少,然而相对于降水的减少而言,蒸散的减少量较小,即水分支出项相对增加,因此,发生土壤水分匮乏的可能性加大,这可能会加剧该地区的暖干化趋势,给作物产量及生态环境带来威胁。     

开垦对黄河三角洲湿地净生态系统CO2交换的影响

近年来, 由于对湿地的不合理利用, 自然湿地被大面积地垦殖为农田, 导致湿地生态系统碳循环的模式发生改变, 从而影响了湿地生态系统碳汇功能。该研究通过涡度相关法, 对山东省东营市黄河三角洲芦苇(Phragmites australis)湿地和开垦多年的棉花(Gossypium spp.)农田的净生态系统CO₂交换(NEE)进行了对比观测, 以探讨该地区典型生态系统NEE的变化规律及其影响因子, 揭示开垦对芦苇湿地NEE和碳汇功能的影响。结果表明: 在生长季, 湿地和农田生态系统NEE的日平均值各月均呈明显的“U”型变化曲线, 非生长季NEE的变幅很小。生长季湿地生态系统日最大净吸收值和释放值分别为16.04 g CO₂·m^-2·d^-1(8月17日)和14.95 g CO₂·m^-2·d^-1(8月9日); 农田生态系统日最大净吸收值和释放值分别为18.99 g CO₂·m^-2·d^-1 (8月22日)和12.23 g CO₂·m^-2·d^-1 (7月29日)。生长季白天两个生态系统NEE与光合有效辐射(PAR)之间呈直角双曲线关系; 非生长季NEE主要受土壤温度(T_s)的影响; 生态系统生长季夜间NEE受T_s和土壤含水量(SWC)的共同影响; 湿地和农田的生态系统呼吸熵(Q₁₀)分别为2.30和3.78。2011年生长季, 黄河三角洲湿地和农田生态系统均表现为CO₂的汇, 总净固碳量分别为780.95和647.35 g CO₂·m^-2, 开垦降低了湿地的碳吸收能力; 而在2011年非生长季, 黄河三角洲湿地和农田生态系统均表现为CO₂的源, CO₂总释放量分别为181.90和111.55 g CO₂·m^-2。全年湿地和农田生态系统总净固碳量分别为599.05和535.80 g CO₂·m^-2。     

西北黑河中游荒漠绿洲农田作物蒸腾与土壤蒸发区分及作物耗水规律

利用中国生态系统研究网络临泽内陆河流域研究站绿洲农田2009年小气候、湍流交换、土壤蒸发和叶片气孔导度等综合观测试验数据,应用Shuttleworth-Wallace(S-W)双源模型以半小时为步长估算了绿洲农田玉米生长季实际蒸散量,并利用涡动相关与微型蒸渗仪实测数据对田间蒸散发量和棵间土壤蒸发量计算结果进行了检验。结果表明:S-W模型较好地估算研究区的蒸散量,并能有效区分农田作物蒸腾和土壤蒸发;全生育期玉米共耗水640 mm,其中作物蒸腾累积量为467 mm,土壤蒸发累积量为173 mm,分别占总量的72.9%和27.1%;日时间尺度上,作物蒸腾和土壤蒸发分别在0—6.3 mm/d和0—4.3 mm/d之间变化,其日平均分别为2.9和1.0 mm/d;田间供水充足,作物蒸腾与土壤蒸发比值明显受作物生长过程影响,播种—出苗期、出苗—拔节期、拔节—抽雄期、抽雄—灌浆期、灌浆—成熟期,其比值分别为0.04、0.8、7.0、5.2和1.4,不同阶段的比值差异主要受叶面积指数影响。     

环境和生物因子对黄河三角洲滨海湿地净生态系统CO2交换的影响

通过涡度相关和微气象观测技术,对黄河三角洲滨海湿地净生态系统CO₂交换(NEE)以及环境、生物因子进行了观测,探究湿地NEE变化规律及环境和生物因子对NEE的影响. 结果表明: 在日尺度上,生长季NEE呈明显“U”型曲线,非生长季变幅较小;在季节尺度上,NEE生长季波动较大,表现为碳汇,非生长季波动较小,表现为碳源;在年尺度上,滨海湿地生态系统表现为碳汇,总净固碳量为-247 g C·m^-2. 白天NEE主要受控于光合有效辐射(PAR),且生态系统表观量子产量(α)与白天生态系统呼吸(R_eco,d)均于8月达到最大值,最大光合速率(A_max)于7月达到最大值;夜间NEE随气温(T_a)呈指数增加趋势,生态系统的温度敏感系数(Q₁₀)为2.5,且土壤含水量(SWC)越高,Q₁₀值越大.非生长季NEE只与净辐射(R_n)呈显著的线性负相关,与其他环境因子无显著相关关系.生长季NEE与R_n、T_a、土壤10 cm温度(T_{s 10})等环境因子以及叶面积指数(LAI)呈显著的线性负相关,但与地上生物量(AGB)无显著相关关系.多元回归分析表明,R_n和LAI对生长季NEE的协同影响达到52%.     

环境和生物因子对黄河三角洲滨海湿地净生态系统CO2交换的影响

通过涡度相关和微气象观测技术,对黄河三角洲滨海湿地净生态系统CO₂交换(NEE)以及环境、生物因子进行了观测,探究湿地NEE变化规律及环境和生物因子对NEE的影响. 结果表明: 在日尺度上,生长季NEE呈明显“U”型曲线,非生长季变幅较小;在季节尺度上,NEE生长季波动较大,表现为碳汇,非生长季波动较小,表现为碳源;在年尺度上,滨海湿地生态系统表现为碳汇,总净固碳量为-247 g C·m^-2. 白天NEE主要受控于光合有效辐射(PAR),且生态系统表观量子产量(α)与白天生态系统呼吸(R_eco,d)均于8月达到最大值,最大光合速率(A_max)于7月达到最大值;夜间NEE随气温(T_a)呈指数增加趋势,生态系统的温度敏感系数(Q₁₀)为2.5,且土壤含水量(SWC)越高,Q₁₀值越大.非生长季NEE只与净辐射(R_n)呈显著的线性负相关,与其他环境因子无显著相关关系.生长季NEE与R_n、T_a、土壤10 cm温度(T_{s 10})等环境因子以及叶面积指数(LAI)呈显著的线性负相关,但与地上生物量(AGB)无显著相关关系.多元回归分析表明,R_n和LAI对生长季NEE的协同影响达到52%.     

森林、草地和农田典型植被蒸散量的差异

利用中国生态系统研究网络(CERN)12个试验站1990—2003年的气象资料,基于大气-土壤-植被传输(SVAT)模型模拟植被的蒸散量,分析了不同气候条件下森林、草地、农田3种植被类型蒸散的变化规律.结果表明:多年平均蒸散量从高到低依次为热带森林、亚热带农田、温带森林、温带农田和草地;对农田生态系统而言,水稻的蒸散量最高,其次为小麦,玉米最低;3种生态系统月蒸散量峰值的年际变化差异显著,森林、草地的月蒸散量峰值年际变化较稳定,农田则较明显.     

毛乌素沙地油蒿-杨柴灌丛生态系统蒸散组分分配及其影响因子

蒸散发(ET)是生态系统水分循环和能量流动的重要组成部分,准确估算ET及其各组分,对认识生态生理过程对水分平衡和植物水分利用策略的影响具有重要意义。本研究于2019年5月20日至9月15日,利用涡度相关技术和微型蒸渗仪对毛乌素沙地油蒿-杨柴灌丛生态系统ET、蒸发(E)和蒸腾(T)进行测定和估算,量化了油蒿-杨柴灌丛生态系统ET组分,并分析ET及其组分的季节特征及影响因素。结果表明: T为毛乌素沙地油蒿-杨柴灌丛生态系统生长季ET的主要组分,T/ET为53.1%。T/ET值随降水减少而升高,E/ET值随降水减少而减少,蒸散组分分配主要受降水调控。在季节尺度上,E与10 cm深处土壤含水量(SWC₁₀)和太阳净辐射(R_n)呈显著正相关,其中,SWC₁₀是E的主要影响因素;T随R_n和叶面积指数(LAI)的升高而升高,随30 cm处土壤含水量(SWC₃₀)的升高呈先升高后降低的单峰趋势,受到SWC₃₀、R_n和LAI的共同影响;水分是ET的主要影响因素。生长季蒸散/降水量(ET/P)为109.2%,5月ET/P为250.5%,表明生长季初期ET耗水部分来自非生长季降水。     

基于Penman-Monteith模型的两个蒸散模型在夏玉米农田的参数修正及性能评价

利用涡度相关系统和小气象系统对2013—2015年夏玉米生长季的蒸散量和气象数据进行实时观测,基于观测数据对以Penman-Monteith模型为基础的FAO-PM模型和KP-PM模型进行分析:首先利用2013和2014年数据对两个模型中的关键参数进行校正,其次利用两个模型对2015年夏玉米农田的日蒸散量进行计算,并与测量值对比,说明两个模型在夏玉米农田的适用性;最后采用分阶段法对KP-PM模型中的经验系数进行修正.结果表明: FAO-PM模型对2015年夏玉米农田日蒸散量的计算值更加接近测量值;利用分阶段法对KP-PM模型进行修正后,模型对日蒸散量的计算效果有了很大提高,且计算值比FAO-PM模型更接近测量值.模型中关键系数与气象条件之间有很大关系,因此利用模型进行蒸散预测时,必须先对模型进行参数校正.该研究可为其他研究人员利用模型估算蒸散量提供方法上的参考.     

太湖流域典型稻麦轮作农田生态系统碳交换及影响因素

利用涡度相关技术观测太湖流域典型稻麦轮作农田生态系统2a净生态系统碳交换(NEE)变化过程,分析其碳交换特征及影响机理,结果表明:太湖流域典型稻麦轮作农田年NEE为-749.49—-785.38 g C m-2a-1,考虑作物籽粒碳和秸秆还田后净吸收88.12 g C m-2a-1,为弱碳汇;稻/麦季日均NEE和白天NEE季节变化直接受作物植被生长影响;麦季夜间NEE与10 cm土壤温度呈显著指数关系,2012/2013年温度敏感系数(Q10)分别为3.03和2.67;当土壤水分低于田间持水量时,麦季夜间NEE主要受土壤温度影响,反之,夜间NEE受土壤温度和水分双重影响;降水对麦季夜间NEE有短时的激发效应;稻季淹水对土壤呼吸产生较明显的阻滞效应,降低了夜间NEE对土壤温度的敏感性,2012和2013年分别为1.88和1.39,稻季淹水与烤田交替变化对土壤呼吸产生明显的抑制或激发的短时效应。     

华北平原冬小麦农田蒸散量

以华北平原冬小麦农田为研究对象,采用涡度相关技术和热红外遥感技术,研究了不同环境条件下土壤含水量与农田蒸散量及作物冠层温度的关系.结果表明,冬小麦在农田郁闭(LAI≥3)、晴天和土壤相对含水量低于田间持水量65%的情况下,蒸发比值日变化正午前后出现相对较低且平稳的变化趋势.在晴天情况下,农田潜热通量与作物冠层温度日变化和季节变化均呈极显著的非线性相关关系,而冠气温差、农田相对蒸散量则与0～100 cm土层的土壤相对含水量密切相关.以13:30～14:00的平均冠层温度值T_c、日最高气温T_{a max}和日净辐射总量Rn_d为统计数据,确立了冬小麦农田日蒸散量ET_d （mm）估算简化模式参数.     

华北山区核桃液流变化特征及对不同时间尺度参考作物蒸散量的响应

核桃是华北山区重要的经济树种,但该地区水资源紧缺,因此研究核桃液流耗水变化及其对气象因子的响应对于加强核桃水分管理、确定水分承载力具有指导意义。2008—2010年连续3个生长季节观测分析了华北山区核桃单株液流变化特征,研究了不同时间尺度(日尺度、月尺度)核桃液流与气象因子、参考作物蒸散量(ET0)的关系,并利用2009—2010年数据建立了拟合方程,在此基础上用拟合方程的估算值与2008年液流实测值进行了对比验证。结果表明:核桃液流具有明显的时间变化特征,最大液流量出现在5、6月份,2008—2010年生长季节(4—9月份)核桃液流量分别为893.23 L、854.88 L和841.77 L,日平均液流量分别为4.96 L、4.75 L和4.68 L。液流年际差异主要是由年降雨量引起。核桃液流与不同时间尺度下气象因子、ET0均有较好的相关性,且随着时间尺度的增加,二者间的相关性变大。利用2009—2010年液流实测值与气象因子、ET0建立的拟合方程,对2008年液流值进行了拟合,并与实测值进行了对比验证,发现不同时间尺度下液流拟合值与实测值均具有较好的一致性。在不能实现对液流进行连续观测的情况下,可以结合当地的气象资料、ET0等因子,利用建立的拟合方程对树木液流进行估算,弥补缺失的数据,具有较高的精确性,可用于指导当地核桃水分管理。     

华北平原参考作物蒸散量变化特征及气候影响因素

参考作物蒸散量是估算作物需水量的关键因子,对指导农田灌溉是有十分重要的现实意义。在气候变化的背景下,利用Penman-Monteith方法,计算华北平原典型站点1961 2007年逐日参考作物蒸散量,并从能量平衡和动力学角度对其分解,分析年际变化和季节变化特征;结合数理统计方法,研究影响参考作物蒸散量及其构成项变化的主次气候因子,为该区农田水分管理提供更有效的科学指导。研究结果表明:在华北平原全区温度显著上升、日照时数,相对湿度,平均风速呈显著下降的背景下,绝大部分站点参考作物蒸散量及构成项呈显著下降趋势。夏季的参考作物蒸散量和辐射项值相对最高,冬季值最低;春季的空气动力学项值相对比例最高。辐射项与空气动力学项年际间呈负相关关系,春夏两季之间呈不显著正相关趋势,秋冬两季呈不显著负相关趋势。辐射项的变化主要受日照时数、风速及温度的影响,其中风速的贡献是负效应;空气动力学项的变化主要受风速、相对湿度及平均温度的影响,相对湿度的贡献是负效应。参考作物蒸散量的变化主要受日照时数、相对湿度、温度日较差和风速的综合影响。此外,降水与其呈显著负相关关系,下降幅度略高于参考作物蒸散量的变化幅度。     

华北平原冬小麦农田生态系统CO2通量特征及其影响因素

利用2014—2015年中国科学院封丘农业生态实验站涡度相关系统观测的冬小麦农田生态系统CO₂通量数据,结合试验地常规气象观测系统的气象数据,分析冬小麦4个生育期(分蘖期、越冬期、拔节期和灌浆期)内CO₂通量的日变化,研究净生态系统碳交换(NEE)的季节变化及其与气象要素的关系.结果表明: 冬小麦整个生育期内NEE为-360.15 g C·m^-2,总初级生产力总量为1920.01 g C·m^-2,冬小麦农田生态系统具有较强的固碳能力.冬小麦农田生态系统CO₂通量具有明显的日变化和季节变化特征,分蘖期表现为碳源,越冬期、拔节期和灌浆期表现为碳汇.表观初始光能利用率平均值为0.03 mg CO₂·μmol^-1,光饱和时的生态系统生产量平均值为1.53 mg CO₂·m^-2·s^-1,月平均生态系统呼吸为193.92 g C·m^-2·month^-1.冬小麦农田生态系统4个生育期NEE与光合有效辐射的相关关系均达到极显著水平.分蘖期、拔节期和灌浆期NEE与饱和水汽压差的相关关系极显著,越冬期达显著水平.冬小麦分蘖期、越冬期和灌浆期NEE日总量与土壤温度呈正相关,拔节期呈负相关关系.     

华北平原冬小麦农田生态系统CO_2通量特征及其影响因素

利用2014—2015年中国科学院封丘农业生态实验站涡度相关系统观测的冬小麦农田生态系统CO_2通量数据,结合试验地常规气象观测系统的气象数据,分析冬小麦4个生育期(分蘖期、越冬期、拔节期和灌浆期)内CO_2通量的日变化,研究净生态系统碳交换(NEE)的季节变化及其与气象要素的关系.结果表明:冬小麦整个生育期内NEE为-360.15g C·m^-2,总初级生产力总量为1920.01 g C·m^-2,冬小麦农田生态系统具有较强的固碳能力.冬小麦农田生态系统CO_2通量具有明显的日变化和季节变化特征,分蘖期表现为碳源,越冬期、拔节期和灌浆期表现为碳汇.表观初始光能利用率平均值为0.03 mg CO_2·μmol^-1,光饱和时的生态系统生产量平均值为1.53 mg CO_2·m^-2·s^-1,月平均生态系统呼吸为193.92g C·m^-2·month^-1.冬小麦农田生态系统4个生育期NEE与光合有效辐射的相关关系均达到极显著水平.分蘖期、拔节期和灌浆期NEE与饱和水汽压差的相关关系极显著,越冬期达显著水平.冬小麦分蘖期、越冬期和灌浆期NEE日总量与土壤温度呈正相关,拔节期呈负相关关系.     

基于改进的双作物系数法估算辽河三角洲芦苇湿地蒸散量

以辽河三角洲湿地芦苇群落为研究对象,利用涡动相关通量、小气候梯度要素、群落内水面蒸发量以及芦苇群落生长参数等数据,基于FAO-56模型的双作物系数法,分别计算作物系数K_c、基础作物系数K_(cb)和水面蒸发系数K_w,分析其日变化动态及主导影响因子,建立基于生物因子和环境因子的小时尺度双作物系数模型。结果如下:(1)芦苇生长初期,K_c和K_(cb)的日变化呈现早晚高、中午略低的多峰波动曲线;在快速生长期和稳定生长期,K_c和K_(cb)白天波动幅度较小,早晚波动幅度较大;生长末期,K_c和K_(cb)夜晚波动幅度较大,白天呈现多峰波动曲线;K_w白天较小、夜晚较高,生长初期白天的数值显著高于其他时期。(2)相关分析表明,气温、相对湿度、风速、株高和叶面积指数是K_c、K_(cb)和K_w的影响因子;基于生物因子和环境因子重新构建双作物系数模型,基于改进的双作物系数法模拟芦苇群落蒸散,决定系数R~2达0.894。(3)利用改进的双作物系数模型和FAO-56模型,对辽河三角洲芦苇群落的蒸发与蒸腾过程进行模拟,实现芦苇群落蒸发过程与蒸腾过程的分离,解决了实际观测无法直接获取芦苇群落蒸腾量的问题,同时提高了芦苇群落蒸散的模拟精度。(4)调整了FAO推荐的芦苇单作物系数常数值,调整后的作物系数更适用于辽河三角洲芦苇湿地。     

Simulating Crop Evapotranspiration Response under Different Planting Scenarios by Modified SWAT Model in an Irrigation District,Northwest China

Modelling crop evapotranspiration (ET) response to different planting scenarios in an irrigation district plays a significant role in optimizing crop planting patterns, resolving agricultural water scarcity and facilitating the sustainable use of water resources. In this study, the SWAT model was improved by transforming the evapotranspiration module. Then, the improved model was applied in Qingyuan Irrigation District of northwest China as a case study. Land use, soil, meteorology, irrigation scheduling and crop coefficient were considered as input data, and the irrigation district was divided into subdivisions based on the DEM and local canal systems. On the basis of model calibration and verification, the improved model showed better simulation efficiency than did the original model. Therefore, the improved model was used to simulate the crop evapotranspiration response under different planting scenarios in the irrigation district. Results indicated that crop evapotranspiration decreased by 2.94% and 6.01% under the scenarios of reducing the planting proportion of spring wheat (scenario 1) and summer maize (scenario 2) by keeping the total cultivated area unchanged. However, the total net output values presented an opposite trend under different scenarios. The values decreased by 3.28% under scenario 1, while it increased by 7.79% under scenario 2, compared with the current situation. This study presents a novel method to estimate crop evapotranspiration response under different planting scenarios using the SWAT model, and makes recommendations for strategic agricultural water management planning for the rational utilization of water resources and development of local economy by studying the impact of planting scenario changes on crop evapotranspiration and output values in the irrigation district of northwest China.     

基于Budyko假设的潮河流域气候和植被变化对实际蒸散发的影响研究

由于参数较少且具有明确的物理学意义,基于水热平衡理论的Budyko假设常用于定量分析以及评价气候变化和植被变化对实际蒸散发的影响,对研究流域水量平衡和能量分配具有重要意义。依据位于我国北方密云水库上游的潮河流域1961—2015年的水文气象数据,选取了4种基于Budyko假设的模型来研究潮河流域水热耦合平衡关系,确定了该流域最适用模型以及模型参数最优值,并且采用情景设置法分析了流域实际蒸散发对气候以及植被变化的响应。结果表明:(1)与经典Budyko模型相比,采用流域下垫面参数修正的Budyko模型计算实际蒸散发的精度更高。其中,傅抱璞模型精度最高,决定系数、相对误差、纳什效率系数和均方根误差分别为0.85、4.30%、0.82和27.66 mm;(2)对傅抱璞模型下垫面参数ω进行优化,确定适用于潮河流域的模型参数取值为2.54,优化后的傅抱璞模型能够更好地反映流域实际蒸散发的变化特征;(3)情景模拟表明,气候变化和植被变化的共同作用导致潮河流域实际蒸散发的上升。其中,气候变化是引起流域蒸散发变化的主要驱动因素。     

Seventeenth century organic agriculture in China: II. Energy flows through an agroecosystem in Jiaxing region

The energy flows in a seventeenth century agroecosystem in Jiaxing region of eastern China were analyzed on the basis of historical data. The agroecosystem included cropping, mulberry-silkworm, livestock, and fishing systems. In terms of energy, the agroecosystem was sustainable. Human labor provided all the power with inputs of about 3700 hr per hectare of farmland. Most or 70% of the labor was expended in the cropping system. Human and animal manure provided most of the nutrients for crop and mulberry production. About two-thirds of the total manure was used in crop production and one-third in the mulberry plantations. The only fossil energy input was a few hand tools. Approximately 55% of the grain was consumed directly by local residents, about one-third of the grain was used to make an alcohol drink and produce distillers'grains, which was fed to pigs, and only 2% of the grains were exported outside the agroecosystem. About two-thirds of the harvested crop residues were used as household fuel, while the remainder was returned to the field as an organic fertilizer. Pork accounted for 85% and silk cocoons 14% of the total animal products produced. Even though the agroecosystem was generally sustainable in terms of energy, the major environmental problem was that two-thirds of the harvested crop residues were used for household fuel. This reduced nutrient cycling in the system. Insufficient land was available to produce fuelwood;thus, crop residues were the primary source of fuel for the people.