Exploring the use of Sentinel-2 datasets and environmental variables to model wheat crop yield in smallholder arid and semi-arid farming systems.

Low levels of agricultural productivity are associated with the persistence of food insecurity, poverty, and other socio-economic stresses. Mapping and monitoring agricultural dynamics and production in real-time at high spatial resolution are essential for ensuring food security and shaping policy interventions. However, an accurate yield estimation might be challenging in some arid and semi-arid regions since input datasets are generally scarce, and access is restricted due to security challenges. This work examines how well Sentinel-2 satellite sensor-derived data, topographic and climatic variables, can be used as covariates to accurately model and predict wheat crop yield at the farm level using statistical models in low data settings of arid and semi-arid regions, using Sulaimani governorate in Iraq as an example. We developed a covariate selection procedure that assessed the correlations between the covariates and their relationships with wheat crop yield. Potential non-linear relationships were investigated in the latter case using regression splines. In the absence of substantial non-linear relationships between the covariates and crop yield, and residual spatial autocorrelation, we fitted a Bayesian multiple linear regression model to model and predict crop yield at 10 m resolution. Out of the covariates tested, our results showed significant relationships between crop yield and mean cumulative NDVI during the growing season, mean elevation, mean end of the season, mean maximum temperature and mean the start of the season at the farm level. For in-sample prediction, we estimated an R2 value of 51 % for the model, whereas for out-of-sample prediction, this was 41 %, both of which indicate reasonable predictive performance. The calculated root-mean-square error for out-of-sample prediction was 69.80, which is less than the standard deviation of 89.23 for crop yield, further showing that the model performed well by reducing prediction variability. Besides crop yield estimates, the model produced uncertainty metrics at 10 m resolution. Overall, this study showed that Sentinel-2 data can be valuable for upscaling field measurement of crop yield in arid and semi-arid regions. In addition, the environmental covariates can strengthen the model predictive power. The method may be applicable in other areas with similar environments, particularly in conflict zones, to increase the availability of agricultural statistics.

Using gridded population and quadtree sampling units to support survey sample design in low-income settings.

BACKGROUND: Household surveys are the main source of demographic, health and socio-economic data in low- and middle-income countries (LMICs). To conduct such a survey, census population information mapped into enumeration areas (EAs) typically serves a sampling frame from which to generate a random sample. However, the use of census information to generate this sample frame can be problematic as in many LMIC contexts, such data are often outdated or incomplete, potentially introducing coverage issues into the sample frame. Increasingly, where census data are outdated or unavailable, modelled population datasets in the gridded form are being used to create household survey sampling frames. METHODS: Previously this process was done by either sampling from a set of the uniform grid cells (UGC) which are then manually subdivided to achieve the desired population size, or by sampling very small grid cells then aggregating cells into larger units to achieve a minimum population per survey cluster. The former approach is time and resource-intensive as well as results in substantial heterogeneity in the output sampling units, while the latter can complicate the calculation of unbiased sampling weights. Using the context of Somalia, which has not had a full census since 1987, we implemented a quadtree algorithm for the first time to create a population sampling frame. The approach uses gridded population estimates and it is based on the idea of a quadtree decomposition in which an area successively subdivided into four equal size quadrants, until the content of each quadrant is homogenous. RESULTS: The quadtree approach used here produced much more homogeneous sampling units than the UGC (1 × 1 km and 3 × 3 km) approach. At the national and pre-war regional scale, the standard deviation and coefficient of variation, as indications of homogeneity, were calculated for the output sampling units using quadtree and UGC 1 × 1 km and 3 × 3 km approaches to create the sampling frame and the results showed outstanding performance for quadtree approach. CONCLUSION: Our approach reduces the manual burden of manually subdividing UGC into highly populated areas, while allowing for correct calculation of sampling weights. The algorithm produces a relatively homogenous population counts within the sampling units, reducing the variation in the weights and improving the precision of the resulting estimates. Furthermore, a protocol of creating approximately equal-sized blocks and using tablets for randomized selection of a household in each block mitigated potential selection bias by enumerators. The approach shows labour, time and cost-saving and points to the potential use in wider contexts.

Forecasting wheat and barley crop production in arid and semi-arid regions using remotely sensed primary productivity and crop phenology: A case study in Iraq.

Crop production and yield estimation using remotely sensed data have been studied widely, but such information is generally scarce in arid and semi-arid regions. In these regions, inter-annual variation in climatic factors (such as rainfall) combined with anthropogenic factors (such as civil war) pose major risks to food security. Thus, an operational crop production estimation and forecasting system is required to help decision-makers to make early estimates of potential food availability. Data from NASA's MODIS with official crop statistics were combined to develop an empirical regression-based model to forecast winter wheat and barley production in Iraq. The study explores remotely sensed indices representing crop productivity over the crop growing season to find the optimal correlation with crop production. The potential of three different remotely sensed indices, and information related to the phenology of crops, for forecasting crop production at the governorate level was tested and their results were validated using the leave-one-year-out approach. Despite testing several methodological approaches, and extensive spatio-temporal analysis, this paper depicts the difficulty in estimating crop yield on an annual base using current satellite low-resolution data. However, more precise estimates of crop production were possible. The result of the current research implies that the date of the maximum vegetation index (VI) offered the most accurate forecast of crop production with an average R2=0.70 compared to the date of MODIS EVI (Avg R2=0.68) and a NPP (Avg R2=0.66). When winter wheat and barley production were forecasted using NDVI, EVI and NPP and compared to official statistics, the relative error ranged from -20 to 20%, -45 to 28% and -48 to 22%, respectively. The research indicated that remotely sensed indices could characterize and forecast crop production more accurately than simple cropping area, which was treated as a null model against which to evaluate the proposed approach.