A multi-resolution global land cover dataset through multisource data aggregation

Recent developments of 30 m global land characterization datasets(e.g., land cover, vegetation continues field) represent the finest spatial resolution inputs for global scale studies. Here, we present results from further improvement to land cover mapping and impact analysis of spatial resolution on area estimation for different land cover types. We proposed a set of methods to aggregate two existing 30 m resolution circa 2010 global land cover maps, namely FROM-GLC(Finer Resolution Observation and Monitoring-Global Land Cover) and FROM-GLC-seg(Segmentation), with two coarser resolution global maps on development, i.e., Nighttime Light Impervious Surface Area(NL-ISA) and MODIS urban extent(MODIS-urban), to produce an improved 30 m global land cover map—FROM-GLC-agg(Aggregation). It was post-processed using additional coarse resolution datasets(i.e., MCD12Q1, GlobCover2009, MOD44 W etc.) to reduce land cover type confusion. Around 98.9% pixels remain 30 m resolution after some post-processing to this dataset. Based on this map, majority aggregation and proportion aggregation approaches were employed to create a multi-resolution hierarchy(i.e., 250 m, 500 m, 1 km, 5 km, 10 km, 25 km, 50 km, 100 km) of land cover maps to meet requirements for different resolutions from different applications. Through accuracy assessment, we found that the best overall accuracies for the post-processed base map(at 30 m) and the three maps subsequently aggregated at 250 m, 500 m, 1 km resolutions are 69.50%, 76.65%, 74.65%, and 73.47%, respectively. Our analysis of area-estimation biases for different land cover types at different resolutions suggests that maps at coarser than 5 km resolution contain at least 5% area estimation error for most land cover types. Proportion layers, which contain precise information on land cover percentage, are suggested for use when coarser resolution land cover data are required.     

Use of multi‐platform,multi‐temporal remote‐sensing data for calibration of a distributed hydrological model: an application in the Arno basin,Italy

Images from satellite platforms are a valid aid in order to obtain distributed information about hydrological surface states and parameters needed in calibration and validation of the water balance and flood forecasting. Remotely sensed data are easily available on large areas and with a frequency compatible with land cover changes. In this paper, remotely sensed images from different types of sensor have been utilized as a support to the calibration of the distributed hydrological model MOBIDIC, currently used in the experimental system of flood forecasting of the Arno River Basin Authority. Six radar images from ERS‐2 synthetic aperture radar (SAR) sensors (three for summer 2002 and three for spring–summer 2003) have been utilized and a relationship between soil saturation indexes and backscatter coefficient from SAR images has been investigated. Analysis has been performed only on pixels with meagre or no vegetation cover, in order to legitimize the assumption that water content of the soil is the main variable that influences the backscatter coefficient. Such pixels have been obtained by considering vegetation indexes (NDVI) and land cover maps produced by optical sensors (Landsat‐ETM). In order to calibrate the soil moisture model based on information provided by SAR images, an optimization algorithm has been utilized to minimize the regression error between saturation indexes from model and SAR data and error between measured and modelled discharge flows. Utilizing this procedure, model parameters that rule soil moisture fluxes have been calibrated, obtaining not only a good match with remotely sensed data, but also an enhancement of model performance in flow prediction with respect to a previous calibration with river discharge data only. Copyright © 2006 John Wiley & Sons, Ltd.     

On the significance of crop‐type information for the simulation of catchment hydrology

Spatially distributed hydrological models require information on the land cover pattern, as it directly influences the generation of run‐off. However, it is not clear which detail of land cover information is suitable for simulating the catchment hydrology. A better understanding of the relationship between the land cover detail and the hydrological processes is therefore required as it would enhance a successful application of the hydrological model. This study investigates the relevance of accurate information about the crop types in the catchment for the simulation of run‐off, baseflow and evapotranspiration. Results reveal that adding knowledge about the crop type in the model simulation is redundant when area‐averaged water flux predictions are intended. On the contrary, when a spatial distribution of water flux predictions is desired, it is meaningful to increase the land cover detail because an increased heterogeneity in the land cover map induces an increased heterogeneity in the hydrological response and locally reduces the uncertainty in the model prediction up to 50%. To assess the land cover uncertainty and to locate the regions for which the reduction in prediction uncertainty is the highest, the thematic uncertainty map (constructed through fuzzy logic) is shown to be a useful tool. Copyright © 2014 John Wiley & Sons, Ltd.     

Estimating landslide susceptibility areas considering the uncertainty inherent in modeling methods

Landslides are one of the most dangerous types of natural disasters, and damage due to landslides has been increasing in certain regions of the world because of increased precipitation. Policy decision makers require reliable information that can be used to establish spatial adaptation plans to protect people from landslide hazards. Researchers presently identify areas susceptible to landslides using various spatial distribution models. However, such data are associated with a high amount of uncertainty. This study focuses on quantifying the uncertainty of several spatial distribution models and identifying the effectiveness of various ensemble methods that can be used to provide reliable information to support policy decisions. The area of study was Inje-gun, Republic of Korea. Ten models were selected to assess landslide susceptibility. Moreover, five ensemble methods were selected for the aggregated results of the 10 models. The uncertainty was quantified using the coefficient of variation and the uncertainty map we developed revealed areas with strongly differing values among single models. A matrix map was created using an ensemble map and a coefficient of variation map. Using matrix analysis, we identified the areas that are most susceptible to landslides according to the ensemble model with a low uncertainty. Thus, the ensemble model can be a useful tool for supporting decision makers. The framework of this study can also be employed to support the establishment of landslide adaptation plans in other areas of the Republic of Korea and in other countries.     

一种综合地学数据的矿产资源潜力制图方法

研究适应信息化时代特征的矿产资源潜力制图新技术、新方法对推动矿产资源评价理论与技术的发展具有重要的意义.笔者把GIS技术、图像分类算法和空间统计学理论进行有机集成,在空间统计学的空间结构分析技术和遥感图像纹理分类算法的基础上,提出了一种以综合地学数据（地质、地球物理、地球化学和遥感图像数据等）为基本数据源的矿产资源潜力自动制图方法.该方法的技术流程为：①数据准备,即对地球物理和地球化学勘探数据进行预处理,生成一个物化遥综合图像文件;②图像空间结构性分析和纹理图像生成,以综合地学图像为研究对象,用空间统计学的结构分析技术研究地学数据综合图像的空间结构性,生成纹理图像;③纹理图像多元分类,用实验变差函数纹理分类方法对研究区进行多元分类,生成分类专题图;④分类后处理,用叠置分析修正空间分类结果,生成区域矿产资源潜力分布图.     

Effects of reduced land cover detail on hydrological model response

Problems related to scale continue to be at the forefront of research in hydrology. Past research into issues of scale has focused mainly on digital elevation model grid size, the appropriate number and size of sub‐areas for subdividing a watershed, parameter transferability between watersheds and appropriate scales for linking hydrological and general circulation models. Much less attention has been given to the effects of scale on the representation of land cover and hydrological model response. Recent studies with respect to changes in land cover and hydrologic response have tended to focus on the issue of land cover maturity and the conversion of land through agricultural and forestry practices. The focus of this study is to examine the impact of the level of detail at which land cover is represented in modelling the hydrological response of Wolf Creek Basin in northwest Canada. A grid‐based land cover map with a spatial resolution of 30 m is coarsened or smoothed using several common grid‐based methods of aggregating categorical data, including: pixel thinning, modal smoothing and modal aggregation. A majority rule method based on polygons is also applied to the 30 m base cover. The SLURP hydrologic model is calibrated for the base cover and used as a reference for comparing simulations for the coarsened or ‘generalized’ land cover maps. Results of the simulations are compared to examine the sensitivity of hydrologic response to generalized land cover information. Comparisons of the SLURP model runs for Wolf Creek suggest that reducing the level of detail of land cover information generally has a limited effect on hydrologic response at the outlet. However, results for averages of water balance components across the basin suggest that the local variability of hydrologic response is affected in general. Copyright © 2007 John Wiley & Sons, Ltd.     

Assessment of calanchi and rill–interrill erosion susceptibilities using terrain analysis and geostochastics: A case study in the Oltrepo Pavese,Northern Apennines,Italy

Soil erosion is one of the most important environmental problems distributed worldwide. In the last decades, numerous studies have been published on the assessment of soil erosion and the related processes and forms using empirical, conceptual and physically based models. For the prediction of the spatial distribution, more and more sophisticated stochastic modelling approaches have been proposed – especially on smaller spatial scales such as river basins. In this work, we apply a maximum entropy model (MaxEnt) to evaluate badlands (calanchi) and rill–interrill (sheet erosion) areas in the Oltrepo Pavese (Northern Apennines, Italy). The aim of the work is to assess the important environmental predictors that influence calanchi and rill–interrill erosion at the regional scale. We used 13 topographic parameters derived from a 12 m digital elevation model (TanDEM-X) and data on the lithology and land use. Additional information about the vegetation is introduced through the normalized difference vegetation index based on remotely sensed data (ASTER images). The results are presented in the form of susceptibility maps showing the spatial distribution of the occurrence probability for calanchi and rill–interrill erosion. For the validation of the MaxEnt model results, a support vector machine approach was applied. The models show reliable results and highlight several locations of the study area that are potentially prone to future soil erosion. Thus, coping and mitigation strategies may be developed to prevent or fight the soil erosion phenomenon under consideration. © 2020 John Wiley & Sons, Ltd.     

Recharge estimation using remotely sensed evapotranspiration in an irrigated catchment in southeast Australia

Reliable estimates of groundwater recharge are required for the sustainable management of surface and ground water resources in semi‐arid regions particularly in irrigated regions. In this study, groundwater recharge was estimated for an irrigated catchment in southeast Australia using a semi‐distributed hydrological model (SWAT). The model was calibrated under the dry climatic conditions for the period from August 2002 to July 2003 using flow and remotely sensed evapotranspiration (ET). The model was able to simulate observed monthly drain flow and spatially distributed remotely sensed ET. Recharge tended to be higher for irrigated land covers, such as perennial pasture, than for non‐irrigated land. On average, the estimated annual catchment recharge ranged between 147 and 289 mm which represented about 40% of the total rainfall and irrigation inputs. The optimized soil parameters indirectly reflected flow bypassing the soil matrix that could be responsible for this substantial amount of recharge. Overall, the estimated recharge was much more than that previously estimated for the wetter years. Copyright © 2011 John Wiley & Sons, Ltd.     

Integration of remotely sensed C factor into SWAT for modelling sediment yield

The C factor, representing the impact of plant and ground cover on soil loss, is one of the important factors of the Modified Universal Soil Loss Equation (MUSLE) in the Soil and Water Assessment Tool (SWAT) to model sediment yield. The daily update of C factors in SWAT was originally determined by land use types and plant growth cycles. This does not reflect the spatial variation of C values that exists within a large land use area. We present a new approach to integrate remotely sensed C factors into SWAT for highlighting the effect of detailed vegetative cover data on soil erosion and sediment yield. First, the C factor was estimated using the abundance of ground components extracted from remote sensing images. Then, the gridding data of the C factor were aggregated to hydrological response units (HRUs), instead of to land use units of SWAT. In the end, the C factor values in HRUs were integrated into SWAT to predict sediment yield by modifying the ysed subroutine. This substitution work not only increases the spatial variation of the C factor in SWAT, but also makes it possible to utilize other sources of C databases rather than those from the United States. The demonstration in the Dage basin shows that the modified SWAT produces reasonable results in water flow simulation and sediment yield prediction using remotely sensed C values. The Nash–Sutcliffe efficiency coefficient (E_NS) and R² for surface runoff range from 0·69 to 0·77 and 0·73 to 0·87, respectively. The coefficients E_NS and R² for sediment yield were generally above 0·70 and 0·60, respectively. The soil erosion risk map based on sediment yield prediction at the HRU level illustrates instructive details on spatial distribution of soil loss. Copyright © 2011 John Wiley & Sons, Ltd.     

Exploiting the use of physical information for the calibration of a lumped hydrological model

In hydrological modelling, the challenge is to identify an optimal strategy to exploit tools and available observations in order to enhance model reliability. The increasing availability of data promotes the use of new calibration techniques able to make use of additional information on river basins. In the present study, a lumped hydrological model—designed with the aim of utilizing remotely sensed data—is introduced and calibrated, adopting four different schemes that adopt, to varying extents, available physical information. The physically consistent conceptualization of the hydrological model used allowed development of a step by step calibration based on a combination of information, such as remotely sensed data describing snow cover, recession curves obtained from streamflow measurements, and time series of surface run‐off obtained with a baseflow mathematical filter applied to the streamflow time‐series. Results suggest that the use of physical information in the calibration procedure tends to increase model reliability with respect to approaches where the parameters are calibrated using an overall statistic based, considerably or exclusively, on streamflow data.     

Macrobenthos habitat mapping in a tidal flat using remotely sensed data and a GIS-based probabilistic model

This paper proposes and tests a method of producing macrofauna habitat potential maps based on a weights-of-evidence model (a probabilistic approach) for the Hwangdo tidal flat, Korea. Samples of macrobenthos were collected during field work, and we considered five mollusca species for habitat mapping. A weights-of-evidence model was used to calculate the relative weights of 10 control factors that affect the macrobenthos habitat. The control factors were compiled as a spatial database from remotely sensed data combined with GIS analysis. The relative weight of each factor was integrated as a species potential index (SPI), which produced habitat potential maps. The maps were compared with the surveyed habitat locations, revealing a strong correlation between the potential maps and species locations. The combination of a GIS-based weights-of-evidence model and remote sensing techniques is an effective method in determining areas of macrobenthos habitat potential in a tidal flat setting.     

On the validation of a coupled water and energy balance model at small catchment scales

Catchment runoff is the most widely used catchment scale measurement in modelling studies, and we have a reasonable degree of confidence in its accuracy. The advent of satellites gives access to a new suite of measurements taken over a defined spatial range. These measurements, principally reflected or emitted radiation, provide hydrologists with new possibilities for quantifying the state of a catchment. Surface temperatures can be readily measured by a satellite on a scale comparable to the size of a small catchment.

In this paper we show that satellite sensed temperatures can provide an important measure of catchment status, which can complement runoff measurements in water balance studies. A one-dimensional model, which couples the land surface energy balance with the soil and surface water balance is tested by comparison with runoff and with remotely sensed surface temperature measurements. Simulations have been run over four years for two small catchments which have a fairly homogeneous vegetation, one being forest and its neighbour pasture. Satellite “surface” temperatures have been interpreted in terms of the energy balance, and used as a test of modelling accuracy. An “effective” surface temperature is calculated as a weighted mean of temperatures of the separate soil and leaf surfaces. This modelled “effective” temperature correlates well with Landsat TM surface temperatures.

When pasture replaces forest, the model predicts a reduction in evapotranspiration of around 30%, a three-fold increase in runoff, and an increase in mean soil moisture status. The change to pasture also results in a rise in mean effective surface temperature of about 4°C, and an increase in summer diurnal temperature range from 10 to 22°C. The winter diurnal temperature range is similar for both vegetation systems.

Inclusion of soil moisture variability in thermal properties results in an increase in mean daily maximum temperature of about 2°C in summer and winter, without much change in daily minima. The daily mean temperature is not significantly affected.     

A CASE STUDY ON RULE-BASED MULTIPLE SOURCE DATA CLASSIFICATION WITH GIS

IINTRODUCTIONLanduse/coverisoneofthemostimportalfactorseffectingsoilandwaterloss.Researchonlanduse/coverandwaterandsoillosswithremotesensinghasbeendonemuchinthepast.Butremotelysensedimageryonlycarriestheinstantaneousandtwo-dimensionalinformationofitsprototypegeographicobjects(ChenandZhao,1990).Thereforemathematicalandphysicalprocessingonremotelysensedimageoftenproduceindefiniteandunreliableresult.Inordertoimprovetheprecisionofclassification,otherdatasetssuchastopographicmaps,thematicmaps…     

Updating of snow depletion curve with remote sensing data

A method for using remotely sensed snow cover information in updating a hydrological model is developed, based on Bayes' theorem. A snow cover mass balance model structure adapted to such use of satellite data is specified, using a parametric snow depletion curve in each spatial unit to describe the subunit variability in snow storage. The snow depletion curve relates the accumulated melt depth to snow‐covered area, accumulated snowmelt runoff volume, and remaining snow water equivalent. The parametric formulation enables updating of the complete snow depletion curve, including mass balance, by satellite data on snow coverage. Each spatial unit (i.e. grid cell) in the model maintains a specific depletion curve state that is updated independently. The uncertainty associated with the variables involved is formulated in terms of a joint distribution, from which the joint expectancy (mean value) represents the model state. The Bayesian updating modifies the prior (pre‐update) joint distribution into a posterior, and the posterior joint expectancy replaces the prior as the current model state. Three updating experiments are run in a 2400 km² mountainous region in Jotunheimen, central Norway (61°N, 9°E) using two Landsat 7 ETM+ images separately and together. At 1 km grid scale in this alpine terrain, three parameters are needed in the snow depletion curve. Despite the small amount of measured information compared with the dimensionality of the updated parameter vector, updating reduces uncertainty substantially for some state variables and parameters. Parameter adjustments resulting from using each image separately differ, but are positively correlated. For all variables, uncertainty reduction is larger with two images used in conjunction than with any single image. Where the observation is in strong conflict with the prior estimate, increased uncertainty may occur, indicating that prior uncertainty may have been underestimated. Copyright © 2006 John Wiley & Sons, Ltd.     

An EcoCity model for regulating urban land cover structure and thermal environment: Taking Beijing as an example

Urban land-use/cover changes and their effects on the eco-environment have long been an active research topic in the urbanization field. However, the characteristics of urban inner spatial heterogeneity and its quantitative relationship with thermal environment are still poorly understood, resulting in ineffective application in urban ecological planning and management.Through the integration of "spatial structure theory" in urban geography and "surface energy balance" in urban climatology, we proposed a new concept of urban surface structure and thermal environment regulation to reveal the mechanism between urban spatial structure and surface thermal environment. We developed the EcoCity model for regulating urban land cover structure and thermal environment, and established the eco-regulation thresholds of urban surface thermal environments. Based on the comprehensive analysis of experimental observation, remotely sensed and meteorological data, we examined the spatial patterns of urban habitation, industrial, infrastructure service, and ecological spaces. We examined the impacts of internal land-cover components(e.g., urban impervious surfaces, greenness, and water) on surface radiation and heat flux. This research indicated that difference of thermal environments among urban functional areas is closely related to the proportions of the land-cover components.The highly dense impervious surface areas in commercial and residential zones significantly increased land surface temperature through increasing sensible heat flux, while greenness and water decrease land surface temperature through increasing latent heat flux. We also found that different functional zones due to various proportions of green spaces have various heat dissipation roles and ecological thresholds. Urban greening projects in highly dense impervious surfaces areas such as commercial, transportation, and residential zones are especially effective in promoting latent heat dissipation efficiency of vegetation, leading to strongly cooling effect of unit vegetation coverage. This research indicates that the EcoCity model provides the fundamentals to understand the coupled mechanism between urban land use structure and surface flux and the analysis of their spatiotemporal characteristics. This model provides a general computational model system for defining urban heat island mitigation, the greening ratio indexes, and their regulating thresholds for different functional zones.     

Conditional simulation of remotely sensed rainfall data using a non-Gaussian v-transformed copula

Quantification of rainfall and its spatial and temporal variability is extremely important for reliable hydrological and meteorological modeling. While rain gauge measurements do not provide reasonable areal representation of rainfall, remotely sensed precipitation estimates offer much higher spatial resolution. However, uncertainties associated with remotely sensed rainfall estimates are not well quantified. This issue is important considering the fact that uncertainties in input rainfall are the main sources of error in hydrologic processes. Using an ensemble of rainfall estimates that resembles multiple realizations of possible true rainfall, one can assess uncertainties associated with remotely sensed rainfall data. In this paper, ensembles are generated by imposing rainfall error fields over remotely sensed rainfall estimates. A non-Gaussian copula-based model is introduced for simulation of rainfall error fields. The v-transformed copula is employed to describe the dependence structure of rainfall error estimates without the influence of the marginal distribution. Simulations using this model can be performed unconditionally or conditioned on ground reference measurements such that rain gauge data are honored at their locations. The presented model is implemented for simulation of rainfall ensembles across the Little Washita watershed, Oklahoma. The results indicate that the model generates rainfall fields with similar spatio-temporal characteristics and stochastic properties to those of observed rainfall data.     

Prediction of vegetation anomalies over an inland river basin in north‐western China

Vegetation cover plays an important role in linking the atmosphere, water, and land and is deemed as a key indicator in the terrestrial ecological system. Therefore, it is of great importance to monitor vegetation dynamics and understand the mechanisms of vegetation change, including that driven by climate change. This study examines (a) the evolution of vegetation dynamics over the Heihe River Basin in the typical arid zone in north‐western China using nonparametric Mann–Kendall test and Thiel Sen's slope; (b) the relationships between remotely sensed vegetation indices (normalized difference vegetation index [NDVI] and enhanced vegetation index [EVI]) and hydroclimatic variables based on correlation analysis; and (c) the prediction of vegetation anomalies using a multiple linear regression model. For the analysis, the Moderate Resolution Imaging Spectroradiometer NDVI/EVI product and the gridded daily meteorological data at a spatial resolution of 0.125° over the period 2001–2010 are considered. The results indicate that vegetation cover improved over a large proportion during 2001–2010, with a significant trend towards warm and wet, characterized by an increase in average annual temperature and precipitation by 0.042 °C/year and 5.8 mm/year, respectively. We test the feasibility of NDVI and EVI in quantifying the responses of vegetation anomaly to climate change and develop a statistical model to predict vegetation dynamics in the basin. The NDVI‐based model is found to be more reliable than the EVI‐based model, partly due to the vegetation characteristics and geomorphologic properties of the study region. The proposed model performs well when there is no lag time between meteorological factors and vegetation indices for grassland and cropland, whereas 1‐month lead time prediction is found to be best for forest. The soil water content is introduced as an extra explanatory variable, which effectively improves the prediction accuracy for different land use types. In general, the predictive ability of the proposed model is stable and satisfactory, and the model can provide useful early warning information for regional water resources management under changing climate.     

Modelling Freshwater Resources at the Global Scale: Challenges and Prospects

Quantification of spatially and temporally resolved water flows and water storage variations for all land areas of the globe is required to assess water resources, water scarcity and flood hazards, and to understand the Earth system. This quantification is done with the help of global hydrological models (GHMs). What are the challenges and prospects in the development and application of GHMs? Seven important challenges are presented. (1) Data scarcity makes quantification of human water use difficult even though significant progress has been achieved in the last decade. (2) Uncertainty of meteorological input data strongly affects model outputs. (3) The reaction of vegetation to changing climate and CO₂ concentrations is uncertain and not taken into account in most GHMs that serve to estimate climate change impacts. (4) Reasons for discrepant responses of GHMs to changing climate have yet to be identified. (5) More accurate estimates of monthly time series of water availability and use are needed to provide good indicators of water scarcity. (6) Integration of gradient-based groundwater modelling into GHMs is necessary for a better simulation of groundwater–surface water interactions and capillary rise. (7) Detection and attribution of human interference with freshwater systems by using GHMs are constrained by data of insufficient quality but also GHM uncertainty itself. Regarding prospects for progress, we propose to decrease the uncertainty of GHM output by making better use of in situ and remotely sensed observations of output variables such as river discharge or total water storage variations by multi-criteria validation, calibration or data assimilation. Finally, we present an initiative that works towards the vision of hyperresolution global hydrological modelling where GHM outputs would be provided at a 1-km resolution with reasonable accuracy.     

The accuracy of drainage network delineation as a function of environmental factors: A case study in Central and Northern Sweden

Drainage networks delineated from Digital Elevation Models (DEMs), are the basis for the modelling of geomorphological and hydrological processes, biogeochemical cycling, and water resources management. Besides providing effective models of water flows, automatically extracted drainage networks based on topography can diverge from reality to varying degrees. The variability of such disagreement within catchments has rarely been examined as a function of the heterogeneity of land cover, soil type, and slope in the catchment of interest. This research gap might not only substantially limit our knowledge of the uncertainty of hydrological prediction, but can also cause problems for users attempting to use the data at a local scale. Using 1:100000 scale land cover maps, Quaternary deposits maps, and 2 m resolution DEMs, it is found that the accuracy of delineated drainage networks tends to be lower in areas with denser vegetation, lower hydraulic conductivity, and higher erodibility. The findings of this study could serve as a guide for the more thoughtful usage of delineated drainage networks in environmental planning, and in the uncertainty analysis of hydrological and biochemical predictions. Therefore, this study makes a first attempt at filling the knowledge gap described above.     

Catchment‐scale sensitivity and uncertainty in water quality modelling

In this paper, we assess the performance of the catchment model SIMulated CATchment model (SIMCAT), to predict nitrate and soluble reactive phosphorus concentrations against four monitoring regimes with different spatial and temporal sampling frequencies. The Generalised Likelihood Uncertainty Estimation (GLUE) uncertainty framework is used, along with a general sensitivity analysis to understand relative parameter sensitivity. Improvements to model calibration are explored by introducing more detailed process representation using the Integrated Catchments model (INCA) water quality model, driven by the European hydrological predictions for the environment model. The results show how targeted sampling of headwater watercourses upstream of point discharges is essential for calibrating diffuse loads and can exert a strong influence on the whole‐catchment model performance. Further downstream, if the point discharges and loads are accurately represented, then the improvement in the catchment‐scale model performance is relatively small as more calibration points are added or frequency is increased. The higher‐order, dynamic model integrated catchments model of phosphorus dynamics, which incorporates sediment and biotic interaction, resulted in improved whole‐catchment performance over SIMCAT, although there are still large epistemic uncertainties from land‐phase export coefficients and runoff. However, the very large sampling errors in routine monitoring make it difficult to invest confidence in the modelling, especially because we know phosphorous transport to be very episodic and driven by high flow conditions for which there are few samples. The environmental modelling community seems to have been stuck in this position for some time, and whilst it is useful to use an uncertainty framework to highlight these issues, it has not widely been adopted, perhaps because there is no clear mechanism to allow uncertainties to influence investment decisions. This raises the question as to whether it might better place a cost on uncertainty and use this to drive more data collection or improved models, before making investment decisions concerning, for example, mitigation strategies. Copyright © 2016 John Wiley & Sons, Ltd.