Search

Accurate and complete digital elevation models (DEMs) play an important fundamental role in geospatial analysis, supporting various engineering applications, human activities, and scientific research. Interferometric synthetic aperture radar (InSAR) plays an increasingly important role in DEM generation. Nonetheless, owing to its inherent characteristics, gaps often appear in regions marked by significant topographical fluctuations, necessitating an extra void-filling process. Traditional void-filling methods have operated directly on preexisting data, succeeding in relatively flat terrain. When facing mountainous regions, there will always be gross errors in elevation values. Regrettably, conventional methods have often disregarded this vital consideration. To this end, this research proposes a DEM void-filling method based on incorporating elevation outlier detection. It accounts for the detection and removal of elevation outliers, thereby mitigating the shortcomings of existing methods and ensuring robust DEM restoration in mountainous terrains. Experiments were conducted to validate the method applicability using TanDEM-X data from Sichuan, China, Hebei, China, and Oregon, America. The results underscore the superiority of the proposed method. Three traditional methods are selected for comparison. The proposed method has different degrees of improvement in filling accuracy, depending on the void status of the local terrain. Compared with the delta surface fill (DSF) method, the root mean squared error (RMSE) of the filling results has improved by 7.87% to 51.87%. The qualitative and quantitative experiments demonstrate that the proposed method is promising for large-scale DEM void-filling tasks. [ABSTRACT FROM AUTHOR]

Satellite-derived aerosol optical depth (AOD) has been extensively utilized for retrieving ground-level PM2.5 distributions. However, the presence of non-random missing data gaps in AOD poses a challenge to directly obtaining the gap-free AOD-derived PM2.5, thereby impeding accurate exposure risk assessment. Here, this study presents a novel and flexible framework that couples stacking and flexible spatiotemporal data fusion (FSDAF) approaches. By integrating multiple models and data sources, this framework aims to generate hourly (24-h) gap-free PM2.5 estimates for the Beijing–Tianjin–Hebei (BTH) region in 2018. This study effectively reconstructed data at least three times more effectively than the original AOD-derived PM2.5, achieving the Pearson coefficient (r), the coefficient determination (R2), root mean squared error (RMSE), and mean absolute error (MAE) values of 0.91, 0.84, 19.38 µg/m3, and 12.17 µg/m3, respectively, based on entire samples. Such strong predictive performance was also exhibited in spatial-based (r: 0.92–0.93, R2: 0.85–0.87, RMSE: 18.13 µg/m3–20.18 µg/m3, and MAE: 11.21 µg/m3–12.52 µg/m3) and temporal-based (r: 0.91–0.98, R2: 0.82–0.96, RMSE: 3.8 µg/m3–21.89 µg/m3, and MAE: 2.71 µg/m3–14.00 µg/m3) validations, indicating the robustness of this framework. Additionally, this framework enables the assessment of annual and seasonal PM2.5 concentrations and distributions, revealing that higher levels are experienced in the southern region, while lower levels prevail in the northern part. Winter exhibits the most severe levels, followed by spring and autumn, with comparatively lower levels in summer. Notably, the proposed framework effectively mitigates bias in calculating population-weighted exposure risk by filling data gaps with calculated values of 51.04 µg/m3, 54.17 µg/m3, 56.24 µg/m3, and 55.00 µg/m3 in Beijing, Tianjin, Hebei, and the BTH region, respectively. [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]

Gravity Recovery and Climate Experiment (GRACE)-derived groundwater storage anomalies (GWSA) have been used to highlight groundwater depletion in regional aquifer systems worldwide. However, the use of GRACE products in smaller areas is limited owing to the coarse spatial resolution of the data product. This study utilized a dynamic downscaling method to improve the GWSA resolution from 1° to 0.05° by constructing a groundwater storage numerical model in the Beijing, Tianjin, and Hebei regions of China (BTH). The results indicate that: (1) the GRACE-derived and calculated GWSA had a good match with an average root mean squared error (RMSE) of 2.61 cm equivalent water height (EWH) and an average Nash–Sutcliffe efficiency coefficient (NSE) of 0.84 for the calibration period. (2) The hydraulic gradient coefficient and specific yield had the highest sensitivity, and transmissivity had the lowest sensitivity; however, different forcing data had no obvious influence on the GWSA. (3) The downscaled results not only exhibited time series variations that were consistent with those of the GRACE-derived solutions but also revealed a finer spatial heterogeneity of the GWSA along with increasing correlation coefficients between the GRACE-derived GWSA and the in situ measurements of groundwater levels by 0.06 and reducing the RMSE by 8.85%. (4) The downscaled results reflected the spatiotemporal change characteristics of groundwater storage in different hydrogeological units and administrative regions well. This study demonstrates the potential applications of the proposed downscaling method for both regional and local groundwater resource management. [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]

Linear regression models are commonly used for estimating ground PM2.5 concentrations, but the global spatial autocorrelation and local spatial heterogeneity of PM2.5 distribution are either ignored or only partially considered in commonly used models for estimating PM2.5 concentrations. Therefore, taking both global spatial autocorrelation and local spatial heterogeneity into consideration, a global-local regression (GLR) model is proposed for estimating ground PM2.5 concentrations in the Yangtze River Delta (YRD) and in the Beijing, Tianjin, Hebei (BTH) regions of China based on the aerosol optical depth data, meteorological data, remote sensing data, and pollution source data. Considering the global spatial autocorrelation, the GLR model extracts global factors by the eigenvector spatial filtering (ESF) method, and combines the fraction of them that passes further filtering with the geographically weighted regression (GWR) method to address the local spatial heterogeneity. Comprehensive results show that the GLR model outperforms the ordinary GWR and ESF models, and the GLR model has the best performance at the monthly, seasonal, and annual levels. The average adjusted R2 of the monthly GLR model in the YRD region (the BTH region) is 0.620 (0.853), which is 8.0% and 7.4% (6.8% and 7.0%) higher than that of the monthly ESF and GWR models, respectively. The average cross-validation root mean square error of the monthly GLR model is 7.024 μg/m3 in the YRD region, and 9.499 μg/m3 in the BTH region, which is lower than that of the ESF and GWR models. The GLR model can effectively address the spatial autocorrelation and spatial heterogeneity, and overcome the shortcoming of the ordinary GWR model that overfocuses on local features and the disadvantage of the poor local performance of the ordinary ESF model. Overall, the GLR model with good spatial and temporal applicability is a promising method for estimating PM2.5 concentrations. [ABSTRACT FROM AUTHOR]

Forest height is an essential parameter used to derive important information about forest ecosystems, such as forest above-ground biomass. In this article, a forest height estimation approach combining P-band and X-band interferometric synthetic aperture radar (InSAR) was introduced. The forest height was estimated using the difference in the penetration of long- and short-wavelength radars to the forest. That is, the P-band and X-band InSAR data were used to extract the digital terrain model (DTM) and digital surface model (DSM), respectively. For the DTM, an improved time-frequency (TF) analysis method was used to reduce the effect of forest scatterers on the extraction of a pure understory terrain phase based on P-band InSAR. For the DSM, a novel compensation algorithm based on a multi-layer model (MLM) was proposed to remove the penetration bias of the X-band. Compared to the existing method based on the infinitely deep uniform volumes (IDUV) model, the MLM-based method is more in line with the characteristics of forest structure and the scattering mechanism for X-band InSAR. The airborne P-band repeat-pass InSAR and spaceborne X-band (TanDEM-X) single-pass InSAR data were used to verify the proposed method over the study area in the Saihanba Forest Farm in Hebei, China. The results demonstrated that the improved TF method can achieve high-precision DTM extraction based on P-band InSAR data, and the root mean square error (RMSE) was 0.94 m. The proposed MLM-based compensation method of the DSM achieved a smaller error (RMSE: 1.67 m) compared to the IDUV-based method (RMSE: 3.01 m). Under the same DTM extracted by P-band InSAR, the estimation accuracy of forest height based on the MLM method was 86.58% (RMSE: 1.81 m), which was 8.49% higher than that of the IDUV-based method (RMSE: 2.98 m). [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]

Grain protein content (GPC) is an important indicator of wheat quality. Earlier estimation of wheat GPC based on remote sensing provided effective decision to adapt optimized strategies for grain harvest, which is of great significance for agricultural production. The objectives of this field study are: (i) To assess the ability of spectral vegetation indices (VIs) of Sentinel 2 data to detect the wheat nitrogen (N) attributes related to the grain quality of winter wheat production, and (ii) to examine the accuracy of wheat N status and GPC estimation models based on different VIs and wheat nitrogen parameters across Analytical Spectra Devices (ASD) and Unmanned Aerial Vehicle (UAV) hyper-spectral data-simulated sentinel data and the real Sentinel-2 data. In this study, four nitrogen parameters at the wheat anthesis stage, including plant nitrogen accumulation (PNA), plant nitrogen content (PNC), leaf nitrogen accumulation (LNA), and leaf nitrogen content (LNC), were evaluated for their relationship between spectral parameters and GPC. Then, a multivariate linear regression method was used to establish the wheat nitrogen and GPC estimation model through simulated Sentinel-2A VIs. The coefficients of determination ( R 2 ) of four nitrogen parameter models were all greater than 0.7. The minimum R 2 of the prediction model of wheat GPC constructed by four nitrogen parameters combined with VIs was 0.428 and the highest R 2 was 0.467. The normalized root mean square error (nRMSE) of the four nitrogen estimation models ranged from 26.333% to 29.530% when verified by the ground-measured data collected from the Beijing suburbs, and the corresponding nRMSE for the GPC-predicted models ranged from 17.457% to 52.518%. The accuracy of the estimated model was verified by UAV hyper-spectral data which had resized to different spatial resolution collected from the National Experimental Station for Precision Agriculture. The normalized root mean square error (nRMSE) of the four nitrogen estimation models ranged from 16.9% to 37.8%, and the corresponding nRMSE for the GPC-predicted models ranged from 12.3% to 13.2%. The relevant models were also verified by Sentinel-2A data collected in 2018 while the minimum nRMSE for GPC invert model based on PNA was 7.89% and the maximum nRMSE of the GPC model based on LNC was 12.46% in Renqiu district, Hebei province. The nRMSE for the wheat nitrogen estimation model ranged from 23.200% to 42.790% for LNC and PNC. These data demonstrate that freely available Sentinel-2 imagery can be used as an important data source for wheat nutrition and grain quality monitoring. [ABSTRACT FROM AUTHOR]