Assessing the performance of a large-scale irrigation system by estimations of actual evapotranspiration obtained by Landsat satellite images resampled with cubic convolution

Remote sensing techniques allow monitoring the Earth surface and acquiring worthwhile information that can be used efficiently in agro-hydrological systems. Satellite images associated to computational models represent reliable resources to estimate actual evapotranspiration fluxes, ET_a, based on surface energy balance. The knowledge of ET_a and its spatial distribution is crucial for a broad range of applications at different scales, from fields to large irrigation districts. In single plots and/or in irrigation districts, linking water volumes delivered to the plots with the estimations of remote sensed ET_a can have a great potential to develop new cost-effective indicators of irrigation performance, as well as to increase water use efficiency. With the aim to assess the irrigation system performance and the opportunities to save irrigation water resources at the “SAT Llano Verde” district in Albacete, Castilla-La Mancha (Spain), the Surface Energy Balance Algorithm for Land (SEBAL) was applied on cloud-free Landsat 5 Thematic Mapper (TM) images, processed by cubic convolution resampling method, for three irrigation seasons (May to September 2006, 2007 and 2008). The model allowed quantifying instantaneous, daily, monthly and seasonal ET_a over the irrigation district. The comparison between monthly irrigation volumes distributed by each hydrant and the corresponding spatially averaged ET_a, obtained by assuming an overall efficiency of irrigation network equal to 85%, allowed the assessment of the irrigation system performance for the area served by each hydrant, as well as for the whole irrigation district. It was observed that in all the investigated years, irrigation volumes applied monthly by farmers resulted generally higher than the corresponding evapotranspiration fluxes retrieved by SEBAL, with the exception of May, in which abundant rainfall occurred. When considering the entire irrigation seasons, it was demonstrated that a considerable amount of water could have been saved in the district, respectively equal to 26.2, 28.0 and 16.4% of the total water consumption evaluated in the three years.     

Atmospheric, oceanic and hydrological contributions to seasonal variations in length of day

 The annual and semiannual residuals derived in the axial angular momentum budget of the solid Earth–atmosphere system reflect significant signals. They must be caused by further excitation sources. Since, in particular, the contribution for the wind term from the atmospheric layer between the 10 and 0.3 hPa levels to the seasonal variations in length of day (LOD) is still missing, it is necessary to extend the top level into the upper stratosphere up to 0.3 hPa. Under the conservation of the total angular momentum of the entire Earth, variations in the oceanic angular momentum (OAM) and the hydrological angular momentum (HAM) are further significant excitation sources at seasonal time scales. Focusing on other contributions to the Earth's axial angular momentum budget, the following data are used in this study: axial atmospheric angular momentum (AAM) data derived for the 10–0.3 hPa layer from 1991 to 1997 for computing the missing wind effects; axial OAM functions as generated by oceanic general circulation models (GCMs), namely for the ECHAM3 and the MICOM models, available from 1975 to 1994 and from 1992 to 1994, respectively, for computing the oceanic contributions to LOD changes, and, concerning the HAM variations, the seasonal estimates of the hydrological contribution as derived by Chao and O'Connor [(1988) Geophys J 94: 263–270]. Using vector representation, it is shown that the vectors achieve a close balance in the global axial angular momentum budget within the estimated uncertainties of the momentum quantities on seasonal time scales. Received: 6 April 2000 / Accepted: 13 December 2000     

青藏高原近地层通量特征的合成分析

利用中日亚洲季风机制合作研究计划 ,设置在青藏高原东部地区的拉萨、日喀则、那曲和林芝 4地 1993年 7月～ 1999年 3月近 7a的自动气象站 (AWS)近地层梯度观测资料 ,确定出分季节的高原地表粗糙度和逐日的地面总体输送系数 ,以此为基础用总体输送公式对地面动量、感热和潜热通量进行了计算 ,并用合成方法分析了 1993～1999年高原近地层通量夏季、冬季的日变化和月变化特征     

洞庭湖镉迁移转化的马尔可夫模型评价

讨论了在现时泥沙、水体、鱼类、藻类的含镉状态下洞庭湖的输入输出通量；并以其为例，在简要介绍马尔可夫链数学方法的基础上，建立了湖泊生态系统中镉迁移转化的马尔可夫链数学模型，初步评价了当洞庭湖不断有镉加入、系统到达状态稳定时各种介质中的镉总量及镉的迁移转化趋势。     

基于通量均衡基准的地壳垂直形变场分析模型

根据全球通量均衡假说建立了地壳垂直形变场的分析模型,在球面近似下求得了该模型的球谐函数级数解并拟合了两个不同区域的地壳垂直形变场,编绘了等值线图,继而在密切平面坐标系下导出了该模型的Fourier级数解。结果表明,该方法不仅具有几何意义,而且具有一定的地球动力学意义。     

青藏高原净生态系统碳交换量的遥感评估模型研究

青藏高原作为世界海拔最高的区域,是全球气候变化的敏感区之一。定量估算这一区域的净生态系统碳交换量（NEE）有利于理解陆地生态系统碳平衡对未来气候变化的响应。本文构建了一个模拟该地区NEE动态变化的净碳收支模型（NCBM）。该模型由来源于MODIS影像的增强型植被指数（EVI）、陆地表面水分指数（LSWI）以及来源于地面观测的空气温度和短波辐射共同驱动,并利用青藏高原地区的3种植被类型（包括高寒灌丛、高寒湿地和高寒草甸）的碳通量长期观测数据对模型进行了校准和验证。结果表明,在模型校准站点年,NCBM模型可以模拟NEE观测值81%的变化,均方根误差（RMSE）为0.03molC/m2/d,模型效率（EF）为0.81。在模型验证站点年,NCBM模型可以预测NEE观测值84%的变化,RMSE为0.03molC/m2/d,EF为0.81。在大多数情况下,NCBM模型可以清晰地模拟各植被类型的NEE季节和年际变化。此外,NCBM模型因为结构简单,模型驱动变量易于获取等优势,具有在区域尺度上模拟NEE时空变化的潜力。但是该模型还需要进一步的改进和发展,特别需要提高对植被非常稀疏地区NEE变化的模拟能力。     

2004年3月18日高纬磁层顶多重通量管事件分析

2004-03-18 23:10～23:50 UT期间,“双星（Double Star）”探测一号卫星（TC-1）在向阳面磁层顶高纬晨侧由内向外穿越磁层顶,其时TC_1的GSM坐标为 (7.5RE, -5.5RE, -5.4RE), RE为地球半径.穿越过程中TC_1观测到了8个通量管和1个磁通量传输事件（FTEs）.在此期间Cluster星簇位于向阳面太阳风内,其GSM坐标为(18.0RE, -3.1RE, -6.2RE),其4颗卫星监测到行星际磁场（IMF）的BZ分量持续南向,BＹ有较大的负值.本文的研究表明：TC_1观测到的前7个通量管具有准周期重现性,周期大约是1～4 min,明显小于以前所观测到的FTEs的平均周期(8～11 min);所有的通量管都具有较强的核心场.本文分别使用最小方差分析法（MVA）和Grad_Shafranov反演方法（GSR）对通量管的轴向进行了分析和对比,发现所有的通量管主轴基本沿晨昏向,结果显示GSR方法在轴向分析上比MVA优越.本文使用GSR方法对通量管的磁场结构进行了分析,恢复出了通量管的磁场在卫星穿越面的结构图;此外,本文还对这次多重通量管事件进行了deHoffmann-Teller(HT)分析,结果表明,所有通量管大致朝南极方向运动,均来源于向日面低纬区域.这说明它们可能起源于向日面低纬区,由该区的磁场分量重联产生.     

风沙流中沙粒浓度分布的实验研究

利用粒子图像测速技术（PIV）,在风洞内对与天然沙接近的石英沙风沙流沙颗粒浓度沿高度分布进行了研究。对 3种不同轴线风速研究结果表明,风沙流中沙粒浓度沿水平方向基本不变,而随高度呈指数衰减,而且其衰减速率与风速密切相关。风速越大,衰减速率越慢;风速越小,衰减速率越快。风沙运动一旦发生时,靠近沙床表面的沙粒浓度随风速变化很小,可以认为趋于稳定,达到最大值。     

Observational estimation of heat budgets on drifting ice and open water over the Arctic Ocean

It is of major scientific interests to determine the parameters of momentum, heat and vapor exchange in the planetary boundary layer in order to study the effects of ocean-ice-atmosphere interactions and their feedback mechanisms on global climate[1]. Lin…     

Comparison of cell- and vertex-centered discretization methods for flow in a two-dimensional discrete-fracture–matrix system

Simulations of flow for a discrete-fracture model in fractured porous rocks have gradually become more practical, as a consequence of increased computer power and improved simulation and characterization techniques. Discrete-fracture models can be formulated in a lower-dimensional framework, where the fractures are modeled in a lower dimension than the matrix, or in an equi-dimensional form, where the fractures and the matrix have the same dimension.