ABSTRACT
Batwa communities inhabiting Uganda's complex topography lack timely access to climate information services (CISs) due to scarce meteorological stations. Tailoring historical satellite data to local spatio-temporal dimensions is central to catalyzing sustainable development goals (SDGs). This study employed satellite data from 1983 to 2023 to inform sustainable and climate-smart decision-making for the Batwa communities in southwestern Uganda. Daily precipitation (DPR) and air temperature data were obtained from the Climate Hazard Group Infrared Precipitation with Station (CHIRPS) and the NASA Power website (NASAPW), respectively. ClimPACT3 model was employed to generate climate indices (CIs) while R-Instat was utilized to perform statistical analysis. Results revealed that the area experiences wetter than normal conditions, with precipitation (PR) ranges of 1,006–1,489 mm. A significant (p-value < 0.05) increase in the annual sum of daily precipitation (ASDPR) > the 95th percentile was observed. The daily minimum temperature (DTMin.) and maximum temperature (DTMax.) ranged from 13.1 to 14.5 °C and 24.0 to 26.4 °C, respectively. The current study connected the implications of the results to the broader SDGs. The synergies between SDG 13 and other SDGs were also established. Therefore, the study's outcomes are imperative in advancing climate resilience for the agricultural, health, and water resources sectors.
HIGHLIGHTS
Satellite datasets were used to assess fluctuations in CIs, while RH-tests were used for DQCD.
ClimPACT3 and R-Instat used for CIs and statistical analyses, respectively.
Implications of the findings connected to the broader SDGs and climate resilience.
Knowledge of variations in extreme CIs could advance SDG 13 (climate action).
Implementing SDG 13 could contribute to the achievement of other SDGs.
INTRODUCTION
A valuable tool for improving sustainable development is historical weather and climate data (HWCD), which, when properly utilized, can enhance the lives of vulnerable communities susceptible to climate change impacts (CCIs) (Dinku et al. 2018). Batwa communities, in the highlands of Kigezi-southwestern Uganda, are particularly vulnerable to climatic threats such as flash floods (FFs), soil erosion, and disease incidence (van Bavel et al. 2020). The frequency of these climate hazards has intensified of late, posing more danger to this Indigenous group (Satyal et al. 2021). The fact that over 70% of the Batwa households live in perpetual deficiency of basic needs such as food, housing, clothing, and education raises concern (Scarpa et al. 2021), especially in the face of a changing climate. The hilly, solitary, and remote terrains in which they live further exacerbate the challenge (Luliro et al. 2022). Harnessing HWCD is essential for these indigenous communities since it facilitates climate-sensitive early-warning system design, planning, research, and the identification of climate risks such as catastrophic weather events. Additionally, the HWCD provides a baseline for future climate projections as well as establishing the climate change (CC) trajectory (Ngoma et al. 2021).
Considering statistics indicating that Africa is the continent most exposed to CCIs (IPCC 2023), harnessing the HWCD becomes even more essential to investigate changes in climate indices (CIs) for early-warning purposes. This should receive the utmost interest since the majority of African economies depend on climate-sensitive industries like agriculture for development. About 80% of agriculture worldwide is rain-fed, according to Dinku (2019), and in Africa, this percentage exceeds 95%. In Africa, agriculture contributes significantly to GDP and provides a living for a larger proportion of the continent's rural population. The application of HWCD in location-specific local contexts would be crucial for African nations to meet their nationally established development plans and the 2030 Agenda SDGs (Dorward et al. 2020). This is due to the fact that successful CC adaptation techniques are frequently context-specific, taking into account local sociocultural, sectoral, and geographic variables that differ from place to place.
Exploitation of the HWCD in Africa is inadequate because the majority of National Meteorological and Hydrological Services (NMHS) gather the data with poor quality checks, thus reducing its relevance to the local community (Dinku 2019). Faniriantsoa & Dinku (2022) reported that the primary challenge to climate data availability and quality in Africa is attributed to a sparse and uneven distribution network of weather stations. The coverage is typically lower in rural areas, where livelihoods are most vulnerable to the adverse impacts of CC. Other problems included low-quality station observations with many missing data points and the inaccessibility of the NMHS's available data. Dinku (2019) claims that insufficient dissemination tools and capacity, scarce resources, and other factors make data accessibility in Uganda a major challenge. In order to overcome these challenges, the present investigation employed rainfall and temperature datasets from CHIRPS and NASAPW, respectively. The selection was based on the datasets' free accessibility, high resolution, and global coverage, which provide a more dependable and contextualized outlook at the local level (Dhanesh et al. 2020; Halimi et al. 2023; Tan et al. 2023; Du et al. 2024).
RESEARCH METHODOLOGY
Study area description
The study focused on Kabale, Kanungu, and Kisoro (KAKAKI) districts of southwestern Uganda (see Figure 1). The geographical latitudes and longitudes of the study area are −1.5° S to −0.3° S and 29.4° E to 30.3° E, respectively. The area is inhabited by a significantly greater number of Batwa communities. The study area is bounded by Ntungamo to the east, the Democratic Republic of Congo to the west, Rukungiri to the north, and Rwanda to the south. According to Satyal et al. (2021), Batwa settlements are scattered throughout the hilly, solitary, and remote parts of these districts, comprising 10–20 households. These areas are frequently susceptible to climatic threats such as FFs, soil erosion, and disease incidence (van Bavel et al. 2020). A study by Scarpa et al. (2021) revealed that about 70% deficiency in basic living standards prevails among the Batwa people. The study location is situated in the tropical Savannah climate zone, with a bimodal (MAM and SOND) rainy season each year. The major crops produced under the rain-fed farming system in the area include Irish potatoes, bananas, coffee, sweet potatoes, beans, and cereals (Scarpa et al. 2021).
Acquisition of the HWCD
Retrieving precipitation data
The traditional methods for measuring rainfall at the local level are rain gauges and weather radars (Shahid et al. 2021). However, because of their limited distribution throughout the research area, it is difficult to record the spatio-temporal variability in precipitation, particularly in the Kigezi region's mountainous areas. The employed daily precipitation data was obtained from CHIRPS, which is accessible at https://app.climateengine.org/climateEngine. The preference for CHIRPS data is attributable to its high resolution (0.05° × 0.05°) and its integration of satellite imagery with ground station data, offering a more comprehensive raster and reliability (Basheer & Elagib 2019; Dhanesh et al. 2020; Du et al. 2024). A study by Shahid et al. (2021) affirms that CHIRPS precipitation data reliably complements ground station data in addition to providing continuous precipitation monitoring over large temporal and spatial scales. Moreover, CHIRPS provides free daily datasets available for extended periods dating back to 1981. To download the precipitation data (PRD), the area of interest was initially zoomed in and a polygon drawn to represent it on the map. Additionally, the CHIRPS 4.8 km dataset was selected under the ‘climate and hydrology’ section, and the ‘variable’ section was modified to ‘precipitation.’ The units were set to ‘milimetres,’ while the ‘statistics’ section was set to ‘total’ to indicate cumulative precipitation. The required data's time period was adjusted to range from January 1, 1983 to December 31, 2023. A graph was generated by clicking on ‘Get Time Series,’ and after clicking on the ‘download button’ in the top right corner of the graph, the PRD was downloaded as a CSV file (Obura et al. 2024c). These steps were repeated for each of the three KAKAKI districts containing the Batwa. Eventually, an average of the three districts was obtained to serve as the representative rainfall dataset.
Obtaining temperature data
The NASAPW was exploited to retrieve the DTMin. and DTMax. datasets at 2 m from https://power.larc.nasa.gov/data-access-viewer/. The NASAPW is a comprehensive, freely accessible, and dependable satellite data source for temperature with high resolution and global coverage, which enables studies to be conducted even on a local scale. NASA implements stringent quality control procedures to guarantee the accuracy of its data. Furthermore, numerous studies have reported on the validity and reliability of NASAPW TED (Aboelkhair et al. 2019; Marzouk 2021; Carrara et al. 2023; Halimi et al. 2023; Tan et al. 2023). The user community ‘Agroclimatology’ was chosen after opening the ‘data access viewer’ on the NASAPW in order to download temperature data (TED). Each of the three KAKAKI districts' lat/long coordinates was input, and the temporal average was set to ‘daily.’ The output file format was set to CSV, and the timeframe was set from January 1, 1983, to December 31, 2023. The minimum and maximum temperatures (MMTs) were nominated for download when the ‘Temperatures/Thermal IR Flux’ section under the ‘Select Parameters’ section was opened. After pressing the ‘Submit’ button, daily native resolution NASA Power (NASAP) CERES/MERRA-2 data was downloaded in a CSV format (Obura et al. 2024c). The MMTs were represented by averaging the three district datasets.
Data quality control diagnostics
The consistency and quality of the daily input data were checked before the indices were computed. A crucial requirement for a robust analysis of climate time series is homogeneity. This study adopted free and user-friendly RH-tests software to perform the data quality control diagnostics (DQCD). Based on the penalised maximal t-test (PMT) or f-test (PMF), the programme can identify and correct multiple change points in a time series, as demonstrated by J. Du et al. (2020) and validated by Liu et al. (2017). Reference stations are necessary for PMT when conducting a homogeneity study, but PMF can be employed as an absolute approach, which means it can be used in isolation or in the absence of neighbouring stations for comparison. Therefore, the PMF method was used in this investigation.
Determining the CIs
The CIs required to comprehend the climatology of the study location were generated utilising ClimPACT v3.1.6 in this investigation. ClimPACT3 is an open-source tool programmed by the Expert Team on Sector-specific Climate Indices (ET-SCI) for the World Meteorological Organisation (WMO). ClimPACT v3.1.6 was chosen because it allows for customization based on local conditions, enhancing the relevance of findings for a specific region (Murara & Mendonça 2019; Dubey et al. 2022). It is also incorporated with diverse internationally recognized CIs for determining precipitation extremes, heat waves, and cold waves, making it a reliable tool for long-term climate data analysis (Nakaegawa & Murazaki 2022). Furthermore, ClimPACT integrates with various climate data formats, promoting wide applicability, supporting research, and climate adaptation planning (WMO 2016). It employs daily parameters against monthly datasets because several significant pieces of information are concealed in monthly climate data, which may be crucial to multiple fields, including agriculture, health, and water resources. Employing daily data to produce CIs provides a comprehensive understanding of multiple variables, including the periodicity and intensity of extreme rainfall and measurements of extremely wet, dry, hot, or cold periods that influence society's economy.
Furthermore, the Mann–Kendall test at 0.05 and Sen's slope estimator were adopted for trend analysis (Agarwal et al. 2021; van der Walt & Fitchett 2021; Şen & Şişman 2023). Ten (10) indices (see Table 1) implemented in this investigation were nominated based on their significance to health (H), agriculture and food security (AFS), as well as water resources and hydrology (WRH). By capturing significant trends and anomalies, CIs helps in summarizing complex climate data to inform decision-making in sectors such as agriculture, health, and water management. In addition, the Standardized Precipitation Evapotranspiration Index (SPEI) was selected in contrast to the Standardized Precipitation Index (SPI) due to its ability to amalgamate the multi-temporal nature of the SPI with its computational simplicity and sensitivity to changes in evaporative demand ignited by temperature trends and fluctuations (Vicente-Serrano et al. 2010). The WMO approved the SPEI as the reference drought index for more effective determination of drought sensitivity and climate risk management of drought hazards on sensitive systems (Mishra & Aadhar 2021).
ET-SCI index . | Description . | Units . | Sector(s) . |
---|---|---|---|
SPEI | A drought metric based on Evaporation and PR. | Unit-less | H, AFS, WRH |
Consecutive Dry Days (CDDs) | Maximum number of CDDs when PR < 1.0 mm. | Days | H, AFS, WRH |
Consecutive Wet Days (CWDs) | Maximum number of CWDs per year when PR ≥ 1.0 mm. | Days | H, AFS, WRH |
Annual Total wet-day rainfall (PRCPTOT) | Sum of daily PR ≥ 1.0 mm | mm | AFS, WRH |
PR20mm | Days on which PR ≥ 20 mm | Days | AFS, WRH |
Max 1-day PR (Rx1day) | Maximum amount of rain that falls in 1-day | mm | H, AFS, WRH |
PR95p | ASDPR > the 95th Percentile. | mm | AFS, WRH |
DTMin. | Daily minimum temperature | °C | H, AFS |
Daily Temperature Range (DTR) | Daily range of maximum and minimum temperature | °C | H, AFS, WRH |
DTMax. | Daily maximum temperature | °C | H, AFS |
ET-SCI index . | Description . | Units . | Sector(s) . |
---|---|---|---|
SPEI | A drought metric based on Evaporation and PR. | Unit-less | H, AFS, WRH |
Consecutive Dry Days (CDDs) | Maximum number of CDDs when PR < 1.0 mm. | Days | H, AFS, WRH |
Consecutive Wet Days (CWDs) | Maximum number of CWDs per year when PR ≥ 1.0 mm. | Days | H, AFS, WRH |
Annual Total wet-day rainfall (PRCPTOT) | Sum of daily PR ≥ 1.0 mm | mm | AFS, WRH |
PR20mm | Days on which PR ≥ 20 mm | Days | AFS, WRH |
Max 1-day PR (Rx1day) | Maximum amount of rain that falls in 1-day | mm | H, AFS, WRH |
PR95p | ASDPR > the 95th Percentile. | mm | AFS, WRH |
DTMin. | Daily minimum temperature | °C | H, AFS |
Daily Temperature Range (DTR) | Daily range of maximum and minimum temperature | °C | H, AFS, WRH |
DTMax. | Daily maximum temperature | °C | H, AFS |
Statistical analysis of the HWCD
The HWCD in CSV format for the study location was imported and analysed in R-Instat 0.7.16.50. The data was checked, reshaped, and arranged using the ‘prepare’ tab. Then, the ‘climatic’ tab was applied to prepare climatic summaries and determine the onset and cessation of the rains and the duration of each season. Furthermore, the statistical analysis was performed for the bimodal rainy seasons using R-Instat. The criterion to establish the rainy season's onset was set to three consecutive days with total precipitation (TPR) ≥ 20 mm and the absence of nine consecutive dry days (CDDs) in the first 21 days. Onset dates were established for the MAM and SOND seasons. The last day of the season to receive the PR > 10 mm was designated as the the season’ cessation. Finally, the output from R-Instat was exported to an Excel spreadsheet for additional modification, where plots for seasonal TPR and the extent of each season were generated. R-Instat was selected for this study because of its easy-to-use interface, which makes complicated statistical processes simpler. It offers quality control tools and also accomodates huge climatic datasets. Advanced statistical techniques like regression, time series analysis, and climatic trend evaluations are integrated into the software, simplifying analysis. R-Instat is an open-source and free software (Fundi et al. 2024), making it easily accessible for use in this study.
RESULTS AND DISCUSSION
Exploration of the SPEI temporal patterns
SPEI range . | Drought description . |
---|---|
SPEI ≤ −2.00 | Extremely dry |
−1.99 ≤ SPEI ≤ −1.50 | Severely dry |
−1.49 ≤ SPEI ≤ −1.00 | Moderately dry |
−0.99 ≤ SPEI ≤0.99 | Near normal |
1.00 ≤ SPEI ≤1.49 | Moderately wet |
1.50 ≤ SPEI ≤1.99 | Severely wet |
SPEI ≥ 2.00 | Extremely wet |
SPEI range . | Drought description . |
---|---|
SPEI ≤ −2.00 | Extremely dry |
−1.99 ≤ SPEI ≤ −1.50 | Severely dry |
−1.49 ≤ SPEI ≤ −1.00 | Moderately dry |
−0.99 ≤ SPEI ≤0.99 | Near normal |
1.00 ≤ SPEI ≤1.49 | Moderately wet |
1.50 ≤ SPEI ≤1.99 | Severely wet |
SPEI ≥ 2.00 | Extremely wet |
A significantly positive SPEI could be valuable for agricultural productivity by increasing soil moisture and water availability, favouring crop growth. This has a positive implication on SDG 2 (Zero Hunger) and especially contributes to achieving SDGTs 2.1, 2.2, and 2.3. In the WRH sector, a positive SPEI trend could improve water supplies but also calls for the improved infrastructural need to manage the excess water. This has direct implications for the efforts to achieve SDGTs 6.1, 6.2, and 6.3. Healthwise, wetter than usual conditions pose significant health risks, including the emergence of water- and vector-borne diseases, communicable diseases, and psychological disorders due to the evacuation of people from flood-vulnerable zones. A study by van Bavel et al. (2020) claims that increased moisture, which supports larger mosquito populations and other insects that spread diseases like dengue fever and the West Nile virus, is a contributing factor in the proliferation of vector-borne diseases in wetter environments. This calls for concerted and prioritized efforts from relevant stakeholders since a laxity directly hampers the achievement of target 3.3 of SDG 3 (Healthy lives and well-being for all).
Temporal variations in PR
Number of consistent dry and wet days
MAM and SOND seasonal onset, cessation, length, and TPR
The MAM's onset window ranged from March 1 to March 31; its cessation window spanned from May 1 to June 30. The SOND season's onset window stretched from September 1 to September 30, and its cessation window was from November 15 to December 15. The MAM and SOND season commencement and cessation dates for the most recent 20 years of research are presented in Table 3, which was generated using R-Instat software.
Year . | MAM season . | SOND season . | ||
---|---|---|---|---|
Onset . | Cessation . | Onset . | Cessation . | |
2004 | 01/03/2004 | 30/05/2004 | 08/09/2004 | 29/11/2004 |
2005 | 03/03/2005 | 30/05/2005 | 05/09/2005 | 22/11/2005 |
2006 | 08/03/2006 | 15/05/2006 | 01/09/2006 | 15/12/2006 |
2007 | 01/03/2007 | 10/06/2007 | 05/09/2007 | 08/12/2007 |
2008 | 05/03/2008 | 06/06/2008 | NA | 04/12/2008 |
2009 | 02/03/2009 | 12/05/2009 | 25/09/2009 | 10/12/2009 |
2010 | 04/03/2010 | 05/06/2010 | 06/09/2010 | 11/12/2010 |
2011 | 11/03/2011 | 28/06/2011 | 03/09/2011 | 06/12/2011 |
2012 | 01/03/2012 | 10/06/2012 | 02/09/2012 | 15/12/2012 |
2013 | 02/03/2013 | 10/05/2013 | 10/09/2013 | 11/12/2013 |
2014 | 11/03/2014 | 12/05/2014 | 05/09/2014 | 08/12/2014 |
2015 | 01/03/2015 | 04/06/2015 | 01/09/2015 | 10/12/2015 |
2016 | 03/03/2016 | 05/05/2016 | 15/09/2016 | 14/12/2016 |
2017 | 22/03/2017 | 06/05/2017 | 02/09/2017 | 30/11/2017 |
2018 | 03/03/2018 | 22/05/2018 | 03/09/2018 | 11/12/2018 |
2019 | 14/03/2019 | 27/06/2019 | 01/09/2019 | 07/12/2019 |
2020 | 01/03/2020 | 06/05/2020 | 25/09/2020 | 07/12/2020 |
2021 | 03/03/2021 | 11/05/2021 | 01/09/2021 | 11/12/2021 |
2022 | 21/03/2022 | 28/05/2022 | 05/09/2022 | 12/12/2022 |
2023 | 01/03/2023 | 04/05/2023 | 09/09/2023 | 30/11/2023 |
Year . | MAM season . | SOND season . | ||
---|---|---|---|---|
Onset . | Cessation . | Onset . | Cessation . | |
2004 | 01/03/2004 | 30/05/2004 | 08/09/2004 | 29/11/2004 |
2005 | 03/03/2005 | 30/05/2005 | 05/09/2005 | 22/11/2005 |
2006 | 08/03/2006 | 15/05/2006 | 01/09/2006 | 15/12/2006 |
2007 | 01/03/2007 | 10/06/2007 | 05/09/2007 | 08/12/2007 |
2008 | 05/03/2008 | 06/06/2008 | NA | 04/12/2008 |
2009 | 02/03/2009 | 12/05/2009 | 25/09/2009 | 10/12/2009 |
2010 | 04/03/2010 | 05/06/2010 | 06/09/2010 | 11/12/2010 |
2011 | 11/03/2011 | 28/06/2011 | 03/09/2011 | 06/12/2011 |
2012 | 01/03/2012 | 10/06/2012 | 02/09/2012 | 15/12/2012 |
2013 | 02/03/2013 | 10/05/2013 | 10/09/2013 | 11/12/2013 |
2014 | 11/03/2014 | 12/05/2014 | 05/09/2014 | 08/12/2014 |
2015 | 01/03/2015 | 04/06/2015 | 01/09/2015 | 10/12/2015 |
2016 | 03/03/2016 | 05/05/2016 | 15/09/2016 | 14/12/2016 |
2017 | 22/03/2017 | 06/05/2017 | 02/09/2017 | 30/11/2017 |
2018 | 03/03/2018 | 22/05/2018 | 03/09/2018 | 11/12/2018 |
2019 | 14/03/2019 | 27/06/2019 | 01/09/2019 | 07/12/2019 |
2020 | 01/03/2020 | 06/05/2020 | 25/09/2020 | 07/12/2020 |
2021 | 03/03/2021 | 11/05/2021 | 01/09/2021 | 11/12/2021 |
2022 | 21/03/2022 | 28/05/2022 | 05/09/2022 | 12/12/2022 |
2023 | 01/03/2023 | 04/05/2023 | 09/09/2023 | 30/11/2023 |
NA = no date fitting the conditions for onset of the season.
From 1983 to 2023, the extent of the MAM season fluctuated more than the SOND season, with ranges between 40–110 days and 73–105 days, respectively (see Figure 5(c)). The shortest MAM season (40 days) occurred in 1984, whereas that for SOND (73 days) was experienced in 2020. The year 1995 recorded the longest MAM season (110 days) from March 1 to June 19, while 2006 experienced the longest SOND season (105 days), with September 1 as the onset and December 15 as the cessation of the season. The increasing MAM seasonal variability poses a risk for climatic hazards like FFs (Ainuddin et al. 2017), which severely impact infrastructure, livelihoods, and food security. These consequently impede sustainable development. Furthermore, the results also depicted more PR inconsistency in the MAM season than in the SOND season, fluctuating between 308–655 mm and 378–623 mm, respectively. The area experiences MAM and SOND seasons' average near-normal rainfall of 404 and 531 mm, respectively. Despite this, the MAM rains were more intense, with an intensity of 135 mm/month compared to SOND's 132 mm/month. Although the MAM season for the year 1984 was the shortest (40 days), it received the highest intensity of rainfall of 253 mm/month. This can trigger climatic hazards for communities. The year 2017 recorded the lowest MAM PR totals (308 mm), whereas the year 2018 received the greatest MAM rains (655 mm). The SOND season's minimum PR totals (378 mm) occurred in 1993, while the largest amounts (623 mm) were received in 2010. Heavy PR over shorter periods exacerbate land degradation (Tofu et al. 2022), destroying terrestrial ecosystems and thus hampering the efforts to achieve SDG 15 (Life on Land) especially targets 15.1–15.6 and 15.8. These findings necessitate the urgency for sustainable land management practices to be sought to preserve ecosystems and promote sustainable agriculture in the face of a changing climate.
Time series evolution of temperature
Implications of SDG 13 (climate action) on other SDGs
CONCLUSIONS
The current study utilized the CHIRPS and NASAPW remote sensing data to investigate fluctuations in extreme CIs. The ClimPACT3 and R-Instat models were applied in the assessment process. The findings show that the region generally receives wetter than normal weather, with PR values from 1,006 to 1,489 mm. There was a significant (p-value < 0.05) increase in the ASDPR > the 95th percentile. Conversely, the DTMin. and DTMax. varied from 13.1 to 14.5 °C and 24.0 to 26.4 °C, respectively. These observations are particularly imperative for the agricultural, health, and water resources sectors. Furthermore, CISs would aid in the achievement of SDG 1, SDG 2, and SDG 13 targets by enhancing food security and climate resilience in Batwa farmer communities. For the reason of making well-informed decisions about water management, the participants in the hydrological and water resources sectors need information on temperature extremes, rainfall variability, and SPEI trends. In order to prepare for an epidemic of food-borne, water-borne, vector-borne, or air-borne diseases, which could help achieve SDG 3, SDG 6, and SDG 11 targets, these outputs are also critical to the health sector. In summary, even though the CHIRPS and NASAPW time series datasets have been modified to improve homogeneity, it is vital to keep in mind that certain parts of these records may still be non-homogeneous when evaluating index changes. To increase reliability, future research should think about using precise station observation data.
AUTHOR CONTRIBUTIONS
M.E. contributed to methodology, formal analysis, writing – original draft; W.J. contributed to data curation, writing – review and editing; and D.O. contributed to conceptualization, visualization, writing – original draft; writing – review and editing.
DATA AVAILABILITY STATEMENT
All relevant data are included in the paper or its Supplementary Information.
CONFLICT OF INTEREST
The authors declare there is no conflict.