Monitoring canopy growth and grain yield of paddy rice in South Korea by using the GRAMI model and high spatial resolution imagery

Monitoring crop conditions and forecasting crop yields are both important for assessing crop production and for determining appropriate agricultural management practices; however, remote sensing is limited by the resolution, timing, and coverage of satellite images, and crop modeling is limited in its application at regional scales. To resolve these issues, the Gramineae (GRAMI)-rice model, which utilizes remote sensing data, was used in an effort to combine the complementary techniques of remote sensing and crop modeling. The model was then investigated for its capability to monitor canopy growth and estimate the grain yield of rice (Oryza sativa), at both the field and the regional scales, by using remote sensing images with high spatial resolution. The field scale investigation was performed using unmanned aerial vehicle (UAV) images, and the regional-scale investigation was performed using RapidEye satellite images. Simulated grain yields at the field scale were not significantly different (p = 0.45, p = 0.27, and p = 0.52) from the corresponding measured grain yields according to paired t-tests (α = 0.05). The model’s projections of grain yield at the regional scale represented the spatial grain yield variation of the corresponding field conditions to within ±1 standard deviation. Therefore, based on mapping the growth and grain yield of rice at both field and regional scales of interest within coverages of a UAV or the RapidEye satellite, our results demonstrate the applicability of the GRAMI-rice model to the monitoring and prediction of rice growth and grain yield at different spatial scales. In addition, the GRAMI-rice model is capable of reproducing seasonal variations in rice growth and grain yield at different spatial scales.     

A workflow for Sustainable Development Goals indicators assessment based on high-resolution satellite data

ABSTRACT

For evaluating the progresses towards achieving the Sustainable Development Goals (SDGs), a global indicator framework was developed by the UN Inter-Agency and Expert Group on Sustainable Development Goals Indicators. In this paper, we propose an improved methodology and a set of workflows for calculating SDGs indicators. The main improvements consist of using moderate and high spatial resolution satellite data and state-of-the-art deep learning methodology for land cover classification and for assessing land productivity. Within the European Network for Observing our Changing Planet (ERA-PLANET), three SDGs indicators are calculated. In this research, harmonized Landsat and Sentinel-2 data are analyzed and used for land productivity analysis and yield assessment, as well as Landsat 8, Sentinel-2 and Sentinel-1 time series are utilized for crop mapping. We calculate for the whole territory of Ukraine SDG indicators: 15.1.1 – ‘Forest area as proportion of total land area’; 15.3.1 – ‘Proportion of land that is degraded over total land area’; and 2.4.1 – ‘Proportion of agricultural area under productive and sustainable agriculture’. Workflows for calculating these indicators were implemented in a Virtual Laboratory Platform. We conclude that newly available high-resolution remote sensing products can significantly improve our capacity to assess several SDGs indicators through dedicated workflows.     

Assessing gaps in irrigated agricultural productivity through satellite earth observations—A case study of the Fergana Valley,Central Asia

Improving crop area and/or crop yields in agricultural regions is one of the foremost scientific challenges for the next decades. This is especially true in irrigated areas because sustainable intensification of irrigated crop production is virtually the sole means to enhance food supply and contribute to meeting food demands of a growing population. Yet, irrigated crop production worldwide is suffering from soil degradation and salinity, reduced soil fertility, and water scarcity rendering the performance of irrigation schemes often below potential. On the other hand, the scope for improving irrigated agricultural productivity remains obscure also due to the lack of spatial data on agricultural production (e.g. crop acreage and yield). To fill this gap, satellite earth observations and a replicable methodology were used to estimate crop yields at the field level for the period 2010/2014 in the Fergana Valley, Central Asia, to understand the response of agricultural productivity to factors related to the irrigation and drainage infrastructure and environment. The results showed that cropping pattern, i.e. the presence or absence of multi-annual crop rotations, and spatial diversity of crops had the most persistent effects on crop yields across observation years suggesting the need for introducing sustainable cropping systems. On the other hand, areas with a lower crop diversity or abundance of crop rotation tended to have lower crop yields, with differences of partly more than one t/ha yield. It is argued that factors related to the infrastructure, for example, the distance of farms to the next settlement or the density of roads, had a persistent effect on crop yield dynamics over time. The improvement potential of cotton and wheat yields were estimated at 5%, compared to crop yields of farms in the direct vicinity of settlements or roads. In this study it is highlighted how remotely sensed estimates of crop production in combination with geospatial technologies provide a unique perspective that, when combined with field surveys, can support planners to identify management priorities for improving regional production and/or reducing environmental impacts.     

Monitoring of water stress in wheat using multispectral indices derived from Landsat-TM

Detection of crop water stress is crucial for efficient irrigation water management. Potential of Satellite data to provide spatial and temporal dynamics of crop growth conditions makes it possible to monitor crop water stress at regional level. This study was conducted in parts of western Uttar Pradesh and Haryana. Multi-temporal Landsat data were used for detecting wheat crop water stress using vegetation indices (VIs), viz. vegetation water stress index (VWSI) and land surface wetness index water stress factor (Ws_LSWI). The estimated water stress from satellite data-based VIs was validated by water stress factor (Ws) derived from flux-tower data. The study observed Ws_LSWI to be better index for water stress detection. The results indicated that Ws_LSWI was superior over other index showing RMSE = 0.12, R² = 0.65, whereas VWSI showed overestimated values with mean RD 4%.     

Geospatial Approach in Assessing Agro-Climatic Suitability of Soybean in Rainfed Agro-Ecosystem

The present study demonstrated the methodology to assess agro-climatic suitability of the soybean crop through integration of crop suitability based on FAO framework of land evaluation and biophysical (water limited) yield potential in the rainfed agro-ecosystem. A long term climatic database (1980–2003) was prepared to compute decadal rainfall and temperature variations of 13 IMD stations in part of Madhya Pradesh state. The climatic database was used in soil water balance software–BUDGET to compute crop specific length of growing period (LGP) and biophysical production potential such as water limited crop yield potential of each soil types for soybean crop. Water limited crop yield potential of soils were found to be varied from 33 to 100 and LGP ranged from 65 to 180 days in the area. FAO based land suitability was analyzed in association with the water limited yield potential for better appraisal of land potential and assess their suitability in rainfed area. FAO based land suitability indicated 2.45 % area as highly suitable and 57.49 % area as moderately suitable. However, integration of water limited crop yield potential with FAO based land suitability lead to agro-climatic suitability analysis indicated 17.60 % and 40.03 % area, respectively as highly suitable and moderately suitable. FAO based land evaluation showed 88.13 % of plains as moderately suitable whereas agro-climatic suitability indicated only 47.79 %. Agro-climatic suitability analysis revealed undulating plateau and undulating plains as most suitable for soybean crop.     

Remote sensing-based global crop monitoring: experiences with China's CropWatch system

Monitoring the production of main agricultural crops is important to predict and prepare for disruptions in food supply and fluctuations in global crop market prices. China's global crop-monitoring system (CropWatch) uses remote sensing data combined with selected field data to determine key crop production indicators: crop acreage, yield and production, crop condition, cropping intensity, crop-planting proportion, total food availability, and the status and severity of droughts. Results are combined to analyze the balance between supply and demand for various food crops and if needed provide early warning about possible food shortages. CropWatch data processing is highly automated and the resulting products provide new kinds of inputs for food security assessments. This paper presents a comprehensive overview of CropWatch as a remote sensing-based system, describing its structure, components, and monitoring approaches. The paper also presents examples of monitoring results and discusses the strengths and limitations of the CropWatch approach, as well as a comparison with other global crop-monitoring systems.     

Assessing the supplementary irrigation for improving crop productivity in water stress region using spatial hydrological model

Water stress during crop cultivation due to inconsistent rainfall is a common phenomenon in maize growing area of Shanmuganadi watershed, located in the semi-arid region of southern peninsular India. The objective is to estimate the supplementary irrigation required to improve the crop productivity during water stress period. Spatial hydrological model, Soil and Water Assessment Tool, has been applied to simulate the watershed hydrology and crop growth for rabi season (October–February) considering the rainfed and irrigated scenarios. The average water stress days of rainfed maize was 60 days with yield of 1.6 t/ha. Irrigated maize with supplementary irrigation of 93–126 mm was resulted in improved yield of 3.8 t/ha with 28 water stress days. The results also suggest that supplemental irrigation can be obtained from groundwater reserves and by adopting early sowing strategy can provide opportunities for improving water productivity in rainfed farming.     

Regional yield estimation for winter wheat with MODIS-NDVI data in Shandong, China

The significance of crop yield estimation is well known in agricultural management and policy development at regional and national levels. The primary objective of this study was to test the suitability of the method, depending on predicted crop production, to estimate crop yield with a MODIS-NDVI-based model on a regional scale. In this paper, MODIS-NDVI data, with a 250 m resolution, was used to estimate the winter wheat (Triticum aestivum L.) yield in one of the main winter-wheat-growing regions. Our study region is located in Jining, Shandong Province. In order to improve the quality of remote sensing data and the accuracy of yield prediction, especially to eliminate the cloud-contaminated data and abnormal data in the MODIS-NDVI series, the Savitzky–Golay filter was applied to smooth the 10-day NDVI data. The spatial accumulation of NDVI at the county level was used to test its relationship with winter wheat production in the study area. A linear regressive relationship between the spatial accumulation of NDVI and the production of winter wheat was established using a stepwise regression method. The average yield was derived from predicted production divided by the growing acreage of winter wheat on a county level. Finally, the results were validated by the ground survey data, and the errors were compared with the errors of agro-climate models. The results showed that the relative errors of the predicted yield using MODIS-NDVI are between −4.62% and 5.40% and that whole RMSE was 214.16 kg ha⁻¹ lower than the RMSE (233.35 kg ha⁻¹) of agro-climate models in this study region. A good predicted yield data of winter wheat could be got about 40 days ahead of harvest time, i.e. at the booting-heading stage of winter wheat. The method suggested in this paper was good for predicting regional winter wheat production and yield estimation.     

Yield Forecasting of Spring Maize Using Remote Sensing and Crop Modeling in Faisalabad-Punjab Pakistan

Real time, accurate and reliable estimation of maize yield is valuable to policy makers in decision making. The current study was planned for yield estimation of spring maize using remote sensing and crop modeling. In crop modeling, the CERES-Maize model was calibrated and evaluated with the field experiment data and after calibration and evaluation, this model was used to forecast maize yield. A Field survey of 64 farm was also conducted in Faisalabad to collect data on initial field conditions and crop management data. These data were used to forecast maize yield using crop model at farmers’ field. While in remote sensing, peak season Landsat 8 images were classified for landcover classification using machine learning algorithm. After classification, time series normalized difference vegetation index (NDVI) and land surface temperature (LST) of the surveyed 64 farms were calculated. Principle component analysis were run to correlate the indicators with maize yield. The selected LSTs and NDVIs were used to develop yield forecasting equations using least absolute shrinkage and selection operator (LASSO) regression. Calibrated and evaluated results of CERES-Maize showed the mean absolute % error (MAPE) of 0.35–6.71% for all recorded variables. In remote sensing all machine learning algorithms showed the accuracy greater the 90%, however support vector machine (SVM-radial basis) showed the higher accuracy of 97%, that was used for classification of maize area. The accuracy of area estimated through SVM-radial basis was 91%, when validated with crop reporting service. Yield forecasting results of crop model were precise with RMSE of 255 kg ha^?1, while remote sensing showed the RMSE of 397 kg ha^?1. Overall strength of relationship between estimated and actual grain yields were good with R² of 0.94 in both techniques. For regional yield forecasting remote sensing could be used due greater advantages of less input dataset and if focus is to assess specific stress, and interaction of plant genetics to soil and environmental conditions than crop model is very useful tool.     

Developing a satellite-based combined drought indicator to monitor agricultural drought: a case study for Ethiopia

Developing a robust drought monitoring tool is vital to mitigate the adverse impacts of drought. A drought monitoring system that integrates multiple agrometeorological variables into a single drought indicator is lacking in areas such as Ethiopia, which is extremely susceptible to this natural hazard. The overarching goal of this study is to develop a combined drought indicator (CDI-E) to monitor the spatial and temporal extents of historic agricultural drought events in Ethiopia. The CDI-E was developed by combining four satellite-based agrometeorological input parameters – the Standardized Precipitation Index (SPI), Land Surface Temperature (LST) anomaly, Standardized Normalized Difference Vegetation Index (stdNDVI) and Soil Moisture (SM) anomaly – for the period from 2001 to 2015. The method used to combine these indices is based on a quantitative approach that assigns a weight to each input parameter using Principal Component Analysis (PCA). The CDI-E results were evaluated using satellite-based gridded rainfall (3-month SPI) and crop yield data for 36 intra-country crop growing zones for a 15-year period (2001 to 2015). The evaluation was carried out for the main rainfall season, Kiremt (June-September), and the short rainfall season, Belg (February-May). The results showed that moderate to severe droughts were detected by the CDI-E across the food insecure regions reported by FEWS NET during Kiremt and Belg rainfall seasons. Relatively higher correlation coefficient values (r > 0.65) were obtained when CDI-E was compared with the 3-month SPI across the majority of Ethiopia. The spatial correlation analyses of CDI-E and cereal crop yields showed relatively good correlations (r > 0.5) in some of the crop growing zones in the northern, eastern and southwestern parts of the country. The CDI-E generally mapped the spatial and temporal patterns of historic drought and non-drought years and hence the CDI-E could potentially be used to develop an agricultural drought monitoring and early warning system in Ethiopia. Moreover, decision makers and donors may potentially use CDI-E to more accurately monitor crop yields across the food-insecure regions in Ethiopia.     

GEE环境下的玉米低温冷害损失快速评估

大范围、及时、准确的灾害损失评估与制图对防灾减灾、农业保险和粮食安全等至关重要。针对传统灾害损失评估方法空间尺度单一、泛化能力差、时效性低,可操作性弱等问题,本文建立了一种遥感产品耦合作物模型的多尺度的灾害损失评估方法MDLA （a Multiscale Disaster Loss Assessment）。该方法利用作物模型的多情景模拟产生大量的灾害样本,结合对应日期的遥感指标构建灾害脆弱性模型,依托Google Earth Engine（GEE）平台将其应用到高分辨率遥感影像和格点灾害指标进行逐象元评估。以鄂伦春自治旗玉米为例,基于精细校准的CERES-Maize模型的模拟,利用两个生长季窗口的LAI和冷积温（CDD）建立统计模型来刻画低温对最终产量的影响,结合Sentinel-2数据逐格点计算完成高精度损失制图。结果显示,校准后的CERES-Maize模拟物候和产量的NRMSE 分别为3.3%和8.9%。冷害情景模拟结果表明不同类型和生育期的低温冷害对玉米产量的影响不尽相同,其中生长峰值期（出苗—吐丝和吐丝—灌浆）最为敏感。回代检验显示,MDLA方法估算精度为11.4%,与历史冷害年份的实际损失相吻合。经评估,鄂伦春2018-08-09的冷害导致玉米减产23.7%,受灾面积1.86×10⁴ ha,其中高海拔地区损失较重（减产率>25%）,低温冷害对该区玉米生产构成了严重的威胁。与现有的统计回归、作物模型模拟以及同化等技术相比,其优势在于：（1）结合遥感观测和作物模型模拟技术能更好地刻画了灾害对产量的影响过程;（2）利用GEE平台快速处理海量遥感数据,提高了灾害损失评估的时效性;（3）不受地面实测数据的限制,易操作,可实现动态、多尺度（象元、田块、村,县等）的损失评估,这为防灾减损、维持粮食丰产稳产提供了保障,也为农业保险的业务化运行提供了思路。     

A simple assessment approach for winter wheat loss risk impacted by water stress

In this study, an empirical assessment approach for the risk of crop loss due to water stress was developed and used to evaluate the risk of winter wheat loss in China, the United States, Germany, France and the United Kingdom. We combined statistical and remote sensing data on crop yields with climate data and cropland distribution to model the effect of water stress from 1982 to 2011. The average value of winter wheat loss due to water stress for the three European countries was about ?931 kg/ha, which was higher than that in China (?570 kg/ha) and the United States (?367 kg/ha). Our study has important implications for the operational assessment of crop loss risk at a country or regional scale. Future studies should focus on using higher spatial resolution remote sensing data, combining actual evapotranspiration to estimate water stress, improving the method for downscaling of statistical crop yield data and establishing more sophisticated zoning methods.     

基于WOFOST模型与UAV数据的玉米生长后期地上生物量估算

地上生物量能够有效反映作物的生长状态,其信息的实时估算对产量预测和农田生产管理都有重要意义。作物生长模型因其详尽的生理生化基础和对生长过程数字化描述能力,成为生物量估算的理想模型。近年来,研究人员利用数据同化算法将时间序列遥感数据同化到作物生长模型中,实现了作物模型由基于气象站的点模拟到区域尺度面模拟的外推,使生物量模拟结果同时具备大范围和机理性两个方面的特点。这一模式下,时间序列的遥感数据质量将对生物量模拟精度产生直接影响,作物生长后期受到光谱饱和的影响,遥感数据的作物冠层信息获取能力会出现明显下降,因此有必要对该阶段遥感数据和作物模型的结合方式进行优化,提升生物量模拟精度。本文针对东北地区春玉米生物量遥感估算存在的问题,提出了利用WOFOST作物模型结合无人机（UAV）遥感数据实现作物生长后期生物量准确估算的新思路。新思路首先利用多光谱遥感数据获取WOFOST模型具备较高空间异质性的土壤速效养分参数以提升模型的空间信息模拟能力,使其能在一定程度上摆脱点尺度模拟的限制。同时,结合集合卡尔曼滤波算法将生长前期无人机（UAV）遥感数据同化到模型中,以缩短模型单独运行时间,减少模型运行过程中的参数误差累积,实现无遥感数据参与下的短期作物生长模拟,并输出生长后期相应的生物量模拟结果。最后,本文利用地面实测数据对新方法的生物量模拟精度进行了评价。结果表明,与全生育期数据同化相比,新方法的生物量估算精度有了明显的提升（全生育期同化：R² = 0.45,RMSE = 4254.30 kg/ha;新方法：R²= 0.86,RMSE = 2216.79 kg/ha）。     

Integration of maximum crop response with machine learning regression model to timely estimate crop yield

Timely and reliable estimation of regional crop yield is a vital component of food security assessment, especially in developing regions. The traditional crop forecasting methods need ample time and labor to collect and process field data to release official yield reports. Satellite remote sensing data is considered a cost-effective and accurate way of predicting crop yield at pixel-level. In this study, maximum Enhanced Vegetation Index (EVI) during the crop-growing season was integrated with Machine Learning Regression (MLR) models to estimate wheat and rice yields in Pakistan's Punjab province. Five MLR models were compared using a fivefold cross-validation method for their predictive accuracy. The study results revealed that the regression model based on the Gaussian process outperformed over other models. The best performing model attained coefficient of determination (R2), Root Mean Square Error (RMSE, t/ ha), and Mean Absolute Error (MAE, t/ha) of 0.75, 0.281, and 0.236 for wheat; 0.68, 0.112, and 0.091 for rice, respectively. The proposed method made it feasible to predict wheat and rice 6– 8 weeks before the harvest. The early prediction of crop yield and its spatial distribution in the region can help formulate efficient agricultural policies for sustainable social, environmental, and economic progress.     

农作物长势综合遥感监测方法

作物收获之前进行大范围作物生长状况评价 ,可以尽早的获得有关作物产量信息。介绍了中国农情遥感监测系统的综合作物长势监测方法。以遥感数据标准化处理、云标识、云污染去除和非耕地去除为基础 ,生成质量一致的遥感数据产品集 ,提取区域作物生长过程。作物长势监测分为实时作物长势监测和作物生长趋势分析。实时的作物长势监测可以定性和定量地在空间上分析作物生长状况 ,分级显示作物生长状况 ,分区域统计水田和旱地中不同长势占的比重。作物生长趋势分析可以进行年际间的生长过程对比 ,从时间轴上反映作物持续生长的差异性 ,统计全国、主产区、省和区划单元 4个尺度的耕地、水田、旱地作物生长过程曲线年际间差异 ,从而为早期的产量预测提供信息。通过处理流程的系统化 ,建设了运行化的作物长势遥感监测分析系统 ,为用户构建了综合的作物实时生长状况 ,苗情的生长趋势分析环境。同时可以依据野外地面实测信息对遥感监测结果进行标定和检验。 1998年以来 ,系统在满足日常运行的前提下 ,技术方法逐渐改进和完善 ,监测结果的精度和可靠性不断得到提高。     

Estimating inter-annual variability in winter wheat sowing dates from satellite time series in Camargue,France

Crop simulation models are commonly used to forecast the performance of cropping systems under different hypotheses of change. Their use on a regional scale is generally constrained, however, by a lack of information on the spatial and temporal variability of environment-related input variables (e.g., soil) and agricultural practices (e.g., sowing dates) that influence crop yields. Satellite remote sensing data can shed light on such variability by providing timely information on crop dynamics and conditions over large areas. This paper proposes a method for analyzing time series of MODIS satellite data in order to estimate the inter-annual variability of winter wheat sowing dates. A rule-based method was developed to automatically identify a reliable sample of winter wheat field time series, and to infer the corresponding sowing dates. The method was designed for a case study in the Camargue region (France), where winter wheat is characterized by vernalization, as in other temperate regions. The detection criteria were chosen on the grounds of agronomic expertise and by analyzing high-confidence time-series vegetation index profiles for winter wheat. This automatic method identified the target crop on more than 56% (four-year average) of the cultivated areas, with low commission errors (11%). It also captured the seasonal variability in sowing dates with errors of ±8 and ±16 days in 46% and 66% of cases, respectively. Extending the analysis to the years 2002–2012 showed that sowing in the Camargue was usually done on or around November 1st (±4 days). Comparing inter-annual sowing date variability with the main local agro-climatic drivers showed that the type of preceding crop and the weather conditions during the summer season before the wheat sowing had a prominent role in influencing winter wheat sowing dates.     

Development of a GIS-based decision support system for assessing land use status

The objectives of this study are to assess land suitability and to predict the spatial and temporal changes in land use types (LUTs) by using GIS-based land use management decision support system. A GIS database with data on climate, topography, soil characteristic, irrigation condition, fertilizer application, and special socioeconomic activities has been developed and used for the evaluation of land productivity for different crops by integrating with a crop growth model—the erosion productivity impact calculator (EPIC). International food policy simulation model (IFPSIM) is also embedded into GIS for the predictions of how crop demands and crop market prices will change under alternative policy scenarios. An inference engine (IE) including land use choice model is developed to illustrate land use choice behavior based on logit models, which allows to analyze how diversified factors ranging from climate changes, crop price changes to land management changes can effect the distribution of agricultural land use. A test for integrated simulation is taken in each 0.1o by 0.1o grid cell to predict the change of agricultural land use types at global level. Global land use changes are simulated from 1992 to 2050.     

Development of a GIS-based decision support system for assessing land use status

IntroductionAgriculturallandusepatternsandtheirchangesaretightlyrelatedwithagriculturepolicyandfoodsecurityissuesundergrowingfooddemand,assess mentofglobalclimatechangeimpactsonagricul ture,environmentalissuesduetotheintensificationofagriculturallandusessuchaswaterpollution,soildegradation,andrecentlywaterscarcityissues.Soasustainableandholisticplanningandmanage mentoflandresourcesshouldcombineallthesere latedinformationwithefficienttoolsforassessmentandevaluationinordertopermitbroad ,interact…     

Yield Estimation Model and Water Productivity of Wheat Crop (Triticum aestivum) in an Irrigation Command Using Remote Sensing and GIS

Crop yield estimation has an important role on economy development and its accuracy and speed influence yield price and helps in deciding the excess or deficit production conditions. The water productivity evaluates the irrigation command through water use efficiency (WUE). Remote sensing (RS) and geographical information system (GIS) techniques were used for crop yield and water productivity estimation of wheat crop (Triticum aestivum) grown in Tarafeni South Main Canal (TSMC) irrigation command of West Bengal State in India. One IRS P6 image and four wide field sensor (WiFS) images for different months of winter season were used to determine the Normalized Difference Vegetation Index (NDVI) and Soil Adjusted Vegetation Index (SAVI) for area under wheat crop. The temporally and spatially distributed spectral growth profile and AREASUM of NDVI (A_NDVI) and SAVI (A_SAVI) with time after sowing of wheat crop were developed and correlated with actual crop yield of wheat (Y_act). The developed relationships between A_SAVI and Y_act resulted high correlation in comparison to that of A_NDVI. Using the developed model the RS based wheat yield (Y_RS) predicted from A_SAVI varied on entire TSMC irrigation command from 22.67 to 33.13 q ha⁻¹ respectively, which gave an average yield of 26.50 q ha⁻¹. The RS generated yield based water use efficiency (WUE_YRS) for water supplied from canal of TSMC irrigation command was found to be 6.69 kg ha⁻¹ mm⁻¹.