Spatial and temporal dynamics of soil moisture after rainfall events along a slope in Regosols of southwest China

High frequency seasonal drought in purple soils (Regosols in FAO taxonomy) of the hilly upland areas of Sichuan basin, China, is one of the key restrictive factors for crop production. In order to manage irrigation and fertilizer application in these soils effectively, the soil water content in a sloped plot with 60 cm soil depth was measured by neutron probe devices to investigate the soil moisture regime during the 1998 rainy season after various amounts of rainfall events. The results showed that variation of soil moisture along the slope positions was highest in the top soil layer during the period of sporadic rainfall that did not induce any runoff. The coefficients of variation of soil moisture at various slope positions (upper, middle, and lower) are 17·36%, 8·95%, 10·25%, 8·58%, 8·05% and 9·21% at the 10 cm, 20 cm, 30 cm, 40 cm, 50 cm and 60 cm soil depths respectively. When surface runoff occurred, the soil moisture dynamics at various positions on the plot were then very different. Soil water content decreased more rapidly on the upper slope than on the middle and lower slope positions. When both surface runoff and throughflow occurred, the soil moisture dynamics in the various layers showed a stable period (soil water content is near constant as time elapses) that lasted about 1 week. Also, the pattern of moisture dynamics is ‘decreasing–stabilization–decreasing’. Thus, irrigation and fertilization management according to the spatial and temporal features of soil moisture dynamics on sloped land can increase the water and fertilizer utilization efficacy by reducing their losses during the stable period. Copyright © 2006 John Wiley & Sons, Ltd.     

Simulating sensitivity of soil moisture and evapotranspiration under heterogeneous soils and land uses

Soil moisture (SM) plays an important role in land surface and atmospheric interactions. It modifies energy balance at the surface and the rate of water cycling between the land and atmosphere. In this paper we provide a sensitivity assessment of SM and ET for heterogeneous soil physical properties and for three land uses including irrigated maize, rainfed maize, and grass at a climatological time-scale by using a water balance model. Not surprisingly, the study finds increased soil water content in the root zone throughout the year under irrigated farming. Soil water depletes to its lowest level under rainfed maize cultivation. We find a ‘land use’ effect as high as 36 percent of annual total evapotranspiration, under irrigated maize compared to rainfed maize and grass, respectively. Sensitivity analyses consisting of comparative simulations using the model show that soil characteristics, like water holding capacity, influence SM in the root zone and affect seasonal total ET estimates at the climatological time-scale. This ‘soils’ effect is smaller than the ‘land use’ effect associated with irrigation but, it is a source of consistent bias for both SM and ET estimates. The ‘climate’ effect basically masks the ‘soils’ effect under wet conditions. These results lead us to conclude that appropriate representation of land use, soils, and climate are necessary to accurately represent the water and energy balance in real landscapes.     

Measurement of soil moisture change in spatially heterogeneous weathered soils using a capacitance probe

A capacitance probe has been employed to measure the spatial and temporal change of soil moisture in a weathered soil profile with heterogeneous texture characteristics. Measurements were made to a depth of 2 m at intervals of 2 cm. The field work was carried out during a wetting-up period in a subtropical monsoonal environment in the Middle Hills of Nepal. Calibrations of the frequency readings against gravimetric and volumetric soil moisture are compared. It was found that the frequency reading correlates much better with the volumetric soil moisture. Although the soil is spatially heterogeneous in its particle size, porosity and bulk density, it was found that a single regression equation can be used to represent the relationship between the frequency reading of the capacitance probe and volumetric soil moisture. The results demonstrate that the capacitance probe enables the investigation of soil moisture variation temporally and spatially at high resolution and can measure a wide moisture range. © 1998 John Wiley & Sons, Ltd.     

Soil moisture dynamics in coastal savanna soils in the tropics under different soil management practices

Abstract

Soil moisture is important for crop cultivation and its adequacy to meet crop-water requirements is determined by the degree of soil management practised and the quantity of water applied to the soil. This study investigates soil moisture dynamics on three plots: Bare (clean, weeds removed), Weedy (kept weedy), and Mulched (cleared of weeds and fully covered with grass mulch) during rainy and dry periods at the Teaching and Research Farm at the University of Cape Coast, in the coastal savanna zone of Ghana. Soil moisture dynamics under different levels of soil compaction were also studied. A Massey Ferguson tractor (MF265) was used to compact the soil at various levels by making 0, 1, 5, 9 and 13 passes. During both the rainy and the dry periods, moisture retention in the soil under bare, weedy and mulched plots increased with depth. During the rainy period, the mean soil moisture retention was in the order: Mulched > Weedy > Bare at both 0–20 cm and 20–40 cm depths. Within a 7–day period, soil moisture measurements from a day after heavy rainfall (intensity > 7 mm h^?1) gave mean moisture losses of 2.7, 4.1 and 3.9% for the Bare, Weedy and Mulched plots, respectively. During the dry period, however, the mean soil moisture retention was of the order: Mulched > Bare > Weedy at both 0–20 cm and 20–40 cm depths. Mean moisture loss during a 7–day dry period was 4.5, 2.9 and 3.4% for the Bare, Weedy and Mulched plots, respectively. Under different levels of soil compaction, the mean moisture retention in the soil increased from 8.3% at 0 pass to 17.8% at 13 passes within the 0–20 cm depth, whilst it decreased from 13.3 to 5.9% from 0 to 13 passes, respectively, within the 20–40 cm depth. It was realized that at less than two passes, the mean soil moisture retention within the 0–20 cm depth was less than the mean moisture retention within the 20–40 cm depth, but the converse happened for more than two passes. The study showed that mulching the soil surface helped to retain enough soil moisture during both the rainy and the dry seasons. Also, soil with high sand content required some sort of soil compaction in order to retain enough moisture at the crop rooting zone.     

地统计学在淡水生态学中的应用

建立在区域化定量理论基础上的地统计学是一门空间统计学分支学科,主要用于研究自然现象的相关性和依赖性。这一理论主要内容包括：半方差图（用以描述研究对象空间相关性）;Ｋｒｉｇｅ空间内插技术（通过空间上抽样点的调查数据对空间上未测点进行估计）;以及通过半方差图求算分形体的分形维数,本文在介绍该理论基本原理的基础上,探讨了其在淡水生态学上的实际应用。     

Determination of bare surface soil moisture from combined temporal evolution of land surface temperature and net surface shortwave radiation

Land surface soil moisture (SSM) is an important variable for hydrological, ecological, and meteorological applications. A multi‐linear model has recently been proposed to determine the SSM content from the combined diurnal evolution of both land surface temperature (LST) and net surface shortwave radiation (NSSR) with the parameters TN (the LST mid‐morning rising rate divided by the NSSR rising rate during the same period) and t_d (the time of daily maximum temperature). However, in addition to the problem that all the coefficients of the multi‐linear model depend on the atmospheric conditions, the model also suffers from the problems of the nonlinearity of TN as a function of the SSM content and the uncertainty of determining the t_d from the diurnal evolution of the LST. To address these problems, a modified multi‐linear model was developed using the logarithm of TN and normalizing t_d by the mid‐morning temperature difference instead of using the TN and t_d. Except for the constant term, the coefficients of all other variables in the modified multi‐linear model proved to be independent of the atmospheric conditions. Using the relevant simulation data, results from the modified multi‐linear model show that the SSM content can be determined with a root mean square error (RMSE) of 0.030m³/m³, provided that the constant term is known or estimated day to day. The validation of the model was conducted using the field measurements at the Langfang site in 2008 in China. A higher correlation is achieved (coefficient of determination: R² = 0.624, RMSE = 0.107m³/m³) between the measured SSM content and the SSM content estimated using the modified multi‐linear model with the coefficients determined from the simulation data. Another experiment is also conducted to estimate the SSM content using the modified model with the constant term calibrated each day by one‐spot measurements at the site. The estimation result has a relatively larger error (RMSE = 0.125m³/m³). Additionally, the uncertainty of the determination of the coefficients is analysed using the field measurements, and the results indicate that the SSM content obtained using the modified model accurately characterizes the surface soil moisture condition. Copyright © 2012 John Wiley & Sons, Ltd.     

Soil water content in southern England derived from a cosmic‐ray soil moisture observing system – COSMOS‐UK

Cosmic‐ray soil moisture sensors have the advantage of a large measurement footprint (approximately 700 m in diameter) and are able to operate continuously to provide area‐averaged near‐surface (top 10–20 cm) volumetric soil moisture content at the field scale. This paper presents the application of this technique at four sites in southern England over almost 3 years. Results show the soil moisture response to contrasting climatic conditions during 2011–2014 and are the first such field‐scale measurements made in the UK. These four sites are prototype stations for a UK COsmic‐ray Soil Moisture Observing System, and particular consideration is given to sensor operating conditions in the UK. Comparison of these soil water content observations with the Joint UK Land Environment Simulator 10‐cm soil moisture layer shows that these data can be used to test and diagnose model performance and indicate the potential for assimilation of these data into hydro‐meteorological models. The application of these large‐area soil water content measurements to evaluate remotely sensed soil moisture products is also demonstrated. Numerous applications and the future development of a national COsmic‐ray Soil Moisture Observing System network are discussed. Copyright © 2016 John Wiley & Sons, Ltd.     

基于土壤墒情抗差的降雨径流预报模型——设计与实验流域应用

本文旨在将实时监测得到的土壤墒情转化为流域水文模型可以直接使用的土壤含水量,论证将实时土壤墒情资料用于实时预报的可行性;利用实时监测土壤墒情,改进传统的模型结构,设计基于实测土壤墒情的降雨径流水文预报模型.采用土壤含水量误差抗差估计技术以抵御观测资料粗差的影响,提高系统的稳定性;并在此基础上提出了土壤含水量系统响应修正方法,以提高模型计算精度.将该模型应用于实验流域——宝盖洞流域进行应用检验,洪水模拟合格率达到92.3%,整体模拟精度达到甲级.     

Potential of soil moisture observations in flood modelling: Estimating initial conditions and correcting rainfall

Rainfall runoff (RR) models are fundamental tools for reducing flood hazards. Although several studies have highlighted the potential of soil moisture (SM) observations to improve flood modelling, much research has still to be done for fully exploiting the evident connection between SM and runoff. As a way of example, improving the quality of forcing data, i.e. rainfall observations, may have a great benefit in flood simulation. Such data are the main hydrological forcing of classical RR models but may suffer from poor quality and record interruption issues. This study explores the potential of using SM observations to improve rainfall observations and set a reliable initial wetness condition of the catchment for improving the capability in flood modelling. In particular, a RR model, which incorporates SM for its initialization, and an algorithm for rainfall estimation from SM observations are coupled using a simple integration method. The study carried out at the Valescure experimental catchment (France) demonstrates the high information content retained by SM for RR transformation, thus giving new possibilities for improving hydrological applications. Results show that an appropriate configuration of the two models allows obtaining improvement in flood simulation up to 15% in mean and 34% in median Nash Sutcliffe performances as well as a reduction of the median error in volume and on peak discharge of about 30% and 15%, respectively.     

Study of water tension differences in heterogeneous sandy soils using surface ERT

Herbaceous vegetation in the Sahel grows almost exclusively on sandy soils which preferentially retain water through infiltration and storage. The hydrological functioning of these sandy soils during rain cycles is unknown. One way to tackle this issue is to spatialize variations in water content but these are difficult to measure in the vadose zone. We investigated the use of Electrical Resistivity Tomography (ERT) as a technique for spatializing resistivity in a non-destructive manner in order to improve our knowledge of relevant hydrological processes. To achieve this, two approaches were examined. First, we focused on a possible link between water tension (which is much easier to measure in the field by point measurements than water content), and resistivity (spatialized with ERT). Second, because ERT is affected by solution non-uniqueness and reconstruction smoothing, we improved the accuracy of ERT inversion by comparing calculated solutions with in-situ resistivity measurements. We studied a natural microdune during a controlled field experiment with artificial sprinkling which reproduced typical rainfall cycles. We recorded temperature, water tension and resistivity within the microdune and applied surface ERT before and after the 3 rainfall cycles. Soil samples were collected after the experiment to determine soil physical characteristics. An experimental relationship between water tension and water content was also investigated. Our results showed that the raw relationship between calculated ERT resistivity and water tension measurements in sand is highly scattered because of significant spatial variations in porosity. An improved correlation was achieved by using resistivity ratio and water tension differences. The slope of the relationship depends on the soil solution conductivity, as predicted by Archie's law when salted water was used for the rain simulation. We found that determining the variations in electrical resistivity is a sensitive method for spatializing the differences in water tension which are directly linked with infiltration and evaporation/drainage processes in the vadose zone. However, three factors complicate the use of this approach. Firstly, the relation between water tension and water content is generally non-linear and dependent on the water content range. This could limit the use of our site-specific relations for spatializing water content with ERT through tension. Secondly, to achieve the necessary optimization of ERT inversion, we used destructive resistivity measurements in the soil, which renders ERT less attractive. Thirdly, we found that the calculated resistivity is not always accurate because of the smoothing involved in surface ERT data inversion. We conclude that further developments are needed into ERT image reconstruction before water tension (and water content) can be spatialized in heterogeneous sandy soils with the accuracy needed to routinely study their hydrological functioning.     

Regional soil moisture simulation for Shaanxi Province using SWAT model validation and trend analysis

The soil moisture in Shaanxi Province,a region with complex topography,is simulated using the distributed hydrological model Soil Water Assessment Tool(SWAT).Comparison and contrast of modeled and observed soil moisture show that the SWAT model can reasonably simulate the long-term trend in soil moisture and the spatiotemporal variability of soil moisture in the region.Comparisons to NCEP/NCAR and ERA40 reanalysis of soil moisture show that the trend of variability in soil moisture simulated by SWAT is more...     

洞庭湖洲滩土壤种子库对土壤水分变化的响应

由于三峡大坝及上游水库群的运行,长江中下游水域水文节律随之发生了改变,导致洞庭湖枯水期提前,进而影响洞庭湖洲滩植被及其土壤种子库的分布格局.本研究在洞庭湖4个自然保护区内选取共11个典型洲滩湿地,沿由水到陆方向根据植被类型将洲滩分为泥沙洲滩、泥沙—湖草洲滩过渡带、湖草洲滩、湖草—南荻洲滩过渡带、南荻洲滩5种洲滩类型.通过样带—样方法调查和采样,并结合湿润和水淹两种条件下的土壤种子库萌发实验,分析了土壤水分变化对洲滩种子库萌发特征的影响及土壤种子库与地表植被的关系.结果显示:①土壤含水量沿水到陆方向由泥沙洲滩向南荻洲滩递减;②不同类型洲滩土壤种子库密度没有显著差异;③温室萌发实验中,水淹条件下土壤种子库物种丰富度和种子库密度显著降低,东洞庭湖自然保护区土壤种子库物种丰富度和种子库密度较高;④地表植被物种丰富度高于土壤种子库,泥沙洲滩土壤种子库与地表植被物种组成的Jaccard相似性指数最低.此外,虉草(Phalaris arundinacea)、芦苇(Phragmites communis)、南荻(Miscanthus sacchariflorus)等只在地表植被中存在,而陌上菜(Lin...     

Information and complexity measures applied to observed and simulated soil moisture time series

Abstract

Time series of soil moisture-related parameters provide important insights into the functioning of soil water systems. Analysis of patterns within such time series has been used in several studies. The objective of this work was to compare patterns in observed and simulated soil moisture contents to understand whether modelling leads to a substantial loss of information or complexity. The time series were observed at four plots in sandy soils within the USDA-ARS OPE3 experimental watershed, for a year; precipitation and evapotranspiration (ET) were measured and estimated, respectively, and used for soil water flow simulation with the HYDRUS-1D software. The information content measures are the metric entropy and the mean information gain, and complexity measures are the fluctuation complexity and the effective measure complexity. These measures were computed based on the binary encoding of soil moisture time series, and used probabilities of patterns, i.e. probabilities of joint or sequential appearance of symbol sequences. The information content of daily soil moisture time series was much smaller than that of rainfall data, and had higher complexity, indicating that soil worked essentially as an information filter. Information content and complexity decreased and increased with depth, respectively, demonstrating the increase in the information filtering action of soil. The information measures of simulated soil moisture content were close to those of the measurements, indicating the successful simulation of patterns in the data. The spatial variability of the information measures for simulated soil moisture content at all depths was less pronounced than the one of measured time series. Compared with precipitation and estimated ET, soil moisture time series had more structure and less randomness in this work. The information measures can provide useful complementary knowledge about model performance and patterns in observation and modelling results.

Citation Pan, F., Pachepsky, Y. A., Guber, A. K., & Hill, R. L. (2011) Information and complexity measures applied to observed and simulated soil moisture time series. Hydrol. Sci. J. 56(6), 1027–1039.     

On the role of soil moisture in daytime evolution of temperatures

This paper explores the relationship between temperature, evaporation and soil moisture using a planetary boundary layer (PBL) model. It focuses on illustrating and quantifying the effect of soil moisture on the evolution of daytime temperatures. A simple convective PBL model coupled with the Penman–Monteith (PM) equation is used to estimate evapotranspiration. Following calibration and sensitivity analysis, the model was used to simulate the relative impact of dry and wet soil moisture conditions on daytime temperatures by changing the surface resistance parameter in the PM equation. It was found that the maximum temperature that can be reached during a day is constrained by the amount of soil moisture and the available net radiation, confirming previously published results. Higher temperatures can be reached with greater net radiation and dry soil moisture conditions. Copyright © 2012 John Wiley & Sons, Ltd.     

Analysis of soil moisture dynamics beneath olive trees

Olive cultivation is a widespread land use in Mediterranean climates. The proper implementation of soil and water conservation practices in groves requires detailed knowledge of the governing hydrological processes. In this work topsoil moisture dynamics under wet and dry conditions and across a small catchment was investigated in the inter row (IR) and directly under the olive tree canopies (UC). We do this using a sensor network (11 stations) and a simple bucket model which was calibrated (June, 2011–2012) and validated (June, 2012–2013). During most of the year the normalized soil moisture contents (s) were greater in the IR than under UC, with an average normalized soil moisture difference of 0.12. The difference between UC and IR normalized soil moisture followed a seasonal pattern, reaching a maximum near 0.30 during spring. An analysis of the normalized soil moisture probability density functions (pdfs) was bimodal, showing characteristic dominant wet and dry soil moisture states, with the highest probability densities for the dry state. Overall the spatial variability of soil moisture was lower UC than in the IR. This was a result of the soil moisture buffering capacity of the canopy with respect to rainfall and evaporation, in addition to observed differences in soil properties. Hourly soil moisture data were successfully modelled (R² > 0.85), both UC and in the IR, yet with the inclusion of a simple formulation for canopy interception for the former. The results provide insight into how olive trees change hydrological processes in their neighbourhood. Copyright © 2016 John Wiley & Sons, Ltd.     

Streamflow predictability in the Saskatchewan/Nelson River basin given macroscale estimates of the initial soil moisture status

Abstract

Soil moisture estimates obtained over large spatial areas will become increasingly available through current and upcoming satellite missions and from numerous land surface parameterization schemes run at global- and continental-scale resolutions. The goal of this research was to evaluate the potential for using macroscale estimates of soil moisture for enhancing streamflow forecasts. Towards this research objective, monthly streamflow estimates were obtained from over 50 gauge locations within the Nelson basin, Canada, for the period 1979–1999. For each streamflow record, multiple linear regression models were used to remove components of the streamflow signal related to previous streamflow, climate teleconnections (e.g. ENSO and AO) and snow water equivalence. Correlations were then assessed between the macroscale soil moisture estimates and the residuals of the multiple linear regression analysis over lead times of one, two and three months. At the one- and two-month lead time, statistically significant relationships between soil moisture and the residuals of streamflow are observed over a large proportion of the gauging locations. The number of catchments with statistically significant relationships decreases significantly after two months and particularly in the months of April—June. This study demonstrates that available macroscale estimates of soil moisture have the potential to enhance streamflow prediction, although further study is suggested to improve upon the soil moisture estimates and their application in a forecast system.     

Spatial variation in solutional denudation and soil moisture over a hillslope hollow

A spatial pattern of relative solutional denudation is described for a hillslope hollow and adjacent spurs at Bicknoller Combe, Somerset. The pattern was obtained from a network of micro-weight loss rock tablets emplaced in the soil. The results show that the hollow is the main locus of solutional denudation. The soil moisture distribution over the hollow indicates that it is a transmission zone for acid soil water percolating from the adjacent spurs to the saturated wedge at the base of the hollow. The wetter acid soils in the hollow are responsible for the relatively higher solutional denudation taking place in the hollow.     

Assessment of drought vulnerability based on the soil moisture PDF

This paper studies the statistics of the soil moisture condition and its monthly variation for the purpose of evaluating drought vulnerability. A zero-dimensional soil moisture dynamics model with the rainfall forcing by the rectangular pulses Poisson process model are used to simulate the soil moisture time series for three sites in Korea: Seoul, Daegu, and Jeonju. These sites are located in the central, south-eastern, and south-western parts of the Korean Peninsular, respectively. The model parameters are estimated on a monthly basis using hourly rainfall data and monthly potential evaporation rates obtained by the Penmann method. The resulting soil moisture simulations are summarized on a monthly basis. In brief, the conclusions of our study are as follows. (1) Strong seasonality is observed in the simulations of soil moisture. The soil moisture mean is less than 0.5 during the dry spring season (March, April, and June), but other months exceed the 0.5 value. (2) The spring season is characterized by a low mean value, a high standard deviation and a positive skewness of the soil moisture content. On the other hand, the wet season is characterized by a high mean value, low standard deviation, and negative skewness of the soil moisture content. Thus, in the spring season, much drier soil moisture conditions are apparent due to the higher variability and positive skewness of the soil moisture probability density function (PDF), which also indicates more vulnerability to severe drought occurrence. (3) Seoul, Daegue, and Jeonju show very similar overall trends of soil moisture variation; however, Daegue shows the least soil moisture contents all through the year, which implies that the south-eastern part of the Korean Peninsula is most vulnerable to drought. On the other hand, the central part and the south-western part of the Korean peninsula are found to be less vulnerable to the risk of drought. The conclusions of the study are in agreement with the climatology of the Korean Peninsula.     

An analysis of infiltration with moisture content distribution in a two‐dimensional discretized water content domain

On the basis of unsaturated Darcy's law, the Talbot–Ogden method provides a fast unconditional mass conservative algorithm to simulate groundwater infiltration in various unsaturated soil textures. Unlike advanced reservoir modelling methods that compute unsaturated flow in space, it only discretizes the moisture content domain into a suitable number of bins so that the vertical water movement is estimated piecewise in each bin. The dimensionality of the moisture content domain is extended from one dimensional to two dimensional in this study, which allows us to distinguish pore shapes within the same moisture content range. The vertical movement of water in the extended model imitates the infiltration phase in the Talbot–Ogden method. However, the difference in this extension is the directional redistribution, which represents the horizontal inter‐bin flow and causes the water content distribution to have an effect on infiltration. Using this extension, we mathematically analyse the general relationship between infiltration and the moisture content distribution associated with wetting front depths in different bins. We show that a more negatively skewed moisture content distribution can produce a longer ponding time, whereas a higher overall flux cannot be guaranteed in this situation. It is proven on the basis of the water content probability distribution independent of soil textures. To illustrate this analysis, we also present numerical examples for both fine and coarse soil textures. Copyright © 2014 John Wiley & Sons, Ltd.     

Measurement and modeling of moisture content above an oscillating water table: implications for beach surface moisture dynamics

This study examined the influence of tidally‐induced oscillations of the beach water table in regulating beach surface moisture dynamics. A series of laboratory experiments were conducted to assess the influence of hysteresis and transient flow effects on surface moisture variability. The experimental apparatus utilized a column of well‐sorted fine sand partially immersed in a reservoir of water. The water level in the reservoir was raised and lowered via a diaphragm‐metering pump to simulate tidally induced fluctuations of the water table, and the moisture content profile within the column was monitored using an array of Delta‐T probes. Moisture contents at specific elevations within the column were utilized as proxies to represent various ‘surface’ elevations (relative to the high water table). Results indicate that surface moisture content behaves in a distinctly hysteretic manner. Examination of water flow scanning curves illustrated that for all surface elevations considered, higher moisture contents for a given pressure head occurred during the drying cycle than during the wetting cycle. This observation is particularly evident with shallow surface elevations (i.e. water table close to the surface) where the Haines Jump phenomenon was found to have a significant influence on moisture content dynamics. Additionally, an assessment of the accuracy of hysteretic and non‐hysteretic models to predict the measured moisture contents demonstrated that hysteretic simulations consistently provide a better representation of the observed moisture contents than non‐hysteretic simulations. A time lag was found between the respective maxima and minima in water table elevation surface moisture content. At the near surface water table positions the time lag ranged between 30 and 100 minutes, and it increased to 240 minutes (four hours) with the high water table at 60 cm below the surface. Copyright © 2013 John Wiley & Sons, Ltd.