A near-global, high resolution land surface parameter dataset for the variable infiltration capacity model

Hydrologic models predict the spatial and temporal distribution of water and energy at the land surface. Currently, parameter availability limits global-scale hydrologic modelling to very coarse resolution, hindering researchers from resolving fine-scale variability. With the aim of addressing this problem, we present a set of globally consistent soil and vegetation parameters for the Variable Infiltration Capacity (VIC) model at 1/16° resolution (approximately 6 km at the equator), with spatial coverage from 60°S to 85°N. Soil parameters derived from interpolated soil profiles and vegetation parameters estimated from space-based MODIS measurements have been compiled into input files for both the Classic and Image drivers of the VIC model, version 5. Geographical subsetting codes are provided, as well. Our dataset provides all necessary land surface parameters to run the VIC model at regional to global scale. We evaluate VICGlobal’s ability to simulate the water balance in the Upper Colorado River basin and 12 smaller basins in the CONUS, and their ability to simulate the radiation budget at six SURFRAD stations in the CONUS.

Radiometric Normalization for Cross-Sensor Optical Gaofen Images with Change Detection and Chi-Square Test

As the number of cross-sensor images increases continuously, the surface reflectance of these images is inconsistent at the same ground objects due to different revisit periods and swaths. The surface reflectance consistency between cross-sensor images determines the accuracy of change detection, classification, and land surface parameter inversion, which is the most widespread application. We proposed a relative radiometric normalization (RRN) method to improve the surface reflectance consistency based on the change detection and chi-square test. The main contribution was that a novel chi-square test automatically extracts the stably unchanged samples between the reference and subject images from the unchanged regions detected by the change-detection method. We used the cross-senor optical images of Gaofen-1 and Gaofen-2 to test this method and four metrics to quantitatively evaluate the RRN performance, including the Root Mean Square Error (RMSE), spectral angle cosine, structural similarity, and CIEDE2000 color difference. Four metrics demonstrate the effectiveness of our proposed RRN method, especially the reduced percentage of RMSE after normalization was more than 80%. Comparing the radiometric differences of five ground features, the surface reflectance curve of two Gaofen images showed more minor differences after normalization, and the RMSE was smaller than 50 with the reduced percentages of about 50–80%. Moreover, the unchanged feature regions are detected by the change-detection method from the bitemporal Sentinel-2 images, which can be used for RRN without detecting changes in subject images. In addition, extracting samples with the chi-square test can effectively improve the surface reflectance consistency.

Dynamic Cosine Method for Normalizing Incidence Angle Effect on C-band Radar Backscattering Coefficient for Maize Canopies Based on NDVI

Wide mode SAR images have an apparent incidence angle effect. The existing incident angle normalization methods assume that the relationship between the incident angle (θ) and the backscattering coefficient (σPQ) does not change with the growth stage of crops, which is in conflict with the real-life situation. Therefore, the normalization results of σPQ based on these existing methods will affect the accuracy of object classification, target recognition, and land surface parameter inversion. Here, the change in θ-σPQ relationship was investigated based on time-series (April to October) σPQ of maize canopies in northeast China, and a dynamic method based on normalized difference vegetation index (NDVI) was developed to normalize the effect of θ on σPQ. Through the accuracy evaluation, the following conclusions are obtained: (1) the dependence (referring to N) of Sentinel 1 C-band σPQ on θ varies with maize NDVI. In addition, the value of N changed from 9.35 to 0.66 at VV polarization from bare soil to biomass peak, and from 6.26 to 0.99 at VH polarization; (2) a dynamic method was proposed to quantify the change of N based on its strong correlation with NDVI, indicated by R2 of 0.82 and 0.80 for VV and VH polarization, respectively; and (3) the overall root mean square error of normalized σPQ based on the newly-developed dynamic method is 0.51 dB, and this accuracy outperforms the original first-order cosine method (1.37 dB) and cosine square law method (1.08 dB) by about 63% and 53% on the whole. This study provides a dynamic framework for normalizing radar backscatter coefficient, improving the retrieval accuracy of land surface parameters from radar remote sensing.

An Efficient Hierarchical Representation Approach of Remote Sensing Application Modeling Based on Distributed Environment

In Earth science, information science, space science, and other disciplines, scientists use the land surface parameter inversion method in their work, applying this to the atmosphere, vegetation, soil, drought, and so on. Multidisciplinary experts sometimes collaborate on a particular application. However, these remote sensing models do not have a unified method of description and management and cannot effectively achieve the sharing of models and data resources. It is also hard to meet user demand for global data and models in the current state, especially in the face of global problems and long-term series problems. In this paper, we examine the scientific questions of the computability and scalability of remote sensing models. This paper adopts a data dependency approach to describe a remote sensing model and implements a hierarchical unified description and management method using modelling based on four layers: a data-processing view, an atomic model view, an on-demand resource package view, and a workflow view. We choose three typical remote sensing models for disaster monitoring as use cases and describe the practical application process of the proposed method. The results demonstrate the advantages and powerful capabilities of this efficient method.

How Does Scale Effect Influence Spring Vegetation Phenology Estimated from Satellite-Derived Vegetation Indexes?

As an important land-surface parameter, vegetation phenology has been estimated from observations by various satellite-borne sensors with substantially different spatial resolutions, ranging from tens of meters to several kilometers. The inconsistency of satellite-derived phenological metrics (e.g., green-up date, GUD, also known as the land-surface spring phenology) among different spatial resolutions, which is referred to as the “scale effect” on GUD, has been recognized in previous studies, but it still needs further efforts to explore the cause of the scale effect on GUD and to quantify the scale effect mechanistically. To address these issues, we performed mathematical analyses and designed up-scaling experiments. We found that the scale effect on GUD is not only related to the heterogeneity of GUD among fine pixels within a coarse pixel, but it is also greatly affected by the covariation between the GUD and vegetation growth speed of fine pixels. GUD of a coarse pixel tends to be closer to that of fine pixels with earlier green-up and higher vegetation growth speed. Therefore, GUD of the coarse pixel is earlier than the average of GUD of fine pixels, if the growth speed is a constant. However, GUD of the coarse pixel could be later than the average from fine pixels, depending on the proportion of fine pixels with later GUD and higher growth speed. Based on those mechanisms, we proposed a model that accounted for the effects of heterogeneity of GUD and its co-variation with growth speed, which explained about 60% of the scale effect, suggesting that the model can help convert GUD estimated at different spatial scales. Our study provides new mechanistic explanations of the scale effect on GUD.

Contribution Analysis of the Spatial-Temporal Changes in Streamflow in a Typical Elevation Transitional Watershed of Southwest China over the Past Six Decades

Attribution analyses on streamflow variation to changing climate and land surface characteristics are critical in studies of watershed hydrology. However, attribution results may differ greatly on different spatial and temporal scales, which has not been extensively studied previously. This study aims to investigate the spatial-temporal contributions of climate change and underlying surface variation to streamflow alteration using Budyko framework. Jiangling River Watershed (JRW), a typical landform transitional watershed in Southwest China, was chosen as the study area. The watershed was firstly divided into eight sub-basins by hydrologic stations, and hydrometeorological series (1954–2015) were divided into sub-intervals to discriminate spatial-temporal features. The results showed that long-term tendencies of hydrometeorological variables, i.e., precipitation (P), potential evapotranspiration (E0), and runoff depth (R), exhibited clear spatial patterns, which were highly related to topographic characteristics. Additionally, sensitivity analysis, which interpreted the effect of one driving factor by unit change, showed that climate factors P and E0, and catchment characteristics (land surface parameter n) played positive, negative, and negative roles in R, according to elastic coefficients (ε), respectively. The spatial distribution of ε illustrated a greater sensitivity and heterogeneity in the plateau and semi-humid regions (upstream). Moreover, the results from attribution analysis showed that the contribution of the land surface factor accounted for approximately 80% of the R change for the entire JRW, with an obvious spatial variation. Furthermore, tendencies of the contribution rates demonstrated regulations across different sub-regions: a decreasing trend of land surface impacts in trunk stream regions and increasing tendencies in tributary regions, and vice versa for climate impacts. Overall, both hydrometeorological variables and contributions of influencing factors presented regularities in long-term tendencies across different sub-regions. More particularly, the impact of the primary influencing factor on all sub-basins exhibited a decreasing trend over time. The evidence that climate and land surface change act on streamflow in a synergistic way, would complicate the attribution analysis and bring a new challenge to attribution analysis.

Land Surface Parameter and State Perturbations in the Global Ensemble Forecast System

The National Centers for Environmental Prediction (NCEP) Global Ensemble Forecast System (GEFS) is underdispersive near the surface, a common characteristic of ensemble prediction systems. Here, several methods for increasing the spread are tested, including perturbing soil initial conditions, soil tendencies, and surface parameters, with physically based perturbations. Perturbations are applied to the soil initial conditions based on empirical orthogonal functions (EOFs) of differences between normalized soil moisture states from two land surface models (LSMs). Perturbations to roughness lengths for heat and momentum, soil hydraulic conductivity, stomatal resistance, vegetation fraction, and albedo are applied, with the amplitude and perturbation scales based on previous research. Soil moisture and temperature tendencies are also perturbed using a stochastic perturbation scheme. The results show that surface perturbations, through their impact on 2-m temperature spread, have a modest positive impact on the skill of short-range ensemble forecasts. However, adjusting the forecasts using an estimate of the systematic bias shows that bias correction has a greater impact on the forecast reliability than surface perturbations, indicating that systematic bias in the model needs to be addressed as well.

Land-surface parameter optimisation using data assimilation techniques: the adJULES system V1.0

Land-surface models (LSMs) are crucial components of the Earth system models (ESMs) that are used to make coupled climate–carbon cycle projections for the 21st century. The Joint UK Land Environment Simulator (JULES) is the land-surface model used in the climate and weather forecast models of the UK Met Office. JULES is also extensively used offline as a land-surface impacts tool, forced with climatologies into the future. In this study, JULES is automatically differentiated with respect to JULES parameters using commercial software from FastOpt, resulting in an analytical gradient, or adjoint, of the model. Using this adjoint, the adJULES parameter estimation system has been developed to search for locally optimum parameters by calibrating against observations. This paper describes adJULES in a data assimilation framework and demonstrates its ability to improve the model–data fit using eddy-covariance measurements of gross primary production (GPP) and latent heat (LE) fluxes. adJULES also has the ability to calibrate over multiple sites simultaneously. This feature is used to define new optimised parameter values for the five plant functional types (PFTs) in JULES. The optimised PFT-specific parameters improve the performance of JULES at over 85 % of the sites used in the study, at both the calibration and evaluation stages. The new improved parameters for JULES are presented along with the associated uncertainties for each parameter.