Modeling Study on the Characteristics of Laminar and Turbulent Argon Plasma Jets Impinging Normally upon a Flat Plate in Ambient Air

Modeling study is performed to compare the flow and heat transfer characteristics of laminar and turbulent argon thermal-plasma jets impinging normally upon a flat plate in ambient air. The combined-diffusion-coefficient method and the turbulence-enhanced combined-diffusion-coefficient method are employed to treat the diffusion of argon in the argon-air mixture for the laminar and the turbulent cases, respectively. Modeling results presented include the flow, temperature and argon concentration fields, the air mass flow-rates entrained into the impinging plasma jets, and the distributions of the heat flux density on the plate surface. It is found that the formation of a radial wall jet on the plate surface appreciably enhances the mass flow rate of the ambient air entrained into the laminar or turbulent plasma impinging-jet. When the plate standoff distance is comparatively small, there exists a significant difference between the laminar and turbulent plasma impinging-jets in their flow fields due to the occurrence of a large closed recirculation vortex in the turbulent plasma impinging-jet, and no appreciable difference is found between the two types of jets in their maximum values and distributions of the heat flux density at the plate surface. At larger plate standoff distances, the effect of the plate on the jet flow fields only appears in the region near the plate, and the axial decaying-rates of the plasma temperature, axial velocity and argon mass fraction along the axis of the laminar plasma impinging-jet become appreciably less than their turbulent counterparts.

Pool Boiling In Microgravity: Recent Results And Perspectives For The Project Depa-Sj10

Two research projects on pool boiling in microgravity have been conducted aboard the Chinese recoverable satellites. Ground-based experiments have also been performed both in normal gravity and in short-term microgravity in the Drop Tower Beijing. Steady boiling of R113 on thin platinum wires was studied with a temperature-controlled heating method, while quasi-steady boiling of FC-72 on a plane plate was investigated with an exponentially increasing heating voltage. In the first case, slight enhancement of heat transfer is observed in microgravity, while diminution is evident for high heat flux in the second one. Lateral motions of bubbles on the heaters are observed before their departure in microgravity. The surface oscillation of the merged bubbles due to lateral coalescence between adjacent bubbles drives it to detach from the heaters. The Marangoni effect on the bubble behavior is also discussed. The perspectives for a new project DEPA-SJ10, which has been planned to be flown aboard the Chinese recoverable satellite SJ-10 in the future, are also presented.

探针移动速度对等离子体射流热流密度测量结果的影响

采用自行研制的中心嵌有铜柱感应件的小尺寸杆状热流探针,在低扰动条件下,对射入大气环境的纯氩层流等离子体射流传向铜探头表面的热流密度进行了动态测量。结果表明,在射流最高温度16500 K、最大轴向速度850 m/s、探针垂直于射流流动方向的移动速度130～260 mm/s的实验参数范围内,随着探针移动速度的提高,测得的热流密度值减小;射流温度和速度越高,探针移动速度对热流密度测量值的影响越大。

高超声速飞行器测热试验研究

在JF-8脉冲风洞中,来流马赫数Ma＝8.0,来流单位长度雷诺数Re/L=1.47×107和2.52×107(1/m)两种试验条件下,对高超声速飞行器1/20缩尺模型进行了表面气动热的测量.模型攻角α=0°,10°,15°,20°,25°和30°.试验给出机身对称面、翼前缘、立尾前缘等处的热流率分布.机头部分最大热流率与由Fay-Riddell公式计算的驻点热流Q0率接近,翼前缘最大热流率在全机身中最大,约为Q0的2倍,因此翼前缘的热环境是最严酷的.

Subcooled Pool Boiling at Microgravity: Preliminary Ground-Based Experimental Results

A temperature-controlled pool boiling (TCPB) device was developed to perform pool boiling heat transfer studies at both normal gravity on Earth and microgravity in the drop tower Beijing and aboard a Chinese recovery satellite. Two platinum wires of 60 ?m in diameter were simultaneously used as heaters and thermometers. The lengths were 30 mm and 40 mm, respectively. The ends of wires were soldered with copper poles to provide low resistance paths for the electric current. The heater resistance, and thus the heater temperature, was kept constant by a feedback circuit similar to that used in constant-temperature hot-wire anemometry. The fluid was R113 at 0.1 Mpa and subcooled by 30 ?C nominally for all cases. The results of the experiments at normal gravity were presented. Four modes, namely single-phase convection, nucleate boiling, transition two-mode boiling, and film boiling were observed. A few data obtained from several preliminary experiments at microgravity in the drop tower Beijing were also presented. A slight increase of the heat flux was obtained.

Experimental Study on Subcooled Pool Boiling in Microgravity Utilizing

Pool Boiling in Microgravity: Recent Results and Perspectives for the Project DEPA-SJ10

Two research projects on pool boiling in microgravity have been conducted aboard the Chinese recoverable satellites. Ground-based experiments have also been performed both in normal gravity and in short-term microgravity in the Drop Tower Beijing. Steady boiling of R113 on thin platinum wires was studied with a temperature-controlled heating method, while quasi-steady boiling of FC-72 on a plane plate was investigated with an exponentially increasing heating voltage. In the first case, slight enhancement of heat transfer is observed in microgravity, while diminution is evident for high heat flux in the second one. Lateral motions of bubbles on the heaters are observed before their departure in microgravity. The surface oscillation of the merged bubbles due to lateral coalescence between adjacent bubbles drives it to detach from the heaters. The Marangoni effect on the bubble behavior is also discussed. The perspectives for a new project DEPA-SJ10, which has been planned to be flown aboard the Chinese recoverable satellite SJ-10 in the future, are also presented.

Daily scale wintertime sea surface temperature and IPC-Navidad variability in the southern Bay of Biscay from 1981 to 2010

The combination of remotely sensed gappy Sea surface temperature (SST) images with the missing data filling DINEOF (data interpolating empirical orthogonal functions) technique, followed by a principal component analysis of the reconstructed data, has been used to identify the time evolution and the daily scale variability of the wintertime surface signal of the Iberian Poleward Current (IPC), or Navidad, during the 1981-2010 period. An exhaustive comparison with the existing bibliography, and the vertical temperature and salinity profiles related to its extremes over the Bay of Biscay area, show that the obtained time series accurately reflect the IPC-Navidad variability. Once a time series for the evolution of the SST signal of the current over the last decades is well established, this time series is used to propose a physical mechanism in relation to the variability of the IPC-Navidad, involving both atmospheric and oceanic variables. According to the proposed mechanism, an atmospheric circulation anomaly observed in both the 500 hPa and the surface levels generates atmospheric surface level pressure, wind-stress and heat-flux anomalies. In turn, those surface level atmospheric anomalies induce mutually coherent SST and sea level anomalies over the North Atlantic area, and locally, in the Bay of Biscay area. These anomalies, both locally over the Bay of Biscay area and over the North Atlantic, are in agreement with several mechanisms that have separately been related to the variability of the IPC-Navidad, i.e. the south-westerly winds, the joint effect of baroclinicity and relief (JEBAR) effect, the topographic beta effect and a weakened North Atlantic gyre.

Planar reflection of gaseous detonations

Pipes containing flammable gaseous mixtures may be subjected to internal detonation. When the detonation normally impinges on a closed end, a reflected shock wave is created to bring the flow back to rest. This study built on the work of Karnesky (2010) and examined deformation of thin-walled stainless steel tubes subjected to internal reflected gaseous detonations. A ripple pattern was observed in the tube wall for certain fill pressures, and a criterion was developed that predicted when the ripple pattern would form. A two-dimensional finite element analysis was performed using Johnson-Cook material properties; the pressure loading created by reflected gaseous detonations was accounted for with a previously developed pressure model. The residual plastic strain between experiments and computations was in good agreement.

During the examination of detonation-driven deformation, discrepancies were discovered in our understanding of reflected gaseous detonation behavior. Previous models did not accurately describe the nature of the reflected shock wave, which motivated further experiments in a detonation tube with optical access. Pressure sensors and schlieren images were used to examine reflected shock behavior, and it was determined that the discrepancies were related to the reaction zone thickness extant behind the detonation front. During these experiments reflected shock bifurcation did not appear to occur, but the unfocused visualization system made certainty impossible. This prompted construction of a focused schlieren system that investigated possible shock wave-boundary layer interaction, and heat-flux gauges analyzed the boundary layer behind the detonation front. Using these data with an analytical boundary layer solution, it was determined that the strong thermal boundary layer present behind the detonation front inhibits the development of reflected shock wave bifurcation.

Simulations and mechanisms of subtropical low-cloud response to climate change

This thesis focuses on improving the simulation skills and the theoretical understanding of the subtropical low cloud response to climate change.

First, an energetically consistent forcing framework is designed and implemented for the large eddy simulation (LES) of the low-cloud response to climate change. The three representative current-day subtropical low cloud regimes of cumulus (Cu), cumulus-over-stratocumulus, and stratocumulus (Sc) are all well simulated with this framework, and results are comparable to the conventional fixed-SST approach. However, the cumulus response to climate warming subject to energetic constraints differs significantly from the conventional approach with fixed SST. Under the energetic constraint, the subtropics warm less than the tropics, since longwave (LW) cooling is more efficient with the drier subtropical free troposphere. The surface latent heat flux (LHF) also increases only weakly subject to the surface energetic constraint. Both factors contribute to an increased estimated inversion strength (EIS), and decreased inversion height. The decreased Cu-depth contributes to a decrease of liquid water path (LWP) and weak positive cloud feedback. The conventional fixed-SST approach instead simulates a strong increase in LHF and deepening of the Cu layer, leading to a weakly negative cloud feedback. This illustrates the importance of energetic constraints to the simulation and understanding of the sign and magnitude of low-cloud feedback.

Second, an extended eddy-diffusivity mass-flux (EDMF) closure for the unified representation of sub-grid scale (SGS) turbulence and convection processes in general circulation models (GCM) is presented. The inclusion of prognostic terms and the elimination of the infinitesimal updraft fraction assumption makes it more flexible for implementation in models across different scales. This framework can be consistently extended to formulate multiple updrafts and downdrafts, as well as variances and covariances. It has been verified with LES in different boundary layer regimes in the current climate, and further development and implementation of this closure may help to improve our simulation skills and understanding of low-cloud feedback through GCMs.

A kinetic theory description for external spherical flows with arbitrary Knudsen number by a moment method

The Maxwell integral equations of transfer are applied to a series of problems involving flows of arbitrary density gases about spheres. As suggested by Lees a two sided Maxwellian-like weighting function containing a number of free parameters is utilized and a sufficient number of partial differential moment equations is used to determine these parameters. Maxwell's inverse fifth-power force law is used to simplify the evaluation of the collision integrals appearing in the moment equations. All flow quantities are then determined by integration of the weighting function which results from the solution of the differential moment system. Three problems are treated: the heat-flux from a slightly heated sphere at rest in an infinite gas; the velocity field and drag of a slowly moving sphere in an unbounded space; the velocity field and drag torque on a slowly rotating sphere. Solutions to the third problem are found to both first and second-order in surface Mach number with the secondary centrifugal fan motion being of particular interest. Singular aspects of the moment method are encountered in the last two problems and an asymptotic study of these difficulties leads to a formal criterion for a "well posed" moment system. The previously unanswered question of just how many moments must be used in a specific problem is now clarified to a great extent.

Stationary Eddies and Zonal Variations of the Global Hydrological Cycle in a Changing Climate

This thesis advances our physical understanding of the sensitivity of the hydrological cycle to global warming. Specifically, it focuses on changes in the longitudinal (zonal) variation of precipitation minus evaporation (P - E), which is predominantly controlled by planetary-scale stationary eddies. By studying idealized general circulation model (GCM) experiments with zonally varying boundary conditions, this thesis examines the mechanisms controlling the strength of stationary-eddy circulations and their role in the hydrological cycle. The overarching goal of this research is to understand the cause of changes in regional P - E with global warming. An understanding of such changes can be useful for impact studies focusing on water availability, ecosystem management, and flood risk.

Based on a moisture-budget analysis of ERA-Interim data, we establish an approximation for zonally anomalous P - E in terms of surface moisture content and stationary-eddy vertical motion in the lower troposphere. Part of the success of this approximation comes from our finding that transient-eddy moisture fluxes partially cancel the effect of stationary-eddy moisture advection, allowing divergent circulations to dominate the moisture budget. The lower-tropospheric vertical motion is related to horizontal motion in stationary eddies by Sverdrup and Ekman balance. These moisture- and vorticity-budget balances also hold in idealized and comprehensive GCM simulations across a range of climates.

By examining climate changes in the idealized and comprehensive GCM simulations, we are able to show the utility of the vertical motion P - E approximation for splitting changes in zonally anomalous P - E into thermodynamic and dynamic components. Shifts in divergent stationary-eddy circulations dominate changes in zonally anomalous P - E. This limits the local utility of the "wet gets wetter, dry gets drier” idea, where existing P - E patterns are amplified with warming by the increase in atmospheric moisture content, with atmospheric circulations held fixed. The increase in atmospheric moisture content manifests instead in an increase in the amplitude of the zonally anomalous hydrological cycle as measured by the zonal variance of P - E. However, dynamic changes, particularly the slowdown of divergent stationary-eddy circulations, limit the strengthening of the zonally anomalous hydrological cycle. In certain idealized cases, dynamic changes are even strong enough to reverse the tendency towards "wet gets wetter, dry gets drier” with warming.

Motivated by the importance of stationary-eddy vertical velocities in the moisture budget analysis, we examine controls on the amplitude of stationary eddies across a wide range of climates in an idealized GCM with simple topographic and ocean-heating zonal asymmetries. An analysis of the thermodynamic equation in the vicinity of topographic forcing reveals the importance of on-slope surface winds, the midlatitude isentropic slope, and latent heating in setting the amplitude of stationary waves. The response of stationary eddies to climate change is determined primarily by the strength of zonal surface winds hitting the mountain. The sensitivity of stationary-eddies to this surface forcing increases with climate change as the slope of midlatitude isentropes decreases. However, latent heating also plays an important role in damping the stationary-eddy response, and this damping becomes stronger with warming as the atmospheric moisture content increases. We find that the response of tropical overturning circulations forced by ocean heat-flux convergence is described by changes in the vertical structure of moist static energy and deep convection. This is used to derive simple scalings for the Walker circulation strength that capture the monotonic decrease with warming found in our idealized simulations.

Through the work of this thesis, the advances made in understanding the amplitude of stationary-waves in a changing climate can be directly applied to better understand and predict changes in the zonally anomalous hydrological cycle.

Aplicação do algoritmo SEBAL para estimativa da evapotranspiração na mesorregião do sul goiano.

As preocupações com o uso da terra têm permeado inúmeros estudos científicos, no âmbito nacional e internacional, voltados para a avaliação dos impactos ambientais causados pelas atividades agropecuárias. Alguns processos do ciclo hidrológico, a exemplo da evapotranspiração apresentam modificações consideráveis, devido às constantes mudanças nos usos dos solos. Desta forma, o presente trabalho busca destacar o problema das rápidas e intensas mudanças no uso do solo oriundas da expansão da atividade agropecuária e seus impactos ao meio ambiente, especialmente sobre o processo da evapotranspiração regional, na mesorregião do Sul Goiano, região típica de cerrado, localizada no centro-oeste brasileiro. A aplicação do algoritmo Surface Energy Balance Algorithm for Land - SEBAL consistiu o cerne da metodologia utilizada, com vista à estimativa dos fluxos de energia e da evapotranspiração em escala regional, obtidos com base no equacionamento do balanço de energia à superfície, complementado por dados de temperatura do ar e velocidade do vento adquiridos em estações meteorológicas (PCDs) instaladas na área de estudo. Foram utilizados dados do sensor MODIS/TERRA dos anos 2006, 2007, 2008, 2009 e 2010. O algoritmo foi testado em sua forma clássica e modificado por alterações nos critérios de seleção dos pixels âncoras, utilizados no procedimento da estimativa do fluxo de calor sensível. Pode-se concluir que a alteração dos critérios influenciou positivamente os resultados obtidos e que os valores da evapotranspiração, na região estudada, indicaram a potencialidade da metodologia empregada para o monitoramento sistemático dos componentes do balanço de energia em escala regional.

内蒙古草原温度梯度上3个典型群落的能量平衡研究

通量测量点的能量收支总是表现出不平衡，即使在地势平坦、植被分布均一、稀疏植被下垫面的情况下也有约30%的能量失衡状况。能量平衡闭合 (EBC) 问题在验证涡度相关系统质量方面，得到了广泛的关注。实验在内蒙古草原3个地点，通过涡度相关系统附近移动能量平衡系统的测定手段，采用能量平衡余项法和最小二乘 (OLS) 线性回归法，研究了土壤热通量、净辐射及可供能量空间变异对能量平衡闭合的影响。 结果显示，EBC 在三个研究点的平均余项为8~19 W m-2，OLS 斜率为0.83~0.96。EBC 在土壤湿润情况的站点要高于干旱站点。 土壤热通量的空间变异三站点平均为白天48 W m-2 (占同时间 Rn 的13%)，夜间15 W m-2 (34%)，平均29 W m-2 (24%)。通过8个工作站的测量，这个变异会造成9% (从0.93到1.01) 的 OLS 斜率差别。夜间的能量平衡不闭合可以由土壤热通量的空间变异解释。如果在本研究的的三个草原站点上忽略了土壤热通量，则会造成较大的余项 (峰值时110 W m-2) 产生，从而使 OLS 斜率增大23%。特别是通量板埋置在地面以下30 mm处时，上层的土壤热储部分占到全部土壤热通量的50%，这不仅影响到 EBC 的大小，更起到调节土壤热通量与“真实的波形” 相一致的作用。如果该部分热储被忽略掉，EBC 余项会增加60 W m-2，OLS 斜率也会变化 (减少) 9%。用大尺度多点测量与涡度塔附近的小尺度测定相比较，后者表现出稍高的闭合率，即 OLS 斜率增加4%。 相对于土壤热通量，净辐射的空间变异较小，三站点平均为白天17 W m-2 (5%)，夜间7 W m-2 (13%)，平均 12 W m-2 (5%)。可以引起3% (从0.88到0.91) 的 OLS 斜率差异。研究结果还表明，风速校正应该在 Q7.1 净辐射仪中应用，校正后的结果与 CNR1 的结果在白天有吻合较好，但在其它时段仍有较大差异，特别是在夜间，风速校正基本不起作用，使得两种仪器间差异达20 W m-2。比较表明，风速校正可以提高白天 Rn 的6%，仅降低夜间0.3%。因此，无论是用余项法还是用 OLS 线性回归法，在比较使用不同仪器的站点间的闭合状况时 (本研究的结果适用于草地 Q7.1 与 CNR1 间的比较)，可以用9：00-15：00 h 时段的数据进行比较，这样可以避免因使用不同仪器的差异所造成的影响。用该时段的数据进行比较，仪器间的差异余项法小于6 W m-2，OLS 法小于3%。 受可供能量空间变异影响，三个站点平均 EBC 的不确定性为白天66 W m-2 (19%)，夜间23 W m-2 (50%)，平均42 W m-2 (36%);或者改变 OLS 斜率11%。用最大值和最小值来衡量，EBC 最大不确定性，正午时在站点 I、II 和 III 中分别为81，114和91 W m-2。故在探讨能量平衡和能量平衡闭合问题时，必须充分考虑到这种不确定性，否则会产生偏差，或者得出错误结论。 研究还表明，即使考虑到白天所有可供能量的最大不确定性，仍然不能使能量平衡闭合。中午 (12：00 h)，站点 I，II 和 III 仍然有14±15，48±12和47±14 W m-2 的失衡不能够归因于可供能量的空间不确定性。因此，其它影响因素也需进行细致的探讨。 在两站点不同测量深度土壤热通量结果的差异性比较实验中，无论在站点 I 还是 III，均表现出一致的结论，即随通量板布置深度加深，其测量结果会越高，与浅层布置的相比，差别可高达150 W m-2。深层土壤热通量的计算仍是个难题，需进一步研究。 在不同植被结构对净辐射测定影响的实验中发现，随刈割强度增加，一天中大部分时间净辐射均减少。正午依次为413，395和388 W m-2。无论正午还是全天合计，重度刈割地点的净辐射均比不刈割对照处理少6%，而且，在整个生长季也少6%，约合40,000 W m-2。测量高度不同，不同处理间对测定结果影响不同：刈割处理中，由于下垫面较均一，结果相差不显著;而对照则表现出较高的差异，用两配对样本T检验表明差异达到极显著 (P<0.000，9:30-15:00 h data)。当使用不同新旧程度的 domes 时，对净辐射结果会产生明显的影响。新 domes 的测量结果白天明显高，晚上明显低，使用了11个月的旧 domes，峰值时，白天低估25 W m-2，晚上高估10 W m-2。说明该差异在进行能量平衡闭合计算时，不能忽略。而全用新的和全用旧的进行比较，晚上仅有2-3 W m-2差异。 考虑生态系统中非生物因子对干扰条件下生物多样性动态和功能的影响，有助于更精确地阐明生物多样性-稳定性功能的关系。为此，设计了一个单因子刈割实验——内蒙古地区一种广泛存在的土地利用方式。主要目的是研究不同强度刈割影响下，微气候变量特别是能量平衡各分量和群落结构的变化及二者的关系。连续4年刈割，占第一位的优势种明显由低矮半灌木冷蒿 (Artemisia frigida) 取代了高大丛生禾草克氏针茅 (Stipa kylovii)。重度刈割下，针茅的盖度、生物量和丛重，群落叶面积、绿色生物量、凋落物量和群落高度一致低于轻度刈割/不刈割处理。微气候由于群落特征的这些变化也呈规律性变化。与对照相比，重度刈割降低了生长季土壤含水量的47.5%，但中午和日均土壤表面温度分别增加了7.4和1.2 °C，并且增加地表下2 cm土壤温度日较差 (日最高与最低温度之差) 4.2 °C。刈割处理由于凋落物少、反射强而表现出较低的净辐射，但土壤热通量显著提高，表现为土层加热和冷却快。因此，重度刈割处理较对照降低了可供能量8%，约合52,000 W m-2 。不同刈割强度间来看，NPP 或 LAI 与土壤热通量和净辐射的比值 (G/Rn) 以及波文比 (H/LE) 间呈负相关。重度刈割处理感热通量显著提高，但潜热通量在处理间差异不显著，表明未刈割处理虽然冠层伸展大，但是并没有导致更大的水分亏缺。未刈割处理增加了抵抗物种改变的能力，而刈割处理在连续一年一割的第四年显著增加了物种数，可能与因刈割影响而导致的群落结构与微气候的改变有关。本研究表明，未刈割处理可以减轻高温干旱季节的高温和干旱胁迫，表现出对环境变化的高抵抗性。未刈割处理的凋落物层和较高的垂直结构所形成的遮荫，可以形成一个阻挡蒸发的篱笆，这是维持其水分的保证。因此，为了恢复退化草原生态系统功能，需要修复能导致微气候变化的植物群落结构，否则难以成功。 本研究立足于原创性的实验研究，在中国特有的自然草原生态系统上开展，结合不同温度梯度的三地区涡度相关系统进行了能量平衡闭合的移动比较实验，以及结合常用土地利用方式的定点能量平衡实验。在翔实的数据基础上，为涡度相关方法的陆地表面能量平衡失衡问题提供了解释。增加了对两个主要能量流——土壤热通量和净辐射空间变异规律的认识，研究对于能量平衡和湍流通量相关研究是有价值的。在三个代表性地区首次利用多个净辐射仪和土壤热通量板的结果与三个标准的涡度相关系统进行了比较，这类量化失衡原因的相关研究应受到高度重视并进一步拓展，以提高对能量失衡的认识，进而推进水热和碳循环研究向更深层次发展。

玉米农田水热碳通量动态及其环境控制机理研究

陆地生态系统与大气之间的水热碳交换是物质、能量循环的关键过程，一直以来都为研究者们所关注。进入20 世纪以来，特别是随着人们对全球气候变暖的逐步认识，气候变化对水热碳交换过程的影响及其对气候变化的响应研究更加备受关注。本研究以2004~2006 年近三年的涡度相关系统连续观测数据为依托，分析了雨养玉米农田水热碳通量的动态及其影响因子。研究表明，玉米农田水热通量(WHF) 呈显著的单峰型日变化， 日最大值出现在正午12:00~13:00，WHF 变化同步。潜热通量(LE)的季节变化规律与日变化相似，冬季小夏季大，年最大值与最小值分别出现在7 月和1 月。显热通量(Hs) 季节变化也呈单峰型，但年最大值出现在5 月，这主要与降水以及作物生长有关。半小时尺度上，WHF 主要受辐射控制，而日峰值受辐射峰值以及植被生长的双重影响；日尺度上，只要有降水过程，Hs 就会随土壤水分的增大而减小，降水停止后逐渐恢复。而降水对LE 的影响受到可用能量(AE)的干扰，表现出复杂的变化趋势。总的来说，降水持续时间越长AE 越少，对LE 的抑制越大；季节尺度上，WHF 受热量与水分的双重制约。Hs 随着天气回暖后第一次较大降水过程的出现呈现明显下降，而LE 则呈现相反的变化趋势。随着雨季到来和作物的生长，Hs 在7 月出现低谷，而LE 呈现相反的趋势随着降水量的增加而增大；年际间WHF 的分布规律大体一致，但因气象条件等的差异，特别是降水的差异造成年际间WHF 略有不同。在不同水文年型下，水分因子的影响作用有显著差异，且WHF 对热量与水分条件变化的敏感程度也不相同。欠水年，水分因子的作用更显著，是制约WHF 变化的主要控制因子，WHF 对水分的变化更敏感；而丰水年，水分因子的影响减弱，热量的盈亏决定着WHF 变化的主要方向。在不同水文年型下，水热碳通量对水热条件的变化表现出不同的响应方式，为研究生态系统对气候变化的响应提供了参考。 净碳(C)吸收期，玉米农田净碳交换(NEE)呈显著的日变化，在日出以后由CO2 释放转变为CO2 吸收，12:30 左右达到一天中的吸收峰值，日落前出现相反的转换。而净C 释放期内，NEE 均为正值且无明显日变化。NEE 季节变化也呈单峰型二次曲线，在7 月下旬或8 月上旬达到年最大吸收率。根据NEE 的正负，一年分为三个阶段：两个C 排放期与一个C 吸收期。一般C 吸收期从6月开始到9 月结束，此前此后均为C 排放期。在半小时、日时间尺度上，光通量密度(PPFD)与NEE 有着相似的变化规律，是控制NEE 的主要因子；在日、季节尺度上，叶面积指数(LAI)和气孔导度(gs)是影响NEE 的主要生物因子，且gs 的影响程度随着发育期的变化而变化，而不同年份间LAI 对NEE 的影响没有显著的差异。几乎在所有时间步长上，土壤温度(Ts)均为生态系统呼吸(Re)的主要控制因子，时间尺度愈短，二者的相关性愈好。总的来说，在较短时间尺度上，高PPFD 与夏季低温将会促进C 的吸收，有利于C 累积。 玉米农田日最大净C 吸收速率(NEEmax, daily)以及吸收释放转换点(NEE=0)均受PPFD 控制。NEEmax, daily 出现时间与PPFDmax, daily 出现时间几乎完全一致，当PPFD 达到1 日内极大值时，净C 吸收也相应达到了日最大值。但NEEmax, daily的量值还受到其它因子的影响。当水分条件充足时，还将受到LAI、gs 等生物因子的控制。NEE 由正转为负的转换点也是由PPFD 决定。当PPFD 稳定大于PPFD*( PPFD*=100 μmol•m-2s -1)时，净C 吸收开始；当PPFD 稳定小于PPFD*时，净C 吸收由此结束。1 日内，PPFD 稳定通过PPFD*之间的时间间隔决定了日净C 吸收的时间长度。日净C 吸收的时间越长，吸收量也越大，且有明显的季节变化，7 月最长9 月最短。 按照热量水分状况将三年分组，分为I 组(水分状况相似，热量条件不同)与II 组(热量条件相似，水分状况不同)。 I 组年际间PPFD 波动是造成C 交换格局变化的关键原因。而II 组年际间C 交换格局不同是由降水量及其不同分布引起的土壤含水量(SWC)变化是造成。SWC 可以解释年际间NEE 变异的97%，而大气水汽压亏缺(VPD)可以解释30.7%；温度因子通过影响C 收支中的呼吸项，间接影响着生态系统的NEE，它可以解释年际间NEE 变异的73.9%，也是造成年际间C 交换格局不同的原因之一；另外，PPFD 和发育期早晚以及净C吸收期长度等也同样影响着C 交换格局的变化。综合两组情况来看，由水分条件年际变化引起的NEE 的波动大于能量年际变化引起的波动。总之，在较长时间尺度上，NEE 对SWC 变化比其对PPFD 变化更敏感，说明在半干旱地区土壤水分条件仍然是决定C 交换格局的主导因子。 NEE 与LE 呈线性相关，它们之间的相关性主要受温度和NEE 的控制，温度越高，二者的相关性越弱，而NEE 越大二者相关性越好。同时，作物蒸腾与土壤蒸发的比例也是影响NEE 与LE 之间关系的主要因素。蒸腾作用所占的比例越大，二者的线性关系越显著，而土壤蒸发比例越大，二者的相关性越弱。总的来说，NEE 与LE 之间的线性关系有明显的季节变化，生长季好于非生长季，夏天好于冬天。 总之，雨养玉米农田水热碳通量既具有其它农田生态系统共有的动态特征，也具有其特有特征。