Understanding the dynamics of soil organic carbon (SOC) is of fundamental importance in land use and management, whether in the current researches or in future scenarios of agriculture systems considering climate change. In order to evaluate SOC stock of the three districts (Delmiro Gouveia, Pariconha, and Inhapi districts) in the semi-arid region of Brazil in rainfed and irrigated agriculture systems under different climate scenarios using the Century model, we obtained RCP4.5 and RCP8.5 climate scenarios derived from the Eta Regional Climate Model (Eta-HadGEM2-ES and Eta-MIROC5) from the National Institute for Space Research, and then input the data of bulk density, pH, soil texture, maximum temperature, minimum temperature, and rainfall into the soil and climate files of the Century model. The results of this study showed that the Eta-HadGEM2-ES model was effective in estimating air temperature in the future period. In rainfed agriculture system, SOC stock under the baseline scenario was lower than that under RCP4.5 and RCP8.5 climate scenarios, while in irrigated agriculture system, SOC stock in the almost all climate scenarios (RCP4.5 and RCP8.5) and models (Eta-HadGEM2-ES and Eta-MIROC5) will increase by 2100. The results of this study will help producers in the semi-arid region of Brazil adopt specific agriculture systems aimed at mitigating greenhouse gas emissions.

Aeolian sandy soil in mining areas exhibits intense evaporation and poor water retention capacity. This study was designed to find a suitable biochar application method to improve soil water infiltration and minimize soil water evaporation for aeolian sand soil. Using the indoor soil column method, we studied the effects of three application patterns (A (0-20 cm was a mixed sample of mixed-based biochar and soil), B (0-10 cm was a mixed sample of mixed-based biochar and soil and 10-20 cm was soil), and C (0-10 cm was soil and 10-20 cm was a mixed sample of mixed-based biochar and soil)), four application amounts (0% (control, CK), 1%, 2%, and 4% of mixed-based biochar in dry soil), and two particle sizes (0.05-0.25 mm (S1) and <0.05 mm (S2)) of mixed-based biochar on water infiltration and evaporation of aeolian sandy soil. We separately used five infiltration models (the Philip, Kostiakov, Horton, USDA-NRCS (United States Department of Agriculture-Natural Resources Conservation Service), and Kostiakov-Lewis models) to fit cumulative infiltration and time. Compared with CK, the application of mixed-based biochar significantly reduced cumulative soil water infiltration. Under application patterns A, B, and C, the higher the application amount and the finer the particle size were, the lower the migration speed of the wetting front. With the same application amount, cumulative soil water infiltration under application pattern A was the lowest. Taking infiltration for 10 min as an example, the reductions of cumulative soil water infiltration under the treatments of A2%_(S2), A4%_(S1), A4%_(S2), A1%_(S1), C2%_(S1), and B1%_(S1) were higher than 30%, which met the requirements of loess soil hydraulic parameters suitable for plant growth. The five infiltration models well fitted the effects of the treatments of application pattern C and S1 particle size (R²>0.980), but the R² values of the Horton model exceeded 0.990 for all treatments (except for the treatment B2%_(S2)). Compared with CK, all other treatments reduced cumulative soil water infiltration, except for B4%_(S2). With the same application amount, cumulative soil water evaporation difference between application patterns A and B was small. Treatments of application pattern C and S1 particle size caused a larger reduction in cumulative soil water evaporation. The reductions in cumulative soil water evaporation under the treatments of C4%_(S1), C4%_(S2), C2%_(S1), and C2%_(S2) were over 15.00%. Therefore, applying 2% of mixed-based biochar with S1 particle size to the underlying layer (10-20 cm) could improve soil water infiltration while minimizing soil water evaporation. Moreover, application pattern was the main factor affecting soil water infiltration and evaporation. Further, there were interactions among the three influencing factors in the infiltration process (application amount×particle size with the most important interaction), while there were no interactions among them in the evaporation process. The results of this study could contribute to the rational application of mixed-based biochar in aeolian sandy soil and the resource utilization of urban and agricultural wastes in mining areas.

Comprehensive assessments of ecosystem services in environments under the influences of human activities and climate change are critical for sustainable regional ecosystem management. Therefore, integrated interdisciplinary modelling has become a major focus of ecosystem service assessment. In this study, we established a model that integrates land use/cover change (LUCC), climate change, and water retention services to evaluate the spatial and temporal variations of water retention services in the Loess Plateau of China in the historical period (2000-2015) and in the future (2020-2050). An improved Markov-Cellular Automata (Markov-CA) model was used to simulate land use/land cover patterns, and ArcGIS 10.2 software was used to simulate and assess water retention services from 2000 to 2050 under six combined scenarios, including three land use/land cover scenarios (historical scenario (HS), ecological protection scenario (EPS), and urban expansion scenario (UES)) and two climate change scenarios (RCP4.5 and RCP8.5, where RCP is the representative concentration pathway). LUCCs in the historical period (2000-2015) and in the future (2020-2050) are dominated by transformations among agricultural land, urban land and grassland. Urban land under UES increased significantly by 0.63×10³ km²/a, which was higher than the increase of urban land under HS and EPS. In the Loess Plateau, water yield decreased by 17.20×10⁶ mm and water retention increased by 0.09×10⁶ mm in the historical period (2000-2015), especially in the Interior drainage zone and its surrounding areas. In the future (2020-2050), the pixel means of water yield is higher under RCP4.5 scenario (96.63 mm) than under RCP8.5 scenario (95.46 mm), and the pixel means of water retention is higher under RCP4.5 scenario (1.95 mm) than under RCP8.5 scenario (1.38 mm). RCP4.5-EPS shows the highest total water retention capacity on the plateau scale among the six combined scenarios, with the value of 1.27×10⁶ mm. Ecological restoration projects in the Loess Plateau have enhanced soil and water retention. However, more attention needs to be paid not only to the simultaneous increase in water retention services and evapotranspiration but also to the type and layout of restored vegetation. Furthermore, urbanization needs to be controlled to prevent uncontrollable LUCCs and climate change. Our findings provide reference data for the regional water and land resources management and the sustainable development of socio-ecological systems in the Loess Plateau under LUCC and climate change scenarios.