Defining Deforestation Patterns Using Satellite Images from 2000 and 2017 : Assessment of Forest Management in Miombo Forests — A Case Study of Huambo Province in Angola

A few studies have recently been published on changes in land use/land cover (LU/LC) of Angolan Miombo forests, however, none have attempted to offer forest management solutions for degraded Miombo forests. Landscapes are witness to past and present natural and social processes influencing the environment, where each period in the past leaves footprints on the landscape’s development, which can be described by a continual decrease in forest area over time. The expansion of degraded areas from 2000 to 20017 began near urban areas where many Miombo forests have been eliminated or highly degraded, particularly in the southwest and northeast of the Huambo province. Large areas of degraded forests were observed along the Benguela railway (Caminho de ferro de Benguela). Our detailed analysis of the landcover map suggests that the impact has been devastating and there is no form of forest protection, which leads to unregulated exploitation. Descriptions of the Miombo forest dynamics are explained using height–diameter curves developed for different vegetation types that provide important insights about forest structures in the management zones. The height–diameter models differed for all vegetation types, and four management zones (MZ) were created based on a set of particular attributes. The vegetation types differed in each management zone, which included agricultural land and bare soil (MZ–E), grassland or savanna (MZ–C), open Miombo forests (MZ–B), and closed Miombo forests (Miombo forests). The four management zones were easily identified on the available maps and the height–diameter models developed represent a fundamental tool for future studies on forest planning.


Introduction
Land use (LU) and land cover (LC) arise from interactions of human activities and global climate change, which can also interact to produce degraded and deforested areas.Historical records tell us that deforestation is not a recent phenomenon; the Mediterranean region was already deforested before the fifth century [1].Identifying and assessing the condition of forests is not easy since people have widely different views of what constitutes degradation.Deforestation occurs when forests are converted to non-forest uses, such as agriculture and road construction.Forest degradation occurs when forest ecosystems lose their capacity to provide important goods and services to people and nature [2].For some people, any forest management activity may cause degradation.For others, a forest is only degraded when it can no longer deliver needed goods and services.There is no globally agreed definition of forest degradation which makes the discussion more complex [3].Deforestation is considered the result of the clearing of land for agriculture [4], both for the large-scale production of global commodities [5,6] and, in Africa, for small-scale production of food and cash crops [7].Deforestation is among the most commonly studied phenomena in the frontier literature, as it is often associated with negative impacts on the global climate and biodiversity [6,8].The scale of the environmental legacy of deforestation is dependent on both the magnitude and timing of historic land cover changes, not merely the snapshot of forest cover that is directly observable today [9][10][11][12][13][14].The most important anthropogenic alterations of natural environments have always been the clearing of forests to establish cropland and pasture, and the exploitation of forests for fuelwood and construction materials.Analysing landscape changes due to past and present natural and social processes at different scales over time constitutes the departure point for landscape management because of close links between forests and physical attributes of the landscape.The mix of LC and LU (landscape composition) usually includes agricultural lands and native vegetation, and human dwellings in villages and urban areas.The spatial arrangement of different land uses and cover types (landscape structure), and the norms and modalities of land governance contribute to the character of a landscape and its management [15,16].
Several published studies have focused on LC and LU changes in tropical forests because of the implications for global carbon cycling, and soil and biodiversity loss [17][18][19].However, little is known about the dynamic of Miombo ecosystems.There is a need for a better understanding of the long-term environmental and socio-economic effects of changes in LC and LU in Miombo ecosystems [20].The quality and quantity of resources, the rates by which these resources change, and the overall distribution of the land cover types should be known, in order to develop management plans and ensure sustainable use of nature resources while preserving biological richness and diversity [21].
In Angola, the over-extraction of wood resources, linked with clearing for agricultural purposes, unregulated burning and, sometimes, overgrazing, creates disorder that influences the health of Miombo forests.Although Angola has approximately 40-60 million hectares of forest (of which 45% is Miombo) largely administered by the government [22][23][24], it lacks the basic data to better understand the forest dynamics.Among sub-Saharan African countries, Angola has one of the highest rates of deforestation and the lowest rate of forest expansion [25].During the military conflict from 1975 to 2002, many people migrated from rural areas to urban centers and often relied on forests as a place to seek sustenance, which contributed to deforestation rates through wood extraction for firewood and charcoal, slash-and-burn cultivation, urban expansion, and logging.
Although illegal logging of valuable timber is often cited as a primary cause of forest degradation in Angola, there is no information to objectively evaluate the problem.Official records estimate that annual deforestation rates were 0.20% between 1990 and 2000, and increased to 2.1% from 1990 to 2005; during the period 1990 to 2010, approximately 5% of Angolan forests were already degraded [26,27].After 27 years of civil war, Angola is now in the process of recovery, and natural resources and resource management play an important role in the rehabilitation and reconstruction process of the country.However, the lack of basic data (e.g., forest types, tree species composition, growth, and height-diameter models) in Angolan Miombo forests is a critical problem for policymakers and managers because they cannot effectively determine rates of allowable harvest or assess the extent of deforestation or degradation.
Regular monitoring and assessment of LU/LC changes are critical for understanding forest dynamics and the extent and impact of anthropogenic and natural changes at regional or global scales.Remotely sensed data have been widely used to classify LC and provide estimates of deforestation patterns in corresponding areas.However, field data collection is also important to formulate models for the most important parameters of interest for forest managers.For example, tree heights and diameters are the most important measures of growth that are often used to evaluate site index, calculate tree volumes, and predict future stand growth [28,29].Tree height and diameter have also been used to estimate merchantable volumes and they are commonly used for planning silvicultural alternatives at the stand level, and for effective forest design and monitoring [30][31][32][33][34][35].Understanding stand structure and the height-diameter relationship of trees provides useful information to establish appropriate countermeasures for sustainable management of endangered forests.Height-diameter relationships have provided adequate estimates of forest site productivity and they are a central element of forest management and planning; site productivity is considered a good measure for even-aged, uneven-aged, and mixed species stands [36,37], such as the Miombo forests of Angola.The estimates provide critical information to allow managers to forecast rates of change and evaluate potential production of wood or biomass under different management prescriptions.Tree growth is not only determined by the physiological and ecological processes of the tree, but also by the site quality.Under normal conditions, there is a strong biological height-diameter relationship which reflects the ecological conditions of the forest and can be a major indicator for evaluating the quality of degraded forest ecosystems [38].The height (H)-diameter (D) relationship, as one of the most important stature characteristics [2], has been recognized as one manifestation of the many ways in which trees adapt to changes in environment, and this process of structural change may last decades or centuries [39].Sustainable forest management plans or forest landscape management plans should focus not only on the spatial distributions of forest stands, but also vertical structures, which could be only possible with the help of several forest models, including height-diameter models.Developing forest stand specific height-diameter models is considered one of the most important tasks in forest planning, inventory design, and monitoring of stand dynamics or stand characteristics, which also determine types and quality of forest landscape management regimes.Forecasting forest change is essential to forest and landscape management [40].Growing interest in issues other than timber production demands a greater use of forest models, including height-diameter and stem volume models.
Several studies have reported the spatial dynamics of deforestation and forest species diversity of the Central Plateau region of Angola [41][42][43][44]; however, no studies have provided adequate tools for the informed management of Angolan Miombo forests, such as height-diameter models, or reported on deforestation patterns that threaten these forests.This paper addresses both of these critical gaps using field data to characterize forest types to create effective management zones.We used remotely sensed data to evaluate deforestation patterns in Miombo forests in the Huambo province of Angola.The topic addressed is of particular importance and warrants consideration mostly due to scarcity of research in this field in Angola.Miombo woodlands are complex landscapes, and the land cover types comprising the Miombo forests have diverse values.They comprise intact and disturbed woodlands, which are patchily distributed amongst fields, wetlands, homesteads, kopjes and termitaria.Few studies of ecosystem values have engaged with this ecological complexity [45].Ecological studies concerned with horizontal components of Miombo forest structure, such as density and basal area, have shown large-scale variations along broad environmental and edaphic gradients.However, variations in the vertical components of the forest structure remain less studied.For example, studies have shown that H declines more sharply with elevation than the diameter.This is despite available evidence suggesting that tree height for a given diameter vary significantly among species and regions.Such variations could hold important implications for carbon storage potential of tropical forests, for instance.This is because tropical tree above-ground biomass and carbon fluxes are usually estimated by applying allometric equations to diameter measurements only.
To examine the importance of different land cover types, in supporting a subset of Miombo woodland values, we identified deforestation patterns of the central plateau of Angola to derive LU strategies that enhance forest management through conservation areas and avoid trade-offs for other ecosystem services.We also modelled some basic forest attributes of the vegetation types to provide suggestions for management zones for future planning efforts within the study area.This link is based on the hypothesis that landscape change is related to variation of the forest structure if H-D relationships differ by vegetation zones and forest types (open forest and closed forest, including zones where forest and savannah overlap).

Study Area
Angola is located on the southwestern coast of Africa, and its climate ranges from tropical humid in the north to extremely arid in the southwest.According to the World Wildlife Fund, Angola has 15 ecoregions, and Miombo forests are the most widespread biogeographic unit [23,46].
Our study focused on Huambo province within the Central Plateau (12 • 29 59.99 S, 15 • 39 59.99 E) of Angola (Figure 1).The dominant soils are ferralsols at higher elevations and alluvial fluvisols at lower elevations.With an area of approximately 34,270 km 2 and a population of approximately 1.9 million inhabitants, Huambo province is divided into 11 municipalities.It is 1700 m a.s.l., with a mean annual temperature of 20 • C, and annual rainfall ranges from 1200 to 1600 mm [47][48][49].The climate is humid meso-thermal with dry winters and warm summers.Most of the precipitation occurs during the warmer summer months (October to April), with heavy precipitation in December and March.Six major rivers originate in Huambo, many of which drain into the Atlantic Ocean.Vegetation in Huambo province consists largely of four types: highland forests (Afromontane forests), Miombo woodlands, swamps, and dry grasslands [50][51][52].The landscape is mainly comprised of Miombo forests, which is subdivided into closed and open forests, and, in lower areas, savanna woodlands intercepted by grasslands are more prominent.These vegetation types are difficult to map because they are scattered and patchy.The purpose of mapping the vegetation is to show their spread and distribution across the whole province.Frequent intense fires, tree cutting for charcoal, and clearing for new agricultural fields affect vegetation distributions and structure.

Study Area
Angola is located on the southwestern coast of Africa, and its climate ranges from tropical humid in the north to extremely arid in the southwest.According to the World Wildlife Fund, Angola has 15 ecoregions, and Miombo forests are the most widespread biogeographic unit [23,46].
Our study focused on Huambo province within the Central Plateau (12° 29′ 59.99″ S, 15° 39′ 59.99″ E) of Angola (Figure 1).The dominant soils are ferralsols at higher elevations and alluvial fluvisols at lower elevations.With an area of approximately 34,270 km 2 and a population of approximately 1.9 million inhabitants, Huambo province is divided into 11 municipalities.It is 1700 m a.s.l., with a mean annual temperature of 20 °C, and annual rainfall ranges from 1200 to 1600 mm [47][48][49].The climate is humid meso-thermal with dry winters and warm summers.Most of the precipitation occurs during the warmer summer months (October to April), with heavy precipitation in December and March.Six major rivers originate in Huambo, many of which drain into the Atlantic Ocean.Vegetation in Huambo province consists largely of four types: highland forests (Afromontane forests), Miombo woodlands, swamps, and dry grasslands [50][51][52].The landscape is mainly comprised of Miombo forests, which is subdivided into closed and open forests, and, in lower areas, savanna woodlands intercepted by grasslands are more prominent.These vegetation types are difficult to map because they are scattered and patchy.The purpose of mapping the vegetation is to show their spread and distribution across the whole province.Frequent intense fires, tree cutting for charcoal, and clearing for new agricultural fields affect vegetation distributions and structure.

Ground Data Collection and Processing
We established 18 circular sample plots (radius: 25.2 m) in open and closed Miombo forests by distributing them in a random fashion to represent the major vegetation types in Huambo Province.Prior to sample plot establishment, we used satellite images of land cover to delineate the potential forest area for sampling.When possible, study sites were identified from the images, and we conducted field surveys at the sites to ensure the areas were suitable for sample plot establishment.The locations of the plots are presented in [53].For each sampled forest patch, the total area of all the plots should comprise at least 2% of the total forest patch area; for example, in a 20-ha forest patch, a total area of 0.4 ha would be sampled.For each sampled tree within a plot, diameter at breast height (DBH) was measured using diameter tape and heights were measured using a hypsometer Vertex IV (Haglof Inc., Långsele, Sweden).Sample plot data allowed us to calculate density (stems•ha −1 ), basal area (BA; m 2 •ha −1 ), and stem volume (m 3 •ha −1 ) for grassland, open and closed Miombo.Analysis of sample plot data was conducted using SAS V.9.4 and MATLAB R2017b professional edition (MathWorks©, Inc., Natick, MA, USA).We determined the parameter values of diameter growth functions using regression analysis [54].

Modelling Height-Diameter Relationships
Data for many tree species were rare in our sampling efforts, thus it was not possible to develop a species-specific height-diameter models for several species.Consequently, we pooled the data from all species to develop a height-diameter model applicable to all species.We used a versatile nonlinear growth function, the Näslund function [55], to fit the data.This function has frequently been used to model the height-diameter relationship for several tree species because of its flexibility [56][57][58][59].The Näslund function is: where H ij and DBH ij are the total height and DBH measurements, respectively, for the jth tree on the ith sample plot, b 1 and b 2 are parameters to be estimated, and ε ij is an error term.A value of 1.3 m, the height at which DBH was measured, was added to avoid the predictions of zero height when DBH approaches zero.
For modelling categorical variables [29,60], we included two dummy variables (V 1 , V 2 ) to describe the effects of three vegetation types on height-diameter relationship (Table 1).

Vegetation Type
The effect was best described when b 1 of the base function (Equation ( 1)) was expressed as a linear function of dummy variables.The base function was extended through redefining its parameter b 1 as a function of tallest tree per sample plot (H max , m) and stand basal area (BA, m 2 ha −1 ), as they significantly contributed to the height-diameter relationship in this setting.Other stand variables, such as stem number per hectare and quadratic mean diameter per sample plot were also subsequently added to the model.However, no significant improvement in the fit statistics was observed.We used an expanded form of the above function with two dummy variables and two stand-level variables: where H maxi is the height of the tallest tree on the ith sample plot; α 1 , α 2 , . . ., α 5 , and b 2 are parameters to be estimated, and V 1 and V 2 are dummy variables that account for the effects of vegetation types.The vectors of errors ε ij and sample plot-level random effects (u i1 , u i2 ) are defined by ε i ∼ N(0, R) and u i ∼ N(0, D), respectively, meaning that vectors ε i and u i were assumed to have normal distributions with zero mean and within-sample plot variance-covariance matrices R and D, respectively [61,62].
The height-diameter model (Equation ( 2)) was estimated using the mixed modelling method.

Modelling Stem Volumes of Miombo Forest Species
We selected an exponential function to most effectively describe the stem volume variations.We modelled stem volume using DBH and total height (HEIGHT) as the main predictors and vegetation types, stand basal area, and DBH-to-quadratic mean diameter ratio (dq) were used as covariate predictors: where SV ij is the stem volume of the jth tree on the ith sample plot; α 1 , α 2 , . . ., α 4 , and b 2 are parameters to be estimated, V 1 and V 2 are dummy variables that account for the effects of vegetation types on the stem volume, and ε ij is an error term that is assumed to be normally distributed with mean zero and constant variance.The volume model was estimated using the nonlinear ordinary least squares regression method.

Remotely Sensed Data and Application
To derive deforestation patterns in Miombo forests, we used Landsat imagery from two different sensors, the Landsat-8 Operational Land Imager (OLI; 2013-2016) and the Landsat-7 Enhanced Thematic Mapper Plus (ETM+; 1999-2016), for the study period from 2000 (August) to 2017 (August).All products were obtained through the United States Geological Survey (USGS) Earth Explorer (EE) tool.

Pre-Processing of Remotely Sensed Data
Due to the size of the study area, we used five frames from Landsat-7 and four frames from Landsat-8, combined the bands (Table 2) for each data frame, and eliminated the null values using the Copy Raster toolset in ArcGIS desktop V.10.5.1 (ESRI Inc., Redlands, CA, USA).Finally, we mosaiced the frames for each year separately using the Mosaic to New Raster tool.Source: [63].
To improve the image quality and eliminate sources of errors, we used the Fast Line-of-sight Atmospheric Analysis of Hypercubes (FLAASH) method, a radiometric correction method, in the Environment for Visualizing Images V.5.1 software (ENVI; Harris Geospatial Solutions, Boulder, CO, USA).We confirmed the results by examining the histogram values produced from the FLAASH procedure.To georeference the images, we used the defined nominal Landsat satellite track (path and row) approach, where the nearest World Reference System (WRS) path and row to the line of sight was the scene's center; the whole process was conducted using the USGS website.In addition, we resampled the image from 2000 from a 15 m spatial resolution to 30 m using the nearest neighbor method.We did this to equalize the cell sizes in either the images of 2000 or 2017.The study of land use classification in Angola is very limited, unfortunately, due to lack of quality satellite information, and for this study area we took maximum advantage of existing spectral sources.Although using different Landsat imagery can cause errors, there are some pre-processing steps that can help users to reduce the errors as much as possible.In this case, we used the value histogram of the same objects in different images to make sure the differences between the digital number DN values were because of land cover changes during the time, not the nature of the image.For example, by comparing dark and artificial objects that did not change during the time (or for which changes were not significant), we were able to verify the accuracy of our methodology.

Data Processing
For the classification process, we separated the work into three stages.First, we decided to base the land use classification on four different classes (vegetation areas, urban areas, water bodies, and others).We then identified certain areas of pixels (training areas) for each class using polygons, and we extracted a signature file from the training samples to supervise the classification process (Table 3).A maximum likelihood classifier was used to classify the study area into the desired classes.Finally, the percentage of area in each class was calculated for 2000 and 2017.The spatial-temporal variability of normalized difference vegetation index NDVI was assessed to study deforestation using harmonic analysis.We first spatially normalized observations to reduce seasonality.Subsequently, we detected deforestation by assessing whether a newly acquired observation (satellite image) in the monitoring period is in an extreme change when compared against spatially normalized values in present time data defined over a reference period.The calculation of the NDVI for multi-date satellite images of Landsat (7,8) was used to perform change detection of the deforestation in Huambo Miombo.The differences between the percentages of each class (2000-2017) were used as a basis for further analyses.The NDVI, as one of the most successfully used vegetation spectral indices, allows comparison between inter-annual and seasonal changes in vegetation.The NDVI measures the amount of green vegetation in an area and is used to distinguish forested from deforested areas.Also, to investigate the impact of railway constriction on deforestation and land-cover change, we compared the land-cover type using a buffer zone with radius of 5 km around the "Caminho de ferro de Benguela" Benguela railway.Using the watershed areas and land cover of the study area, we devised a method to delineate management zones (MZ).First, topographic attributes were calculated from elevation data using TAPES (Terrain Analysis Programs for the Environmental Sciences), a useful program designed to handle the complexity of individual landscape processes and patterns encountered in delineating the appropriate spatial and temporal scales [64].
We then removed imperfections or sinks of the digital elevation model (DEM) with the fill tool in the ArcGIS spatial analyst toolbox.Fill layer was used to calculate the hydrologic flow direction layer using the direction toolset in the spatial analyst toolbox, hydrology.The hydrologic flow direction raster was analysed as an input of the basin layer to find all sets of connected pixels that belong to the same drainage basin.Basin tool was used to create a raster delineating all drainage basins, hereafter referred to as watersheds.A watershed is the upslope area that contributes flow to a common outlet as concentrated drainage.After transforming the results of the basin tool to a vector polygon based on visual assessment of land cover type, we specified a MZ for each watershed area.Finally, for simplicity in calculating the percentage of each MZ, we merged similar MZs together.

Height-Diameter Relationships
The mixed-effects model described a large proportion of the total variations in the heightdiameter relationships without significant residual trends; all parameter estimates were significant (p < 0.05) (Table 4).Including the sample plot-level random effects in the model significantly increased the fit with reference to the model fitted without random effects; the R 2 value increased by 14% in comparison to the ordinary least squares model.Sample plot-specific height-diameter curves simulated using only the fixed portion of the mixed model (mean model) and mixed-effects model (both fixed and random portions included) seemed different.The latter curves passed through the middle of the data for each plot, which indicated that calibration using the height measurement of at least one randomly selected tree worked well for the sample plot-specific model.2).AIC: Akaike's information criterion; BIC: Schwarz's information criterion; α1, …, α5, b2: fixed parameters; ui1, ui2 = random effect parameters; σ 2 ui1: variance of ui1; σ 2 ui2: variance of ui2; σ 2 : residual variance.

Height-Diameter Relationships
The mixed-effects model described a large proportion of the total variations in the height-diameter relationships without significant residual trends; all parameter estimates were significant (p < 0.05) (Table 4).Including the sample plot-level random effects in the model significantly increased the fit with reference to the model fitted without random effects; the R 2 value increased by 14% in comparison to the ordinary least squares model.Sample plot-specific height-diameter curves simulated using only the fixed portion of the mixed model (mean model) and mixed-effects model (both fixed and random portions included) seemed different.The latter curves passed through the middle of the data for each plot, which indicated that calibration using the height measurement of at least one randomly selected tree worked well for the sample plot-specific model.
The selected model (Näslund function) is based on the growth theory (faster increase of height in earlier stages and slower increase in later stages) [56,65].Näslund function fitted best to our data due to its flexibility in fitting height-diameter relationships for a number of tree species (Figure 3) [58].The selected model (Näslund function) is based on the growth theory (faster increase of height in earlier stages and slower increase in later stages) [56,65].Näslund function fitted best to our data due to its flexibility in fitting height-diameter relationships for a number of tree species (Figure 3) [58].

Volume Assessment
The exponential function described a large proportion of the stem volume variations [R 2 = 0.846; RMSE = 0.0276 (m 3 ha −1 ); AIC = 4106] without significant residual trends, and all parameter estimates were highly significant (p < 0.0001; Table 5).

Volume Assessment
The exponential function described a large proportion of the stem volume variations [R 2 = 0.846; RMSE = 0.0276 (m 3 ha −1 ); AIC = 4106] without significant residual trends, and all parameter estimates were highly significant (p < 0.0001; Table 5).The vegetation type-specific model behaved significantly differently for different vegetation types, which suggested that vegetation types could also affect stem volumes (Figure 4).The vegetation type-specific model behaved significantly differently for different vegetation types, which suggested that vegetation types could also affect stem volumes (Figure 4).5), mean values of covariate predictors, and height at five-meter intervals.

Deforestation Patterns
The loss of natural vegetation and plantation is the most visible evidence of deforestation pattern in Huambo.The main argument for this assertion is the ongoing vegetation loss (17% from 2000 to 2017) and agricultural expansion (confirmed by other studies) which is also verified by the remote sensing analysis (Tables 6 and 7).This pattern was accentuated along the railway line (Caminho de ferro de Benguela) where natural forest was exploited to supply wood for the steam locomotives.Figure 5 of forested areas in Angola in 2000 and 2017 shows that many forests have suffered from deforestation.From our detailed analysis of the LC map, the impact is devastating and natural forests in Huambo appear to not be managed at all.  the curves were produced using parameter estimates (Table 5), mean values of covariate predictors, and height at five-meter intervals.

Deforestation Patterns
The loss of natural vegetation and plantation is the most visible evidence of deforestation pattern in Huambo.The main argument for this assertion is the ongoing vegetation loss (17% from 2000 to 2017) and agricultural expansion (confirmed by other studies) which is also verified by the remote sensing analysis (Tables 6 and 7).This pattern was accentuated along the railway line (Caminho de ferro de Benguela) where natural forest was exploited to supply wood for the steam locomotives.Figure 5 of forested areas in Angola in 2000 and 2017 shows that many forests have suffered from deforestation.From our detailed analysis of the LC map, the impact is devastating and natural forests in Huambo appear to not be managed at all.The change of forest cover is evident in Figure 6 where the most patterns are shown in conversion of natural forest to agriculture land and disappearance of forest plantation (Figure 6).This is a result of the logging and slash-and-burn agriculture of subsistence, where farmers typically plant their crops, deplete the soil, and move on after a couple of years.Figure 7 shows an area which remained unchanged between 2000 and 2017.The change of forest cover is evident in Figure 6 where the most patterns are shown in conversion of natural forest to agriculture land and disappearance of forest plantation (Figure 6).This is a result of the logging and slash-and-burn agriculture of subsistence, where farmers typically plant their crops, deplete the soil, and move on after a couple of years.Figure 7 shows an area which remained unchanged between 2000 and 2017.

Delineated Management Zones
Four management zones (MZ) were aligned with elevation attributes of the landscape, whereby MZ-A and MZ-B represent lower elevations in the north-western portion of the study area (Figure 8).MZ-E comprises upland areas (from 1484 m to 1897 m) in the southern portion of the study area; MZ-C can be interpreted as a combination of elevation areas from 1275 m to 1609 m above sea level

Delineated Management Zones
Four management zones (MZ) were aligned with elevation attributes of the landscape, whereby MZ-A and MZ-B represent lower elevations in the north-western portion of the study area (Figure 8).MZ-E comprises upland areas (from 1484 m to 1897 m) in the southern portion of the study area; MZ-C can be interpreted as a combination of elevation areas from 1275 m to 1609 m above sea level

Delineated Management Zones
Four management zones (MZ) were aligned with elevation attributes of the landscape, whereby MZ-A and MZ-B represent lower elevations in the north-western portion of the study area (Figure 8).MZ-E comprises upland areas (from 1484 m to 1897 m) in the southern portion of the study area; MZ-C can be interpreted as a combination of elevation areas from 1275 m to 1609 m above sea level (Figure 8).For a thorough description of the management zones, see [54].Vegetation types also differ in each management zone, whereby MZ-E represents agricultural land and bare soil, MZ-C consists of grasslands or savannas that include MZ-D, as described in [54], MZ-B contains open Miombo forests, and MZ-A is comprised of closed Miombo forests.It was not possible to represent the plantations by a single management zone because they were found distributed across other zones with similar elevation patterns.
(Figure 8).For a thorough description of the management zones, see [54].Vegetation types also differ in each management zone, whereby MZ-E represents agricultural land and bare soil, MZ-C consists of grasslands or savannas that include MZ-D, as described in [54], MZ-B contains open Miombo forests, and MZ-A is comprised of closed Miombo forests.It was not possible to represent the plantations by a single management zone because they were found distributed across other zones with similar elevation patterns.

Tree Height-Diameter Model and Description of the Management Zones
Tree height and DBH are important variables in forest inventories and timber management as these variables are necessary for estimating stand volume and biomass.Tree height data modelled from DBH can be employed as an input variable in various forest models, together with DBH, such as growth and yield models, site productivity models, crown models, biomass models, and carbon budget models [61].These models serve as important tools in forest management decision-making.Realistic description of stand growth and stem volume, and simulation of stand structure and dynamics, is possible only with use of the stand-specific height-diameter models [66].The heightdiameter models may be used to evaluate site productivity [36,45,67,68].Developing precise standspecific height-diameter models is considered one of the most important tasks in forest planning, inventory design, and monitoring of stand dynamics [56].Precise height-diameter models developed with subject-specific (e.g., sample plot-specific) variations are included in the model [57,59,66].The presented sample plot-specific mixed effects height-diameter model can be used for characterization

Tree Height-Diameter Model and Description of the Management Zones
Tree height and DBH are important variables in forest inventories and timber management as these variables are necessary for estimating stand volume and biomass.Tree height data modelled from DBH can be employed as an input variable in various forest models, together with DBH, such as growth and yield models, site productivity models, crown models, biomass models, and carbon budget models [61].These models serve as important tools in forest management decision-making.Realistic description of stand growth and stem volume, and simulation of stand structure and dynamics, is possible only with use of the stand-specific height-diameter models [66].The height-diameter models may be used to evaluate site productivity [36,45,67,68].Developing precise stand-specific height-diameter models is considered one of the most important tasks in forest planning, inventory design, and monitoring of stand dynamics [56].Precise height-diameter models developed with subject-specific (e.g., sample plot-specific) variations are included in the model [57,59,66].The presented sample plot-specific mixed effects height-diameter model can be used for characterization of the vertical stand structures, estimation of volume and biomass, and simulation of stand dynamics for each vegetation zone of the studied forest.
The major challenge in our study is that a large number of tree species existed with relatively few observations for the majority [69].The task of developing a species-specific height-diameter model for several species seems intractable, at least initially, and the mixed species approach used is advantageous when the species is treated as a random effect to circumvent the need for grouping the data [70].However, we did not consider this because of computational complexity, but species-specific effects on the height-diameter relationships were included using the dummy variable modelling approach.This approach is considered suitable for modelling height-diameter relationships using the data originated from multi-layered and multi-species stands [57,71,72].We chose the best functional form (Näslund function) to adequately describe the variations in the height-diameter relationship of individual trees for each vegetation type (Figure 3).Näslund function is based on the growth theory (faster increase of height in earlier stages and slower increase in the later stages) [56,61,65,73] and is best suited to our data due to its pronounced flexibility for fitting data acquired from multiple tree species [58,59].By looking at the growth curve patterns, we hoped to discover some similar patterns or principles that govern the growth of the species in the respective vegetation types.The height-diameter model is significantly different for each vegetation type, which indicates that the vegetation types significantly influence height-diameter relationships (Figure 3).The grassland zone was found inside MZ-E (1,484 m to 1897 m a.s.l.), which is the most heavily deforested area, followed by open Miombo forests of MZ-C and MZ-B (Figure 8).

Deforestation Patterns of Miombo Forests in Huambo Province
One of the objectives of the study was to identify deforestation patterns and assess important management criteria in degraded Miombo forest areas.By analysing spatial and temporal information from maps, satellite data and field observations, it was possible to map, quantify, analyze and explain the natural and human influences on the landscape over time.The landscape patterns represent the key to studying correlation between woodland vegetation and landforms.It provided the opportunity to focus on issues related to vegetation variability, forest zoning and forest fragmentation.
From a historical perspective, deforestation in Angola accelerated after the civil war and not during the colonial era.Many believe that with the arrival of the Portuguese in Angola, intensive exploitation of forest resources began due to the high timber demands for ships at that time.Large quantities of precious timber were exported to Portugal [74], and timber was needed to supply steam trains on the Benguela railway (Caminho de ferro de Benguela), both of which increased pressure on forests.Over time, natural forests were unable to recover and regrow, thus many turned to new introduced species, such as Eucalyptus spp), Thuja spp., Pinus spp.and Grevillea robusta, which were planted in very large areas along the railway [48,75].Within five kilometers of the railway, deforestation was particularly evident between 2000 and 2017 (Table 6).Local populations use the plantations for firewood.The eucalypts plantations alongside the Benguela Railway in Huambo are some of the largest in all of Africa [76]; they were originally planted to provide firewood for the steam locomotives.The plantations were quickly deforested, although natural regeneration can help forest growth and recovery (Figures 6 and 7).In some areas with poor access, some plantations were largely untouched and conserved (Figure 8).Most of the plantations belong to private enterprises (mainly the railway and pulp companies), while the remaining plantations are state property managed by the Ministry of Agriculture and Rural Development.
During the Angolan civil war, large portions of the population relocated to the cities for security reasons; however, with the Luena Peace Agreement in April 2002 [77] and the end of the civil war, millions of people were forced to leave the cities and return to their homes, and this initiated a period of rapid expansion of agricultural fields.In addition, trees were rapidly harvested by the returning populations who needed firewood for cooking and heating (Figure 5).
From 2000 to 2017, approximately 17% of forest vegetation was lost (Figures 5 and 6; Table 5).The decreased areas of vegetation and the increase of urban areas illustrate the primary changes in vegetation cover throughout the entire study area.Within a radius of 30 to 50 km around Huambo Province, many areas appear to be savanna-like because they contain only a few scattered trees.The expansion of degraded areas initiated near urban areas, and those patterns were accentuated in the south-western and north-eastern portions of the province (Figure 5); many of the areas lost were grasslands or savannas.
From a global perspective, managed forests increase over time, which is largely a result of forest management planning [78,79].In Angola, however, forested area has decreased over time, which may be explained by the lack of forest management plans and intensive charcoal and firewood production for cooking and heating, as suggested by several other studies [27,42,78].Schneibel et al suggested that the expansion of agricultural land was a primary contributor to deforestation rates, and they concluded that even during the conflict period and the post-ceasefire period (2002), the expansion of agriculture land was not reduced [42].
Evidence of deforested areas is well documented in Angola, however, no management actions have been described or assessed for these degraded areas.Romeiras et al documented the biogeographic patterns of timber species, and determined that Angolan timber species are under increasing risk [23]; they recommended conservation and sustainable management measures be developed to ensure future forests.Our results addressed this need through the creation of management zones.This creation of management zones is an important first step for future forest planning and highlighting management suggestions considering the current policies in the country.
In the absence of baseline information on which to base assessments of rates of deforestation, multi-temporal analysis of remote sensed data seems to be the most readily available option.
The forest management zoning system presented in our article is based on topographic attributes and is linked to the management zones proposed by [54].The zones can help in the sustainable management of native forests in meeting current and future community needs.The zones are easily distinguishable in the map, as they largely follow the topography of the area.
We recommend that each management zone should be examined to determine the appropriate operations and long-term strategies to ensure that these needs are met.For example, in MZ-E, a predominantly agricultural area, grassland and savanna vegetation consists of shorter trees and a mean DBH of 4.9 cm.In MZ-E, management objectives should seek to maximize the protection of cleared (non-forested) areas and forbid all activities that require complete vegetation removal (e.g., timber harvest, removal of forest products, grazing by domestic stock) to allow forests to regenerate naturally, or even artificially in certain areas.For MZ-C, the objective should be to conserve forest ecosystems while also allowing some production activities, however, with limited tree removal.For MZ-A and B, they could be reserved for the management of native forests for the production of timber, charcoal, and firewood production utilizing the full range of silviculture options, as described in [54].Conservation should address the environmental values that are dependent on the forest structure.Moreover, activities such as grazing and timber supply, for lease fencing and firewood, respectively, should be allowed only on MZ-B.To determine final forest management zoning in MZ-A and MZ-B, more field investigations and data collection are recommended.Plantation zones should be managed to maximize sustainable timber production based on reliable estimates of annual allowable harvest.However, more studies are needed to adequately evaluate the present state of Angolan forest plantations.
Since growth models (diameter and height growth models), height-diameter models, site productivity models, and stem volume models are important tools for decision-making in sustainable forest management planning, they were considered in this study.However, growth models and site productivity models (dominant height-age models) were not able to be developed due to lack of time series data (repeatedly measured data).The land cover assessment helped in developing management zones where, possibly in the future, we will be able to do this when we collect data from the same sample plots.All these models help determine the annual increment and allowable harvest in every year from any forest, and therefore our height-diameter models and volume models will be helpful in formulating better management plans for Miombo forests.

Conclusions
The primary focus of this study was to define patterns of deforestation and propose a management framework for the degraded areas based on our intensive characterization of the forest attributes, including the use of basic height-diameter models.The expansion of agricultural land has been the primary driver of deforestation.The conclusion is that even during the conflict period and afterward, the deforestation rate never decreased.In all cases, informed creation of management zones that are consistent over time and transparent for the purposes of systematic monitoring is critical to the success of sustainable forestry in Angola.Assessments of forest degradation using high-resolution imagery can be useful to map tree cover loss in support of this goal.
Based on our results (Figures 5 and 6), it appears that it is not simply the number of trees removed that threatens the stability of Miombo forests, but also the patterns of deforestation associated with the demand for new agricultural lands and charcoal and firewood production for cooking and heating.The results demonstrate that deforestation is a serious environmental problem in Angola, where the ongoing vegetation loss (17% from year 2000 to 2017), including woodland forest, in the past few decades is accentuated.However, the use of remote-sensing technology in land use classification in Angola is limited, particularly in forest or deforestation change analysis.
The four management zones are distinguishable on the map and growth height-diameter models developed represent a fundamental tool for future studies on forest planning for Miombo forests in Huambo Province, Angola.Using a combination of LC and LU data with topographic data, we concluded that certain types of changes are influenced by terrain morphology.

21 Figure 2 .
Figure 2. Frequency distribution of the tree species in Miombo forests in Huambo.

Figure 2 .
Figure 2. Frequency distribution of the tree species in Miombo forests in Huambo.

Figure 3 .
Figure 3. Sample plot-specific height-diameter curves overlaid on the measured data.Curves were produced with stand basal area (BA) variation at 2 m 2 •ha −1 interval within the observed data range.(a) curves produced with ordinary least square height-diameter model, (b) curve produced with mixed-effect mode for which one randomly selected tree per sample plot was used to estimate random effect parameters of Equation (2) and calibrate (localize) the model.

Figure 3 .
Figure 3. Sample plot-specific height-diameter curves overlaid on the measured data.Curves were produced with stand basal area (BA) variation at 2 m 2 •ha −1 interval within the observed data range.(a) curves produced with ordinary least square height-diameter model, (b) curve produced with mixed-effect mode for which one randomly selected tree per sample plot was used to estimate random effect parameters of Equation (2) and calibrate (localize) the model.

Figure 4 .
Figure 4. Stem volume (m 3 ha −1 ) curves by diameter at breast height (DBH) for each vegetation type; the curves were produced using parameter estimates (Table5), mean values of covariate predictors, and height at five-meter intervals.

Figure 4 .
Figure 4. Stem volume (m 3 ha −1 ) curves by diameter at breast height (DBH) for each vegetation type; the curves were produced using parameter estimates (Table5), mean values of covariate predictors, and height at five-meter intervals.

Figure 5 .
Figure 5. Deforestation trends between 2000 and 2017 in a satellite-derived NDVI.

Figure 5 .
Figure 5. Deforestation trends between 2000 and 2017 in a satellite-derived NDVI.

Figure 6 .
Figure 6.Example of forest degraded areas in different landscape categories between 2000 (top row) and 2017 (bottom row): natural forest (left), developed areas (middle), and forest plantations (top right).

Figure 6 .
Figure 6.Example of forest degraded areas in different landscape categories between 2000 (top row) and 2017 (bottom row): natural forest (left), developed areas (middle), and forest plantations (top right).

Figure 6 .
Figure 6.Example of forest degraded areas in different landscape categories between 2000 (top row) and 2017 (bottom row): natural forest (left), developed areas (middle), and forest plantations (top right).

Figure 8 .
Figure 8. Classification of the study area into four management zones using watershed areas and land cover types.

Figure 8 .
Figure 8. Classification of the study area into four management zones using watershed areas and land cover types.

Table 1 .
Representation of matrix calculation.

Table 3 .
Training data used in supervised classification process.

Table 4 .
Parameter estimates, variance components, and fit statistics of the mixed effects model (Equation

Table 5 .
Parameter estimates and fit statistics of the volume model (Equation3).

Table 5 .
Parameter estimates and fit statistics of the volume model (Equation3).

Table 6 .
The percentages of land use (LU) and land cover (LC) changes from 2000 to 2017 for the entire study area.The value for other land cover is quite large and represents bare land and agriculture land together.

Table 6 .
The percentages of land use (LU) and land cover (LC) changes from 2000 to 2017 for the entire study area.The value for other land cover is quite large and represents bare land and agriculture land together.

Table 7 .
The percentages of LU/LC changes for the "Caminho de ferro de Benguela" Benguela railway line (Figure1).

Table 7 .
The percentages of LU/LC changes for the "Caminho de ferro de Benguela" Benguela railway line (Figure1).