# Mapping Climate Parameters over the Territory of Botswana Using GMT and Gridded Surface Data from TerraClimate

Polina Lemenkova

► To cite this version:

Polina Lemenkova. Mapping Climate Parameters over the Territory of Botswana Using GMT and Gridded Surface Data from TerraClimate. ISPRS International Journal of Geo-Information, 2022, 11 (9), pp.473. 10.3390/ijgi11090473. hal-03765979

HAL Id: hal-03765979

https://hal.science/hal-03765979v1

Submitted on 31 Aug 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
Article

# Mapping Climate Parameters over the Territory of Botswana Using GMT and Gridded Surface Data from TerraClimate

Polina Lemenkova ID

Laboratory of Image Synthesis and Analysis (LISA), École polytechnique de Bruxelles, Brussels Faculty of Engineering, Campus de Solbosch CP131/3, Université Libre de Bruxelles (ULB), Building L, Avenue Franklin D. Roosevelt 50, B-1050 Brussels, Belgium; polina.lemenkova@ulb.be or pauline.lemenkova@gmail.com; Tel.: +32-471860459

**Abstract:** This articles presents a new series of maps showing the climate and environmental variability of Botswana. Situated in southern Africa, Botswana has an arid to semi-arid climate, which significantly varies in its different regions: Kalahari Desert, Makgadikgadi Pan and Okavango Delta. While desert regions are prone to droughts and periods of extreme heat during the summer months, other regions experience heavy downpours, as well as episodic and unpredictable rains that affect agricultural activities. Such climatic variations affect social and economic aspects of life in Botswana. This study aimed to visualise the non-linear correlations between the topography and climate setting at the country's scale. Variables included T °C min, T °C max, precipitation, soil moisture, evaporation, transpiration (PET and AET), downward surface shortwave radiation, vapour pressure and vapour pressure deficit (VPD), wind speed and Palmer Drought Severity Index (PDSI). The dataset was taken from the TerraClimate source and GEBCO for topographic mapping. The mapping approach included the use of Generic Mapping Tools (GMT), a console-based scripting toolset, which enables the use of a scripting method of automated mapping. Several GMT modules were used to derive a set of climate parameters for Botswana. The data were supplemented with the adjusted cartographic elements and inspected by the Geospatial Data Abstraction Library (GDAL). The PDSI in Botswana in 2018 shows stepwise variation with seven areas of drought: (1) -3.7 to -2.2. (extreme); (2) -2.2 to -0.8 (strong, southern Kalahari); (3) -0.8 to 0.7 (significant, central Kalahari); (4) 0.7 to 2.1 (moderate); (5) 2.1 to 3.5 (lesser); (6) 3.5 to 4.9 (low); (7) 4.9 to 6.4 (least). The VPD has a general trend towards the south-western region (Kalahari Desert, up to 3.3), while it is lower in the north-eastern region of Botswana (up to 1.4). Other values vary respectively, as demonstrated in the presented 12 maps of climate and environmental inventory in Botswana.

**Keywords:** Africa; Botswana; cartography; climate; computer science; data science; drought; GMT; mapping; programming

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Lemenkova, P. Mapping Climate Parameters over the Territory of Botswana Using GMT and Gridded Surface Data from TerraClimate. *ISPRS Int. J. Geo-Inf.* **2022**, *11*, 473. https://doi.org/10.3390/ijgi11090473

Academic Editor: Godwin Yeboah
and Wolfgang Kainz

Received: 26 June 2022

Accepted: 28 August 2022

Published: 31 August 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2022 by the author. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

## 1. Introduction

Mapping climate parameters using advanced cartographic methods is presently the subject of intense research. Approaches to image processing [1,2], climate data modelling [3] and landscape studies [4] have largely converged on the requirements for the importance of cartographic visualisation for a proper environmental analysis for feature extraction and evaluation. While standard methods of handling spatial data such as Geographic Information Systems (GIS) are popular for climate and environmental mapping [5–7], superior performance has been demonstrated by incorporating advanced scripting solutions for complex cartographic tasks. A large variety of spatial data processing paradigms have been applied to the problem of cartographic visualisation, including spatial analysis, data management and simulation, data retrieval from the remote sensing sources, pattern recognition and data-driven modelling.

One of the most popular GIS-based modelling tools for cartographic data handling is presented by the commercial software ArcGIS and its open source analogue QGIS. These
software structure spatial data as visible layers in the Graphical User Interface (GUI) where layers represent thematically structured and semantically close data types. Such structure enables connections between spatial variables through the analysis of attributive tables where metadata are stored in GIS. Thus, the cartographic variables are kept within each layer so that the parameters of the data are conditionally independent given the semantic connectivity and identical projection. The factorial nature of these data enables GIS-based mapping, which forms the basis for the majority of existing maps.

However, both as a data visualisation and mapping tool and especially as a method for spatial data modelling, GIS-based methods have shown certain drawbacks. The data processing in GIS is mostly performed using active interaction of the user for traditional mapping workflow where many handmade operations needed to be performed manually. Besides a time-consuming workflow for mapping, another drawback of the GIS methods is that, unlike scripting methods, the complete workflow processing is not automated. At the same time, automated methods in cartography significantly increase the accuracy and precision of plotting. As a result, the usual strategies of mapping largely benefit from applying the programming approach in cartography.

This paper presents the use of scripting cartographic toolset Generic Mapping Tools (GMT) [8] for a thematic mapping of Botswana. The GMT differs principally from the traditional GIS approaches. It is defined as a set of commands controlled using embedded programming language, which defines each cartographic parameters using a line of code. In such a way, the cartographic processing is being performed by writing and executing a script, without the need for GUI and the layer-based paradigm of GIS. Such a performance enables us to finely tune the selected parameters of cartographic plot.

As a technical instrument for mapping, the GMT is interesting in that, like programming languages, e.g., Python or R, its syntax encodes the conditional parameters of visualised objects using built-in language with commands and options while simultaneously maintaining the simple independent structure of map which does not require creating layers and a GIS project as such. This underlies the efficient performance of GMT as against the GIS workflow. Moreover, the GMT effectiveness is accomplished by the flexibility of data formats, integration with a variety of tools, including Geospatial Data Abstraction Library (GDAL) and Python and exploiting multi-format data through straightforward importing and reformatting.

In this study, the GMT was used for mapping Botswana using a set of climate and environmental variables. Botswana is a country situated in southern Africa (Figure 1), where drought is a reoccurring natural hazard [9]. The interest in the occurrence and distribution of droughts and climatic extremes in Botswana has been driven by the need to analyse possible consequences from droughts that affect both social and environmental aspects of the country. Such issues include, for instance, range degradation [10] and the evolution and adaptation of agro-pastoral livelihoods to recurrent droughts [11]. Moreover, land changes results in exposed dryland agricultural landscapes, which in turn involves expanding bare lands and rocky surfaces [12]. Biochemical analysis includes the evaluation of the effects from evapotranspiration on the evolution of dissolved inorganic carbon in the north-western Botswana [13]. Other case includes climate assessment and its effects on vegetation [14].

A variety of the environmental studies in Botswana are focused on the analysis of rainfall repeatability, intensity, and consequences in different regions of the country. Such and similar studies often include the modelling of the meteorological data, which enables both retrospective assessment and prognostic modelling of climate change variations to be performed, aiming to predict and mitigate possible effects from climatic extremes. For instance, previous study [15] performed a retrospective analysis of climate effects in the arid environment of Botswana to determine changes in rainfall series and to identify correlations between rainfall and temperature. A prognosis estimating the beginning and end of summer rainfalls may highlight regional uncertainties in rainfalls [16]. Rainfall
variability and trends for analysis of their effects on vulnerable dryland agriculture through
climate variation are important for climate risk management in Botswana [17].

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Topographic map of Botswana. Mapping: GMT. Source: author.</figcaption></figure>

Mapping climate parameters is a challenging task with many applications, such as
environmental risk assessment, agricultural monitoring, and natural resources manage-
ment. A key ingredient to cartographic data visualisation is the design of an effective
workflow and the selection of software tools. A traditional approach to this task is to apply
Geographic Information System (GIS) for spatial data processing. Mapping is then per-
formed by conventional GIS programs with existing Graphical User Interface (GUI). GIS is
a commonly accepted conceptual tool in cartographic data processing, widely documented
in the existing works [18–20]. Existing cartographic methods following the GIS approach
may make use of well-known methods of data modelling and visualisation, as reflected in
a variety of reports [21–25].

To progress beyond the traditional methods of GIS, the current investigation presents
an alternative approach to cartographic data processing using console-based scripts of
GMT [8]. The GMT utilises a modular system of GMT which automatically decomposes
the complete cartographic workflow into a series of continuous piecewise components
of scripts, i.e., lines of code written in GMT syntax. In such a way, the script handles the
complete process of spatial data from data formatting and modelling to generating a final
layout, presenting a map, and converting the output data. In view of this, in this paper we study the benefits of exploiting the GMT scripting toolset for mapping environmental and climate parameters in Botswana.

The goal of the presented research is to demonstrate the functionality of the console-based scripting in GMT for climate and environmental mapping of Botswana. To contribute to the existing climate studies of Botswana, this study presented a series of the new maps based on the open datasets. The maps are made by command line which increased the speed and accuracy of mapping and represented spatial data of various complexity.

The general objective was to map the variability of environmental and climate parameters of the country using the TerraClimate dataset. The variables in this dataset are designed such that one can evaluate the GMT approach with respect to various data dimensions and various parameters over the study area. Specifically for the territory of Botswana, while mapping the TerraClimate dataset, we address handling and visualising the diverse extent of values in climate parameters, from small to large ones. This is important, since difficulties in spatial data processing arise when data have contrasting dimensions, e.g., temperature, topography, and PDSI values differ and various colour scales for proper mapping and effective cartographic visualisation are required, which was managed using the GMT solutions.

The specific objective was to support the environmental and climate analysis in Botswana through the advanced cartographic visualisation. It involves the analysis of the distribution of temperature extremes in Botswana, occurrence of drought events, variation in the potential and actual evapotranspiration visualised on the maps. The correlations between the soil moisture and precipitation, climate parameters and topographic landforms can be considered for land management [26–29], crop and vegetation monitoring [30,31], meteorological evaluation [32], nature conservation, and sustainable development in Botswana [33,34]. In light of the above, the focus of the presented research is on the cartographic application of novel programming methods for climate and environmental assessment of Botswana, which includes the following objectives:

i To produce a series of climatic maps of Botswana and their associated topographic, environmental, and climate attributes, referenced to the regional and local geospatial geographic scale of the country;

ii To perform script-based mapping techniques, which present a new cartographic methodology. This includes the demonstrating of the programming approach of GMT, which differs from the traditional GIS, because it integrates the cartographic methods and the programming syntax;

iii To apply the open source datasets of GEBCO and TerraClimate with a spatial extent of Botswana, aimed at addressing and visualising the occurrence, trends and strength of different climate variables associated with environmental factors;

iv To visualise the climate parameters of Botswana: temperature extremes $T_{min}/T_{max}$, precipitation, soil moisture, potential evapotranspiration (PET), actual evapotranspiration (AET), downward surface shortwave radiation, vapour pressure (VP), vapour pressure deficit (VPD), wind speed, Palmer Drought Severity Index (PDSI).

The automated methods and independent modules of GMT were used to process climate data showing mapping climate and environmental variables of Botswana. This included reference to the cartographic workflow of GMT, described in the existing technical literature with examples of using the program [35–39]. Using a scripting approach, this paper presents an alternative cartographic algorithm for the stable and accurate mapping and visualisation of the climate parameters over Botswana in the presence of various nonlinear effects from topography and climate on the environmental setting of the country, following the examples in the existing studies [40,41].
## 2. Study Area

### 2.1. Relief

The topography of Botswana has its highest points along the south-east regions of the country, in the Kalahari Basin. In the rest of the regions, the landscapes are levelled, which creates unique climate setting. The territory of the country is covered by a diverse range of landscapes, including a complex mosaic of grass, sands, and savanna [42]. The Okavango Delta is presented by the alluvial fan of uniform topographic gradient with complex patterns of vegetation, hydrology, and sedimentation [43].

A notable topographic feature of Botswana is the Kalahari Desert, which occupies up to 70% of the territory of country. The Kalahari Desert is mostly covered by red sands and sandstones with quartzose sediments and sandstones. Such biochemical composition makes it prone to chemical weathering, producing high volumes of clay minerals [44]. Together with the extreme degrees of temperature and occasional fires, the geologic and soil setting of the Kalahari Desert limits the distribution of vegetation and controls the types of dominating plants. As a result, the Kalahari is mostly covered by species adapted to such environmental conditions [45,46].

A contrasting example of landscape is the Makgadikgadi Pan, which is a salt pan with specific dune, barchan, and crescentic sand landforms, located over the dry savanna in the north-eastern area of the country and formed as a tectonic surface expression in the Makgadikgadi Rift Zone [47,48]. Being one of the largest salt flats in the world [49], it impacts the neighbouring lands through the soil system. Thus, the eolian salts from the evaporite-covered Makgadikgadi Depression, such as chlorides, sodium, and bicarbonates, are being transported to the soil in the adjacent regions. These soluble salts are then further transported down to the ground soil water, which significantly degrades water quality [50].

### 2.2. Climate

Botswana is a country with dominating semi-arid climate, vast areas lacking surface water, and high extreme temperatures, as seen in Figures 2 and 3.

Due to the extreme temperatures during the summer period and droughts [51,52], Botswana is one of the world's least populated countries. The population of Botswana needs to adapt to special living setting and adjust their social activities as they are constrained by the environmental factors: livestock diseases, human-wildlife conflicts, drought, water supply deficit, resource scarcity, fragmented landscapes, and bush encroachment of savannah ecosystems [53–55]. Furthermore, specific environmental and climate conditions of the country and periods of severe drought contribute to the development of diseases, such as malaria or cholera. This is proved in several existing studies that analysed a correlation between the distribution of epidemiological diseases and environmental variations in Botswana: summer rainfalls, mean annual temperatures, and topographic altitude [56–59]. At the same time, extreme temperature varies regionally, according to the distribution of landforms, which is also influenced by water areas and rivers, see Figure 3.

The links between the topographic structure of the country and climate processes can be revealed through a regional analysis of the variation of climate parameters with regard to the topography. For example, a notable landform of Botswana is presented by the Okavango Delta, which is located in the north-western part of the country (Figure 1). It is one of the world's largest inland deltas, losing about 98% of its water through evapotranspiration [60]. The wetlands of the Okavango Delta are notable for the variety of ecosystems present, with a large biodiversity level in fauna and flora. The region of Okavango provides the livelihood of the local communities as well as the a tourism destination that substantially contributes to the revenue of Botswana [61,62].
<figure><img src="image_5.png" /><figcaption>**Figure 2.** T °C minimum (*T*<sub>min</sub>) in Botswana. Mapping: GMT. Source: author.</figcaption></figure>

In such regions, special measures and actions were introduced as state subsidies for drought relief management aimed at food and nutrition security [63]. The assessment of social vulnerability to droughts in Botswana acts as a measure for the reduction of climate and environmental disaster risks [64]. Prevoius study [65] recognised that hydrological droughts in Botswana may have several months of lag after the climatological droughts, which shows a complex and non-linear correlation between the hydrological, topographic and environmental factors associated with climate variability, which depend on precipitation pattern, Figure 4.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** T °C maximum (*T*<sub>max</sub>) in Botswana. Mapping: GMT. Source: author.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Precipitation in Botswana (2018). Mapping: GMT. Source: author.</figcaption></figure>

### 3. Data

This study used a TerraClimate, a free global climate dataset [66], with publicly available data: https://climate.northwestknowledge.net/TERRACLIMATE (accessed date: 21 August 2022). The TerraClimate includes several variables showing climate and environmental parameters, based on regular measurements of relevant data. These include the following key parameters: extreme air temperatures (minimum/maximum), precipitation, climatic water deficit, soil water equivalent, runoff, soil moisture, Palmer Drought Severity Index (PDSI), reference (potential) precipitation based on theoretical computations, actual evapotranspiration based on real measurements using remote sensing data, downward shortwave radiation at the surface (SRAD), vapour pressure (vap), vapour pressure deficit (vpd) and wind speed.

The TerraClimate dataset presents reference data for investigation of hydrological and climatic changes, which can be applied for agricultural modelling, such as crop water demand estimates, agricultural plant zoning, and natural resource management. To analyse the behaviour of these variables and variability over the territory of Botswana. The hydrological applications of TerraClimate include, for instance, surface runoff simulation and estimation of basin budget. Further applications of TerraClimate can be applied to the
modelling vegetation health and forest monitoring using correlations between the severity of droughts, hydrological and temperature regime and mapping soil moisture, Figure 5.

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Soil moisture in Botswana. Mapping: GMT. Source: author.</figcaption></figure>

The important part of terraClimate variables is presented by the Palmer Drought Severity Index (PDSI) which was used for visualising the strength of drought over the territory of Botswana. Besides the PDSI, there are other types of indices used for meteorological and climate research, for instance: Standardized Precipitation Index (SPI), Standardized Precipitation Evapotranspiration Index (SPEI) [67,68] or drought models [69,70]. However, the PDSI was chosen due to its availability in TerraClimate dataset. Additionally, it is a reliable and reputable index, known since its development [71]. The PDSI is based on the measurement of dryness using data on precipitation and temperature.

The topographic mapping was based on the GEBCO dataset covering the Earth's surface, https://www.gebco.net/ [72,73]. GEBCO has two significant advantages as a topographic data source: it is open source and publicly available and has an unprecedentedly high resolution (15 arc-second). Its high topographic precision is inherently influenced by the algorithms of data capture, processing and modelling [74–76]. The GEBCO data is well suited to reliably visualise the topography and geomorphic structures in the processed maps, due to its high resolution, accuracy, and precision [77–81].
## 4. Methodology

The methodology of recent environmental studies is mostly focused on utilising the traditional approaches of data modelling associated with the GUI-based GIS software [82–92]. In contrast, this study is inspired by the approach of GMT [8], which performs illustrative mapping through running the scripts via the command line from the console, e.g., Figure 6. This yields print-quality maps based on the raster and vector layers processed using a combination of various GMT modules, as required in scripting.

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Potential evapotranspiration (PET). Data: WorldClim. Mapping: GMT. Source: author.</figcaption></figure>

GMT has gained popularity in topographic, marine geologic and geophysical mapping where highly detailed and accurate data for modelling the phenomena of the Earth processes and structures are essential [93–99].

The 12 GMT scripts used for plotting maps presented in this manuscript are publicly available in the open GitHub repository of the author with a full access to the GMT codes: https://github.com/paulinelemenkova/Mapping_Botswana GMT_Scripts. In contrast to classical GIS where only the GUI-based data processing with pre-defined functions can be performed, in a GMT scenario much more geo-information can be adjusted using flexible functionality of the code syntax: layout, design, inserted maps, added cartographic elements, scale, etc. In particular, the screen position of the legend map in the final layout
10. `readx0y0wh < tmp'`: here, the temporary file with the data is created.

11. `gmtpsxy - R - J - O - K - T - X - $x0 - Y - $y0 >> $ps'`: here, the plotting is finished, this part of the script is ready, and `-O -K` signifies the overlay and code continuation.

Using the same approach of scripting presented above, all the other climate maps have been plotted using the TerraClimate open source dataset for the year 2018, Figures 7 and 8.

<figure><img src="image_10.png" /><figcaption>**Figure 7.** Actual evapotranspiration (AET) in Botswana. Mapping: GMT. Source: author.</figcaption></figure>

The most essential code snippets used for plotting climate maps are as follows:

* Extracting a subset of AET for Botswana: `gmtrgrdcutTerraClimate_aet2018.nc - R19.5/29.5/-27/-17.5 - Gbw_aet.nc`.

* Adding grids for all the maps: `- Bpxg1f0.5a1 - Bpyg1f0.5a1 - Bsxg1 - Bsyg1'`. Here, the 'f' stands for the frequency of minor ticks, 'g' for major, and 'a' for annotations.

* Converting PostScript file to the image JPG file with 720 dpi resolution using GhostScript (PET of Botswana): `gmtpsconvertBWpET.ps - A0.5c - E720 - Tj - Z'`.

* Generating colour palette and adjusting it for the extent of values in PDSI map of Botswana: `gmtdmakecpt - Cqual - mixed - 12.cpt - V - T - 5.1/12.1 > pauline.cpt`.

* Using the Geospatial Data Abstraction Library (GDAL) library for analysis of the statistics for $T_{min}$ raster grid of Botswana: `gdalinfo - statsbw_tim.nc`.
The general climate maps included following parameters plotted on the maps: $T_{min}/T_{max}$, soil moisture, precipitation, downward surface shortwave radiation, vapour pressure, vapour pressure deficit, and wind speed, and the derived reference indices (PDSI, PET, AET). All the maps were plotted using Mercator projection for compatibility reasons, in order to achieve relevant matching coordinate reference grid at the identical spatial extent and congruent resolution in each map for the compatibility and agreement of the map series. The maps were adjusted for the extent of the territory of Botswana and the surrounding regions (neighbouring countries) by subtracting the grids from the corresponding NetCDF files using `-R` module, e.g., `gmtgrdcutTerraClimate_tmin_2018.nc -R19.5/29.5/-27/-17.5 -Gbw_tim.nc'`. Here, `-R19.5/29.5/-27/-17.5` stands for the coordinates in WESN (West-East-South-North) convention with coordinates of the southern hemisphere as negative values (here, the example is provided for the downward surface shortwave radiation map in Figure 8.

<figure><img src="image_11.png" /><figcaption>**Figure 8.** Downward surface shortwave radiation. Mapping: GMT. Source: author.</figcaption></figure>

The implementation of the GMT algorithm for cartographic data processing is straightforward. In addition to the parameters of modules needed to plot the maps, the optimising setup of general GMT parameters was done. This included, for instance, the annotation offset, which is required to tightly plot the auxiliary texts, using the following expression: `MAP_ANNOT_OFFSET0.2c`. Based on the preceding observations of the layouts of maps,
the width and colour of map countours was defined using following expression and pa-
rameters: MAP TICK PEN PRIMARY thinner, dimgray, which overwrote the previous
parameters while disregarding the earlier settings, see Figure 9.

Effectively, we applied similar parameters for the cartographic grid using the ex-
pression MAP GRID PEN thinnest, dimgray, in order to keep the grid visualised yet
not disturb the main content of maps. Furthermore, we introduce the template of the
output geographical coordinate of the maps by formatting each layout using the fol-
lowing command: gmtsetFORMAT_GEO_MAP = dddF. In these parameters, the goal
is to adjust the same parameters of plotting for all the maps from the map series, to
make them consistently designed for comparative analysis. Using similar approach, the
script initialised using the expression ps = Topo_BW.ps where all the commands exe-
cuted by modules were inserted. In the end of script, it was followed the by command
gmpscconvert TopoB W.ps - A0.2c - E720 - Tj - Z for data conversion.

<figure><img src="image_12.png" /><figcaption>**Figure 9.** Vapor pressure in Botswana (2018). Mapping: GMT. Source: author.</figcaption></figure>

**5. Results**

The analysis of the distribution of the temperature extremes over Botswana (Figures 2 and 3) is important to understand the distribution of PDSI values since these two parameters are linked. So, as can be seen in Figure 2 (Tmin in Botswana in 2018), the minimal values over the country are recorded in the south-eatsern and southern regions of the country, in
the proximity of the Lobatse city and the Limpopo river valley, where values are lower and reach +16.5 °C and further to the east (region of South Africa). The maximal temperature (*T*<sub>max</sub>) in Botswana in 2018 is shown in Figure 3.

An important aspect of the *T*<sub>max</sub> distribution over Botswana depicted in Figure 3 is that the highest values of the *T*<sub>max</sub> are in the southern region of Kalahari (up to +38 °C, bright red colour in Figure 3). This is partially caused by the typical hydrological and meteorological process of the evaporation due to the lack of water as that all of the absorbed light goes into raising the surface temperature in the Kalahari Desert. The precipitation map (Figure 4) exhibited a certain degree of spatially dependent distribution of rainfalls, which generally increases to the north. However, the isolines demonstrate the rather wedge-shaped sharp decrease of that values that reach the maximum between the 21° and 22° E in the centre of the Kalahari Desert.

The Okavango Delta shows the transitional values between 80 and 111 mm (lilac and red colours in Figure 4). The Ngotwane and Marico basins have values of 38 to 48 mm and the Makgadikgadi Pan has dominating values of 80 to 95 mm. If we compare these values to the different regions of Botswana, it becomes clear that the precipitation is higher in the northern regions of the country, and in the southern along the river of Limpopo. The study of soil properties is important for environmental and landscape research [101–105], as soil systems directly affects the plant health through linked interactions. Thus, the variations in soil moisture indicate the environmental health of the ecosystems and dynamics in land cover types. Figure 5 illustrates the variations in soil moisture which show the least moisture in the Kalahari Desert (0–1 mm/m, i.e., extreme dry, yellow colour in Figure 5) with a slight increase in the NW and selected central districts of Kalahari (1–2 mm/m, beige), Figure 5.

Further gradual increase in soil moisture (2 to 5 mm/m, green colours) is notable in the north-western region of the country, specifically in the Okawanga Delta and on the border with Namibia. Moderate variations (up to 10 mm/m) are also notable around the Sehithwa village (Figure 5). The north-western bias in the trend of data distribution is notable in the north of the country owing to the precipitation, which increases northwards. At the same time, the temperature decrease is associated with climate effects in Botswana and direct impact from the Kalahari Desert in the southern region of the country, Figure 10.

The potential evapotranspiration (PET) (Figure 6) values are visualised continuously from the data extracted for Botswana. Here the data are ranging from a minimum of 133 mm to a maximum of 250 mm. Basic statistical inspection shows the mean of 196.048 mm and the standard deviation of 26.166 mm. The maximal is clearly visible in the lower region of Kalahari region (over 240 mm). Using the GDAL inspection by 'gdalinfo', the actual evapotranspiration (AET) Figure 7) was analysed and the data variations is shown as follows: a minimum of 13 mm, a maximum of 148 mm, a mean of 73.225 mm, and the standard deviation of 28.295 mm. Here, the minimal values correspond to the location of the Kalahari Desert in southern Botswana (values 0 to 50 mm). The gradients of isolines decrease (become less curved) in the northward direction towards the Okavango Delta, where the values reach up to 120 mm.

The curves of the continuous fields showing the downward surface shortwave radia-
tion are illustrated in Figure 8, which reflects the total amount of shortwave radiation, both
direct and diffuse, that reaches the Earth's surface over the area in Botswana. The overall
increase in the values can be noticed as increasing towards the southern region of Kalahari
where the values reach up to 326 Wm⁻² (bright red colours in Figure 8) from the lowest in
the northern region with the minimal values of 234 Wm⁻² (light beige colours in Figure 8).
The vapour pressure (Figure 9) and vapour pressure deficit (Figure 10) correlate in general
trend increasing towards the southern regions of Kalahari Desert and reaching values of
2.56 atm from the minimal of 1.17 atm in the north of the country for vapour pressure.
<figure><img src="image_13.png" /><figcaption>**Figure 10.** Vapor pressure deficit in Botswana. Mapping: GMT. Source: author.</figcaption></figure>

As shown in Figure 10 (Vapour Pressure Deficit, VPD of Botswana), there is a general trend for vapour pressure to be higher in the SW (Kalahari Desert with the highest values of 3.3, bright yellow colour) and lower in the NE of Botswana (up to 1.4, purple colours in Figure 10). The VPD shows the deficit between the practical (real) value of the moisture in the air and the theoretical amount of moisture that the air can contain when saturated. The measured wind speed (Figure 11) reaches the maximal values in the area of Jwaneng eastward of Kalahari. The inspection of data gives the following data range: minimum at 1.8 m/s, maximum at 3.4 m/s, mean at 2.587 m/s, and standard deviation of values at 0.352 m/s (Figure 11).
<figure><img src="image_14.png" /><figcaption>**Figure 11.** Wind speed in Botswana. Mapping: GMT. Source: author.</figcaption></figure>

In general, the PDSI drought levels in Botswana in 2018 (Figure 12) show stepwise variation of index with roughly seven areas: (1) −3.7 to −2.2. (extreme drought); (2) −2.2 to −0.8 (strong drought, southern Kalahari); (3) −0.8 to 0.7 (significant drought): central Kalahari; (4) 0.7 to 2.1 (moderate drought); (5) 2.1 to 3.5 (lesser drought); (6) 3.5 to 4.9 (low drought); (7) 4.9 to 6.4 (least drought), Figure 12. The zero level is highlighted by the thick white isoline for the convenience of the classification (Figure 12). One possibility to explain the PDSI values (Figure 12) showing the minimal values, i.e., extreme drought with values −2.2 to −3.8), is that the minimal temperatures are recorded as the deepest area around the Lobaste, which is situated in the southern region of the country (Figure 2). Furthermore, this also corresponds to the low precipitation level in this region, compare with Figure 4, 35 to 61 mm, green colours: region of Jwaneng, Gaborone, area southwards of Kalahari.

In this region, the Downward Shortwave Surface Radiation (DSSR) is the highest, with values exceeding 318 Wm⁻² (compare Figure 12 with Figure 8). The next (lesser) level of drought severity with values of −0.8 to −2.2 (light plum colour in Figure 12) was examined with map of the vapour pressure having consecutively decreasing values southwards of Kalahari (1.6 to 1.2) correlating to the second level of drought of Botswana (−0.8 to −2.2).
<figure><img src="image_15.png" /><figcaption>**Figure 12.** PDSI in Botswana. Mapping: GMT. Source: author.</figcaption></figure>

**6. Discussion**

The proposed framework of the GMT-based mapping produces a series of accurate maps through the efficient cartographic workflow. The speed of running the scripts is only several seconds for even complex data handling, which ensures the effectiveness of the cartographic workflow and reduces the human-based routine. Another important difference between the GIS and GMT is that each script is technically defined in the lines of code. This guarantees plotted map to be resulted as designed, in contrast to the hidden flaws and occasional errors of GIS interface. In other words, both the data processing and data visualisation in GMT are defined from a fixed commands of code lines that handle the data in a straightforward way and results in the programmed output. Using the presented map of PDSI, the results can be further evaluated for the analysis of four types of drought classification according to the regional setting of Botswana:

- Meteorological drought (dominating dry weather conditions);

- Hydrological drought (noticeable low water supply);

- Agricultural drought (crops are being affected);

- Socioeconomic drought.
Here, the fourth type affects the social system and well-being of the population the most. These results continue the existing investigations on the complex interrelationships between climate, environmental, and topographic factors in Botswana [106-108].

The maps can be used as a background for a more complex analysis of climate and environmental effects on social well-being [109–112]. This includes the Kalahari Desert, Makgadikgadi Pan, and Okavango Delta as the most important regions and landforms of Botswana, for determination of the possibility of droughts appearance in Botswana, using several multi-temporal climate datasets for different years and decades.

## 7. Conclusions

With ever-expanding geospatial data collection, the analysis of the geo-information is becoming a fundamental research issue for many applications in Earth sciences. These include a wide variety of cases starting from data visualisation for detecting environmental trends, or data processing for climate change analysis, to more complex tasks of climate modelling, including prognosis and forecasting. Mapping spatial data aims to visualise environmental parameters in the most effective way in order to highlight variations in parameters of topography, temperature, soil moisture, land use dynamic, and many more [113–117]. To this end, cartographic visualisation has to bridge the gap between technical issues of spatial data processing and semantic interpretation of geo-information as a source of knowledge in environmental and Earth science.

The conceptual approach of data handling in GIS includes data semantics, such as points, lines, and polygons in vector types, and pixels in raster types, which can be read and processed by a variety of functions and modules existing in various GIS [118–120]. With enough functionality and robust data, such methods can be used for data visualisation and modelling. However, obtaining higher automation and speed of data processing becomes more difficult in traditional GIS. This is caused due to subjective misclassification and errors originating from manual data handling. Processing large spatial datasets increases the need for automation to avoid repetitive cartographic operations and ensure quick and seamless data formatting. Moreover, while GIS methods are more or less capable of handling spatial data as objects, including surface modelling, they are often more sensitive to visualisation and mapping. In many cases, preparing final graphical map layouts in GIS remains a handmade process that requires optimisation and automation.

Currently, most geospatial data processing tasks and the analysis of geo-information require that spatial features and objects plotted on the maps are extracted from the raw datasets effectively and rapidly, to produce the reliable and readable maps. The GMT method works in a unified framework which can jointly optimise technical data processing and produce the aesthetically adjusted maps. Compared to the traditional GIS, applying the scripting approach of GMT benefits from being a simple and straightforward approach. It is based on the programming approach of GMT efficient optimisation of geodata processing where each cartographic operation is processed by a single module with various finely adjusted options through the extended functionality. Once written, scripts are easily amenable and adjustable to repetitive data processing that require similar tasks. However, it also requires certain skills in programming for writing and executing the scripts.

This study presented the new series of climate and environmental inventory maps of Botswana for 2018 using GMT scripting based on the TerraClimate and GEBCO datasets. The main contribution of this paper is to demonstrate that a straightforward programming method of modern cartography can be used to represent climate and environmental parameters accurately and effectively without requiring additional commercial GIS. The presented maps can be used as a background for further regional climate studies of Botswana. Usually based on processing a variety of thematic datasets, environmental and applied geological studies depict spatio-temporal trends in ecosystems related to the Earth's surface and processes using topographic, climate, and environmental datasets [121–127]. Practical application of the presented maps includes environmental assessment and agricultural analysis in Botswana, as it is strongly affected by climate variability and characteristics,
such as periods of drought, extreme temperatures and vapour pressure, as well as seasonal changes in rainfalls [128–131].

Advances in the data processing can be recommended for future studies, including the assessment of correlation between the selected environmental and climate variables and their location and spatio-temporal variations. Among others, climate parameters can be used as datasets: PDSI, AET, PET, $T_{max}/T_{min}$, soil moisture, precipitation intensity and repeatability with respect to local topography. Furthermore, as a direction for future research, one may also consider statistical data analysis, which is possible in GMT, besides mapping. For instance, plotting histograms of data distribution for climatic and topographic indicators can be performed using module 'pshistogram'.

**Funding:** This research received no external funding.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The GitHub repository of the author with available GMT scripts: https://github.com/paulinelemenkova/Mapping_Botswana_GMT_Scripts (accessed on 25 June 2022).

**Acknowledgments:** The author cordially thanks the three anonymous reviewers for their comments, corrections and suggestions, which improved the initial version of the manuscript.

**Conflicts of Interest:** The author declares no conflict of interest.

## Abbreviations

The following abbreviations are used in this manuscript:

<table><tr><td>AET</td><td>Actual Evapotranspiration</td></tr><tr><td>DSSR</td><td>Downward Shortwave Surface Radiation</td></tr><tr><td>GDAL</td><td>Geospatial Data Abstraction Library</td></tr><tr><td>GEBCO</td><td>General Bathymetric Chart of the Oceans</td></tr><tr><td>GIS</td><td>Geographic Information System</td></tr><tr><td>GMT</td><td>Generic Mapping Tools</td></tr><tr><td>GUI</td><td>Graphical User Interface</td></tr><tr><td>NetCDF</td><td>Network Common Data Form</td></tr><tr><td>PDSI</td><td>Palmer Drought Severity Index</td></tr><tr><td>PET</td><td>Potential Evapotranspiration</td></tr><tr><td>SRTM</td><td>Shuttle Radar Topography Mission</td></tr><tr><td>SPI</td><td>Standardized Precipitation Index</td></tr><tr><td>SPEI</td><td>Standardized Precipitation Evapotranspiration Index</td></tr><tr><td>VPD</td><td>Vapor Pressure Deficit</td></tr><tr><td>VP</td><td>Vapor Pressure</td></tr></table>

## References

1. Turner, M.D. Methodological Reflections on the Use of Remote Sensing and Geographic Information Science in Human Ecological Research. *Hum. Ecol.* **2003**, *31*, 255–279. [CrossRef]

2. Hamandawana, H. The impacts of herbivory on vegetation in Moremi Game Reserve, Botswana: 1967–2001. *Reg. Environ. Chang.* **2012**, *12*, 1–15. [CrossRef]

3. Maoyi, M.L.; Abiodun, B.J. Investigating the response of the Botswana High to El Niño Southern Oscillation using a variable resolution global climate model. *Theor. Appl. Climatol.* **2022**, *147*, 1601–1615. [CrossRef]

4. Segobye, A. Divided commons: The political economy of Southern Africa's cultural heritage landscapes—Observations of the central Kalahari game reserve, Botswana. *Archaeologies* **2006**, *2*, 52–72. [CrossRef]

5. Ceccato, P.; Fernandes, K.; Ruiz, D.; Allis, E. Climate and environmental monitoring for decision making. *Earth Perspect.* **2014**, *1*, 16. [CrossRef]

6. de Sherbinin, A. Climate change hotspots mapping: What have we learned? *Clim. Chang.* **2014**, *123*, 23–37. [CrossRef]

7. Herslund, L.B.; Jalayer, F.; Jean-Baptiste, N.; Jørgensen, G.; Kabisch, S.; Kombe, W.; Lindley, S.; Nyed, P.K.; Pauleit, S.; Printz, A.; et al. A multi-dimensional assessment of urban vulnerability to climate change in Sub-Saharan Africa. *Nat. Hazards* **2016**, *82*, 149–172. [CrossRef]
8. Wessel, P.; Luis, J.F.; Uieda, L.; Scharroo, R.; Wobbe, F.; Smith, W.H.F.; Tian, D. The Generic Mapping Tools version 6. *Geochem. Geophys. Geosyst.* **2019**, *20*, 5556–5564. [CrossRef]

9. Batisani, N. The Spatio-Temporal-Severity Dynamics of Drought in Botswana. *J. Environ. Prot.* **2011**, *2*, 803–816. [CrossRef]

10. Vanderpost, C.; Ringrose, S.; Matheson, W.; Arntzen, J. Satellite based long-term assessment of rangeland condition in semi-arid areas: An example from Botswana. *J. Arid. Environ.* **2011**, *75*, 383–389. [CrossRef]

11. Mogotsi, K.; Nyangito, M.M.; Nyariki, D.M. The role of drought among agro-pastoral communities in a semi-arid environment: The case of Botswana. *J. Arid. Environ.* **2013**, *91*, 38–44. [CrossRef]

12. Akinyemi, F.O.; Mashame, G. Analysis of land change in the dryland agricultural landscapes of eastern Botswana. *Land Use Policy* **2018**, *76*, 798–811. [CrossRef]

13. Atekwana, E.A.; Molwalefhe, L.; Kgaodi, O.; Cruse, A.M. Effect of evapotranspiration on dissolved inorganic carbon and stable carbon isotopic evolution in rivers in semi-arid climates: The Okavango Delta in North West Botswana. *J. Hydrol. Reg. Stud.* **2016**, *7*, 1–13. [CrossRef]

14. Ringrose, S.; Chipanshi, A.; Matheson, W.; Chanda, R.; Motoma, L.; Magole, I.; Jellema, A. Climate- and Human-Induced Woody Vegetation Changes in Botswana and Their Implications for Human Adaptation. *Environ. Manag.* **2002**, *30*, 98–109. [CrossRef]

15. Kenabatho, P.K.; Parida, B.P.; Moalafhi, D.B. The value of large-scale climate variables in climate change assessment: The case of Botswana's rainfall. *Phys. Chem. Earth Parts A/B/C* **2012**, *50–52*, 64–71. [CrossRef]

16. Byakatonda, J.; Parida, B.P.; Kenabatho, P.K.; Moalafhi, D.B. Prediction of onset and cessation of austral summer rainfall and dry spell frequency analysis in semiarid Botswana. *Theor. Appl. Climatol.* **2019**, *135*, 101–117. [CrossRef]

17. Batisani, N.; Yarnal, B. Rainfall variability and trends in semi-arid Botswana: Implications for climate change adaptation policy. *Appl. Geogr.* **2010**, *30*, 483–489. [CrossRef]

18. Longley, P.A.; Goodchild, M.F.; Maguire, D.J.; Rhind, D.W. *Geographic Information Science and Systems*, 4th ed.; John Wiley & Sons: Hoboken, NJ, USA, 2015.

19. Chrisman, N. *Exploring Geographic Information Systems*, 2nd ed.; John Wiley & Sons: Hoboken, NJ, USA, 2001.

20. Vankova, L.; Krejza, Z.; Kocourkova, G.; Laciga, J. Geographic Information System Usage Options in Facility Management. *Procedia Comput. Sci.* **2022**, *196*, 708–716. doi: 10.1016/j.procs.2021.12.067. [CrossRef]

21. Rahman, M.M.; Szabó, G. Sustainable Urban Land-Use Optimization Using GIS-Based Multicriteria Decision-Making (GIS-MCDM) Approach. *ISPRS Int. J. Geo-Inf.* **2022**, *11*, 313. [CrossRef]

22. Visvalingam, M. Cartography, GIS and Maps in Perspective. *Cartogr. J.* **1989**, *26*, 26–32. [CrossRef]

23. Nilsson, A.; Mentis, D.; Korkovelos, A.; Otwani, J. A GIS-Based Approach to Estimate Electricity Requirements for Small-Scale Groundwater Irrigation. *ISPRS Int. J. Geo-Inf.* **2021**, *10*, 780. [CrossRef]

24. Lemenkova, P. The visualization of geophysical and geomorphologic data from the area of Weddell Sea by the Generic Mapping Tools. *Stud. Quat.* **2021**, *38*, 19–32. [CrossRef]

25. Xing, Z.; Guo, W. A New Urban Space Analysis Method Based on Space Syntax and Geographic Information System Using Multisource Data. *ISPRS Int. J. Geo-Inf.* **2022**, *11*, 297. [CrossRef]

26. Reed, M.; Dougill, A. Linking degradation assessment to sustainable land management: A decision support system for Kalahari pastoralists. *J. Arid. Environ.* **2010**, *74*, 149–155. [CrossRef]

27. Jovanovic, N.Z.; Annandale, J.G.; Claassens, A.S.; Lorentz, S.A.; Tanner, P.D. Modeling Irrigation with Gypsiferous Mine Water: A Case Study in Botswana. *Mine Water Environ.* **2001**, *20*, 65–72. [CrossRef]

28. Basupi, L.V.; Quinn, C.H.; Dougill, A.J. Historical perspectives on pastoralism and land tenure transformation in Ngamiland, Botswana: What are the policy and institutional lessons? *Pastoralism* **2017**, *7*, 24. [CrossRef]

29. Mfundisi, K.B. Overview of an integrated management plan for the Okavango Delta Ramsar Site, Botswana. *Wetlands* **2008**, *28*, 538–543. [CrossRef]

30. Osupile, K.; Yahya, A.; Samikannu, R. A Review on Agriculture Monitoring Systems using Internet of Things (IoT). In Proceedings of the 2022 International Conference on Applied Artificial Intelligence and Computing (ICAAIC), Kuala Lumpur, Malaysia, 9–10 December 2022; pp. 1565–1572. [CrossRef]

31. Tsheboeng, G.; Murray-Hudson, M.; Kashe, K. Regeneration status of riparian tree species in two sites that differ in land-use in the Okavango Delta, Botswana. *J. For. Res.* **2017**, *28*, 1073–1082. [CrossRef]

32. Moses, O. Weather systems influencing Botswana rainfall: The case of 9 December 2018 storm in Mahalapye, Botswana. *Model. Earth Syst. Environ.* **2019**, *5*, 1473–1480. [CrossRef]

33. Mogome-Ntsatsi, K.; Adeola, O.A. Promoting environmental awareness in Botswana: The role of community education. *Environmentalist* **1995**, *15*, 281–292. [CrossRef]

34. Lange, G. Wealth, Natural Capital, and Sustainable Development: Contrasting Examples from Botswana and Namibia. *Environ. Resour. Econ.* **2004**, *29*, 257–283. [CrossRef]

35. Kasalica, V.; Lamprecht, A.L. Automated composition of scientific workflows: A case study on geographic data manipulation. In Proceedings of the 2018 IEEE 14th International Conference on e-Science (e-Science), Amsterdam, The Netherlands, 29 October–1 November 2018; pp. 362–363. [CrossRef]

36. Lemenkova, P. Automatic Data Processing for Visualising Yap and Palau Trenches by Generic Mapping Tools. *Cartogr. Lett.* **2019**, *27*, 72–89. [CrossRef]
37. Lemenkova, P. Seismicity in the Afar Depression and Great Rift Valley, Ethiopia. *Environ. Res. Eng. Manag.* **2022**, *78*, 83–96. [CrossRef]

38. Spinellis, D. Drawing Tools. *IEEE Softw.* **2009**, *26*, 12–13. [CrossRef]

39. Lemenkova, P. GMT Based Comparative Analysis and Geomorphological Mapping of the Kermadec and Tonga Trenches, Southwest Pacific Ocean. *Geogr. Tech.* **2019**, *14*, 39–48. [CrossRef]

40. Yajima, T.; Gwandu, K.R.; Sigweni, M.A. Mapping mafic and ultramafic units using remote sensing data in north east Botswana. In *Proceedings of the 2011 IEEE International Geoscience and Remote Sensing Symposium, Vancouver, BC, Canada, 24–29 July 2011*; pp. 2218–2220. [CrossRef]

41. Ceccato, P.; Vancutsem, C.; Temimi, M. Monitoring air and Land Surface Temperatures from remotely sensed data for climate-human health applications. In *Proceedings of the 2010 IEEE International Geoscience and Remote Sensing Symposium, Honolulu, HI, USA, 25–30 July 2010*; pp. 178–180. [CrossRef]

42. Campbell, A. The use of wild food plants, and drought in Botswana. *J. Arid. Environ.* **1986**, *11*, 81–91. [CrossRef]

43. Ellery, W.N.; McCarthy, T.S.; Smith, N.D. Vegetation, hydrology, and sedimentation patterns on the major distributary system of the Okavango fan, Botswana. *Wetlands* **2003**, *23*, 357–375. [CrossRef]

44. Franchi, F.; Kelepile, T.; Di Capua, A.; De Wit, M.C.J.; Kemiso, O.; Lasarwe, R.; Catuneanu, O. Lithostratigraphy, sedimentary petrography and geochemistry of the Upper Karoo Supergroup in the Central Kalahari Karoo Sub-Basin, Botswana. *J. Afr. Earth Sci.* **2021**, *173*, 104025. [CrossRef]

45. Masunga, G.S.; Moe, S.R.; Pelekekae, B. Fire and Grazing Change Herbaceous Species Composition and Reduce Beta Diversity in the Kalahari Sand System. *Ecosystems* **2013**, *16*, 252–268. [CrossRef]

46. Elliott, D.R.; Thomas, A.D.; Hoon, S.R.; Sen, R. Niche partitioning of bacterial communities in biological crusts and soils under grasses, shrubs and trees in the Kalahari. *Biodivers. Conserv.* **2014**, *23*, 1709–1733. [CrossRef]

47. Eckardt, F.D.; Cotterill, F.P.; Flügel, T.J.; Kahle, B.; McFarlane, M.; Rowe, C. Mapping the surface geomorphology of the Makgadikgadi Rift Zone (MRZ). *Quat. Int.* **2016**, *404*, 115–120. [CrossRef]

48. McFarlane, M.; Long, C. Pan floor 'barchan' mounds, Ntwetwe Pan, Makgadikgadi, Botswana: Their origin and palaeoclimatic implications. *Quat. Int.* **2015**, *372*, 108–119. [CrossRef]

49. Burrough, S.L.; Thomas, D.S.G.; Bailey, R.M.; Davies, L. From landform to process: Morphology and formation of lake-bed barchan dunes, Makgadikgadi, Botswana. *Geomorphology* **2012**, *161–162*, 1–14. [CrossRef]

50. Wood, W.W.; Eckardt, F.D.; Kraemer, T.F.; Eng, K., Quantitative Eolian Transport of Evaporite Salts from the Makgadikgadi Depression (Ntwetwe and Sua Pans) in Northeastern Botswana: Implications for Regional Ground-Water Quality. In *Sabhka Ecosystems. Tasks for Vegetation Science*; Öztürk, M., Böer, B., Barth, H.J., Clüsener-Godt, M., Khan, M., Breckle, S.W., Eds.; Springer: Dordrecht, The Netherlands, 2010; Voume 46; pp. 27–37. [CrossRef]

51. Hemming, C.F. Drought in Botswana. *J. Arid. Environ.* **1980**, *3*, 181. [CrossRef]

52. Karikari, S.K.; Tabona, T.T. Constitutive traits and selective indices of Bambara groundnut (*Vigna subterranea* (L.) Verdc) landraces for drought tolerance under Botswana conditions. *Phys. Chem. Earth Parts A/B/C* **2004**, *29*, 1029–1034. [CrossRef]

53. Badimo, D.; Mogomotsi, P.K.; Mogomotsi, G.E.J. Participatory mapping of water collection by Sexaxa community in Okavango, Botswana. *Phys. Chem. Earth* **2020**, *124*, 102963. [CrossRef]

54. Basupi, L.V.; Quinn, C.H.; Dougill, A.J. Adaptation strategies to environmental and policy change in semi-arid pastoral landscapes: Evidence from Ngamiland, Botswana. *J. Arid. Environ.* **2019**, *166*, 17–27. [CrossRef]

55. Charis, G.; Danha, G.; Muzenda, E. Waste valorisation opportunities for bush encroacher biomass in savannah ecosystems: A comparative case analysis of Botswana and Namibia. *Procedia Manuf.* **2019**, *35*, 974–979. [CrossRef]

56. Craig, M.H.; Sharp, B.L.; Mabaso, M.L.H.; Kleinschmidt, I. Developing a spatial-statistical model and map of historical malaria prevalence in Botswana using a staged variable selection procedure. *Int. J. Health Geogr.* **2007**, *6*. [CrossRef]

57. Tonnang, H.E.; Tchouassi, D.P.; Juarez, H.S.; Igweta, L.K.; Djouaka, R.F. Zoom in at African country level: Potential climate induced changes in areas of suitability for survival of malaria vectors. *Int. J. Health Geogr.* **2014**, *13*. [CrossRef]

58. Motsholapheko, M.R.; Ngwenya, B.N. Access to Water Resources and Household Vulnerability to Malaria in the Okavango Delta, Botswana. In *African Handbook of Climate Change Adaptation*; Oguge, N., Ayal, D., Adeleke, L., da Silva, I., Eds.; Springer: Cham, Switzerland, 2021; pp. 1227–1246. [CrossRef]

59. Charnley, G.E.C.; Kelman, I.; Green, N.; Hinsley, W.; Gaythorpe, K.A.M.; Murray, K.A. Exploring relationships between drought and epidemic cholera in Africa using generalised linear models. *BMC Infect. Dis.* **2021**, *21*, 1177. [CrossRef]

60. Moses, O.; Hambira, W.L. Effects of climate change on evapotranspiration over the Okavango Delta water resources. *Phys. Chem. Earth Parts A/B/C* **2018**, *105*, 98–103. [CrossRef]

61. Milzow, C.; Kgotlhang, L.; Bauer-Gottwein, P.; Meier, P.; Kinzelbach, W. Regional review: The hydrology of the Okavango Delta, Botswana—Processes, data and modelling. *Hydrogeol. J.* **2009**, *17*, 1297–1328. [CrossRef]

62. Wolski, P.; Murray-Hudson, M. Flooding dynamics in a large low-gradient alluvial fan, the Okavango Delta, Botswana, from analysis and interpretation of a 30-year hydrometric record. *Hydrol. Earth Syst. Sci.* **2005**, *2*, 1865–1892. [CrossRef]

63. Belbase, K.; Morgan, R. Food security and nutrition monitoring for drought relief management: The case of Botswana. *Food Policy* **1994**, *19*, 285–300. [CrossRef]

64. Dintwa, K.F.; Letamo, G.; Navaneetham, K. Quantifying social vulnerability to natural hazards in Botswana: An application of cutter model. *Int. J. Disaster Risk Reduct.* **2019**, *37*, 101189. [CrossRef]
65. Byakatonda, J.; Parida, B.P.; Kenabatho, P.K. Relating the dynamics of climatological and hydrological droughts in semiarid Botswana. *Phys. Chem. Earth. Parts A/B/C* **2018**, *105*, 12–24. [CrossRef]

66. Abatzoglou, J.; Dobrowski, S.; Parks, S.; Hegewisch, K.C. TerraClimate, a high-resolution global dataset of monthly climate and climatic water balance from 1958–2015. *Sci. Data* **2018**, *5*, 170191. [CrossRef]

67. Byakatonda, J.; Parida, B.P.; Kenabatho, P.K.; Moalafhi, D.B. Influence of climate variability and length of rainy season on crop yields in semiarid Botswana. *Agric. For. Meteorol.* **2018**, *248*, 130–144. [CrossRef]

68. Byakatonda, J.; Parida, B.P.; Moalafhi, D.B.; Kenabatho, P.K. Analysis of long term drought severity characteristics and trends across semiarid Botswana using two drought indices. *Atmos. Res.* **2018**, *213*, 492–508. [CrossRef]

69. Ujeneza, E.L.; Abiodun, B.J. Drought regimes in Southern Africa and how well GCMs simulate them. *Clim. Dyn.* **2015**, *44*, 1595–1609. [CrossRef]

70. Akinyemi, F.O.; Abiodun, B.J. Potential impacts of global warming levels 1.5 °C and above on climate extremes in Botswana. *Clim. Chang.* **2019**, *154*, 387–400. [CrossRef]

71. Palmer, W. *Meteorological Drought*; Research Paper no. 45; U.S. Department of Commerce Weather Bureau: Washington, DC, USA, 1965.

72. Ermel, H. Der deutsche Beitrag zur Neuherstellung der General Bathymetric Chart of the Oceans (GEBCO). *Dtsch. Hydrogr. Z.* **1966**, *19*, 49–57. [CrossRef]

73. Schenke, H. General Bathymetric Chart of the Oceans (GEBCO). In *Encyclopedia of Marine Geosciences*; Springer: Dordrecht, The Netherlands, 2016; pp. 268–269. [CrossRef]

74. Bracke, S.; Sochacki, S.; Radetzky, M. Multivariate Data Analytics in Surface Topography Assessments: Case Study High Precision Fine Grinding Processes. *Procedia Manuf.* **2019**, *39*, 1752–1761. [CrossRef]

75. Smith, W.H.F. On the accuracy of digital bathymetric data. *J. Geophys. Res.* **1993**, *98*, 9591–9603. [CrossRef]

76. Lemenkova, P. Mapping Ghana by GMT and R scripting: Advanced cartographic approaches to visualize correlations between the topography, climate and environmental setting. *Adv. Geod. Geoinf.* **2022**, *71*, e16. [CrossRef]

77. Jakobsson, M.; Macnab, R. A Comparison Between GEBCO Sheet 5.17 and the International Bathymetric Chart of the Arctic Ocean (IBCAO) Version 1.0. *Mar. Geophys. Res.* **2006**, *27*, 35–48. [CrossRef]

78. Lemenkova, P. Topography of the Aleutian Trench south-east off Bowers Ridge, Bering Sea, in the context of the geological development of North Pacific Ocean. *Baltica* **2021**, *34*, 27–46. [CrossRef]

79. Hall, J.K. GEBCO Centennial Special Issue—Charting the Secret World of the Ocean Floor: The GEBCO Project 1903–2003. *Mar. Geophys. Res.* **2006**, *27*, 1–5. [CrossRef]

80. Lemenkova, P. GEBCO Gridded Bathymetric Datasets for Mapping Japan Trench Geomorphology by Means of GMT Scripting Toolset. *Geod. Cartogr.* **2020**, *46*, 98–112. [CrossRef]

81. Marks, K.M.; Smith, W.H.F. An Evaluation of Publicly Available Global Bathymetry Grids. *Mar. Geophys. Res.* **2006**, *27*, 19–34. [CrossRef]

82. Kolawole, O.D.; Wolski, P.; Ngwenya, B.; Mmopelwa, G. Ethno-meteorology and scientific weather forecasting: Small farmers and scientists' perspectives on climate variability in the Okavango Delta, Botswana. *Clim. Risk Manag.* **2014**, *4–5*, 43–58. [CrossRef]

83. Lemenkova, P. Tanzania Craton, Serengeti Plain and Eastern Rift Valley: Mapping of geospatial data by scripting techniques. *Est. J. Earth Sci.* **2022**, *71*, 61–79. [CrossRef]

84. Ringrose, S.; Chanda, R.; Nkambwe, M.; Sefe, F. Environmental change in the mid-Boteti area of north-central Botswana: Biophysical processes and human perceptions. *Environ. Manag.* **1996**, *20*, 397–410. [CrossRef]

85. Lindh, P.; Lemenkova, P. Evaluation of Different Binder Combinations of Cement, Slag and CKD for S/S Treatment of TBT Contaminated Sediments. *Acta Mech. Autom.* **2021**, *15*, 236–248. [CrossRef]

86. Manisa, M.; Das, R.K.; Segobye, M.; Maphale, L. Developing local geoid model to assess accuracy of orthometric heights from GPS-based ellipsoidal heights in Botswana. *Spat. Inf. Res.* **2016**, *24*, 607–616. [CrossRef]

87. Lemenkova, P. Robust Vegetation Detection Using RGB Colour Composites and Isoclust Classification of the Landsat TM Image. *Geomat. Landmanag. Landsc.* **2021**, *4*, 147–167. [CrossRef]

88. Ringrose, S.; Musisi-Nkambwe, S.; Coleman, T.; Nellis, D.; Bussing, C. Use of Landsat Thematic Mapper Data to Assess Seasonal Rangeland Changes in the Southeast Kalahari, Botswana. *Environ. Manag.* **1999**, *23*, 125–138. [CrossRef]

89. Lemenkova, P. Dataset compilation by GRASS GIS for thematic mapping of Antarctica: Topographic surface, ice thickness, subglacial bed elevation and sediment thickness. *Czech Polar Rep.* **2021**, *11*, 67–85. [CrossRef]

90. Pachka, H.; Annelise, T.; Alan, K.; Tshikae, P.; Kgori, P.; Chevalier, V.; Paweska, J.; Jori, F. Rift Valley fever vector diversity and impact of meteorological and environmental factors on Culex pipiens dynamics in the Okavango Delta, Botswana. *Parasites Vectors* **2016**, *9*, 434. [CrossRef]

91. Pricope, N.G. Variable-source flood pulsing in a semi-arid transboundary watershed: The Chobe River, Botswana and Namibia. *Environ. Monit. Assess.* **2013**, *185*, 1883–1906. [CrossRef] [PubMed]

92. Lemenkova, P. Sentinel-2 for High Resolution Mapping of Slope-Based Vegetation Indices Using Machine Learning by SAGA GIS. *Transylv. Rev. Syst. Ecol. Res.* **2020**, *22*, 17–34. [CrossRef]

93. Virden, W.; Habermann, T.; Glover, G.; Divins, D.; Sharman, G.; Fox, C. Multibeam bathymetric data at NOAA/NGDC. In Proceedings of the Oceans '04 MTS/IEEE Techno-Ocean '04 (IEEE Cat. No.04CH37600), Kobe, Japan, 9–12 November 2004; Volume 2, pp. 1159–1162. [CrossRef]
94. Lemenkova, P. Submarine tectonic geomorphology of the Pliny and Hellenic Trenches reflecting geologic evolution of the southern Greece. *Rud.-Geološko-Naft. Zb.* **2021**, *36*, 33–48. [CrossRef]

95. Senturk, S.; Cakir, Z.; Berk Ustundag, B. The potential of Sentinel-IA interferometric SAR data in monitoring of surface subsidence caused by overdrafting groundwater in agricultural areas. In Proceedings of the 2016 Fifth International Conference on Agro-Geoinformatics (Agro-Geoinformatics), Tianjin, China, 18–20 July 2016; pp. 1–4. [CrossRef]

96. Lemenkova, P. Using GMT for 2D and 3D Modeling of the Ryukyu Trench Topography, Pacific Ocean. *Misc. Geogr.* **2021**, *25*, 213–225. [CrossRef]

97. Lindh, P.; Lemenkova, P. Geochemical tests to study the effects of cement ratio on potassium and TBT leaching and the pH of the marine sediments from the Kattegat Strait, Port of Gothenburg, Sweden. *Baltica* **2022**, *35*, 47–59. [CrossRef]

98. Lemenkova, P. Geomorphology of the Puerto Rico Trench and Cayman Trough in the Context of the Geological Evolution of the Caribbean Sea. *Ann. Univ. Mariae-Curie-Skłodowska, Sect. B*—Geogr. Geol. Mineral. Petrogr. **2020**, *75*, 115–141. [CrossRef]

99. Shi, H.; Du, Z.; Lu, Y.; Hu, X.; Ke, X. Amery Ice Shelf Digital Elevation Model from GLAS and GMT. In Proceedings of the 2009 Third International Symposium on Intelligent Information Technology Application, Nanchang, China, 21–22 November 2009; Volume 2; pp. 129–133. [CrossRef]

100. Lemenkova, P. Mapping submarine geomorphology of the Philippine and Mariana trenches by an automated approach using GMT scripts. *Proc. Latvian Acad. Sci. Sect. B* **2022**, *76*, 258–266. [CrossRef]

101. Totolo, O.; Mosweu, S. Spatial Variability of Selected Soil Properties in Relation to Different Land Uses in Northern Kgalagadi (Matsheng), Botswana. *Int. J. Geosci.* **2012**, *3*, 659–663. [CrossRef]

102. Lindh, P.; Lemenkova, P. Seismic velocity of P-waves to evaluate strength of stabilized soil for Svenska Cellulosa Aktiebolaget Biorefinery Östrand AB, Timrå. *Bull. Pol. Acad. Sci. Tech. Sci.* **2022**, *70*, 1–9. [CrossRef]

103. Mosweu, S.; Munyati, C.; Kabanda, T. Modification of soil properties by Prosopis L. in the Kalahari Desert, South-Western Botswana. *Open J. Ecol.* **2013**, *3*, 145–150. [CrossRef]

104. Lemenkov, V.; Lemenkova, P. Using TeX Markup Language for 3D and 2D Geological Plotting. *Found. Comput. Decis. Sci.* **2021**, *46*, 43–69. [CrossRef]

105. Zeilhofer, P. Soil Mapping in the Pantanal of Mato Grosso, Brazil, using Multitemporal Landsat TM data. *Wetl. Ecol. Manag.* **2006**, *14*, 445–461. [CrossRef]

106. McGill, B.M.; Altchenko, Y.; Hamilton, S.K.; Kenabatho, P.K.; Sylvester, S.R.; Villholth, K.G. Complex interactions between climate change, sanitation, and groundwater quality: A case study from Ramotswa, Botswana. *Hydrogeol. J.* **2019**, *27*, 997–1015. [CrossRef]

107. Maoyi, M.L.; Abiodun, B.J. How well does MPAS-atmosphere simulate the characteristics of the Botswana High? *Clim. Dyn.* **2021**, *57*, 2109–2128. [CrossRef]

108. Ramberg, L.; Hancock, P.; Lindholm, M.; Meyer, T.; Ringrose, S.; Sliva, J.; Van As, J.; Vander Post, C. Species diversity of the Okavango Delta, Botswana. *Aquat. Sci.* **2006**, *68*, 310–337. [CrossRef]

109. Asefa, S. Enhancing food access in Africa: The Botswana experience. *Stud. Comp. Int. Dev.* **1991**, *26*, 59–83. [CrossRef]

110. Vanderpost, C. Population change and environmental problems in the Mid-Boteti region of Botswana. *GeoJournal* **1995**, *35*, 521–529. [CrossRef]

111. Rahm, D.; Swatuk, L.; Matheny, E. Water Resource Management in Botswana: Balancing Sustainability and Economic Development. *Environ. Dev. Sustain.* **2006**, *8*, 157–183. [CrossRef]

112. Mogomotsi, P.K.; Sekelemani, A.; Mogomotsi, G.E.J. Climate change adaptation strategies of small-scale farmers in Ngamiland East, Botswana. *Clim. Chang.* **2020**, *159*, 441–460. [CrossRef]

113. Moreri, K.K.; Maphale, L.; Nkhwanana, N. Optimizing dispatch and home delivery services utilizing GIS in Botswana: Botswana Post case study. *Spat. Inf. Res.* **2017**, *25*, 565–573. [CrossRef]

114. Mpofu, G.; Darkoh, M.K.; Gwebu, T. Peri-urbanization landuse dynamics: An analysis of evolving patterns and their impacts on Gabane Village, Botswana. *GeoJournal* **2018**, *83*, 725–741. [CrossRef]

115. Lemenkova, P. Console-Based Mapping of Mongolia Using GMT Cartographic Scripting Toolset for Processing TerraClimate Data. *Geosciences* **2022**, *12*, 140. [CrossRef]

116. Ramberg, L.; Wolski, P.; Krah, M. Water balance and infiltration in a seasonal floodplain in the Okavango Delta, Botswana. *Wetlands* **2006**, *26*, 677–690. [CrossRef]

117. Lancaster, N. Response of Aeolian Processes and Landforms to Climate Change and Variability. In *Treatise on Geomorphology*, 2nd ed.; Schroder, J.J.F., Ed.; Academic Press: Oxford, UK, 2022; pp. 318–339. [CrossRef]

118. Radulović, M.; Brdar, S.; Mesaroš, M.; Lukić, T.; Savić, S.; Basarin, B.; Crnojević, V.; Pavić, D. Assessment of Groundwater Potential Zones Using GIS and Fuzzy AHP Techniques—A Case Study of the Titel Municipality (Northern Serbia). *ISPRS Int. J. Geo-Inf.* **2022**, *11*, 257. [CrossRef]

119. Lemenkova, P. Geophysical Mapping of Ghana Using Advanced Cartographic Tool GMT. *Kartogr. Geoinf.* **2021**, *20*, 16–37. [CrossRef]

120. Karydas, C.; Jiang, B. Scale Optimization in Topographic and Hydrographic Feature Mapping Using Fractal Analysis. *ISPRS Int. J. Geo-Inf.* **2020**, *9*, 631. [CrossRef]

121. Akinyemi, F.O.; Kgomo, M.O. Vegetation dynamics in African drylands: An assessment based on the Vegetation Degradation Index in an agro-pastoral region of Botswana. *Reg. Environ. Chang.* **2019**, *19*, 2027–2039. [CrossRef]
122. Ishimoto, Y.; Yabuta, S.; Kgokong, S.; Motsepe, M.; Tominaga, J.; Teramoto, S.; Konaka, T.; Mmopelwa, G.; Kawamitsu, Y.; Akashi, K.; et al. Environmental evaluation with greenhouse gas emissions and absorption based on life cycle assessment for a Jatropha cultivation system in frost- and drought-prone regions of Botswana. *Biomass Bioenergy* **2018**, *110*, 33–40. [CrossRef]

123. Klaučo, M.; Gregorová, B.; Stankov, U.; Marković, V.; Lemenkova, P. Determination of ecological significance based on geostatistical assessment: A case study from the Slovak Natura 2000 protected area. *Open Geosci.* **2013**, *5*, 28–42. [CrossRef]

124. Eze, P.N.; Madani, N.; Adoko, A.C. Multivariate Mapping of Heavy Metals Spatial Contamination in a Cu–Ni Exploration Field (Botswana) Using Turning Bands Co-simulation Algorithm. *Nat. Resour. Res.* **2019**, *28*, 109–124. [CrossRef]

125. Lemenkova, P. Handling Dataset with Geophysical and Geological Variables on the Bolivian Andes by the GMT Scripts. *Data* **2022**, *7*, 74. [CrossRef]

126. Chipanshi, A.C.; Chanda, R.; Totolo, O. Vulnerability Assessment of the Maize and Sorghum Crops to Climate Change in Botswana. *Clim. Chang.* **2003**, *61*, 339–360. [CrossRef]

127. Klaučo, M.; Gregorová, B.; Koleda, P.; Stankov, U.; Marković, V.; Lemenkova, P. Land planning as a support for sustainable development based on tourism: A case study of Slovak Rural Region. *Environ. Eng. Manag. J.* **2017**, *16*, 449–458. [CrossRef]

128. Maruatona, P.B.; Moses, O. Assessment of the onset, cessation, and duration of rainfall season over Botswana. *Model. Earth Syst. Environ.* **2022**, *8*, 1657–1668. [CrossRef]

129. Parida, B.; Moalafhi, D. Regional rainfall frequency analysis for Botswana using L-Moments and radial basis function network. *Phys. Chem. Earth Parts A/B/C* **2008**, *33*, 614–620. [CrossRef]

130. Moses, O.; Ramotonto, S. Assessing forecasting models on prediction of the tropical cyclone Dineo and the associated rainfall over Botswana. *Weather. Clim. Extrem.* **2018**, *21*, 102–109. [CrossRef]

131. Dahlberg, A.C. Vegetation diversity and change in relation to land use, soil and rainfall—A case study from North-East District, Botswana. *J. Arid. Environ.* **2000**, *44*, 19–40. [CrossRef]