Local Perception of Drivers of Land-Use and LandCover Change Dynamics across Dedza District , Central Malawi Region

Research on Land Use and Land Cover (LULC) dynamics, and an understanding of the drivers responsible for these changes, are very crucial for modelling future LULC changes and the formulation of sustainable and robust land-management strategies and policy decisions. This study adopted a mixed method consisting of remote sensing and Geographic Information System (GIS)-based analysis, focus-group discussions, key informant interviews, and semi-structured interviews covering 586 households to assess LULC dynamics and associated LULC change drivers across the Dedza district, a central region of Malawi. GIS-based analysis of remotely sensed data revealed that barren land and built-up areas extensively increased at the expense of agricultural and forest land between 1991 and 2015. Analysis of the household-survey results revealed that the perceptions of respondents tended to validate the observed patterns during the remotely sensed data-analysis phase of the research, with 57.3% (n = 586) of the respondents reporting a decline in agricultural land use, and 87.4% (n = 586) observing a decline in forest areas in the district. Furthermore, firewood collection, charcoal production, population growth, and poverty were identified as the key drivers of these observed LULC changes in the study area. Undoubtedly, education has emerged as a significant factor influencing respondents’ perceptions of these drivers of LULC changes. However, unsustainable LULC changes observed in this study have negative implications on rural livelihoods and natural-resource management. Owing to the critical role that LULC dynamics play to rural livelihoods and the ecosystem, this study recommends further research to establish the consequences of these changes. The present study and future research will support decision makers and planners in the design of tenable and coherent land-management strategies.


Introduction
Land-use and land-cover (LULC) change has become a key research-priority area, attracting much interest from the global scientific community since the 1970s [1][2][3]. Particularly, the attention on LULC dynamics occurring at the local scale has arisen due to an inherent ecosystem, and socioeconomic impact at the national, regional, and even global level [4,5]. Natural causes and anthropogenic activities are responsible for LULC dynamics changes globally, with the latter overriding natural causes [6,7]. These changes are described by complex multitemporal and scale interactions of social, demographic, economic, institutional, and environmental factors [8][9][10][11]. These changes have serious socioeconomic and environmental impact on rural livelihoods in many regions of Sub-Saharan Africa (SSA) [12]. In some parts of the SSA region, population growth, high poverty levels, settlements, fuelwood, charcoal production, and agricultural expansion were reported as contributory factors for LULC changes [13][14][15][16][17][18]. More research with regard to location, nature, magnitude, extent, and rate of land-use and land-cover dynamics is still required in the context of SSA, where high population growth coupled with infertile land and overexploitation of other natural resources such water and forests is prevailing [19].
Malawi's economy is entirely dependent on agriculture and other related sectors, especially forests and fisheries. Due to its reliance on rain-fed agriculture and exposure to floods and droughts, Malawi is among southern Africa's most climate-change-vulnerable countries [20]. Almost 85% of Malawi's population live in rural and marginalized areas, and approximately 80% of this population entirely depend on natural-resource endowments for their subsistence, household income, and livelihoods [21][22][23][24]. The high dependence on natural resources such as land, forests, and water puts pressure on these resources, leading to overexploitation, forest degradation, and deforestation [25,26]. Recent studies have revealed that deforestation and forest degradation in Malawi are due to uncontrolled firewood collection, infrastructure development, agriculture expansion, illegal charcoal production, shifting cultivation, urbanization, high population, and tobacco-curing by smallholder farmers and estate owners [23,26] Like any other country in the Sub-Saharan Africa (SSA) region, Malawi's LULC has experienced rapid and extensive changes over the past decades due to significant transformations caused by human-environment interactions [27]. Despite the fact that few studies on LULC changes have been done in Malawi, research on the factors contributing to these changes at the national and even local level remains scant. Thus, few studies have explained LULC change dynamics at the national level [28][29][30][31][32]. Studies on LULC dynamics and the associated drivers on the local scale are vital for seeking viable, feasible, appropriate, and coherent natural-resource management strategies. Several researchers have emphasized that understanding LULC drivers is a perplexing question in global science, and these drivers are still a contentious issue; further research is indispensable [33][34][35]. The causes of LULC changes are intricate and dynamic, and they vary from one place to another [36]. In other words, globally identified drivers of LULC changes are location-specific, varying from region to region depending on the socioeconomic and biophysical factors prevailing that location. It is worth noting that LULC change drivers are also time-specific. For instance, a driver identified 10 years ago may not be valid in recent times if remedial solutions are put in place by the actors. It is, therefore, impossible to generalize that LULC trends/changes occurring on a broader spatial scale and the drivers influencing these changes are inherent landscapes [35,37]. Examination of LULC driver dynamics is a requisite as far as resolving environmental and socioeconomic challenges, biodiversity conservation ,reduction and management of LUCC changes impacts and consequences at local, national, regional and global level is concerned [38,39].
It is worthwhile noting that inclusive research on the drivers and impacts of LULC dynamics in Dedza is beneficial to readily comprehend the inter-relationships between locals and natural resources. Any management intervention strategies to properly address the drivers of LULC changes and the development of sustainable land-use systems in the study area should begin with local empirical evidence and understanding the underlying drivers of changing LULC. A profound understanding of the complex interdependence between LULC changes and rural livelihoods, together with the coping strategies that local communities use to address such changes, are fundamental for decision-making by policymakers, planners, and other stakeholders [13]. Estimating the rate, nature, type, and pattern of LULC changes in any landscape, as well as understanding factors that influence these changes, are also essential for projecting future changes [40,41].
Remote-sensing (RS) and GIS technologies only identify the nature, extent, and rate of LULC changes on the landscape; however, they do not provide an explanation about the underlying causes of LULC dynamics on the landscape [42,43]. Despite this, RS has demonstrated its effectiveness and applicability in investigating the relationship that exists between people and the environment in which they live [44]. Therefore, this study aims at quantifying LULC changes and assessing the local perceptions of drivers of LULC change between 1991 and 2015 in Dedza. Thus, the study captured local communities' perceptions of LULC change trends and the drivers of these changes in the study area. Some researchers have reported that observed LULC dynamics on any landscape is a reflection of aggregated decisions at the household level in response to policy and an institutional environment over a period of time [45][46][47]. The findings of this study are envisioned to form the basis for a robust understanding of the LULC change dynamics that planners, environmentalists, decision-makers, and other stakeholders could use in formulating sound management and environmental planning strategies, or guidelines for the maintenance of ecosystem services, and conservation and utilization of natural resources in Dedza or alternative districts with similar settings.

Study Area
The study was conducted in Dedza, located in the central region of Malawi, bordering Lilongwe district, Ntcheu to the south, Mangochi to the east, Salima to the northeast, and Mozambique to the west ( Figure 1). The district covers a geographical area of about 362,400 ha [48,49]. Physiography is characterized by uplands and lowlands with uneven terrain. The district is divided into three topographic zones, namely, the Lilongwe plain (altitude, 1100-1300 m), the Dedza highlands (1200-2200 m), and the Dedza escarpments (1000-1500 m). The district has a subtropical highland climate [50]. Mean annual temperatures are relatively low and fluctuate between 14 and 21 °C, with an average temperature of 15.5 °C (the coldest months are June and July, while November is the hottest month). Rainfall occurs between the months of November and March, with a mean annual rainfall ranging from 800 to 1200 mm. The district has experienced climate-related disasters and extreme events such as floods and droughts [51]. The district is characterized by generally ferruginous soils that are deep and brown to reddish in color [52]. Clay and sandy loam soils are predominant in the study area [49,53]. Agriculture is the major land use in Dedza, with major crops grown in the area being maize (Zea mays), Irish potatoes (Solanum tuberosum), sweet potatoes (Ipomoea batatas), groundnuts (Arachis hypogaea L.), beans (Phaseolus vulgaris L), and soybeans (Glycine max). Rice and cotton are also grown along the lakeshore and valleys. People in the district also keep livestock comprising of cattle, goats, pigs, sheep, and poultry. The economy and livelihoods of the majority of the communities of the study area are primarily based on natural resources, especially land, forests, and water [49,51]. Other economic activities and sources of livelihood strategies include small and medium enterprises (SMEs), arts and crafts, quarrying, and fishing. The district has three land-tenure systems, namely, government land, customary land, and private-leasehold land. Dedza has an estimated population of 624,445, with an annual population growth rate of 2.6% [51]. It is one of the most densely populated districts in Malawi, with a population density of 172 persons per km 2 compared to the national average of 139 persons per km 2 . The average family size in the studied landscape is 6 persons against the national average of 4.4 persons per household.

Data Acquisition and Preprocessing
Geospatial and remote-sensing data are reliable sources for understanding and ascertaining the drivers of LULC changes of any landscape [54]. In this study, change-detection analysis using multiple sets of spatiotemporal Landsat images for 1991, 2001, and 2015 was used to establish LULC changes in Dedza. Table 1 summarizes the characteristics of the multitemporal satellite data used in this research. ArcGIS 10.6 and ERDAS IMAGINE 9.3 software were used to perform standard imageprocessing techniques, including extraction, geometric correction or georeferencing, atmospheric correction, topographic correction, layer stacking (band selection and combination), image enhancement, and subsetting (clipping). The three images were also registered to a common Universal Transverse Mercator (UTM) co-ordinate system, Zone 36S, with World Geocoded System (UTM WGS 84) projection parameters.

Image Classsification and Land-Use and Land-Cover Dynamics
Images were classified using hybrid classification that combines supervised and unsupervised classification algorithms. The two methods were used to reduce spectral reflectance noise, especially singling out agricultural land from built-up areas and bare land. A Maximum Likelihood Classification (MLC) algorithm was performed for each image (Equations 1 and 2). Studies have shown that MLC is the most common, successful, and widely adopted classification algorithm [55][56][57][58]. A classification scheme of 6 classes was developed based on physiographical knowledge of the study area, supporting ancillary data, researchers' prior local knowledge, and visual interpretation using the historical function of Google Earth. The 6 LULC classes were categorized as water bodies, wetlands, agricultural land, forest, built-up areas and barren land ( Table 2). A stratified random sampling method was employed to collect 221 points for accuracy assessment. Google Earth images were used to extract the reference data. Accuracy assessment was determined using the kappa coefficient, overall accuracy, producer and user accuracy, which were derived from the error (confusion) matrix as discussed in References [59] and [60]. In order to continue with LULC analysis, the 2015 LULC map was subjected to a minimum of 85% overall accuracy as recommended by References [61] and [62]. The classified 2015 images were used as reference to classify historical images. In this case, the used signatures for the 2015 images were superimposed on older images. Considerations were made to ensure that the images were captured at comparable phenological dates during the study period. In addition, historical images (1991 and 2001) were further visually interpreted, taking into account image tone, texture, shape, and class patterns. Table 2. Land-use land-cover (LULC) classes used in this study.

Water bodies
Rivers, permanent open water, lakes, ponds, reservoirs.
Wetland Permanent and seasonal grasslands along lake, river, and streams, marshy land and swamps. Agricultural land All cultivated and uncultivated agricultural lands areas, such as farmlands, crop fields including fallow lands/plots, and horticultural lands.

Forest
Protected forests, plantations, deciduous forests, mixed forest lands, and forests on customary land. Built-up areas Residential, commercial and service, industrial, socioeconomic infrastructure, and mixed urban and other urban, transportation, roads, and airports.
Barren land Areas around and within forest-protected areas with no or very little vegetation cover, including exposed soils, stock quarry, rocks, landfill sites, and areas of active excavation.
LULC change analysis was determined using a post-classification comparison (PCC) technique, and this resulted in a cross-tabulation (transition) matrix. The LULC change-transition matrix was computed using the overlay procedure in ArcGIS in order to quantify the area converted from a particular LULC class to another LULC category during the study period. The annual rate of change was also determined using the procedure by References [63][64][65]. Equation (1) provides a benchmark for comparing LULC changes that are not sensitive to differing periods between study periods.
where is the annual rate of change for each class, and and are areas of each LULC class at and , respectively.

Household Surveys
Face-to-face interviews in the form of key informant interviews, focus-group discussions guided by a checklist, and semi-structured household questionnaires were used in this study. The questionnaires were comprised of both open-and closed-ended questions to gather information about the perceptions of the local communities on LULC changes, and the drivers of these changes in Dedza during the studied period (1991 to 2015). A questionnaire was preferred for this study as it provides insight into the drivers of LULC changes [66]. The study employed a random sampling method to select respondents for the household interviews. The structured questionnaire was first pretested in 20 households in the Traditional Authority (TA) of Kaphuka (but not included in the sampled households for this study); then, modifications were made before the actual interviews of the sampled households. The questionnaire was administered to 586 households from 23 October 2017 to 10 November 2017 from 4 TAs, namely, Senior Chief Kachindamoto, Inkosi Kaphuka, Senior Chief Kachere, and TA Kasumbu. Additionally, the questionnaire was administered to respondents who (i) were aged 20 years and above, (ii) had lived in the respective area for at least 10 years, and (iii) were implicit decision-makers in the household, and/or, in the absence of a family head, it was made with appropriate representative and knowledgeable member of the household. The questionnaire had 7 sections covering the socioeconomic characteristics of the household, perceptions of local communities on LULC changes, and their causes (Appendix A). Each household interview lasted between 30 and 60 minutes.

Focus-Group Discussions and Key Informant Interviews
Focus-group discussions (FGDs) and key informant interviews were carried out to triangulate the obtained information from the household interviews and gain an in-depth and detailed understanding of local people's perceptions on LULC changes that had taken place in the studied landscape, and the associated underlying causes perceived to have contributed to the changes. A total of 4 FGDs were carried out in 4 TAs targeting the Area Development Committees (ADCs) where household interviews were conducted in the same period. FGDs facilitated by the researcher were carried out according to the procedure proposed by Reference [67], and were guided by a checklist of questions related to LULC changes and their driving forces. Each FGD consisted of 10-15 people and lasted between 120 and 180 minutes. A purposive sampling method was used to identify key informants based on their knowledge on the study area. In this study, key informants were exclusively technical members from the Dedza district council that were familiar with the issues in the study area. These technical members included the district commissioner, and researchers and officers from agriculture, natural-resource, and environmental institutions and organizations.

Other Datasets
Other data used in this study were climate (temperature and rainfall) data from 1991 to 2015, which were obtained from the Malawi Department of Climate Change and Meteorological Services (DCCMS) under the Ministry of Natural Resources, Energy, and Mining. Population data were obtained from the National Statistical Office of Malawi (NSO). Population estimations before 1991 and after 2008 were calculated by extrapolating the closest census data and annual growth rates using the formula adopted by Reference [14]: where and are total populations at Times 1 and 2, respectively; = exponential population constant; = number of years between two census enumerations; and = annual population growth rate.

Statistical Analysis
The study used a combination of data-analytical approaches and techniques including GISbased processing, descriptive statistics, and regression analysis. LULC change analyses were done using ArcGIS, QGIS, and ERDAS Imagine software. The socioeconomic data derived from the questionnaire were entered, processed, coded, and analyzed using Statistical Package for Social Sciences (SPSS) version 20 and subsequently subjected to further analysis. Descriptive-statistics analysis was used to describe socioeconomic variables of the households and summarize their responses and ranking of drivers of LULC changes. Ranking the drivers of LULC changes perceived by respondents (household surveys) was computed with the principle of weighted average using the ranking index adopted by References [68] and [69]: where = value given for the least-ranked level (for example, if the least rank is the 10th, then = 10, = 9, = 1; = counts of the least ranked level (in the above example, the count of the 10th rank = , and the count of the 1st rank = ). Data collected through FGDs and key informant interviews were qualitatively analyzed [70]. A nonparametric test (Pearson's chi-square) was used to ascertain the differences/associations between socioeconomic variables and respondent perceptions on drivers of LULC changes. Logisticregression analysis was performed to identify the key drivers of LULC changes in Dedza at the household level (Equation 6). By determining the drivers of LULC changes at the household level, the dependent variable was local people's perception of drivers for LULC changes and/or the perceived drivers identified, while independent variables included socioeconomic characteristics, such as age, gender, family size, education, and land-holding size. Logistic analysis at the household level estimated the probability of the effects of the independent variables on the dependent variables [71]: where = dependent variable indicating the likelihood that = 1, = the intercept, … n = coefficients of associated independent variables, and … = independent variables.

Accuracy Assessment
Accuracy assessment based on error (confusion matrices) showed an overall accuracy of 91.86%, with a kappa coefficient of 0.866 (Table 3). There were slight differences in user and producer accuracies of individual classes but the results of the datasets showed higher overall accuracy. These results provided a fundamental platform for subsequent analysis of LULC changes. Overall accuracy = 91.86%, Kkappa coefficient = 0.866. Figure 2 shows the spatial representation of LULC types from 1991 to 2015. The proportionate coverage area of each of the six classes extracted in Dedza from 1991 to 2015 of LULC change trends are summarized in Table 4 and Figure 3. At the beginning of the study period (1991), agricultural land was the most dominant LULC, covering 71.3% of the total studied area, followed by barren land (24.53%), forest (2.64%), wetlands (0.96%), water (0.37%), and built-up areas (0.2%) ( Table 4). The trend continued up to 2015 except for built-up areas. During the studied period (1991-2015), built-up areas substantially expanded almost tenfold (i.e., 950%) and barren land slightly increased, from 24.53% to 25.85%. Conversely, agriculture land, forest, wetlands, and water bodies drastically decreased in the same period ( Figure 5). The highest net loss was in agricultural land, followed by forest land. Despite these transformations, changes did not occur at equal rates. Results revealed that the area occupied by water bodies decreased by 34.8%, wetlands by 26.1%, forests by 37.2%, and agricultural land by 2.6% between 1991 and 2015. Built-up areas and barren land increased at an annual rate of 9.8% and 0.22% yr −1 . On the other hand, forests experienced strong loss at an annual rate of 1.94% yr −1; followed by agricultural land, wetlands, and water declining at a corresponding rate of change of 0.11%, 1.26%, and 1.78% yr −1 , respectively.    Table 5 shows the cross-tabulation change matrix for the changed areas and their corresponding percentages from one LULC class to another in comparison with the total area of each LULC class from 1991 to 2015. Despite the fact that all LULC classes have undergone changes in the study area, the degree of these changes was inherently different. Conversions occurred across the whole study area. During the study period, 96.03% of agricultural land remained unchanged, followed by barren land (93.72%), built-up areas (86.20%), water bodies (64.39%), wetlands (50.8%), and forest (30.23%). This clearly indicates that forest experienced the highest conversion with almost 70% of its total area converted to barren land (61.48%) and the rest to other LULC classes. The majority of agricultural land was converted to built-up areas (7244.91 ha) and barren land (2,960.01 ha), while the majority of barren land was converted to forest (2,803.86 ha) and agricultural land (2,162.61 ha). Even though built-up areas did not change much, almost 7244 ha were gained from agricultural land (7244.91 ha). Note: Bold numbers on the diagonal represent unchanged LULC proportions from 1991 to 2015 and their corresponding percentages, while others are the areas changed from one class to another.

Socioeconomic and Demographic Characteristics of Sampled Households
The socioeconomic and demographic attributes of the sampled households are presented in Table 6. The results revealed that the age of the respondents ranged from 20 to 97 years, with an average of 39.2 years. About 93.3% of the interviewees lived in the study area throughout the studied period. The majority (78.7%) of the respondents were married, about 63.3% of the sampled households were female, and 71.7% of the households were male-headed. The results also indicated that household size ranged from one person to 13 people, with an average of 5.6 persons. It is also worth noting that a larger proportion (96.1%) of the interviewees owned land, with 5.9% being landless. The farm size of the respondents varied from 0.25 to 13 acres, with an average of 2.32 acres. With respect to their education status, 77.8% of the respondents were literate (64.3% and 13.5% attended primary and secondary school, respectively), and 22.2% had never attended school. Approximately 82% of the sampled households were engaged in farming activities, and a small portion of the respondents (18%) were involved in on-farm activities, such as businesses, professional work, and craft work. The mean household income of the respondents was USD721.30 (MK 286,843.26) per year. Farming was ranked as the most important source of income in Dedza. Income from self-employment opportunities, such as businesses, handcraft, and trade, were ranked second, followed by piece works or occasional jobs, Village Loan Savings (VLS), full-time private/government employment, sale of forest produce, and renting out land.

Local-Community Perceptions on Observed Trends of LULC Changes and Proximity to Infrastructure
Significant differences were found among the interviewed households in perceptions regarding LULC changes and distance to different infrastructures such as main roads, health centers, schools, and towns (p < 0.001). Respondents perceived that agricultural land and forest cover significantly declined (p < 0.001) in the studied landscape. Results shows that 57.3% and 87.4% of local communities correctly perceived that agricultural land and forest, respectively, had declined ( Figure  4). Almost half of the respondents (53.4%) perceived that distance from water bodies remained the same over the studied period. Conversely, distance to infrastructures such as main roads, health centers, bus stops, and towns remained unchanged except for distance to markets and schools, which significantly declined (p < 0.001). Key informants from different institutions and FGDs also correctly perceived that agricultural land and forest cover drastically declined from 1991 to 2015.

Ranked Drivers of LULC Changes
The respondents identified 24 factors (12 proximate drivers and 12 underlying drivers) as important drivers contributing to LULC changes in Dedza, especially during the period under review (Tables 7 and 8). Fuelwood collection, charcoal production, timber, construction, and agriculture expansion were the top five ranked proximate drivers of LULC changes in the study area, with fire collection and charcoal production ranked first and second, respectively (Table 7). Similar results were also revealed during key informant interviews and FGDs in which firewood collection, charcoal production, settlements, and agricultural expansion were identified as the main causes of LULC in the study area. With respect to underlying causes of LULC drivers in the study area, the interviewed households identified population growth as the most important underlying driver contributing to LULC, followed by poverty, lack of financial resources, lack of law enforcement, and demand for timber (  Figure 5). Observed rainfall data between 1991 and 2015 were consistent with the local communities' perceptions, as indicated by declining unreliable rainfall ( Figure 6).

Household-Level Logistic Regression of Perceived Drivers of LULC Changes
Results revealed that education level negatively and significantly affected (p < 0.05) high perceptions of local communities on firewood collection, agricultural expansion, poverty, and population growth as LULC drivers in Dedza (Table 9). Charcoal production and settlements were not significantly influenced by age, gender, education level, land-holding size, and household size.

Land-Use and Land-Cover Change Dynamics
The post-classification comparison results for change-detection analysis and the change matrix from 1991 to 2015 revealed the extent of LULC changes occurring in different LULC classes throughout the study period. Dedza experienced substantial and increased rates of LULC changes between 1991 and 2015. Agricultural and barren land are the major LULC classes accounting for almost 96% of the total landscape in both 1991 and 2015. Most agricultural land, forest land, and water bodies from 1991 were intensively converted to built-up areas, barren land, and agricultural land, respectively. Recently, agricultural land in Dedza was developed for residential, commercial, and business purposes. The expansion rate of built-up areas on other LULC categories increased following the development of residential areas for commercial, academic, and business purposes. Barren land expanded at the expense of forest land and wetlands. The presence of major roads in the study area accelerated the expansion of built-up areas and exploitation of resources. Communities in the study area also correctly perceived that built-up areas and barren land had increased over the past years, with a decline in agricultural land, rivers, wetlands, and forest land. Additionally, as observed during field visits, demand for agricultural land and wetlands to be converted to residential land, and also land prices for these lands, had increased over the past years. Additionally, the use of older respondents (≥ 20 years) provided an accurate historical narrative of LULC changes in the study area, confirming the results of the observed LULC changes interpreted from remotely sensed data in the period of 1991-2015. Similar findings of other researchers showed that LULC changes occurred in related settings. For example, woodlands declined by 88.5%, while urban areas increased by 143% between 1984 and 2013 in the Likangala River catchment in Malawi [32]. Increased built-up areas and reduction in forest land and fresh water of the Upper Shire River Catchment of Malawi was also reported [28]. Contrary to the findings in this study, both authors found an increase in agricultural land in their study areas. It was reported that 20,747 hectares of forest land were lost between 1990 and 2008 in Malawi's Dzalanyama Forest Reserve, of which 64% of forest land was lost between 2000 and 2008 [29]. A recent study revealed that built-up areas increased by about tenfold at the expense of grasslands, shrub-bush land, and woodlands in the Central Rift Valley of Ethiopia between 1973 and 2014 [72]. Similar observations of the expansion in built-up areas, accompanied by a decline in forest land and agricultural land, were also made by other studies [15,29,40,43,69,73]

Drivers of LULC Changes
The research findings, based on the household surveys, FGDs, and key informant interviews, pointed to local communities perceiving firewood collection, charcoal production, agricultural expansion, settlements, and timber as the important proximate drivers of LULC changes in Dedza. These proximate drivers were triggered by high poverty levels, population growth, unreliable rainfall, lack of law enforcement by government, poor access to an alternative-energy supply, and high cost of agricultural input.
The majority of the local communities felt that population growth increased during the study period. Indeed, the population of Dedza has increased by 28% since 1998. This is also confirmed by the results of the population model used in this study, which simulated an increase in population for the studied period in Dedza from 1991 to 2015 [48]. Household surveys, FGDs, and key informants perceived that the rapid increase of the population in the study area was largely due to high fertility rates, early marriages, high birth rates, reduced mortality, polygamy, immigration, and illiteracy. Dedza shares its border with Mozambique, and during the war, economic instability, and the drought crisis, people from Mozambique would migrate to Dedza to survive. Earlier studies in Malawi also found population pressure as one of the drivers of LULC changes [32,74]. In other parts of the world, population growth was also reported as the main driver of LULC changes [14,[16][17][18]75] Firewood collection and charcoal production are the top two important proximate drivers of LULC changes in Dedza between 1991 and 2015. This is also directly associated with the use of threestone open-fire stoves by 88.2% of the interviewees, while the rest use charcoal stoves for cooking. This kind of domestic cooking stove enables households to use more firewood, thereby exacerbating deforestation and forest degradation. The use of three-stone open-fire stoves results in indoor-air pollution, which severely impacts human health, particularly the vulnerable populace, such as children and women. These results are also directly connected with the wide use of biomass as the main source of energy for the majority of the Malawi population. The use of charcoal and fuelwood for energy in the district is triggered by high poverty levels and low coverage of electricity and alternative sources of energy. Approximately 90% of Malawi's population relies on charcoal and firewood for energy [48,76]. This explains the forest-cover loss in the study area between 1991 and 2015. Proximity of Dedza to Lilongwe, the capital of Malawi, offers a market for forest products, and this exacerbates the collection of illegal firewood and the charcoal produced for harvested poles and timber for construction from government forest reserves in Dedza. The persistence of electricity blackouts (load shedding 8 to 24 hours) in Malawi (evidenced in Appendix B) also encourages the overdependence of local communities and urban dwellers on charcoal and firewood in order to meet increased demand in urban and rural areas. The inefficient production and unsustainable use of biomass energy sources in Malawi adversely contributes to environmental degradation, such as high deforestation, desertification, and soil erosion.
Among the perceived important drivers indirectly contributing to LULC changes in Dedza is poverty. Local communities are unable to buy agricultural inputs due to high poverty levels, high cost of agricultural inputs, and lack of financial resources. The majority of the local communities in the district are characterized by high levels of poverty and lack of alternative livelihood sources. Harvesting and selling of forest produce and products such as poles, timber, firewood, and charcoal are among the sources of income for most of the communities in the study area. Local communities living in Dedza and the surrounding districts are also forced to clear forests for additional cultivated land or to sustain their livelihoods as an immediate and quick source of income. As perceived by key informants and through focus-group discussions, Dedza rainfall has been very variable. The rural communities in Dedza depend on the sales of forest produce as a common survival strategy in the case of land degradation, decline or failure of crop production, soil infertility, frequent and prolonged droughts, and unreliable rainfall. Overdependence and unsustainable extraction of natural resources without alternative economic strategies, such as forests, land, and water, results in serious environmental problems including soil erosion, biodiversity loss and disintegration, natural-resource depletion, water and air pollution, deforestation, and forest degradation. The results of this study resonate with other similar studies in Africa where high poverty levels were reported as the contributory factors for LULC changes [14,15,77,78]. This study has further revealed that, among main socioeconomic determinants, the education level of rural communities significantly affected their perceptions toward LULC drivers in the study area.

Conclusions
The study has examined LULC changes using multitemporal remotely sensed images in conjunction with household surveys, FGDs, and key informant interviews to establish their drivers in Dedza during the period of 1991-2015. There was a substantial decline in forest land, agricultural land, wetlands, and water, while built-up areas and barren land drastically increased over the studied period. Firewood collection, charcoal production, population growth, and poverty were ranked as the important drivers perceived by local communities to be responsible for LULC dynamics in the studied area. The findings also depict that education level significantly affected interviewees' perceptions toward some of the drivers of LULC changes. The drivers identified in this study can be used as a tool for land-use planning, as well as input for modelling future LULC changes for the development of effective land-management strategies, guidelines, and policies for informed decisionmaking in Dedza and other districts with similar settings in Malawi. Appropriately tenable strategies and policies are urgently needed in the study area to address or avert undesirable LULC changes taking place in Dedza. Based on these results, the study recommends further studies to investigate the impact and consequences of these LULC changes on the rural livelihoods of the studied area so that landscape-management decisions and strategies are made based on scientific findings.  Acknowledgments: This research study was supported by the Organization for Women in Science for the Developing World (OWSD), University of Pretoria, the National Research Foundation (NRF) of South Africa, and the Schlumberger Foundation through its Faculty for the Future Program. We would also like to thank the Malawi Department of Forestry, the Malawi Department of Climate Change and Meteorological Services, and the NSO for the ancillary data used in this study. Vote of thanks to technical members of Dedza district council, chiefs and local communities who participated in this household survey and enumerators who tirelessly assisted with data collection.

Conflicts of Interest:
The authors declare no conflict of interest.   11. Indicate the number of farms you have, and their size, purpose, and distance from home.

Farm land Size (acres)
Purpose/Use (consumption, sale, or both) Distance from home