Search

Forest biomass estimation is undoubtedly one of the most pressing research subjects at present. Combining multi-source remote sensing information can give full play to the advantages of different remote sensing technologies, providing more comprehensive and rich information for aboveground biomass (AGB) estimation research. Based on Landsat 8, Sentinel-2A, and ALOS2 PALSAR data, this paper takes the artificial coniferous forests in the Saihanba Forest of Hebei Province as the object of study, fully explores and establishes remote sensing factors and information related to forest structure, gives full play to the advantages of spectral signals in detecting the horizontal structure and multi-dimensional synthetic aperture radar (SAR) data in detecting the vertical structure, and combines environmental factors to carry out multivariate synergistic methods of estimating the AGB. This paper uses three variable selection methods (Pearson correlation coefficient, random forest significance, and the least absolute shrinkage and selection operator (LASSO)) to establish the variable sets, combining them with three typical non-parametric models to estimate AGB, namely, random forest (RF), support vector regression (SVR), and artificial neural network (ANN), to analyze the effect of forest structure on biomass estimation, explore the suitable AGB of artificial coniferous forests estimation of machine learning models, and develop the method of quantifying saturation value of the combined variables. The results show that the horizontal structure is more capable of explaining the AGB compared to the vertical structure information, and that combining the multi-structure information can improve the model results and the saturation value to a great extent. In this study, different sets of variables can produce relatively superior results in different models. The variable set selected using LASSO gives the best results in the SVR model, with an R 2 values of 0.9998 and 0.8792 for the training and the test set, respectively, and the highest saturation value obtained is 185.73 t/ha, which is beyond the range of the measured data. The problem of saturation in biomass estimation in boreal medium- and high-density forests was overcome to a certain extent, and the AGB of the Saihanba area was better estimated. [ABSTRACT FROM AUTHOR]

This paper considers extinction coefficient changes with height caused by the inhomogeneous distribution of scatterers in heterogeneous forests and uses the InSAR phase center height histogram and Gaussian function to fit the normalized extinction coefficient curve so as to reflect the vertical structure of the heterogeneous forest. Combining polarization decomposition based on the physical model and the PolInSAR parameter inversion method, the ground and volume coherence matrices can be separated based on the polarization characteristics and interference coherence diversity. By combining the new abovementioned parameters, the semi-empirical improved RVoG inversion model can be used to both quantify the effects of temporal decorrelation on coherence and phase errors and avoid the effects of small vertical wavenumbers on the large temporal baseline of spaceborne data. The model provided robust inversion for the height of the coniferous forest and enhanced the parameter estimation of the forest structure. This study addressed the influence of vertical structure differences on the extinction coefficient, though the coherence of the ground and volume in sparse vegetation areas could not be accurately estimated, and the oversensitivity of temporal decorrelation caused by inappropriate vertical wavenumbers. According to this method we used spaceborne L-band ALOS-2 PALSAR data on the Saihanba forest in Hebei Province acquired in 2020 for the purpose of height inversion, with a temporal baseline range of 14–70 days and the vertical wavenumber range of 0.01–0.03 rad/m. The results are further validated using sample data, with R 2 reaching 0.67. [ABSTRACT FROM AUTHOR]

In terrestrial ecosystems, vegetation is sensitive to climate change and human activities. Its spatial-temporal changes also affect the ecological and social environment. In this paper, we considered the Beijing–Tianjin–Hebei region to study the spatio-temporal vegetation patterns. The detailed analysis of a moderate-resolution imaging spectroradiometer (MODIS) data were carried out through the Google Earth Engine (GEE) platform. Our results show a slow and tortuous upward trend in the average leaf area index (LAI) in the study region for the periods 2001–2020. Specifically, Beijing had the highest LAI value, with an average of 1.64 over twenty years, followed by Hebei (1.30) and Tianjin (1.04). Among different vegetation types, forests had the highest normalized difference vegetation index (NDVI) with the range of 0.62–0.78, followed by shrubland (0.58–0.75), grassland (0.34–0.66), and cropland (0.38–0.54) over the years. Spatially, compared to the whole study area, index value in the northwestern part of the Beijing–Tianjin–Hebei region increased greatly in many areas, such as northwest Beijing, Chengde, and Zhangjiakou, indicating a significant ecological optimization. Meanwhile, there was ecological degradation in the middle and southeast regions, from Tangshan southeastward to Handan, crossing Tianjin, Langfang, the east part of Baoding, Shijiazhuang, and the west part of Cangzhou. Air temperature and precipitation were positively and significantly correlated with net primary production (NPP) and precipitation stood out as a key driver. Additionally, an intensification of the urbanization rate will negatively impact the vegetation NPP, with the shrubland and forest being affected most relative to the cropland. [ABSTRACT FROM AUTHOR]

High-resolution leaf area index (LAI) maps from remote sensing data largely depend on empirical models, which link field LAI measurements to the vegetation index. The existing empirical methods often require the field measurements to be sufficient for constructing a reliable model. However, in many regions of the world, there are limited field measurements available. This paper presents a prior knowledge-based (PKB) method to derivate LAI with limited field measurements, in an effort to improve the accuracy of empirical model. Based on the assumption that the experimental sites with the same vegetation type can be represented by similar models, a priori knowledge for crops was extracted from the published models in various cropland sites. The knowledge, composed of an initial guess of each model parameter with the associated uncertainty, was then combined with the local field measurements to determine a semi-empirical model using the Bayesian inversion method. The proposed method was evaluated at a cropland site in the Huailai region of Hebei Province, China. Compared with the regression method, the proposed PKB method can effectively improve the accuracy of empirical model and LAI estimation, when the field measurements were limited. The results demonstrate that a priori knowledge extracted from the universal sites can provide important auxiliary information to improve the representativeness of the empirical model in a given study area. [ABSTRACT FROM AUTHOR]

The aeromagnetic survey is a common remote sensing tool for detecting iron deposits. The local wavenumber of a magnetic anomaly is used to interpret the edges or positions of the sources, and can involve first- or second-order local wavenumbers. In this paper, we derived a linear equation between the second-order local wavenumber and the source location; therefore, we propose a constraint of first and second-order local wavenumbers. Tests on synthetic data show that the source parameters, computed using a combination of equations that involved different-order local wavenumbers, are closer to the true values and show a smaller spread in estimated values. For gridded data, we proved that the different-order combination allowed us to accurately estimate the source position. When applied to the aeromagnetic data from Hebei province, China, we refined the location of most magnetic features, which we interpreted as possible iron deposits. [ABSTRACT FROM AUTHOR]

Global climate change and rapid urbanization have caused increases in urban floods. Urban flood risk assessment is a vital method for preventing and controlling such disasters. This paper takes the central region of Cangzhou city in Hebei Province as an example. Detailed topographical information, such as the buildings and roads in the study area, was extracted from GF-2 data. By coupling the two models, the SWMM and MIKE21, the spatial distribution of the inundation region, and the water depth in the study area under different return periods, were simulated in detail. The results showed that, for the different return periods, the inundation region was generally consistent. However, there was a large increase in the mean inundation depth within a 10-to-30-year return period, and the increase in the maximum inundation depth and inundation area remained steady. The comprehensive runoff coefficient in all of the scenarios exceeded 0.8, indicating that the drainage system in the study area is insufficient and has a higher flood risk. The flood risk of the study area was evaluated based on the damage curve, which was obtained from field investigations. The results demonstrate that the loss per unit area was less than CNY 250/m2 in each return period in the majority of the damaged areas. Additionally, the total loss was mainly influenced by the damaged area, but, in commercial areas, the total loss was highly sensitive to the inundation depth. [ABSTRACT FROM AUTHOR]

Beijing-Tianjin-Hebei (BTH) has been suffering from severe groundwater storage (GWS) consumption and land subsidence (LS) for a long period. The overexploitation of groundwater brings about severe land subsidence, which affects the safety and development of BTH. In this paper, we utilized multi-frame synthetic aperture radar datasets obtained by the Rardarsat-2 satellite to monitor land subsidence's temporal and spatial distribution in the BTH from 2012 to 2016 based on multi-temporal interferometric synthetic aperture radar (MT-InSAR). In addition, we also employed the Gravity Recovery and Climate Experiment (GRACE) mascon datasets acquired by the Center for Space Research (CSR) and Jet Propulsion Laboratory (JPL) to obtain the GWS anomalies (GWSA) of BTH from 2003 to 2016. Then we evaluate the accuracy of the results obtained. Furthermore, we explored the relationship between the regional GWSA and the average cumulative subsidence in the BTH. The total volume change of subsidence is 59.46% of the total volume change of groundwater storage. Moreover, the long-term decreasing trend of the GWSA (14.221 mm/year) and average cumulative subsidence (17.382 mm/year) show a relatively high consistency. Finally, we analyze the heterogeneity of GWS change (GWSC) and LS change (LSC) in the four typical areas by the Lorenz curve model. The implementation of the South-to-North Water Diversion Project (MSWDP) affects the heterogeneity of GWSC and LSC. It can be seen that the largest heterogeneity of LSC lags behind the GWSC in the Tianjin-Langfang-Hengshui-Baoding area. The largest uneven subsidence in Beijing and Tianjin occurred in 2015, and the largest uneven subsidence in Hengshui-Baoding occurred in 2014. After that, the heterogeneity of subsidence gradually tends to stable. [ABSTRACT FROM AUTHOR]

ZY-1 02D is China's first civil hyperspectral (HS) operational satellite, developed independently and successfully launched in 2019. It can collect HS data with a spatial resolution of 30 m, 166 spectral bands, a spectral range of 400~2500 nm, and a swath width of 60 km. Its competitive advantages over other on-orbit or planned satellites are its high spectral resolution and large swath width. Unfortunately, the relatively low spatial resolution may limit its applications. As a result, fusing ZY-1 02D HS data with high-spatial-resolution multispectral (MS) data is required to improve spatial resolution while maintaining spectral fidelity. This paper conducted a comprehensive evaluation study on the fusion of ZY-1 02D HS data with ZY-1 02D MS data (10-m spatial resolution), based on visual interpretation and quantitative metrics. Datasets from Hebei, China, were used in this experiment, and the performances of six common data fusion methods, namely Gram-Schmidt (GS), High Pass Filter (HPF), Nearest-Neighbor Diffusion (NND), Modified Intensity-Hue-Saturation (IHS), Wavelet Transform (Wavelet), and Color Normalized Sharping (Brovey), were compared. The experimental results show that: (1) HPF and GS methods are better suited for the fusion of ZY-1 02D HS Data and MS Data, (2) IHS and Brovey methods can well improve the spatial resolution of ZY-1 02D HS data but introduce spectral distortion, and (3) Wavelet and NND results have high spectral fidelity but poor spatial detail representation. The findings of this study could serve as a good reference for the practical application of ZY-1 02D HS data fusion. [ABSTRACT FROM AUTHOR]

Accurate determination of phenological information of crops is essential for field management and decision-making. Remote sensing time-series data are widely used for extracting phenological phases. Existing methods mainly extract phenological phases directly from individual remote sensing time-series, which are easily affected by clouds, noise, and mixed pixels. This paper proposes a novel method of phenological phase extraction based on the time-weighted dynamic time warping (TWDTW) algorithm using MODIS Normalized Difference Vegetation Index (NDVI) 5-day time-series data with a spatial resolution of 500 m. Firstly, based on the phenological differences between winter wheat and other land cover types, winter wheat distribution is extracted using the TWDTW classification method, and the results show that the overall classification accuracy and Kappa coefficient reach 94.74% and 0.90, respectively. Then, we extract the pure winter-wheat pixels using a method based on the coefficient of variation, and use these pixels to generate the average phenological curve. Next, the difference between each winter-wheat phenological curve and the average winter-wheat phenological curve is quantitatively calculated using the TWDTW algorithm. Finally, the key phenological phases of winter wheat in the study area, namely, the green-up date (GUD), heading date (HD), and maturity date (MD), are determined. The results show that the phenological phase extraction using the TWDTW algorithm has high accuracy. By verification using phenological station data from the Meteorological Data Sharing Service System of China, the root mean square errors (RMSEs) of the GUD, HD, and MD are found to be 9.76, 5.72, and 6.98 days, respectively. Additionally, the method proposed in this article is shown to have a better extraction performance compared with several other methods. Furthermore, it is shown that, in Hebei Province, the GUD, HD, and MD are mainly affected by latitude and accumulated temperature. As the latitude increases from south to north, the GUD, HD, and MD are delayed, and for each 1° increment in latitude, the GUD, HD, and MD are delayed by 4.84, 5.79, and 6.61 days, respectively. The higher the accumulated temperature, the earlier the phenological phases occur. However, latitude and accumulated temperature have little effect on the length of the phenological phases. Additionally, the lengths of time between GUD and HD, HD and MD, and GUD and MD are stable at 46, 41, and 87 days, respectively. Overall, the proposed TWDTW method can accurately determine the key phenological phases of winter wheat at a regional scale using remote sensing time-series data. [ABSTRACT FROM AUTHOR]

The forest canopy height (FCH) plays a critical role in forest quality evaluation and resource management. The accurate and rapid estimation and mapping of the regional forest canopy height is crucial for understanding vegetation growth processes and the internal structure of the ecosystem. A stacking algorithm consisting of multiple linear regression (MLR), support vector machine (SVM), k-nearest neighbor (kNN), and random forest (RF) was used in this paper and demonstrated optimal performance in predicting the forest canopy height by synergizing Sentinel-2 images acquired from the cloud-based computation platform Google Earth Engine (GEE) with data from ICESat-2 (Ice, Cloud, and Land Elevation Satellite-2). This research was conducted to achieve continuous mapping of the canopy height of plantations in Saihanba Mechanical Forest Plantation, which is located in Chengde City, northern Hebei province, China. The results show that stacking achieved the best prediction accuracy for the forest canopy height, with an R2 of 0.77 and a root mean square error (RMSE) of 1.96 m. Compared with MLR, SVM, kNN, and RF, the RMSE obtained by stacking was reduced by 25.2%, 24.9%, 22.8%, and 18.7%, respectively. Since Sentinel-2 images and ICESat-2 data are publicly available, this opens the door for the accurate mapping of the continuous distribution of the forest canopy height globally in the future. [ABSTRACT FROM AUTHOR]

The data assimilation technique is an effective method for reducing initial condition errors in numerical weather prediction (NWP) models. This paper evaluated the potential of the weather research and forecasting (WRF) model and its three-dimensional data assimilation (3DVar) module in improving the accuracy of rainfall-runoff prediction through coupled atmospheric-hydrologic systems. The WRF model with the assimilation of radar reflectivity and conventional surface and upper-air observations provided the improved initial and boundary conditions for the hydrological process; subsequently, three atmospheric-hydrological systems with variable complexity were established by coupling WRF with a lumped, a grid-based Hebei model, and the WRF-Hydro modeling system. Four storm events with different spatial and temporal rainfall distribution from mountainous catchments of northern China were chosen as the study objects. The assimilation results showed a general improvement in the accuracy of rainfall accumulation, with low root mean square error and high correlation coefficients compared to the results without assimilation. The coupled atmospheric-hydrologic systems also provide more accurate flood forecasts, which depend upon the complexity of the coupled hydrological models. The grid-based Hebei system provided the most stable forecasts regardless of whether homogeneous or inhomogeneous rainfall was considered. Flood peaks before assimilation were underestimated more in the lumped Hebei model relative to the other coupling systems considered, and the model seems more applicable for homogeneous temporal and spatial events. WRF-Hydro did not exhibit desirable predictions of rapid flood process recession. This may reflect increasing infiltration due to the interaction of atmospheric and land surface hydrology at each integration, resulting in mismatched solutions for local runoff generation and confluence. [ABSTRACT FROM AUTHOR]

The long-term overexploitation of groundwater leads to serious land subsidence and threatens the safety of Beijing-Tianjin-Hebei (BTH). In this paper, an interferometric point target analysis (IPTA) with small baseline subset InSAR (SBAS-InSAR) technique was used to derive the land subsidence in a typical BTH area from 2012 to 2018 with 126 Radarsat-2 and 184 Sentinel-1 images. The analysis reveals that the average subsidence rate reached 118 mm/year from 2012 to 2018. Eleven subsidence features were identified: Shangzhuang, Beijing Airport, Jinzhan and Heizhuanghu in Beijing, Guangyang and Shengfang in Langfang, Wangqingtuo in Tianjin, Dongguang in Cangzhou, Jingxian and Zaoqiang in Hengshui and Julu in Xingtai. Comparing the different types of land use in subsidence feature areas, the results show that when the land-use type is relatively more complex and superimposed with residential, industrial and agricultural land, the land subsidence is relatively more significant. Moreover, the land subsidence development patterns are different in the BTH areas because of the different methods adopted for their water resource development and utilization, with an imbalance in their economic development levels. Finally, we found that the subsidence changes are consistent with groundwater level changes and there is a lag period between land subsidence and groundwater level changes of approximately two months in Beijing Airport, Jinzhan, Jingxian and Zaoqiang, of three months in Shangzhuang, Heizhuanghu, Guangyang, Wangqingtuo and Dongguang and of four months in Shengfang. [ABSTRACT FROM AUTHOR]

The southern part of the Hebei Province is one of China's major crop-producing regions. Due to the continuous decline in groundwater level, agricultural water use is facing significant challenges. Precision agricultural irrigation management is undoubtedly an effective way to solve this problem. Based on multisource data (time series soil moisture active passive (SMAP) data, Moderate Resolution Imaging Spectroradiometer (MODIS) normalized difference vegetation index (NDVI) and evapotranspiration (ET), and meteorological station precipitation), the irrigation signal (frequency, timing and area) is detected in the southern part of the Hebei Province. The SMAP data was processed by the 5-point moving average method to reduce the error caused by the uncertainty of the microwave data derived SM. Irrigation signals can be detected by removing the precipitation effect and setting the SM change threshold. Based on the validation results, the overall accuracy of the irrigation signal detection is 77.08%. Simultaneously, considering the spatial resolution limitation of SMAP pixels, the SMAP irrigation area was downscaled using the winter wheat area extracted from MODIS NDVI. The analytical results of 55 winter wheat samples (5 samples in a group) showed that winter wheat covered by one SMAP pixel had an 82.72% growth consistency in surface water irrigation period, which can indicate a downscaling effectiveness. According to the above statistical analysis, this paper considers that although the spatial resolution of SMAP data is insufficient, it can reflect the change of SM more sensitively. In areas where the crop pattern is relatively uniform, the introduction of high-resolution crop pattern distribution can be used not only to detect irrigation signals but also to validate the effectiveness of irrigation signal detection by analyzing crop growth consistency. Therefore, the downscaling results can indicate the true winter wheat irrigation timing, area and frequency in the study area. [ABSTRACT FROM AUTHOR]