陕北黄土高原植被覆盖度动态变化分析研究

以Landsat TM归一化植被指数(NDVI)为数据源,运用像元二分模型提取陕北黄土高原1990、2000、2010年夏季的植被覆盖度,分析陕北黄土高原植被覆盖度的空间变化情况。结果显示:研究区植被覆盖度呈波动上升趋势。不同等级植被覆盖度在数量和空间位置上的转移较为活跃。大于等于60%的植被覆盖度和小于等于40%的植被覆盖度在空间上呈西南—东北两个方向扩张的分布趋势。受气候和人为等因素的影响,陕北黄土高原植被改善良好。     

河南省植被覆盖度时空变化及其原因分析

基于河南省2010—2019年MODIS NDVI数据,利用像元二分模型估算植被覆盖度,利用一元线性回归获取变化趋势,并利用Pearson相关系数法探究河南省的植被指数和植被覆盖度的时空变化原因,可为河南省生态现状、生态环境建设规划和布局提供参考。结果表明：(1)2010—2019年河南省年度最大NDVI均值存在波动,呈微弱下降趋势;(2)河南省较高和中等植被覆盖度占比较大,且存在中等植被覆盖度向高植被覆盖度和低或无植被覆盖度的转换,社会因素是中等分化到低或无植被覆盖度的主要因素;(3)不同类型植被对降水变化的敏感性存在差异,森林、降水丰富地区和以灌溉为主的农业区受降水变化影响较弱,城镇、乡村和以降水为主的农业区受降水影响较强。     

中国北方地区植被覆盖度遥感估算及其变化分析

为了分析中国北方地区2000年之后植被覆盖度的时空分布及其变化,利用MODIS光谱反射率数据计算归一化植被指数,采用像元二分模型对中国北方地区2000—2012年植被覆盖度进行定量估算,分析研究区13 a间植被覆盖度的时空变化特征。研究结果表明:植被覆盖度年内变化特征体现在最大植被覆盖度一般出现在7和8月份,与中国北方地区植被的生长季相一致;整个中国北方地区年最大植被覆盖度呈现缓慢增长的趋势,其增长速率为每年0.2%;年最大植被覆盖度变化的空间分布具有较大差异,其中东北、华北和黄土高原等三北防护林工程建设区的年最大植被覆盖度有较明显的增长。     

喀什市城市植被及其覆盖度研究

针对国内外对城市植被及其覆盖度方面的研究不够深入的现状,该文以喀什市1990年、2000年及2010年3期Landsat TM/ETM+遥感影像为数据源,结合归一化植被指数(NDVI)与植被覆盖度估算方法,对喀什市植被覆盖区土地类型、面积变化等方面进行研究。结果表明:1990年至2010年,高覆盖区所占比例由1990年到2000年逐步增加,2010年稍有减少;中高覆盖区面积出现先减小后增加的变化趋势;中覆盖区面积逐年增加;低覆盖区呈先减小后增加的趋势;极低覆盖区变化最慢,2010年比1990年稍有增加。     

阜新地区植被覆盖度变化提取及分析

植被覆盖度是反应地区生态环境的重要指标,利用1995,2007年的两期TM遥感数据,以归一化植被指数(NDVI)像元二分法为植被覆盖度估算模型,计算阜新地区不同时期的植被覆盖度并得出阜新地区植被覆盖度等级图以及阜新地区植被覆盖度变化等级图。得出如下结论:1995年到2007年阜新地区植被覆盖度退化面积为64.817%,好转面积为6.547%,基本无变化区域为28.636%,阜新地区植被覆盖度退化严重。     

鄱阳湖生态经济区成立前后植被覆盖变化

针对鄱阳湖生态经济区成立后,国家加大区域生态保护力度,区域植被覆盖度发生变化的问题,该文采用像元二分法估算植被覆盖度,并分析其在生态经济区成立前后变化特征。结果表明:鄱阳湖生态经济区植被覆盖情况较好,中高植被覆盖度和高植被覆盖度区域面积所占比例较高。生态经济区成立前后,植被覆盖变化呈现差异,成立之前(1995—2010年)植被生长发生负向转移,退化现象较为严重,发生退化的区域主要分布在鄱阳湖周边和南部;成立之后(2010—2015年)植被生长发生正向转移,在省会和城市周围呈现退化趋势,其他区域呈现增长趋势。农业耕作效率和林地资源保护是生态经济区成立前后植被覆盖变化差异的主要原因。     

南方丘陵区植被覆盖度遥感估算的地形效应评估

植被覆盖变化是生态环境领域的核心研究内容之一,但其估算精度常受到地形效应、土壤背景、大气效应等各种因素影响。以Landsat 8 OLI为遥感数据源,基于像元二分模型,分别利用归一化差值植被指数(NDVI)、经Cosine-C校正的归一化差值植被指数(NDVI)和归一化差值山地植被指数(NDMVI)建立植被覆盖度估算模型,以评估南方丘陵区植被覆盖度的地形效应。结果表明,3种植被覆盖度估算模型均能削弱地形效应,但消除或抑制地形效应影响的能力不同。比较而言,基于NDMVI指数构建的植被覆盖度估算模型的地形效应最小,更适合地形复杂区域的植被覆盖度遥感估算;基于Cosine-C校正的NDVI植被指数构建的植被覆盖度估算模型的地形效应次之,但存在一定的过度校正现象;基于NDVI植被指数构建的植被覆盖度估算模型的地形效应最大,尤其当坡度≥10°时,阴坡植被覆盖度比阳坡明显偏低。     

基于MOD13Q1数据的大湄公河次区域植被覆盖时空变化分析

基于MOD13Q1-NDVI数据,采用最大值合成法提取2000—2017年的月植被指数,分别从月平均NDVI和年平均NDVI两个角度分析了大湄公河次区域植被覆盖的时空变化特征.结果表明:2000年以来,大湄公河次区域年平均NDVI总体上呈波动增长;但不同时段植被覆盖度变化不同,2000—2005年植被覆盖度总体略有降低...     

海南岛植被覆盖变化驱动因子及相对作用评价

为了解海南岛植被覆盖变化影响因子,利用MODIS归一化植被指数(normalized difference vegetation index,NDVI)数据,采用残差分析法、趋势分析法和相对作用分析法对海南岛2001—2015年间植被覆盖变化特征进行研究,定量分析了气候变化和人类活动对海南岛植被覆盖变化的贡献。研究表明,2001—2015年间,海南岛植被总体呈现增加趋势,其增长速率为0.024/10 a,植被覆盖增加区域占全岛的77.77%,植被退化区域所占比例仅为22. 23%。在植被改善区域,气候变化和人类活动的相对作用分别为31.04%和68.96%,其中相对作用大于50%的区域(即以人类活动为主)占整个植被改善区域的80.79%;植被退化区域,气候变化和人类活动的相对作用分别为35.03%和64.97%,其中相对作用大于50%的区域(即以人类活动为主)占整个植被退化区域的75. 59%。整体来看,气候变化对海南岛植被生长有促进作用,人类活动对海南岛地区植被的建设作用大于破坏作用。     

Landsat时序数据支持下的藏东南植被时空变化监测

本文利用GEE平台和1990—2019年巴宜区Landsat遥感影像,采用像元二分模型、相关性分析等方法分析了巴宜区植被覆盖度的时空变化特征与驱动力。研究结果表明：①1990—2019年巴宜区植被覆盖度总体呈稳中有增的趋势,其中,河谷区域增加明显,而高海拔区域相对稳定;②1990—2019年巴宜区气温呈显著升高,降水略有下降,总体呈“暖干化”,气温较降水量对植被覆盖变化更明显,但气候变化对植被覆盖变化影响总体不明显;③1990—2019年巴宜区植被覆盖变化与人类活动有很好的相关性,其中,低、中低、中、中高植被覆盖区域,呈显著的负相关,而高植被覆盖区域呈正相关。本文基于遥感大数据和地理云计算的植被覆盖监测动态监测和定量分析方法,能对高山峡谷区生态评估和演替分析提供一定的技术支撑和科学数据。     

Estimate of heavy metals in soil and streams using combined geochemistry and field spectroscopy in Wan-sheng mining area,Chongqing, China

Heavy metals contaminated soils and water will become a major environmental issue in the mining areas. This paper intends to use field hyper-spectra to estimate the heavy metals in the soil and water in Wan-sheng mining area in Chongqing. With analyzing the spectra of soil and water, the spectral features deriving from the spectral of the soils and water can be found to build the models between these features and the contents of Al, Cu and Cr in the soil and water by using the Stepwise Multiple Linear Regression (SMLR). The spectral features of Al are: 480 nm, 500 nm, 565 nm, 610 nm, 680 nm, 750 nm, 1000 nm, 1430 nm, 1755 nm, 1887 nm, 1920 nm, 1950 nm, 2210 nm, 2260 nm; The spectral features of Cu are: 480 nm, 500 nm, 610 nm, 750 nm, 860 nm, 1300 nm, 1430 nm, 1920 nm, 2150 nm, 2260 nm; And the spectral features of Cr are: 480 nm, 500 nm, 610 nm, 715 nm, 750 nm, 860 nm, 1300 nm, 1430 nm, 1755 nm, 1920 nm, 1950 nm. With these features, the best models to estimate the heavy metals in the study area were built according to the maximal R². The R² of the models of estimating Al, Cu and Cr in the soil and water are 0.813, 0.638, 0.604 and 0.742, 0.584, 0.513 respectively. And the gradient maps of these three types of heavy metals’ concentrations can be created by using the Inverse distance weighted (IDW).The gradient maps indicate that the heavy metals in the soil have similar patterns, but in the North-west of the streams in the study area, the contents are of great differences. These results show that it is feasible to predict contaminated heavy metals in the soils and streams due to mining activities by using the rapid and cost-effective field spectroscopy.     

Determination of the factors governing soil erodibility using hyperspectral visible and near-infrared reflectance spectroscopy

Soil erodibility, which is difficult to estimate and upscaling, was determined in this study using multiple spectral models of soil properties (soil organic matter (SOM), water-stable aggregates (WSA) > 0.25 mm, the geometric mean radius (D_g)). Herein, the soil erodibility indicators were calculated, and soil properties were quantitatively analyzed based on laboratory simulation experiments involving two selected contrasting soils. In addition, continuous wavelet transformation was applied to the reflectance spectra (350–2500 nm) of 65 soil samples from the study area. To build the relationship, the soil properties that control erodibility were identified prior to the spectral analysis. In this study, the SOM, D_g and WSA >0.25 mm were selected to represent the most significant soil properties controlling erodibility and describe the erodibility indicator based on a logarithmic regression model as a function of SOM or WSA > 0.25 mm. Five, six and three wavelet features were observed to calibrate the estimated soil properties model, and the best performance was obtained with a combination feature regression model for SOM (R² = 0.86, p < 0.01), D_g (R² = 0.79, p < 0.01) and WSA >0.25 mm (R² = 0.61, p < 0.01), respectively. One part of the wavelet features captured amplitude variations in the broad shape of the reflectance spectra, and another part captured variations in the shape and depth of the soil dry substances. The wavelet features for the validated dataset used to predict the SOM, WSA >0.25 mm and D_g were not significantly different compared with the calibrated dataset. The synthesized spectral models of soil properties, and the formation of a new equation for soil erodibility transformed from the spectral models of soil properties are presented in this study. These results show that a spectral analytical approach can be applied to complex datasets and provide new insights into emerging dynamic variation with erodibility estimation.     

Determination of the crop row orientations from Formosat-2 multi-temporal and panchromatic images

This paper presents a technique developed for the retrieval of the orientation of crop rows, over anthropic lands dedicated to agriculture in order to further improve estimate of crop production and soil erosion management. Five crop types are considered: wheat, barley, rapeseed, sunflower, corn and hemp. The study is part of the multi-sensor crop-monitoring experiment, conducted in 2010 throughout the agricultural season (MCM’10) over an area located in southwestern France, near Toulouse. The proposed methodology is based on the use of satellite images acquired by Formosat-2, at high spatial resolution in panchromatic and multispectral modes (with spatial resolution of 2 and 8 m, respectively). Orientations are derived and evaluated for each image and for each plot, using directional spatial filters (45° and 135°) and mathematical morphology algorithms. “Single-date” and “multi-temporal” approaches are considered. The single-date analyses confirm the good performances of the proposed method, but emphasize the limitation of the approach for estimating the crop row orientation over the whole landscape with only one date. The multi-date analyses allow (1) determining the most suitable agricultural period for the detection of the row orientations, and (2) extending the estimation to the entire footprint of the study area. For the winter crops (wheat, barley and rapeseed), best results are obtained with images acquired just after harvest, when surfaces are covered by stubbles or during the period of deep tillage (0.27 > R² > 0.99 and 7.15° > RMSE > 43.02°). For the summer crops (sunflower, corn and hemp), results are strongly crop and date dependents (0 > R² > 0.96, 10.22° > RMSE > 80°), with a well-marked impact of flowering, irrigation equipment and/or maximum crop development. Last, the extent of the method to the whole studied zone allows mapping 90% of the crop row orientations (more than 45,000 ha) with an error inferior to 40°, associated to a confidence index ranging from 1 to 5 for each agricultural plot.     

Exploring hyperspectral bands and estimation indices for leaf nitrogen accumulation in wheat

Hyperspectral sensing can provide an effective means for fast and non-destructive estimation of leaf nitrogen (N) status in crop plants. The objectives of this study were to design a new method to extract hyperspectral spectrum information, to explore sensitive spectral bands, suitable bandwidth and best vegetation indices based on precise analysis of ground-based hyperspectral information, and to develop regression models for estimating leaf N accumulation per unit soil area (LNA, g N m⁻²) in winter wheat (Triticum aestivum L.). Three field experiments were conducted with different N rates and cultivar types in three consecutive growing seasons, and time-course measurements were taken on canopy hyperspectral reflectance and LNA under the various treatments. Then, normalized difference spectral indices (NDSI) and ratio spectral indices (RSI) based on the original spectrum and the first derivative spectrum were constructed within the range of 350–2500 nm, and their relationships with LNA were quantified. The results showed that both LNA and canopy hyperspectral reflectance in wheat changed with varied N rates, with consistent patterns across different cultivars and seasons. The sensitive spectral bands for LNA existed mainly within visible and near infrared regions. The best spectral indices for estimating LNA in wheat were found to be NDSI (R₈₆₀, R₇₂₀), RSI (R₉₉₀, R₇₂₀), NDSI (FD₇₃₆, FD₅₂₆) and RSI (FD₇₂₅, FD₅₁₆), and the regression models based on the above four spectral indices were formulated as Y = 26.34x^1.887, Y = 5.095x − 6.040, Y = 0.609 e^3.008x and Y = 0.388x^1.260, respectively, with R² greater than 0.81. Furthermore, expanding the bandwidth of NDSI (R₈₆₀, R₇₂₀) and RSI (R₉₉₀, R₇₂₀) from 1 nm to 100 nm at 1 nm interval produced the LNA monitoring models with similar performance within about 33 nm and 23 nm bandwidth, respectively, over which the statistical parameters of the models became less stable. From testing of the derived equations, the model for LNA estimation on NDSI (R₈₆₀, R₇₂₀), RSI (R₉₉₀, R₇₂₀), NDSI (FD₇₃₆, FD₅₂₆) and RSI (FD₇₂₅, FD₅₁₆) gave R² over 0.79 with more satisfactory performance than previously reported models and physical models in wheat. It can be concluded that the present hyperspectral parameters of NDSI (R₈₆₀, R₇₂₀), RSI (R₉₉₀, R₇₂₀), NDSI (FD₇₃₆, FD₅₂₆) and RSI (FD₇₂₅, FD₅₁₆) can be reliably used for estimating LNA in winter wheat.     

Nonlinear regional economic dynamics: continuous-time specification,estimation and stability analysis

This paper presents an innovative approach to the study of regional economic dynamics within a nonlinear continuous-time econometric framework—a generalized specification of the Lotka–Volterra system of equations. This specification, which accounts for interdependent behavior of three industrial sectors and spillover effects of activities in neighboring regions, is employed in an analysis of five Italian regions between 1980 and 2003. For these regions, we report estimation results, characterize the varying systems dynamics, analyze the models’ local and global stability properties, and determine via sensitivity analyses which structural features appear to exert the greatest influence on these properties.

Study on Recovering the Earth＇s Potential Field Based on GOCE Gradiometry

Given the second radial derivative Vrr（P） ｜δs of the Earth＇s gravitational potential V（P） on the surface δS corresponding to the satellite altitude, by using the fictitious compress recovery method, a fictitious regular harmonic field rrVrr（P）^＊ and a fictitious second radial gradient field V：（P） in the domain outside an inner sphere Ki can be determined, which coincides with the real field V（P） in the domain outside the Earth. Vrr^＊（P）could be further expressed as a uniformly convergent expansion series in the domain outside the inner sphere, because rrV（P）^＊ could be expressed as a uniformly convergent spherical harmonic expansion series due to its regularity and harmony in that domain. In another aspect, the fictitious field V^＊（P） defined in the domain outside the inner sphere, which coincides with the real field V（P） in the domain outside the Earth, could be also expressed as a spherical harmonic expansion series. Then, the harmonic coefficients contained in the series expressing V^＊（P） can be determined, and consequently the real field V（P） is recovered. Preliminary simulation calculations show that the second radial gradient field Vrr（P） could be recovered based only on the second radial derivative V（P）｜δs given on the satellite boundary. Concerning the final recovery of the potential field V（P） based only on the boundary value Vrr （P）｜δs, the simulation tests are still in process.     

An assessment of chlorophyll-a concentration spatio-temporal variation using Landsat satellite data,in a small tropical reservoir

Traditional approaches to monitoring aquatic systems are often limited by the need for data collection which often is time-consuming, expensive and non-continuous. The aim of the study was to map the spatio-temporal chlorophyll-a concentration changes in Malilangwe Reservoir, Zimbabwe as an indicator of phytoplankton biomass and trophic state when the reservoir was full (year 2000) and at its lowest capacity (year 2011), using readily available Landsat multispectral images. Medium-spatial resolution (30 m) Landsat multispectral Thematic Mapper TM 5 and ETM+ images for May to December 1999–2000 and 2010–2011 were used to derive chlorophyll-a concentrations. In situ measured chlorophyll-a and total suspended solids (TSS) concentrations for 2011 were employed to validate the Landsat chlorophyll-a and TSS estimates. The study results indicate that Landsat-derived chlorophyll-a and TSS estimates were comparable with field measurements. There was a considerable wet vs. dry season differences in total chlorophyll-a concentration, Secchi disc depth, TSS and turbidity within the reservoir. Using Permutational multivariate analyses of variance (PERMANOVA) analysis, there were significant differences (p < 0.0001) for chlorophyll-a concentration among sites, months and years whereas TSS was significant during the study months (p < 0.05). A strong positive significant correlation among both predicted TSS vs. chlorophyll-a and measured vs. predicted chlorophyll-a and TSS concentrations as well as an inverse relationship between reservoir chlorophyll-a concentrations and water level were found (p < 0.001 in all cases). In conclusion, total chlorophyll-a concentration in Malilangwe Reservoir was successfully derived from Landsat remote sensing data suggesting that the Landsat sensor is suitable for real-time monitoring over relatively short timescales and for small reservoirs. Satellite data can allow for surveying of chlorophyll-a concentration in aquatic ecosystems, thus, providing invaluable data in data scarce (limited on site ground measurements) environments.     

On the multivariate total least-squares approach to empirical coordinate transformations. Three algorithms

The multivariate total least-squares (MTLS) approach aims at estimating a matrix of parameters, Ξ, from a linear model (Y−E _Y = (X−E _X ) · Ξ) that includes an observation matrix, Y, another observation matrix, X, and matrices of randomly distributed errors, E _Y and E _X . Two special cases of the MTLS approach include the standard multivariate least-squares approach where only the observation matrix, Y, is perturbed by random errors and, on the other hand, the data least-squares approach where only the coefficient matrix X is affected by random errors. In a previous contribution, the authors derived an iterative algorithm to solve the MTLS problem by using the nonlinear Euler–Lagrange conditions. In this contribution, new lemmas are developed to analyze the iterative algorithm, modify it, and compare it with a new ‘closed form’ solution that is based on the singular-value decomposition. For an application, the total least-squares approach is used to estimate the affine transformation parameters that convert cadastral data from the old to the new Israeli datum. Technical aspects of this approach, such as scaling the data and fixing the columns in the coefficient matrix are investigated. This case study illuminates the issue of “symmetry” in the treatment of two sets of coordinates for identical point fields, a topic that had already been emphasized by Teunissen (1989, Festschrift to Torben Krarup, Geodetic Institute Bull no. 58, Copenhagen, Denmark, pp 335–342). The differences between the standard least-squares and the TLS approach are analyzed in terms of the estimated variance component and a first-order approximation of the dispersion matrix of the estimated parameters.     

Validation of the SEBS-derived sensible heat for FY3A/VIRR and TERRA/MODIS over an alpine grass region using LAS measurements

In this study, sensible heat (H) calculation using remote sensing data over an alpine grass landscape is conducted from May to September 2010, and the calculation is validated using LAS (large aperture scintillometers) measurements. Data from two remote sensing sensors (FY3A-VIRR and TERRA-MODIS) are analysed. Remote sensing data, combined with the ground meteorological observations (pressure, temperature, wind speed, humidity) are fed into the SEBS (Surface Energy Balance System) model. Then the VIRR-derived sensible heat (VIRR_SEBS_H) and MODIS-derived sensible heat (MODIS_SEBS_H) are compared with the LAS-estimated H, which are obtained at the respective satellite overpass time. Furthermore, the similarities and differences between the VIRR_SEBS_H and MODIS_SEBS_H values are investigated. The results indicate that VIRR data quality is as good as MODIS data for the purpose of H estimation. The root mean square errors (rmse) of the VIRR_SEBS_H and MODIS_SEBS_H values are 45.1098 W/m² (n = 64) and 58.4654 W/m² (n = 71), respectively. The monthly means of the MODIS_SEBS_H are marginally higher than those of VIRR_SEBS_H because the satellite overpass time of the TERRA satellite lags by 25 min to that of the FT3A satellite. Relative evaporation (EFr), which is more time-independent, shows a higher agreement between MODIS and VIRR. Many common features are shared by the VIRR_SEBS_H and the MODIS_SEBS_H, which can be attributed to the SEBS model performance. In May–June, H is over-estimated with more fluctuations and larger rmse, whereas in July–September, H is under-estimated with fewer fluctuations and smaller rmse. Sensitivity analysis shows that potential temperature gradient (delta_T) plays a dominant role in determining the magnitude and fluctuation of H. The largest rmse and over-estimation in H occur in June, which could most likely be attributed to high delta_T, high wind speed, and the complicated thermodynamic state during the transitional period when bare land transforms to dense vegetation cover.     

Plant phenolics and absorption features in vegetation reflectance spectra near 1.66 μm

Past laboratory and field studies have quantified phenolic substances in vegetative matter from reflectance measurements for understanding plant response to herbivores and insect predation. Past remote sensing studies on phenolics have evaluated crop quality and vegetation patterns caused by bedrock geology and associated variations in soil geochemistry. We examined spectra of pure phenolic compounds, common plant biochemical constituents, dry leaves, fresh leaves, and plant canopies for direct evidence of absorption features attributable to plant phenolics. Using spectral feature analysis with continuum removal, we observed that a narrow feature at 1.66 μm is persistent in spectra of manzanita, sumac, red maple, sugar maple, tea, and other species. This feature was consistent with absorption caused by aromatic CH bonds in the chemical structure of phenolic compounds and non-hydroxylated aromatics. Because of overlapping absorption by water, the feature was weaker in fresh leaf and canopy spectra compared to dry leaf measurements. Simple linear regressions of feature depth and feature area with polyphenol concentration in tea resulted in high correlations and low errors (% phenol by dry weight) at the dry leaf (r² = 0.95, RMSE = 1.0%, n = 56), fresh leaf (r² = 0.79, RMSE = 2.1%, n = 56), and canopy (r² = 0.78, RMSE = 1.0%, n = 13) levels of measurement. Spectra of leaves, needles, and canopies of big sagebrush and evergreens exhibited a weak absorption feature centered near 1.63 μm, short ward of the phenolic compounds, possibly consistent with terpenes. This study demonstrates that subtle variation in vegetation spectra in the shortwave infrared can directly indicate biochemical constituents and be used to quantify them. Phenolics are of lesser abundance compared to the major plant constituents but, nonetheless, have important plant functions and ecological significance. Additional research is needed to advance our understanding of the spectral influences of plant phenolics and terpenes relative to dominant leaf biochemistry (water, chlorophyll, protein/nitrogen, cellulose, and lignin).