Relating the depth of the water table to the depth of weathering

Weathering of bedrock creates and occludes permeability, affecting subsurface water flow. Often, weathering intensifies above the water table. On the contrary, weathering can also commence below the water table. To explore relationships between weathering and the water table, a simplified weathering model for an eroding hillslope was formulated that takes into account both saturated and unsaturated subsurface water flow (but does not fully account for changes in dissolved gas chemistry). The phreatic line was calculated using solutions to mathematical treatments for both zones. In the model, the infiltration rate at the hill surface sets both the original and the eventual steady-state position of the water table with respect to the weathering reaction front. Depending on parameters, the weathering front can locate either above or below the water table at steady state. Erosion also affects the water table position by changing porosity and permeability even when other hydrological conditions (e.g. hydraulic conductivity of parent material, infiltration rate at the surface) do not change. The total porosity in a hill (water storage capacity) was found to increase with infiltration rate (all else held constant). This effect was diminished by increasing the erosion rate. We also show examples of how the infiltration rate affects the position of the water table and how infiltration rate affects weathering advance. Published 2020. This article is a U.S. Government work and is in the public domain in the USA     

Insights on stream temperature processes through development of a coupled hydrologic and stream temperature model for forested coastal headwater catchments

Stream temperature controls a number of biological, chemical, and physical processes occurring in aquatic environments. Transient snow cover and advection associated with lateral throughflow inputs can have a dominant influence on stream thermal regimes for headwater catchments in the rain‐on‐snow zone. Most existing stream temperature models lack the ability to properly simulate these processes. We developed and evaluated a conceptual‐parametric catchment‐scale stream temperature model that includes the role of transient snow cover and lateral advection associated with throughflow. The model consists of routines for simulating canopy interception, snow accumulation and melt, hillslope throughflow runoff and temperature, and stream channel energy exchange processes. The model was used to predict discharge and stream temperature for a small forested headwater catchment near Vancouver, Canada, using long‐term (1963–2013) weather data to compute model forcing variables. The model was evaluated against 4 years of observed stream temperature. The model generally predicted daily mean stream temperature accurately (annual RMSE between 0.57 and 1.24 °C) although it overpredicted daily summer stream temperatures by up to 3 °C during extended low streamflow conditions. Model development and testing provided insights on the roles of advection associated with lateral throughflow, channel interception of snow, and surface–subsurface water interactions on stream thermal regimes. This study shows that a relatively simple but process‐based model can provide reasonable stream temperature predictions for forested headwater catchments located in the rain‐on‐snow zone.     

Constraining water table depth simulations in a land surface model using estimated baseflow

Several recent studies have shown the significance of representing groundwater in land surface hydrologic simulations. However, optimal methods for model parameter calibration in order to realistically simulate baseflow and groundwater depth have received little attention. Most studies still use globally constant groundwater parameters due to the lack of available datasets for calibration. Moreover, when models are calibrated, various parameter combinations are found to exhibit equifinality in simulated total runoff due to model parameter interactions. In this study, a simple lumped groundwater model is incorporated into the Community Land Model (CLM), in which the water table is interactively coupled to soil moisture through the groundwater recharge fluxes. The coupled model (CLMGW) is successfully validated in Illinois using a 22-year (1984–2005) monthly observational dataset. Baseflow estimates from the digital recursive filter technique are used to calibrate the CLMGW parameters. The advantage obtained from incorporating baseflow calibration in addition to traditional calibration based on measured streamflow alone is demonstrated by a Monte Carlo-type simulation analysis. Using the optimal parameter sets identified from baseflow calibration, flow partitioning and water table depth simulations using CLMGW are improved, and the equifinality problem is alleviated. For other regions that lack observations of water table depth, the baseflow calibration approach can be used to enhance parameter estimation and constrain water table depth simulations.     

Temporal variation of chemical and mechanical weathering in NE Iceland: Evaluation of a steady-state model of erosion

This study critically assesses the temporal sensitivity of the steady-state model of erosion that has been applied to chemical and mechanical weathering studies of volcanic islands and the continents, using only one sample from each catchment. The model assumes a geochemical mass balance between the initially unweathered rock of a drainage basin and the dissolved and solid loads of the river.Chemical composition of 178 samples of suspended and dissolved inorganic river constituents, collected in 1998–2002, were studied from five basaltic river catchments in NE Iceland. The Hydrological Service in Iceland has monitored the discharge and the total suspended inorganic matter concentration (SIM) of the glacial rivers for ~ four decades, making it possible to compare modelled and measured SIM fluxes.Concentration of SIM and grain size increased with discharge. As proportion of clay size particles in the SIM samples increased, concentrations of insoluble elements increased and of soluble decreased. The highest proportion of altered basaltic glass was in the clay size particles.The concentration ratio of insoluble elements in the SIM was used along with data on chemical composition of unweathered rocks (high-Mg basalts, tholeiites, rhyolites) to calculate the pristine composition of the original catchment rocks. The calculated rhyolite proportions compare nicely with area-weighted average proportions, from geological maps of these catchments.The calculated composition of the unweathered bedrock was used in the steady-state model, together with the chemical composition of the suspended and dissolved constituents of the river. Seasonal changes in dissolved constituent concentrations resulted in too low modelled concentrations of SIM_mod at high discharge (and too high SIM_mod at low discharge). Samples collected at annual average river dissolved load yielded SIM_mod concentrations close to the measured ones. According to the model, the studied rivers had specific mechanical denudation rates of 1.3–3.0 kg/m²/yr whereas the average measured rates were 0.8–3.5 kg/m²/yr which are among the highest on Earth.This study validates the use of a steady-state model of erosion to estimate mechanical weathering rates at the scale of a river catchment when the collected riverine dissolved load represents the average chemical composition over a mean hydrological year.     

MHYDAS‐Erosion: a distributed single‐storm water erosion model for agricultural catchments

In this paper, we present MHYDAS‐Erosion, a dynamic and distributed single‐storm water erosion model developed as a module of the existing hydrological MHYDAS model. As with many catchment erosion models, MHYDAS‐Erosion is able to simulate sediment transport, erosion and deposition by rill and interrill processes. Its originality stems from its capacity to integrate the impact of land management practices (LMP) as key elements controlling the sedimentological connectivity in agricultural catchments. To this end, the water‐sediment pathways are first determined by a specific process‐oriented procedure defined and controlled by the user, which makes the integration of LMP easier. The LMP dynamic behaviours are then integrated into the model as a time‐dependent function of hydrological variables and LMP characteristics. The first version of the model was implemented for vegetative filters and tested using water and sediment discharge measurements at three nested scales of a densely instrumented catchment (Roujan, OMERE Observatory, southern France). The results of discharge and soil loss for simulated rainfall events have been found to acceptably compare with available data. The average R² values for water and sediment discharge are 0·82 and 0·83, respectively. The sensitivity of the model to changes in the proportion of LMP was assessed for a single rain event by considering three scenarios of the Roujan catchment management with vegetative filters: 0% (Scenario 1), 18% (Scenario 2, real case) and 100% (Scenario 3). Compared to Scenario 2 (real case), soil losses decreased for Scenario 3 by 65% on the agricultural plot scale, 62% on the sub‐catchment scale and 45% at the outlet of the catchment and increased for Scenario 1 by 0% on the plot scale, 26% on the sub‐catchment scale and 18% at the outlet of the catchment. Copyright © 2010 John Wiley & Sons, Ltd.     

Analysis of the contributions of topographic,soil, and vegetation features on the spatial distributions of surface soil moisture in a steep natural forested headwater catchment

Surface soil moisture has been extensively studied for various land uses and landforms. Although many studies have reported potential factors that control surface soil moisture over space or time, the findings have not always been consistent, indicating a need for identification of the main factors. This study focused on the static controls of topographic, soil, and vegetation features on surface soil moisture in a steep natural forested headwater catchment consisting of three hillslope units of a gully area, side slope, and valley‐head slope. Using a simple correlation analysis to investigate the effects of the static factors on surface soil moisture at depths of 0–20 cm at 470 points in 13 surveys, we addressed the characteristics of surface soil moisture and its main controlling factors. The results indicated that the mean of surface soil moisture was in the decreasing order of gully area > valley‐head slope > side slope. The relationship between the mean and standard deviation of surface soil moisture showed a convex‐upward shape in the headwater catchment, a negative curvilinear shape in the gully area, and positive curvilinear shapes at the side and valley‐head slopes. At the headwater catchment and valley‐head slope, positive contributions of soil porosity and negative contributions of slope gradient and saturated hydraulic conductivity were the main controlling factors of surface soil moisture under wetter conditions, whereas positive contributions of topographic wetness index and negative contributions of vegetation density were the main controlling factors of surface soil moisture under drier conditions. At the side slope underlain by fractured bedrocks, only saturated hydraulic conductivity and vegetation density were observed to be the controlling factors. Surface soil moisture in the gully area was mainly affected by runoff rather than were static features. Thus, using hillslope units is effective for approximately estimating the hydrological behaviours of surface moisture on a larger scale, whereas dependency between the main static factors and moisture conditions is helpful for estimating the spatial distributions of surface moisture on a smaller scale.     

Combination of CFCs and stable isotopes to characterize the mechanism of groundwater–surface water interactions in a headwater basin of the North China Plain

Mountainous areas are characterized by steep slopes and rocky landforms, with hydrological conditions varying rapidly from upstream to downstream, creating variable interactions between groundwater and surface water. In this study, mechanisms of groundwater–surface water interactions within a headwater catchment of the North China Plain were assessed along the stream length and during different seasons, using hydrochemical and stable isotope data, and groundwater residence times estimated using chlorofluorocarbons. These tracers indicate that the river is gaining, due to groundwater discharge in the headwater catchment both in the dry and rainy seasons. Residence time estimation of groundwater using chlorofluorocarbons data reveals that groundwater flow in the shallow sedimentary aquifer is dominated by the binary mixing of water approximating a piston flow model along 2 flow paths: old water, carried by a regional flow system along the direction of river flow, along with young water, which enters the river through local flow systems from hilly areas adjacent to the river valley (particularly during the rainy season). The larger mixing ratio of young water from lateral groundwater recharge and return flow of irrigation during the rainy season result in higher ion concentrations in groundwater than in the dry season. The binary mixing model showed that the ratio of young water versus total groundwater ranged from 0.88 to 0.22 and 1.0 to 0.74 in the upper and lower reaches, respectively. In the middle reach, meandering stream morphology allows some loss of river water back into the aquifer, leading to increasing estimates of the ratio of young water (from 0.22 to 1). This is also explained by declining groundwater levels near the river, due to groundwater extraction for agricultural irrigation. The switch from a greater predominance of regional flow in the dry season, to more localized groundwater flow paths in the wet season is an important groundwater–surface water interactions mechanism, with important catchment management implications.     

Base flow separation for soil erosion simulation in a granitic forested headwater catchment using a process-based model,GeoWEPP

Distributed erosion models, which simulate the physical processes of water flow and soil erosion, are effective for predicting soil erosion in forested catchments. Although subsurface flow through multiple pathways is dominant for runoff generation in forested headwater catchments, the process-based erosion model, Geo-spatial interface for Water Erosion Prediction Project(Geo WEPP), does not have an adequate subsurface component for the simulation of hillslope water flow. In the current study, t...     

Travel times for snowmelt-dominated headwater catchments: Influences of wetlands and forest harvesting,and linkages to stream water quality

The time it takes water to travel through a catchment, from when it enters as rain and snow to when it leaves as streamflow, may influence stream water quality and catchment sensitivity to environmental change. Most studies that estimate travel times do so for only a few, often rain-dominated, catchments in a region and use relatively short data records (<10 years). A better understanding of how catchment travel times vary across a landscape may help diagnose inter-catchment differences in water quality and response to environmental change. We used comprehensive and long-term observations from the Turkey Lakes Watershed Study in central Ontario to estimate water travel times for 12 snowmelt-dominated headwater catchments, three of which were impacted by forest harvesting. Chloride, a commonly used water tracer, was measured in streams, rain, snowfall and as dry atmospheric deposition over a 31 year period. These data were used with a lumped convolution integral approach to estimate mean water travel times. We explored relationships between travel times and catchment characteristics such as catchment area, slope angle, flowpath length, runoff ratio and wetland coverage, as well as the impact of harvesting. Travel time estimates were then used to compare differences in stream water quality between catchments. Our results show that mean travel times can be variable for small geographic areas and are related to catchment characteristics, in particular flowpath length and wetland cover. In addition, forest harvesting appeared to decrease mean travel times. Estimated mean travel times had complex relationships with water quality patterns. Results suggest that biogeochemical processes, particularly those present in wetlands, may have a greater influence on water quality than catchment travel times.     

Simulating daily, monthly and annual water balances in a land surface model using alternative root water uptake schemes

Hydrological simulations at multi-temporal time scales by a widely used land surface model (LSM) are investigated under contrasting vegetation and meteorological conditions. Our investigation focuses particularly on the effects of two different representations of root water uptake and root profile on simulated evapotranspiration (ET) and soil moisture by the Integrated BIosphere Simulator (IBIS). For this purpose, multi-year eddy covariance measurements, collected at four flux-tower sites across North America, were used to gauge IBIS simulations with: (a) its standard version (IBIS2.1), in which static root water uptake (RWU) and root profile schemes are incorporated; and (b) a modified version in which dynamic RWU and root profile schemes replaces the static schemes used in the standard version. Overall, our results suggest that the modified version of the model performs more realistically than the standard version, particularly when high atmospheric demand for evaporation is combined with high atmospheric vapour pressure deficit and low soil water availability. The overall correlation between simulated and measured monthly ET rates at the simulated sites reached 0.87 and 0.91 for the standard and the modified versions, respectively. Our results also show that the incorporation of the dynamic RWU in IBIS yields improved simulations of ET under very dry conditions, when soil moisture falls down to very low levels. This suggests that adequate representations of vegetation responses to drought are needed in LSMs as many state of the art climate models projections of future climate indicate more frequent and/or more intense drought events occurring in some regions of the globe. Our analysis also highlighted the urgent need for adequate methodologies to correct field measurements that exhibit energy imbalances in order to provide rigorous assessments of land surface model simulations of heat and mass exchanges between the land surface and the atmosphere.     

Combining a rainfall–runoff model and a regionalization approach for flood and water resource assessment in the western Po Valley,Italy

ABSTRACT

Selecting the best structure and parameterization of rainfall–runoff models is not straightforward and depends on a broad number of factors. In this study, the “Modello Idrologico Semi-Distribuito in continuo” (MISDc) was tested on 63 mountainous catchments in the western Po Valley (Italy) and the optimal model parameters were regionalized using different strategies. The model performance was evaluated through several indexes analysing hydrological regime, high-flow condition and flow–duration curve (FDC). In general, MISDc provides a good fit behaviour with a Kling-Gupta Efficiency index greater than 0.5 for 100% and 84% of cases for calibration and validation, respectively. Concerning the regionalization, spatial proximity approach is the most accurate solution obtaining satisfactory performance. Lastly, the predicted FDCs showed an excellent similarity with the observed ones. Results encourage to apply MISDc over the study area for flood forecasting and for assessing water resources availability thanks to the modest computational efforts and data requirements.     

Intercomparison across scales between remotely-sensed land surface temperature and representative equilibrium temperature from a distributed energy water balance model

Abstract

Quantifying the reliability of distributed hydrological models is an important task in hydrology to understand their ability to estimate energy and water fluxes at the agricultural district scale as well the basin scale for water resources management in drought monitoring and flood forecasting. In this context, the paper presents an intercomparison of simulated representative equilibrium temperature (RET) derived from a distributed energy water balance model and remotely-sensed land surface temperature (LST) at spatial scales from the agricultural field to the river basin. The main objective of the study is to evaluate the use of LST retrieved from operational remote sensing data at different spatial and temporal resolutions for the internal validation of a distributed hydrological model to control its mass balance accuracy as a complementary method to traditional calibration with discharge measurements at control river cross-sections. Modelled and observed LST from different radiometric sensors located on the ground surface, on an aeroplane and a satellite are compared for a maize field in Landriano (Italy), the agricultural district of Barrax (Spain) and the Upper Po River basin (Italy). A good ability of the model in reproducing the observed LST values in terms of mean bias error, root mean square error, relative error and Nash-Sutcliffe index is shown.

Editor Z.W. Kundzewicz; Associate editor D. Gerten     

Improvement of the Noah land surface model for warm season processes: evaluation of water and energy flux simulation

The Noah model is a land surface model of the National Centers for Environmental Prediction. It has been widely used in regional coupled weather and climate models (i.e. Weather Research and Forecasting Model, Eta Mesoscale Model) and global coupled weather and climate models (i.e. National Centers for Environmental Prediction Global Forecast System, Climate Forecast System). Therefore, its continued improvement and development are keys to enhancing our weather and climate forecast ability and water and energy flux simulation accuracy. North American Land Data Assimilation System phase 1 (NLDAS‐1) experiments indicated that the Noah model exhibited substantial bias in latent heat flux, total runoff and land skin temperature during the warm season, and such bias can significantly affect coupled weather and climate models. This paper presents a study to improve the Noah model by adding model parameterization processes such as including seasonal factor on leaf area index and root distribution and selecting optimal model parameters. We compared simulated latent heat flux, mean annual runoff and land skin temperature from the Noah control and test versions with measured latent heat flux, land surface skin temperature, mean annual runoff and satellite‐retrieved land surface skin temperature. The results show that the test version significantly reduces biases in latent heat, total runoff and land skin temperature simulation. The test version has been used for the NLDAS phase 2 (NLDAS‐2) to produce 30‐year water flux, energy flux and state variable products to support the US drought monitor of National Integrated Drought Information System. Copyright © 2012 John Wiley & Sons, Ltd.     

Soil loss,water ponding and sediment deposition variations as a consequence of rainfall intensity and land use: a multi‐criteria analysis

Prediction of areas prone to land degradation in agricultural catchments is a complex task. This is due to the difficulties encountered in data gathering over wide regions and in the translation of existing scientific knowledge to a quantitative and spatially explicit risk assessment system. This paper incorporates the use of remotely sensed data, terrain analysis and a multi‐criteria mechanism for evaluating risks of soil loss, water ponding, and sediment deposition in a mid‐size agricultural Mediterranean catchment, under 80 years of intensive cultivation. The research uses simulations to study the effect of topographic attributes, soil characteristics, vegetation cover, rainfall intensity and human activities on the three above‐mentioned processes. The results show that, from the methodological point of view, the integration of knowledge from several experts yields better predictive results than relying on a single expert, even the one found to be most consistent. Also, the use of a simple weighted linear combination was more useful than the more sophisticated computerized programming technique. From the phenomenological point of view, the increase in rainfall intensity and land‐use transformation from orchard to field‐crops has led to a significant increase in soil loss and sediment yield, while extreme changes in tillage direction have only yielded minor changes in water ponding. The developed system's predictive capabilities also show that the outcomes can be used as a basis for decisions on catchment management in regions of high environmental sensitivity. Copyright © 2009 John Wiley & Sons, Ltd.     

Parameter transferability under changing climate: case study with a land surface model in the Durance watershed,France

Abstract

In physically-based land surface models, the parameters can all be prescribed a priori but calibration can be used to enhance the realism of the simulations in well instrumented domains. In such a case, the transferability of calibrated parameters under non-stationary conditions needs to be addressed, especially in the context of climate change. To this end, we used the Catchment Land Surface Model (CLSM) in the Upper Durance watershed located in the French Alps, which experienced a significant increase in temperature over the last century. The CLSM is forced by a 50-year meteorological dataset of good quality. Four parameters of the CLSM (one related to snow processes and three to soil properties) are calibrated against discharge observations with a multi-objective algorithm. First, the robustness of the CLSM parameterizations is tested by the Differential Split Sample Test (DSST). The simulations show good performances over a wide range of retrospective climatic conditions, except when the parameters are calibrated over a period with a large contribution of snowmelt to annual mean discharge. Then, the use of a climate change scenario reveals that the parameterizations of soil moisture processes in the CLSM are responsible for an increasing dispersion among simulations when facing dry and warm conditions. However, the differences between the simulated changes of river discharge remain very small. This work shows that calibration conveys some uncertainties, but they are moderate in the studied case, and pertain to the most conceptual parameterizations of this physically-based model.     

Reproduction of Pechora runoff hydrographs with the help of a model of heat and water exchange between the land surface and the atmosphere (SWAP)

The potentialities of a procedure for calculating the Pechora River runoff from the pan-Arctic river basin are studied. The procedure is based on the use of a model describing heat and water exchange between the land surface and the atmosphere and two variants of input data sets relying on global databases on meteorological characteristics and land surface parameters and data of standard measurements of meteorological characteristics in combination with parameters of the land surface of the basin, taken from global databases. In both cases, use was made of the method for optimizing part of the most important model parameters, including both land surface parameters and correction factors for some meteorological elements.     

Evaluation of multiple surface soil moisture for Korean regional flux monitoring network sites: Advanced Microwave Scanning Radiometer E,land surface model,and ground measurements

In the past few decades, there have been great developments in remotely sensed soil moisture, with validation efforts using land surface models (LSMs) and ground‐based measurements, because soil moisture information is essential to understanding complex land surface–atmosphere interactions. However, the validation of remotely sensed soil moisture has been very limited because of the scarcity of the ground measurements in Korea. This study validated Advanced Microwave Scanning Radiometer E (AMSR‐E) soil moisture data with the Common Land Model (CLM), one of the most widely used LSMs, and ground‐based measurements at two Korean regional flux monitoring network sites. There was reasonable agreement regarding the different soil moisture products for monitoring temporal trends except National Snow and Ice Data Centre (NSIDC) AMSR‐E soil moisture, albeit there were essential comparison limitations by different spatial scales and soil depths. The AMSR‐E soil moisture data published by the National Aeronautics and Space Administration and Vrije Universiteit Amsterdam (VUA) showed potential to replicate temporal variability patterns (root‐mean‐square errors = 0·10–0·14 m³ m^?3 and wet BIAS = 0·09 ? 0·04 m³ m^?3) with the CLM and ground‐based measurements. However, the NSIDC AMSR‐E soil moisture was problematic because of the extremely low temporal variability and the VUA AMSR‐E soil moisture was relatively inaccurate in Gwangneung site characterized by complex geophysical conditions. Additional evaluations should be required to facilitate the use of recent and forthcoming remotely sensed soil moisture data from Soil Moisture and Ocean Salinity and Soil Moisture Active and Passive missions at representative future validation sites. Copyright © 2011 John Wiley & Sons, Ltd.     

Predicting dry‐season flows with a monthly rainfall–runoff model: Performance for gauged and ungauged catchments

Hydrologic models are useful to understand the effects of climate and land‐use changes on dry‐season flows. In practice, there is often a trade‐off between simplicity and accuracy, especially when resources for catchment management are scarce. Here, we evaluated the performance of a monthly rainfall–runoff model (dynamic water balance model, DWBM) for dry‐season flow prediction under climate and land‐use change. Using different methods with decreasing amounts of catchment information to set the four model parameters, we predicted dry‐season flow for 89 Australian catchments and verified model performance with an independent dataset of 641 catchments in the United States. For the Australian catchments, model performance without catchment information (other than climate forcing) was fair; it increased significantly as the information to infer the four model parameters increased. Regressions to infer model parameters from catchment characteristics did not hold for catchments in the United States, meaning that a new calibration effort was needed to increase model performance there. Recognizing the interest in relative change for practical applications, we also examined how DWBM could be used to simulate a change in dry‐season flow following land‐use change. We compared results with and without calibration data and showed that predictions of changes in dry‐season flow were robust with respect to uncertainty in model parameters. Our analyses confirm that climate is a strong driver of dry‐season flow and that parsimonious models such as DWBM have useful management applications: predicting seasonal flow under various climate forcings when calibration data are available and providing estimates of the relative effect of land use on seasonal flow for ungauged catchments.     

The effects of coal gangue and fly ash on the hydraulic properties and water content distribution in reconstructed soil profiles of coal‐mined land with a high groundwater table

Soil water systems have been severely degraded in coal‐mined and subsiding land, where a shallow groundwater table is also present. The present paper discussed the effects of fly ash (FA) and coal gangue (CG) as filling materials on the hydraulic properties and water content distribution in a profile for the purpose of rehabilitating subsided lands. The saturated water content, water characteristic curve, and water diffusivity of local soil, FA, CG, and a mixture of FA and CG (“mixed filling”) were characterized. A column experiment was conducted to investigate the changes of water content in profiles reconstructed from the combination of soil and filling materials, including soil only, FA, CG, and a mixture of FA and CG, which were used to fill the lower part of the reconstructed profile. The mixture of FA and CG was found to possess similar hydraulic properties to those of the soil, particularly high water‐holding capacity and permeability. Moreover, the volumetric water contents in the whole profile containing the mixture of FA and CG were consistent with those of the profile reconstructed with soil only. As a result, it is recommended to adopt the mixture of FA and CG for reconstructing the lower profile of the land to alleviate or rehabilitate subsided land in coal mines.