Spatiotemporal estimation of snow depth using point data from snow stakes,digital terrain models,and satellite data

Snow availability in Alpine catchments plays an important role in water resources management. In this paper, we propose a method for an optimal estimation of snow depth (areal extension and thickness) in Alpine systems from point data and satellite observations by using significant explanatory variables deduced from a digital terrain model. It is intended to be a parsimonious approach that may complement physical‐based methodologies. Different techniques (multiple regression, multicriteria analysis, and kriging) are integrated to address the following issues: We identify the explanatory variables that could be helpful on the basis of a critical review of the scientific literature. We study the relationship between ground observations and explanatory variables using a systematic procedure for a complete multiple regression analysis. Multiple regression models are calibrated combining all suggested model structures and explanatory variables. We also propose an evaluation of the models (using indices to analyze the goodness of fit) and select the best approaches (models and variables) on the basis of multicriteria analysis. Estimation of the snow depth is performed with the selected regression models. The residual estimation is improved by applying kriging in cases with spatial correlation. The final estimate is obtained by combining regression and kriging results, and constraining the snow domain in accordance with satellite data. The method is illustrated using the case study of the Sierra Nevada mountain range (Southern Spain). A cross‐validation experiment has confirmed the efficiency of the proposed procedure. Finally, although it is not the scope of this work, the snow depth is used to asses a first estimation of snow water equivalent resources.     

基于垂直剖面函数的天顶对流层延迟插值算法

通过分析由ERA-Interim气象再分析资料积分方法得到的天顶对流层总延迟随高程变化的规律,提出一种基于垂直剖面函数的天顶对流层延迟(ZTD)插值算法。该算法以ZTD的垂直分布规律为基础,通过垂直剖面函数实现ZTD在高程方向上的精准投影延拓,可以避免因高差较大造成的空间内插结构畸形。采用IGS站提供的高精度对流层产品进行实验验证表明,该算法相对于传统算法能够有效提高ZTD改正值的精度,尤其在高差超过1 km的情况下,相对于反距离加权法精度提升了96%,相对于空间回归法精度提升了79%。     

Inference and uncertainty of snow depth spatial distribution at the kilometre scale in the Colorado Rocky Mountains: the effects of sample size,random sampling,predictor quality,and validation procedures

Historically, observing snow depth over large areas has been difficult. When snow depth observations are sparse, regression models can be used to infer the snow depth over a given area. Data sparsity has also left many important questions about such inference unexamined. Improved inference, or estimation, of snow depth and its spatial distribution from a given set of observations can benefit a wide range of applications from water resource management, to ecological studies, to validation of satellite estimates of snow pack. The development of Light Detection and Ranging (LiDAR) technology has provided non‐sparse snow depth measurements, which we use in this study, to address fundamental questions about snow depth inference using both sparse and non‐sparse observations. For example, when are more data needed and when are data redundant? Results apply to both traditional and manual snow depth measurements and to LiDAR observations. Through sampling experiments on high‐resolution LiDAR snow depth observations at six separate 1.17‐km² sites in the Colorado Rocky Mountains, we provide novel perspectives on a variety of issues affecting the regression estimation of snow depth from sparse observations. We measure the effects of observation count, random selection of observations, quality of predictor variables, and cross‐validation procedures using three skill metrics: percent error in total snow volume, root mean squared error (RMSE), and R². Extremes of predictor quality are used to understand the range of its effect; how do predictors downloaded from internet perform against more accurate predictors measured by LiDAR? Whereas cross validation remains the only option for validating inference from sparse observations, in our experiments, the full set of LiDAR‐measured snow depths can be considered the ‘true’ spatial distribution and used to understand cross‐validation bias at the spatial scale of inference. We model at the 30‐m resolution of readily available predictors, which is a popular spatial resolution in the literature. Three regression models are also compared, and we briefly examine how sampling design affects model skill. Results quantify the primary dependence of each skill metric on observation count that ranges over three orders of magnitude, doubling at each step from 25 up to 3200. Whereas uncertainty (resulting from random selection of observations) in percent error of true total snow volume is typically well constrained by 100–200 observations, there is considerable uncertainty in the inferred spatial distribution (R²) even at medium observation counts (200–800). We show that percent error in total snow volume is not sensitive to predictor quality, although RMSE and R² (measures of spatial distribution) often depend critically on it. Inaccuracies of downloaded predictors (most often the vegetation predictors) can easily require a quadrupling of observation count to match RMSE and R² scores obtained by LiDAR‐measured predictors. Under cross validation, the RMSE and R² skill measures are consistently biased towards poorer results than their true validations. This is primarily a result of greater variance at the spatial scales of point observations used for cross validation than at the 30‐m resolution of the model. The magnitude of this bias depends on individual site characteristics, observation count (for our experimental design), and sampling design. Sampling designs that maximize independent information maximize cross‐validation bias but also maximize true R². The bagging tree model is found to generally outperform the other regression models in the study on several criteria. Finally, we discuss and recommend use of LiDAR in conjunction with regression modelling to advance understanding of snow depth spatial distribution at spatial scales of thousands of square kilometres. Copyright © 2012 John Wiley & Sons, Ltd.     

Assessing regression‐based statistical approaches for downscaling precipitation over North America

This paper assesses linear regression‐based methods in downscaling daily precipitation from the general circulation model (GCM) scale to a regional climate model (RCM) scale (45‐ and 15‐km grids) and down to a station scale across North America. Traditional downscaling experiments (linking reanalysis/dynamical model predictors to station precipitation) as well as nontraditional experiments such as predicting dynamic model precipitation from larger‐scale dynamic model predictors or downscaling dynamic model precipitation from predictors at the same scale are conducted. The latter experiments were performed to address predictability limit and scale issues. The results showed that the downscaling of daily precipitation occurrence was rarely successful at all scales, although results did constantly improve with the increased resolution of climate models. The explained variances for downscaled precipitation amounts at the station scales were low, and they became progressively better when using predictors from a higher‐resolution climate model, thus showing a clear advantage in using predictors from RCMs driven by reanalysis at its boundaries, instead of directly using reanalysis data. The low percentage of explained variances resulted in considerable underestimation of daily precipitation mean and standard deviation. Although downscaling GCM precipitation from GCM predictors (or RCM precipitation from RCM predictors) cannot really be considered downscaling, as there is no change in scale, the exercise yields interesting information as to the limit in predictive ability at the station scale. This was especially clear at the GCM scale, where the inability of downscaling GCM precipitation from GCM predictors demonstrates that GCM precipitation‐generating processes are largely at the subgrid scale (especially so for convective events), thus indicating that downscaling precipitation at the station scale from GCM scale is unlikely to be successful. Although results became better at the RCM scale, the results indicate that, overall, regression‐based approaches did not perform well in downscaling precipitation over North America. Copyright © 2013 John Wiley & Sons, Ltd.     

网络大数据下的中国现代食甜习惯空间分布特征及其影响因素研究

饮食地理文化作为地域文化中最具地方特色的重要元素,在现代人口大规模流动背景下呈现出全新的多样化局面,而基于传统认知的“南甜北咸”的地域分异已然不能代表中国现代食甜分布的空间特征。因此,本文采用网络爬虫技术,获取我国大陆31个省会城市共计约2000万条美食消费数据,从传统类菜品、主食类菜品、饮料类和甜品类菜品4个方面计算城市食甜度,在ArcGIS、MySQL软件支持下,借助GIS空间分析和数理统计方法探究我国现代食甜习惯的空间分布特征,分析影响食甜分布的因素。研究发现：① 中国食甜在空间分布上存在显著的地域分异特征,聚类分析评价参数R ²高达0.88,现代食甜习惯总体呈现“东高北中,西微内低”的包围式格局;② 从整体抑或局部角度,在1%显著性水平上莫兰指数均为正,中国食甜分布呈现显著的空间正相关关系,形成特色鲜明的3个地理集聚区,即以苏浙沪闽为主的东南沿海高甜集聚区,以渝黔川为主的西南内陆低甜集聚区和以陕宁为主的西北内陆低甜集聚区;③ 构建了中国现代食甜习惯分布影响因素模型,其拟合精度为0.82,分析结果显示降水、湿度、气温等气象要素及地理位置是影响现代我国食甜空间分布的重要因素。     

A smoothed cap model for variably saturated soils and its robust numerical implementation

This paper deals with the numerical implementation of a cap model for unsaturated soils. It provides a brief review of existing cap model approaches, based on which an improved model formulated in terms of generalised effective stress and matric suction is derived and described in detail. Although the proposed model is a multisurface plasticity model, it can efficiently be implemented using only single‐surface projections because of the smoothness of the model, which is obtained by construction. Numerical algorithms are provided for these single‐surface stress projections, using a single‐equation approach whenever possible. The robustness of the utilised single‐equation approaches is enhanced by proposing problem‐fitted start‐up procedures based on investigations of the nonlinear projection equations. A comparison of the model response with extensive material test data is used to validate the model and to demonstrate the robust application of the approach to silty sands and low to medium plasticity clays. Copyright © 2015 John Wiley & Sons, Ltd.     

Calibration of paired watersheds: Utility of moving sums in presence of externalities

Historically, paired watershed studies have been used to quantify the hydrological effects of land use and management practices by concurrently monitoring 2 similar watersheds during calibration (pretreatment) and post‐treatment periods. This study characterizes seasonal water table and flow response to rainfall during the calibration period and tests a change detection technique of moving sums of recursive residuals (MOSUM) to select calibration periods for each control–treatment watershed pair when the regression coefficients for daily water table elevation were most stable to minimize regression model uncertainty. The control and treatment watersheds were 1 watershed of 3–4‐year‐old intensely managed loblolly pine (Pinus taeda L.) with natural understory, 1 watershed of 3–4‐year‐old loblolly pine intercropped with switchgrass (Panicum virgatum), 1 watershed of 14–15‐year‐old thinned loblolly pine with natural understory (control), and 1 watershed of switchgrass only. The study period spanned from 2009 to 2012. Silvicultural operational practices during this period acted as external factors, potentially shifting hydrologic calibration relationships between control and treatment watersheds. MOSUM results indicated significant changes in regression parameters due to silvicultural operations and were used to identify stable relationships for water table elevation. None of the calibration relationships developed using this method were significantly different from the classical calibration relationship based on published historical data. We attribute that to the similarity of historical and 2010–2012 leaf area index on control and treatment watersheds as moderated by the emergent vegetation. Although the MOSUM approach does not eliminate the need for true calibration data or replace the classic paired watershed approach, our results show that it may be an effective alternative approach when true data are unavailable, as it minimizes the impacts of external disturbances other than the treatment of interest.     

Concentration–discharge relationships in the coal mined region of the New River basin and Indian Fork sub‐basin,Tennessee, USA

For many basins, identifying changes to water quality over time and understanding current hydrologic processes are hindered by fragmented and discontinuous water‐quality and hydrology data. In the coal mined region of the New River basin and Indian Fork sub‐basin, muted and pronounced changes, respectively, to concentration–discharge (C–Q) relationships were identified using linear regression on log‐transformed historical (1970s–1980s) and recent (2000s) water‐quality and streamflow data. Changes to C–Q relationships were related to coal mining histories and shifts in land use. Hysteresis plots of individual storms from 2007 (New River) and the fall of 2009 (Indian Fork) were used to understand current hydrologic processes in the basins. In the New River, storm magnitude was found to be closely related to the reversal of loop rotation in hysteresis plots; a peak‐flow threshold of 25 cubic meters per second (m³/s) segregates hysteresis patterns into clockwise and counterclockwise rotational groups. Small storms with peak flow less than 25 m³/s often resulted in dilution of constituent concentrations in headwater tributaries like Indian Fork and concentration of constituents downstream in the mainstem of the New River. Conceptual two or three component mixing models for the basins were used to infer the influence of water derived from spoil material on water quality. Copyright © 2012 John Wiley & Sons, Ltd.     

Inferring air–water temperature relationships from river and catchment properties

Water temperature (Tw) is a key determinant of freshwater ecosystem status and cause for concern under a changing climate. Hence, there is growing interest in the feasibility of moderating rising Tw through management of riparian shade. The Loughborough University Temperature Network (LUTEN) is an array of 36 water and air temperature (Ta) monitoring sites in the English Peak District set‐up to explore the predictability of local Tw, given Ta, river reach, and catchment properties. Year 1 of monitoring shows that 84%–94% of variance in daily Tw is explained by Ta. However, site‐specific logistic regression parameters exhibit marked variation and dependency on upstream riparian shade. Perennial spring flows in the lower River Dove also affect regression model parameters and strongly buffer daily and seasonal mean Tw. The asymptote of the models (i.e. maximum expected Tw) is particularly sensitive to groundwater inputs. We conclude that reaches with spring flows potentially offer important thermal refuges for aquatic organisms against expected long‐term warming of rivers and should be afforded special protection. Copyright © 2013 John Wiley & Sons, Ltd.     

基于Google Earth土地承包经营权工作底图的制作

对利用Google Earth影像制作农村土地承包经营权工作底图的方法和技术流程进行了全面探讨。首先在CASS 7.0软件下对集体土地所有权行政界线进行坐标转换,再使用自编程序提取界线拐点的WGS-84大地坐标,通过"地图下载器"下载Google Earth影像并拼接输出TIF影像,利用Arc Map软件的投影变换功能,将墨卡托投影的TIF影像转换为高斯投影影像,最后,在CASS 7.0软件下插入变换后的正射影像,叠加行政界线,形成完整的承包经营权外业调查工作底图。