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.     

Varying water utilization of Haloxylon ammodendron plantations in a desert‐oasis ecotone

Haloxylon ammodendron is a desert shrub used extensively in China for restoring degraded dry lands. An understanding of the water source used by H. ammodendron plantations is critical achieving sustainable vegetation restoration. We measured mortality, shoot size, and rooting depth in 5‐, 10‐, 20‐, and 40‐year‐old H. ammodendron plantations. We examined stable isotopic ratios of oxygen (δ¹⁸O) in precipitation, groundwater, and soil water in different soil layers and seasons, and in plant stem water to determine water sources at different shrub ages. We found that water acquisition patterns in H. ammodendron plantations differed with plantation age and season. Thus, the main water source for 5‐year‐old shrubs was shallow soil water. Water sources of 10‐year‐old shrubs shifted depending on the soil water conditions during the season. Although their tap roots could absorb deep soil water, the plantation main water sources were from soil water, and about 50% of water originated from shallow and mid soil. This pattern might occur because main water sources in these plantations were changeable over time. The 20‐ and 40‐year‐old shrubs acquired water mainly from permanent groundwater. We conclude that the main water source of a young H. ammodendron plantation was soil water recharged by precipitation. However, when roots reached sufficient depth, water originated mainly from the deep soil water, especially in the dry season. The deeply rooted 20‐ and 40‐year‐old shrubs have the ability to exploit a deep and reliable water source. To achieve sustainability in these plantations, we recommend a reduction in the initial density of H. ammodendron in the desert‐oasis ecotone to decelerate the consumption of shallow soil water during plantation establishment.     

Quantitative estimation of submarine groundwater discharge using airborne thermal infrared data acquired at two different tidal heights

Submarine groundwater discharge (SGD) plays an important role in coastal biogeochemical processes and hydrological cycles, particularly off volcanic islands in oligotrophic oceans. However, the spatial and temporal variations of SGD are still poorly understood owing to difficulty in taking rapid SGD measurements over a large scale. In this study, we used four airborne thermal infrared surveys (twice each during high and low tides) to quantify the spatiotemporal variations of SGD over the entire coast of Jeju Island, Korea. On the basis of an analytical model, we found a linear positive correlation between the thermal anomaly and squares of the groundwater discharge velocity and a negative exponential correlation between the anomaly and water depth (including tide height and bathymetry). We then derived a new equation for quantitatively estimating the SGD flow rates from thermal anomalies acquired at two different tide heights. The proposed method was validated with the measured SGD flow rates using a current meter at Gongcheonpo Beach. We believe that the method can be effectively applied for rapid estimation of SGD over coastal areas, where fresh groundwater discharge is significant, using airborne thermal infrared surveys.     

The exponential decline in saturated hydraulic conductivity with depth: a novel method for exploring its effect on water flow paths and transit time distribution

The strong vertical gradient in soil and subsoil saturated hydraulic conductivity is characteristic feature of the hydrology of catchments. Despite the potential importance of these strong gradients, they have proven difficult to model using robust physically based schemes. This has hampered the testing of hypotheses about the implications of such vertical gradients for subsurface flow paths, residence times and transit time distribution. Here we present a general semi‐analytical solution for the simulation of 2D steady‐state saturated‐unsaturated flow in hillslopes with saturated hydraulic conductivity that declines exponentially with depth. The grid‐free solution satisfies mass balance exactly over the entire saturated and unsaturated zones. The new method provides continuous solutions for head, flow and velocity in both saturated and unsaturated zones without any interpolation process as is common in discrete numerical schemes. This solution efficiently generates flow pathlines and transit time distributions in hillslopes with the assumption of depth‐varying saturated hydraulic conductivity. The model outputs reveal the pronounced effect that changing the strength of the exponential decline in saturated hydraulic conductivity has on the flow pathlines, residence time and transit time distribution. This new steady‐state model may be useful to others for posing hypotheses about how different depth functions for hydraulic conductivity influence catchment hydrological response. Copyright © 2016 John Wiley & Sons, Ltd.     

Revised runoff curve number for runoff prediction in the Loess Plateau of China

The soil conservation service (now Natural Resources Conservation Service) Curve Number (SCS-CN), one of the most commonly used methods for surface runoff prediction. The runoff calculated by this method was very sensitive to CN values. In this study, CN values were calculated by both arithmetic mean (CN_C) and least square fit method (CN_F) using observed rainfall-runoff data from 43 sites in the Loess Plateau region, which are considerably different from the CN₂ values obtained from the USDA-SCS handbook table (CN_T). The results showed that using CN_C instead of CN_T for each watershed produce little improvement, while replacing CN_T with CN_F improves the performance of the original SCS-CN method, but still performs poorly in most study sites. This is mainly due to the SCS-CN method using a constant CN value and discounting of the temporal variation in rainfall-runoff process. Therefore, three factors—soil moisture, rainfall depth and intensity—affecting the surface runoff variability are considered to reflect the variation of CN in each watershed, and a new CN value was developed. The reliability of the proposed method was tested with data from 38 watersheds, and then applied to the remaining five typical watersheds using the optimized parameters. The results indicated that the proposed method, which boosted the model efficiencies to 81.83% and 74.23% during calibration and validation cases, respectively, performed better than the original SCS-CN and the Shi and Wang (2020b) method, a modified SCS-CN method based on tabulated CN value. Thus, the proposed method incorporating the influence of the temporal variability of soil moisture, rainfall depth, and intensity factors suggests an accurate runoff prediction for general applications under different hydrological and climatic conditions on the Loess Plateau region.     

An examination of snow redistribution over smooth land‐fast sea ice

An understanding of temporal evolution of snow on sea ice at different spatial scales is essential for improvement of snow parameterization in sea ice models. One of the problems we face, however, is that long‐term climate data are routinely available for land and not for sea ice. In this paper, we examine the temporal evolution of snow over smooth land‐fast first‐year sea ice using observational and modelled data. Changes in probability density functions indicate that depositional and drifting events control the evolution of snow distribution. Geostatistical analysis suggests that snowdrifts increased over the study period, and the orientation was related to the meteorological conditions. At the microscale, the temporal evolution of the snowdrifts was a product of infilling in the valleys between drifts. Results using two shore‐based climate reporting stations (Paulatuk and Tuktoyuktuk, NWT) suggest that on‐ice air temperature and relative humidity can be estimated using air temperature recorded at either station. Wind speed, direction and precipitation on ice cannot be accurately estimated using meteorological data from either station. The temporal evolution of snow distribution over smooth land‐fast sea ice was modelled using SnowModel and four different forcing regimes. The results from these model runs indicate a lack of agreement between observed distribution and model outputs. The reasons for these results are lack of meteorological measurements prior to the end of January, lack of spatially adequate surface topography and discrepancies between meteorological variables on land and ice. Copyright © 2009 John Wiley & Sons, Ltd.     

面向WebGL的矢量数据三维绘制技术

矢量数据的叠加显示能够提高三维虚拟地球的表达效果与分析能力。受限于GPU的计算精度,在三维虚拟地球中矢量数据绘制普遍存在抖动现象和深度冲突现象。对基于WebGL的矢量数据三维绘制中计算精度问题进行了分析,提出了使用CPU RTC技术和GPU RTE技术提高顶点变换的精度,使用多视锥渲染算法和深度平面技术解决深度缓存精度问题。实验证明,这几种技术和算法可以有效缓解抖动现象和深度冲突现象导致的视觉干扰,改善了各种尺度和范围的矢量数据在三维地形上的叠加显示效果。     

基于FY3A/MERSI资料分析广东省气溶胶光学厚度分布

利用国产极轨气象卫星FY3A的MERSI AOD产品分析2010—2013年广东省气溶胶光学厚度的分布规律。结果表明:MERSI AOD产品与地面太阳光度计实测数据的相关系数为0.72,其平均绝对值误差为0.12,均方根误差为0.15,数据精度可满足研究需要;从AOD的空间分布看,珠三角西翼东翼山区五市,其中佛山市、东莞市、中山市为广东省AOD均值最高的地区,梅州市、河源市为广东省AOD均值最低的地区;从AOD的时间分布看,2010—2013年间,AOD呈现先升高后降低的趋势,2011年为拐点,与此同时,AOD还表现出明显的季节变化特征,春季为AOD高值期,夏季、秋季次之,冬季最低。     

太湖地区HJ-1卫星CCD数据反演气溶胶及大气校正

针对HJ-1卫星CCD数据,利用改进的暗像元法反演气溶胶光学厚度(AOD),再利用反演的AOD对其进行大气校正。将反演的气溶胶与地基太阳光度计数据进行对比验证,发现当反演的AOD值大于0.2时,反演值与地基观测值的相关系数为0.964,符合MODIS业务化反演AOD的精度要求。再将反演得到的气溶胶带入6S辐射传输模型中,对HJ-1卫星CCD数据进行大气校正实验。结果表明,该方法能有效提高HJ-1卫星CCD数据大气校正的精度,更好地复原地物的真实光谱信息。     

Stabilisation/solidification and bioaugmentation treatment of petroleum drill cuttings

Petroleum drill cuttings are usually treated by techniques suitable for particular contaminant groups. The significance of this study consists in the development of a treatment technology that can simultaneously handle the hydrocarbon and metal constituents of drill cuttings. Bioaugmentation is combined with stabilisation/solidification (S/S), within S/S monoliths and in granulated S/S monoliths. Portland cement was used for S/S treatment at 30% binder dosage. Bioaugmentation treatment involved two bacterial densities of a mixed culture bio-preparation. The effects of inclusion of compost, fertiliser and activated carbon were also evaluated. After 28 days, the combined S/S and bioaugmentation treatments recorded up to 15% higher total petroleum hydrocarbon (TPH) loss than control S/S treatment without bioaugmentation. Embedding fertiliser, activated carbon and higher bacterial density within S/S monoliths resulted in the highest (99%) TPH reduction but higher concentrations of metals. The addition of compost and lower bacterial density to granulated S/S monoliths led to similar (98%) TPH degradation and lower amounts of metals. The results suggest that with better mixture optimisation, combining S/S and bioaugmentation could engender more sustainable treatment of drill cuttings.