A parameterisation for the threshold shear velocity to initiate deflation of dry and wet sediment

A new parameterisation for the threshold shear velocity to initiate deflation of dry and wet particles is presented. It is based on the balance of moments acting on particles at the instant of particle motion. The model hence includes a term for the aerodynamic forces, including the drag force, the lift force and the aerodynamic-moment force, and a term for the interparticle forces. The effect of gravitation is incorporated in both terms. Rather than using an implicit function for the effect of the aerodynamic forces as reported earlier in literature, a constant aerodynamic coefficient was introduced. From consideration of the van der Waals force between two particles, it was further shown that the effect of the interparticle cohesion force between two dry particles on the deflation threshold should be inversely proportional to the particle diameter squared. The interparticle force was further extended to include wet bonding forces. The latter were considered as the sum of capillary forces and adhesive forces. A model that expresses the capillary force as a function of particle diameter squared and the inverse of capillary potential was deduced from consideration of the well-known model of Fisher and the Young–Laplace equation. The adhesive force was assumed to be equal to tensile strength, and a function which is proportional to particle diameter squared and the inverse of the potential due to adhesive forces was derived. By combining the capillary-force model and the adhesive force model, the interparticle force due to wet bonding was simplified and written as a function of particle diameter squared and the inverse of matric potential. The latter was loglinearly related to the gravimetric moisture content, a relationship that is valid in the low-moisture content range that is important in the light of deflation of sediment by wind. By introducing a correction to force the relationship to converge to zero moisture content at oven dryness, the matric potential–moisture content relationship contained only one unknown model parameter, viz. moisture content at −1.5 MPa. Working out the model led to a rather simple parameterisation containing only three coefficients. Two parameters were incorporated in the term that applies to dry sediment and were determined by using experimental data as reported by Iversen and White [Sedimentology 29 (1982) 111]. The third parameter for the wet-sediment part of the model was determined from wind-tunnel experiments on prewetted sand and sandy loam aggregates. The model was validated using data from wind-tunnel experiments on the same but dry sediment, and on data obtained from simulations with the model of Chepil [Soil Sci. Soc. Am. Proc. 20 (1956) 288]. The experiments showed that soil aggregates should be treated as individual particles with a density equal to their bulk density. Furthermore, it was shown that the surface had to dry to a moisture content of about 75% of the moisture content at −1.5 MPa before deflation became sustained. The threshold shear velocities simulated with our model were found to be in good agreement with own observations and with simulations using Chepil's model.     

Leaf area index retrieval based on canopy reflectance and vegetation index in eastern China

The aim of this paper is to investigate the feasibility of using Landsat TM data to retrieve leaf area index (LAI). To get a LAI retrieval model based ground reflectance and vegetation index, detailed field data were collected in the study area of eastern China, dominated by bamboo, tea plant and greengage. Plant canopy reflectance of Landsat TM wavelength bands has been inversed using software of 6S. LAI is an important ecological parameter. In this paper, atmospheric corrected Landsat TM imagery was utilized to calculate different vegetation indices (VI), such as simple ratio vegetation index (SR), shortwave infrared modified simple ratio (MSR), and normalized difference vegetation index (NDVI). Data of 53 samples of LAI were measured by LAI-2000 (LI-COR) in the study area. LAI was modeled based on different reflectances of bands and different vegetation indices from Landsat TM and LAI samples data. There are certainly correlations between LAI and the reflectance of Tm3, TM4, TM5 and TM7. The best model through analyzing the results is LAI = 1.2097*MSR + 0.4741 using the method of regression analysis. The result shows that the correlation coefficient R2 is 0.5157, and average accuracy is 85.75%. However, whether the model of this paper is suitable for application in subtropics needs to be verified in the future.     

The tectono‐metamorphic evolution of the very low‐grade hangingwall constrains two‐stage gneiss dome formation in the Montagne Noire (Southern France)

The Montagne Noire in the southernmost French Massif Central is made of an ENE‐elongated gneiss dome flanked by Palaeozoic sedimentary rocks. The tectonic evolution of the gneiss dome has generated controversy for more than half a century. As a result, a multitude of models have been proposed that invoke various tectonic regimes and exhumation mechanisms. Most of these models are based on data from the gneiss dome itself. Here, new constraints on the dome evolution are provided based on a combination of very low‐grade petrology, K–Ar geochronology, field mapping and structural analysis of the Palaeozoic western Mont Peyroux and Faugères units, which constitute part of the southern hangingwall of the dome. It is shown that southward‐directed Variscan nappe‐thrusting (D₁) and a related medium‐P metamorphism (M₁) are only preserved in the area furthest away from the gneiss dome. The regionally dominant pervasive tectono‐metamorphic event D₂/M₂ largely transposes D₁ structures, comprises a higher metamorphic thermal gradient than M₁ (transition low‐P and medium‐P metamorphic facies series) and affected the rocks between c. 309 and 300 Ma, post‐dating D₁/M₁ by more than 20 Ma. D₂‐related fabrics are refolded by D₃, which in its turn, is followed by dextral‐normal shearing along the basal shear zone of both units at c. 297 Ma. In the western Mont Peyroux and Faugères units, D₂/M₂ is largely synchronous with shearing along the southern dome margin between c. 311 and 303 Ma, facilitating the emplacement of the gneiss dome into the upper crust. D₂/M₂ also overlaps in time with granitic magmatism and migmatization in the Zone Axiale between c. 314 and 306 Ma, and a related low‐P/high‐T metamorphism at c. 308 Ma. The shearing that accompanied the exhumation of the dome therefore was synchronous with a peak in temperature expressed by migmatization and intrusion of melts within the dome, and also with the peak of metamorphism in the hangingwall. Both, the intensity of D₂ fabrics and the M₂ metamorphic grade within the hangingwall, decrease away from the gneiss dome, with grades ranging from the anchizone–epizone boundary to the diagenetic zone. The related zonation of the pre‐D₃ metamorphic field gradients paralleled the dome. These observations indicate that D₂/M₂ is controlled by the exhumation of the Zone Axiale, and suggest a coherent kinematic between the different crustal levels at some time during D₂/M₂. Based on integration of these findings with regional geological constraints, a two‐stage exhumation of the gneiss dome is proposed: during a first stage between c. 316 and 300 Ma dome emplacement into the upper crust was controlled by dextral shear zones arranged in a pull‐apart‐like geometry. The second stage from 300 Ma onwards was characterized by northeast to northward extension, with exhumation accommodated by north‐dipping detachments and hangingwall basin formation along the northeastern dome margin.     

基于温度植被干旱指数的电网输电线路走廊干旱分析

以湖北省输电线路走廊地区作为研究区,利用2013年1~9月MODIS卫星影像数据,处理得到月尺度的归一化植被指数(Normalized Differential Vegetation Index,NDVI)与地表温度(Land Surface Temperature,LST)数据,构建NDVI-Ts特征空间,计算得到温度植被干旱指数(Temperature Vegetation Dryness Index,TVDI),用TVDI监测结果分析湖北省输电线路走廊区域2013年干旱时空分布情况。结果表明,湖北省输电线路走廊地区TVDI和土壤含水量之间存在显著的负相关,相关系数达到0.525(p0.05),由MODIS卫星影像计算得到TVDI影像可以有效表明湖北省输电线路走廊地区的土壤含水情况。     

承德初霜冻预报指标分析及其预报方法

应用常规观测资料与ncep1°×1°再分析资料,对承德市各县区初霜的时空分布特征及南北县区初霜出现时间的稳定性进行了分析,在对117个初霜环流形势分析的基础上,建立了初霜预报的天气概念模型;筛选因子,分县区确立了预报指标;采用概率、指标、数值预报产品相结合的方法及概率区间取值法建立了分县区初霜预报模型;检验表明,该模型对于中重度初霜冻预报准确率可达95%以上,对轻度初霜的预报准确率可达68%以上,无漏报出现。     

Change Analysis of Three Major Antarctic Ice Shelves Based on Multi-source Remote Sensing Data

In recent years, melting and calving happen in the Antarctic ice shelves. In this paper, seven periods of coastlines were applied to provide an analysis of the ice front changes of Ross, Filchner-Ronne and Amery ice shelves with the inclusion of 1997 and 2000 Radarsat products, 2003/2004 and 2008/2009 MODIS products and 2006, 2012, 2015 coastline which were extracted from MODIS images. Change area, SCE (Shoreline Change Envelope) and NSM (Net Shoreline Movement) were applied to analyze the variation of the ice shelf front. The results shows that, the ice front of Amery ice shelf has advanced since 1997 and the total outward extension distance of the ice front was about 20 to 25 km while the advance area reached 3.03×10³ km². Ross ice shelf and Filchner-Ronne ice shelf continued to advance after ice calving events under the driver action of glacier. However, the advance area was less than the retreat area and the net change area is respectively -9.39×10³ km² and -5.86×10³ km². The retreat distance of the collapse area were up to 53 km and 39 km in the two biggest ice shelves.     

Interactive discovery of sequential patterns in time series of wind data

Wind speed and direction vary over space and time due to the interactions between different pressures and temperature gradients within the atmospheric layers. Near the earth’s surface, these interactions are modulated by topography and artificial structures. Hence, characterizing wind behaviour over large areas and long periods is a complex but essential task for various energy-related applications. In this study, we present a novel approach to discover wind patterns by integrating sequential pattern mining and interactive visualization techniques. The approach relies on the use of the Linear time Closed pattern Miner sequence algorithm in conjunction with a time sliding window that allows the discovery of all sequential patterns present in the data. These patterns are then visualized using integrated 2D and 3D coordinated multiple views and visually explored to gain insight into the characteristics of the wind from a spatial, temporal and attribute (type of wind pattern) point of view. This proposed approach is used to analyse 10 years of hourly wind speed and direction data for 29 weather stations in the Netherlands. The results show that there are 15 main sequential patterns in the data. The spatial task shows that weather stations located in the same region do not necessarily experience similar wind pattern. For within the selected time interval, similar wind patterns can be observed in different stations and in the same station at different times of occurrence. The attribute task discovered that the repetitive occurrences of chosen pattern indicate as regular wind behaviour at different weather stations that persisted continuously over time. The results of these tasks show that the proposed interactive discovery facilitates the understanding of wind dynamics in space and time.     

中国中部典型高山气象站风速的均一性检验及其变化特点

采用二相回归方法并结合台站历史沿革信息,在对中国中部典型高山站南岳和庐山1960-2017年平均风速资料进行均一性检验和订正的基础上,分析其变化特征及其与周边低海拔台站的差异,并利用NCEP/NCAR再分析风速资料对其差异进行验证。结果表明：南岳站平均风速序列存在一个由测风仪器变更而导致的非均一点,而庐山站不存在非均一点;南岳和庐山年及四季平均风速均显著高于周边台站,且高山站以春季和夏季风速最大,而低海拔台站各季节风速差异较小;近58 a高山站及周边低海拔台站的年及四季平均风速均呈显著的减小趋势,但高山站的减小速率显著高于低海拔台站;同区域NCEP/NCAR的1000 hPa和850 hPa平均风速变化的差异与高山站和低海拔台站的差异基本一致,说明中低空和地面风速的这种差异在中国中部地区具有一定的普遍性。     

冷空气强风在大型城市中的精细结构和形成机制

地面强风可直接造成人员伤亡和经济损失,严重影响出行安全、工农业生产等社会秩序。强风的发生与天气系统和复杂下垫面的共同作用密切相关,在城市区域尤为明显。受数值模拟技术和计算资源的限制,对实际天气条件下大范围城区的强风现象进行建筑物分辨率的大规模数值模拟研究仍是一个挑战。本研究利用中尺度气象模式嵌套流体计算动力模式的超高分辨率局地气象预报系统,对强冷空气过程造成广州市区的一次强风事件进行数值模拟,深入探讨强风的精细结构和形成机制。结果表明,伴随着强冷空气入侵,广州市区的平均风速和风场高频扰动均明显增强,且在城市冠层顶尤为明显,呈现区域不均匀的三维结构,数值模拟与地面观测相一致。较大范围的强风速和阵风主要出现在建筑物较为低矮的老城区上空,并持续影响下游河道等开阔区域。在高层建筑密集的新城区,虽然整体风速明显减弱,但能在平行风向的街道狭管和下游区域形成局地强风。特别是,超高层建筑群引起显著的垂直环流,导致强风扰动向下传播,造成最大风速达8 m s?1的地面局地强风,阵风指数接近2。上游建筑群引起的风场扰动呈现大尺度湍流结构,能沿着平均气流传播影响数公里之远的下游地区。当风场扰动经过广州塔等单体超高层建筑时,可在其两侧绕流区再次加强,形成局地强风。局地强风和阵风还出现在垂直于风向排列的沿江高层建筑群的侧边,与建筑屏风的阻挡效应和缺口溢出有关。研究结果促进认识城市强风的时空特征和物理机制,有助于提升城市气象的精细化预报水平。     

延庆-张家口地区复杂地形冬季山谷风特征分析

基于2016年12月—2017年2月和2017年12月—2018年2月两年冬季的近地面自动气象站逐时观测数据以及张家口探空数据分析延庆-张家口一带（包括张家口崇礼、赤城、海坨、小五台山区,延怀、怀涿、洋河、蔚县盆地以及北京延庆、昌平、怀柔部分平原地区）复杂地形的风场精细化时、空分布特征,揭示不同复杂地形下局地风场的时、空变化规律,加深对复杂地形动力、热力作用对近地面风场影响的认识,为冬季山区风场预报以及复杂地形数值模式改进提供参考。结果表明:晴朗小风天风持续性作为矢量平均风速和标量平均风速的比值,可以作为研究风场变化规律的重要参数。根据风持续性的日变化特征,可以将研究区域内所有站点分为10种类型,分别代表不同局地地形特征的影响,风持续与风向变化的相关也很强。研究区域主要有3种类型的地形风:斜坡风、峡谷风以及较大尺度的山区平原风。不同地形特征下的风场、风持续性存在明显不同的日变化特征,山风和谷风相互转化的时间也不同,山区最早,盆地次之,平原区最晚;山风时段持续时间较谷风时段长,风速小;晴朗小风天实测风反映了实际风场的特征,而排除环境背景风场,弱化地形动力作用后整个冬季的局地风作为理论山谷风,更能反映热力作用下的山谷风特征。