Generation of continuous surface soil moisture dataset using combined optical and thermal infrared images

Surface soil moisture (SSM) is a critical variable for understanding water and energy flux between the atmosphere and the Earth's surface. An easy to apply algorithm for deriving SSM time series that primarily uses temporal parameters derived from simulated and in situ datasets has recently been reported. This algorithm must be assessed for different biophysical and atmospheric conditions by using actual geostationary satellite images. In this study, two currently available coarse‐scale SSM datasets (microwave and reanalysis product) and aggregated in situ SSM measurements were implemented to calibrate the time‐invariable coefficients of the SSM retrieval algorithm for conditions in which conventional observations are rare. These coefficients were subsequently used to obtain SSM time series directly from Meteosat Second Generation (MSG) images over the study area of a well‐organized soil moisture network named REMEDHUS in Spain. The results show a high degree of consistency between the estimated and actual SSM time series values when using the three SSM dataset‐calibrated time‐invariable coefficients to retrieve SSM, with coefficients of determination (R²) varying from 0.304 to 0.534 and root mean square errors ranging from 0.020 m³/m³ to 0.029 m³/m³. Further evaluation with different land use types results in acceptable debiased root mean square errors between 0.021 m³/m³ and 0.048 m³/m³ when comparing the estimated MSG pixel‐scale SSM with in situ measurements. These results indicate that the investigated method is practical for deriving time‐invariable coefficients when using publicly accessed coarse‐scale SSM datasets, which is beneficial for generating continuous SSM dataset at the MSG pixel scale.     

Temporal stability of soil moisture in an oasis of northwestern China

Temporal stability of soil moisture has been widely used in hydrological monitoring since it emerged. However, the spatial analysis of temporal stability at the landscape scale is often limited because of insufficient sampling numbers. This work made an effort to investigate the spatial variations of temporal stability of soil moisture in an oasis landscape. The specific objectives of the study were to explore the spatial patterns of temporal stability and to determine the controlling factors of temporal stability in the desert oasis. A time series of soil moisture measurements were gathered on 23 occasions at 118 locations over 3 years in a rectangular transect of approximately 100 km². The nonparametric Spearman's rank correlation coefficient, standard deviation of relative difference (SDRD), and mean absolute bias error (MABE) were used to quantify the temporal stability of soil moisture. Results showed that the temporal stability of soil moisture was depth dependent and season dependent. The spatial pattern of soil moisture in a deep soil layer and between two same seasons generally had a high temporal stability. SDRD and MABE were spatially autocorrelated and exhibited strong spatial structures in the geographic space. The concept of temporal stability can be extended to describe the time‐stable areas of soil moisture with geostatistics. There were great differences between SDRD and MABE in describing the temporal stability of soil moisture and in identifying the controlling factors of temporal stability. In this case, MABE was a better alternative to estimate the areal mean soil moisture using representative locations than SDRD. Land use type, soil moisture condition, and soil particle composition were the dominant controls of temporal stability in the oasis. These insights could help to better understand the essence of temporal stability of soil moisture in arid regions.     

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.     

Tropical precipitation anomalies and d‐excess evolution during El Niño 2014‐16

The last 2014‐16 El Niño event was among the three strongest episodes on record. El Niño considerably changes annual and seasonal precipitation across the tropics. Here, we present a unique stable isotope data set of daily precipitation collected in Costa Rica prior to, during, and after El Niño 2014‐16, in combination with Lagrangian moisture source and precipitation anomaly diagnostics. δ²H composition ranged from ‐129.4 to +18.1 (‰) while δ¹⁸O ranged from ‐17.3 to +1.0 (‰). No significant difference was observed among δ¹⁸O (P=0.186) and δ²H (P=0.664) mean annual compositions. However, mean annual d‐excess showed a significant decreasing trend (from +13.3 to +8.7 ‰) (P<0.001) with values ranging from +26.6 to ‐13.9 ‰ prior to and during the El Niño evolution. The latter decrease in d‐excess can be partly explained by an enhanced moisture flux convergence across the southeastern Caribbean Sea coupled with moisture transport from northern South America by means of an increased Caribbean Low Level Jet regime. During 2014‐15, precipitation deficit across the Pacific domain averaged 46% resulting in a very severe drought; while a 94% precipitation surplus was observed in the Caribbean domain. Understanding these regional moisture transport mechanisms during a strong El Niño event may contribute to a) better understanding of precipitation anomalies in the tropics and b) re‐evaluate past stable isotope interpretations of ENSO events in paleoclimatic archives within the Central America region.     

Monitoring over‐winter soil water dynamics in a freezing and snow‐covered environment using a thermally insulated tensiometer

Water potential below a frozen soil layer was continuously monitored over an entire winter period (using thermally insulated tensiometers sheltered in a heated chamber) along with other soil, snow and atmospheric variables. In early winter, the freezing front advanced under a thin snow cover, inducing upward soil water flow in the underlying unfrozen soil. The freezing front started to retreat when the snow cover became thick enough to insulate the soil, resulting in the reversal of the flow direction in the unfrozen zone. These data provide a clear illustration of soil water dynamics, which have rarely been monitored with a tensiometer. Copyright © 2005 John Wiley & Sons, Ltd.     

基于温度植被干旱指数的电网输电线路走廊干旱分析

以湖北省输电线路走廊地区作为研究区,利用2013年1~9月MODIS卫星影像数据,处理得到月尺度的归一化植被指数(Normalized Differential Vegetation Index,NDVI)与地表温度(Land Surface Temperature,LST)数据,构建NDVI-Ts特征空间,计算得到温度植被干旱指数(Temperature Vegetation Dryness Index,TVDI),用TVDI监测结果分析湖北省输电线路走廊区域2013年干旱时空分布情况。结果表明,湖北省输电线路走廊地区TVDI和土壤含水量之间存在显著的负相关,相关系数达到0.525(p0.05),由MODIS卫星影像计算得到TVDI影像可以有效表明湖北省输电线路走廊地区的土壤含水情况。     

Variability in bedform morphology and kinematics with transport stage

Bedform geometry is widely recognized to be a function of transport stage. Bedform aspect ratio (height/length) increases with transport stage, reaches a maximum, then decreases as bedforms washout to a plane bed. Bedform migration rates are also linked to bedform geometry, in so far as smaller bedforms in coarser sediment tend to migrate faster than larger bedforms in finer sediment. However, how bedform morphology (height, length and shape) and kinematics (translation and deformation) change with transport stage and suspension have not been examined. A series of experiments is presented where initial flow depth and grain size were held constant and the transport stage was varied to produce bedload dominated, mixed‐load dominated and suspended‐load dominated conditions. The results show that the commonly observed pattern in bedform aspect ratio occurs because bedform height increases then decreases with transport stage, against a continuously increasing bedform length. Bedform size variability increased with transport stage, leading to less uniform bedform fields at higher transport stage. Total translation‐related and deformation‐related sediment fluxes all increased with transport stage. However, the relative contribution to the total flux changed. At the bedload dominated stage, translation‐related and deformation‐related flux contributed equally to the total flux. As the transport stage increased, the fraction of the total load contributed by translation increased and the fraction contributed by deformation declined because the bedforms got bigger and moved faster. At the suspended‐load dominated transport stage, the deformation flux increased and the translation flux decreased as a fraction of the total load, approaching one and zero, respectively, as bedforms washed out to a plane bed.     

Upper-Tropospheric Environment–Tropical Cyclone Interactions over the Western North Pacific: A Statistical Study

Based on 25-year(1987–2011) tropical cyclone(TC) best track data, a statistical study was carried out to investigate the basic features of upper-tropospheric TC–environment interactions over the western North Pacific. Interaction was defined as the absolute value of eddy momentum flux convergence(EFC) exceeding 10 m s~(-1)d~(-1). Based on this definition, it was found that 18% of all six-hourly TC samples experienced interaction. Extreme interaction cases showed that EFC can reach~120 m s~(-1)d~(-1) during the extratropical-cyclone(EC) stage, an order of magnitude larger than reported in previous studies.Composite analysis showed that positive interactions are characterized by a double-jet flow pattern, rather than the traditional trough pattern, because it is the jets that bring in large EFC from the upper-level environment to the TC center. The role of the outflow jet is also enhanced by relatively low inertial stability, as compared to the inflow jet. Among several environmental factors, it was found that extremely large EFC is usually accompanied by high inertial stability, low SST and strong vertical wind shear(VWS). Thus, the positive effect of EFC is cancelled by their negative effects. Only those samples during the EC stage, whose intensities were less dependent on VWS and the underlying SST, could survive in extremely large EFC environments, or even re-intensify. For classical TCs(not in the EC stage), it was found that environments with a moderate EFC value generally below ~25 m s~(-1)d~(-1) are more favorable for a TC's intensification than those with extremely large EFC.     

Effects of surface flux parameterization on the numerically simulated intensity and structure of Typhoon Morakot (2009)

The effects of surface flux parameterizations on tropical cyclone(TC) intensity and structure are investigated using the Advanced Research Weather Research and Forecasting(WRF-ARW) modeling system with high-resolution simulations of Typhoon Morakot(2009).Numerical experiments are designed to simulate Typhoon Morakot(2009) with different formulations of surface exchange coefficients for enthalpy(C_K) and momentum(C_D) transfers,including those from recent observational studies based on in situ aircraft data collected in Atlantic hurricanes.The results show that the simulated intensity and structure are sensitive to C_K and C_D,but the simulated track is not.Consistent with previous studies,the simulated storm intensity is found to be more sensitive to the ratio of C_K/C_D than to C_K or C_D alone.The pressure-wind relationship is also found to be influenced by the exchange coefficients,consistent with recent numerical studies.This paper emphasizes the importance of C_D and C_K on TC structure simulations.The results suggest that C_D and C_K have a large impact on surface wind and flux distributions,boundary layer heights,the warm core,and precipitation.Compared to available observations,the experiment with observed C_D and C_K generally simulated better intensity and structure than the other experiments,especially over the ocean.The reasons for the structural differences among the experiments with different C_D and C_K setups are discussed in the context of TC dynamics and thermodynamics.     

Long-Term Variability of the Wind Field over the Indian Ocean Based on ERA-Interim Reanalysis

A detailed study of long-term variability of winds using 30 years of data from the European Centre for Medium-range Weather Forecasts global reanalysis (ERA-Interim) over the Indian Ocean has been carried out by partitioning the Indian Ocean into six zones based on local wind extrema. The trend of mean annual wind speed averaged over each zone shows a significant increase in the equatorial region, the Southern Ocean, and the southern part of the trade winds. This indicates that the Southern Ocean winds and the southeast trade winds are becoming stronger. However, the trend for the Bay of Bengal is negative, which might be caused by a weakening of the monsoon winds and northeast trade winds. Maximum interannual variability occurs in the Arabian Sea due to monsoon activity; a minimum is observed in the subtropical region because of the divergence of winds. Wind speed variations in all zones are weakly correlated with the Dipole Mode Index (DMI). However, the equatorial Indian Ocean, the southern part of the trade winds, and subtropical zones show a relatively strong positive correlation with the Southern Oscillation Index (SOI), indicating that the SOI has a zonal influence on wind speed in the Indian Ocean. Monsoon winds have a decreasing trend in the northern Indian Ocean, indicating monsoon weakening, and an increasing trend in the equatorial region because of enhancement of the westerlies. The negative trend observed during the non-monsoon period could be a result of weakening of the northeast trade winds over the past few decades. The mean flux of kinetic energy of wind (FKEW) reaches a minimum of about 100?W?m^?2 in the equatorial region and a maximum of about 1500?W?m^?2 in the Southern Ocean. The seasonal variability of FKEW is large, about 1600?W?m^?2, along the coast of Somalia in the northern Indian Ocean. The maximum monthly variability of the FKEW field averaged over each zone occurs during boreal summer. During the onset and withdrawal of monsoon, FKEW is as low as 50?W?m^?2. The Southern Ocean has a large variation of about 1280?W?m^?2 because of strong westerlies throughout the year.