1. Introduction
Decision support systems (DSS) originated as important tools to aid decision-makers in business administration and industrial management in the 1960s [
1]. DSS technologies were subsequently widely adopted in many other disciplinary areas, including natural resource management, in which DSS had become widespread by the mid-1980s [
2]. Spatial DSS in particular began evolving in the late 1990s, coinciding with significant advances in geographic information systems (GIS) and computing power in this timeframe [
3]. The Ecosystem Management Decision Support (EMDS) system was one of the earliest examples of a spatial DSS for environmental analysis and planning [
4,
5], but the system has steadily evolved since that time, and remains in wide use up to the present [
6]. The EMDS system is an application framework for constructing DSS for integrated, multiscale analysis and planning [
7]. It supports landscape-level analyses through logic and decision engines integrated with the ArcGIS
® geographic information system (GIS, ArcGIS 10.5, 2017, Environmental Systems Research Institute, Redlands, CA, USA), and the open source GIS systems QGIS, MapWindow, and DotSpatial.
The logic engine evaluates landscape data against a logic model designed in the NetWeaver Developer
® system (NetWeaver2 Developer, 2017, Rules of Thumb, PA, USA) [
8] to derive logic-based interpretations of complex ecosystem conditions, such as watershed conditions, wildfire potential, and landscape integrity. Logic modeling can be defined as a form of knowledge-based reasoning in which expert or tacit knowledge is used to construct a representation of how to think about and solve a problem [
9]. In the EMDS context, in particular, logic is used to reason about the meaning of environmental states (e.g., given the data from the prescriptions on a management area, what can we infer about the score of each stand with respect to the criteria values for each nondominated solution?). EMDS decision models are built with Criterium DecisionPlus
® (CDP, version 4.0, 2014, InfoHarvest, Seattle, WA, USA), which implements the Analytical Hierarchy Process (AHP; [
10]), the Simple Multi-Attribute Rating Technique (SMART; [
11]), or their combination. The AHP belongs to the general class of multicriteria decision models (MCDMs; [
12]), while SMART derives from the closely related field of multiattribute utility theory (MAUT; [
13]). Engines that implement Multi-Criteria Decision Methods (MCDM) solutions in particular have become increasingly common in contemporary DSS architectures for environmental decision support, beginning in about 1995, either as stand-alone MCDM solutions or as MCDM components of larger systems such as EMDS. In contrast, the use of engines that implement formal logic in environmental decision support is less common in contemporary systems. The latter have often been built by encapsulating growth and yield models and management planning methods rather than logic approaches [
14]. However, logic processing has remained at the core of EMDS functionality since 1995 because of its capability to model the relatively large, complex, and often abstract problems associated with decision support for ecosystem management and adaptive management [
6].
On the other hand, the web-based DSS SADfLOR builds from the experience of developing decision support systems to enhance forestry and natural resource management planning [
15,
16,
17]. Its modular structure encompasses an information management module that encapsulates a web GIS and that is instrumental to organize administrative, geographic, and inventory forestry data. A simulation module projects the forest over time according to management options defined by stakeholders. It generates the management problem decision space. The latter is translated into a resource capability model by a matrix generator submodule. Decision-makers and stakeholders may input information that is translated into policy models (e.g., constraints on ecosystem services flows). An optimization module brings together the resource capability and the policy models. This module may be used to generate a solution to the management-planning problem by an a priori preference modeling multiple-criteria approach (e.g., goal programming). It may be used further to provide stakeholders and decision-makers with information about trade-offs between ecosystem services and trigger the application of an a posteriori preference modeling a multiple-criteria method (e.g., Pareto frontier). In this case, decision-makers can interact with the Pareto frontier submodule, analyze trade-offs, and select various nondominated solutions according to their preferences. In contexts characterized by multiple decision-makers, this submodule may be used to facilitate negotiation and consensus building (e.g., as seen in Reference [
18]). Nevertheless, diverse and conflicting preferences and objectives suggest the potential of combining the multiple-criteria functionality of both the web-based SADfLOR and CDP and NetWeaver to further enhance the negotiation of solutions.
A multicriteria decision problem involves choosing one among various nondominated alternatives on how well those alternatives rate against a chosen set of structured criteria [
19]. Thus, CDP and NetWeaver as multicriteria DSS can be used to compare the performance of multiple alternative solutions following a predefined model with criteria values resulting from the optimization process of SADfLOR. Both knowledge-based DSS are useful to solve this problem. NetWeaver is based on fuzzy logic networks that are well suited as logic frameworks for integrating model results from various analytical systems, such as simulators, linear programs, Bayesian belief networks, and production-rule systems [
20]. On the other hand, CDP provides extensive support for implementing this decision process, including setting priorities, explaining the derivation of conclusions, checking model robustness, and providing trade-off analyses. The criteria are weighted in terms of importance to the decision-maker and, in the end, after the evaluation by the model, the alternatives are ordered by their decision scores as a prioritized ranking of alternatives by preference [
19].
In the present study, we describe the use of a combination of decision support approaches encapsulated in the EMDS and the web-based SADfLOR systems to recommend optimal selections of alternative management combinations in the study area. Firstly, we use a spatially explicit multiple-criteria method within SADfLOR to generate both the decision space and the criteria space of the management-planning problem. This decision approach presents ecosystem services trade-off information in the form of Pareto frontiers to support the negotiation of alternative bundles of ecosystem services In the first aspatial phase of the analysis, we use the CDP component of EMDS to support further this negotiation by performing a multicriteria decision analysis (MCDA) that assesses the aggregate performance of selected points on the Pareto frontier for the study area as a whole. In the second spatial phase of the analysis, the NetWeaver component of EMDS is used to recommend selection of specific points on the Pareto frontier in specific landscape units of the study area. In the discussion, we conclude by considering the advantages and complexities of employing a hybrid DSS such as EMDS/SADfLOR to realize the potential of the combination of the decision approaches to support collaborative forest ecosystem management planning.
2. Materials and Methods
2.1. Study Area
For this study, we considered the Zonas de Intervenção Florestal (ZIF) Paiva and Entre-Douro e Sousa (ZIF_VS), located in northwestern Portugal (
Figure 1). The maximum elevation reaches about 700 m, yearly average temperatures range between 10 °C and 15 °C, and the main soil types are anthrosols, leptosols, and regosols. It spans three counties (Castelo de Paiva, Penafiel, and Paredes), covering an area of about 14,388 ha, and it is classified into 1976 stands for forest management purposes. The area is dominated by pure eucalypt (
Eucalyptus globulus Labill) stands and mixed stands of eucalypt and maritime pine (
Pinus pinaster Ait.)—about 66% and 33% of the ZIF_VS area, respectively, and hardwoods (mainly chestnut in 1% of the case study area). Community (local parish) property accounts for 35% of the ZIF_VS area. Medium and large private properties (area greater than 5 ha) extend over 60% of the ZIF_VS area. The remaining 5% is owned by small or very small forest owners [
19]. The stands are managed by each landowner with the support of the Local Forest Owners Association that is responsible for the development of the joint landscape-level management plan. The most important products are eucalypt pulpwood and maritime pine sawlogs, but the range of ecosystem services includes further hardwood (chestnut) sawlogs, carbon sequestration, and biodiversity. ZIF Paiva stands belong to 185 associated landowners, while ZIF_VS has 191 associated landowners. The web-based SADfLOR DSS has evolved with the involvement and participation of several stakeholders, ranging from the Instituto da Conservação da Natureza e das Florestas (ICNF), to the forest industry and nongovernmental organizations (e.g., forest owners associations) [
15,
21]. It was used to simulate all ZIF_VS stand-level prescriptions (59,133) over a 90-year planning horizon.
2.2. Overview of Decision Tools and Analytical Steps
In this work, we focused on the integrated application of three decision support approaches encapsulated in two DSS. The decision support approaches were: a multiple-criteria Pareto frontier method available in the web-based SADfLOR DSS, an aspatial strategic MCDA method available in the Criterium DecisionPlus (CDP) component of the EMDS DSS, and a spatial logic-based approach available in the NetWeaver component of EMDS.
The web-based tool encapsulates a multicriteria mathematical programming model and displays the trade-offs between ecosystem services in the form of Pareto frontiers. This analysis supports the selection of packages of ecosystem services according to the preferences of decision-makers. For that purpose, the ecosystem services are classified by the tool into two sets of criteria. The first includes wood products (e.g., total wood), while the second includes nonwood products as well as nonmarket services (e.g., biodiversity). The interface of the web-based tool asks the decision maker whether the provision of each ecosystem service (criteria) is to be maximized or minimized. After the selection by decision-makers of the list of ecosystem services to maximize or minimize in the bundle, the tool launches the process for generating the Pareto frontiers (e.g., as in References [
22,
23]). This information provides a basis for the negotiation of the packages of ecosystem services to be supplied by the decision-makers [
18]. These packages correspond to nondominated points in the Pareto frontier. After its selection, the web-based tool triggers the solution of the multicriteria mathematical programming problem associated to each point, i.e., it retrieves the management plan that may convey the corresponding package of services.
In this study, an aspatial strategic MCDA in CDP and a spatial logic-based approach in NetWeaver were both used to choose the best nondominated solution among a set of alternatives generated by the Pareto frontier tool in SADfLOR. CDP is an MCDA DSS that is useful for ranking alternatives based on how well each rates against a set of decision criteria [
24]. CDP implements both the analytic hierarchy process (AHP; [
10,
25]) and SMART [
11] techniques. We use the AHP method to derive weights on criteria and subcriteria by the pair-wise comparison technique, but we use SMART to score the performance of the attributes of alternatives. The combined use of AHP and SMART offers the advantage of eliminating, or at least greatly reducing, the well-known problem of rank reversal of alternatives in the classic AHP approach, in which scores on attributes of alternatives are derived by pair-wise comparisons [
24]. NetWeaver provides a logic-based framework that is useful to evaluate the conflicting alternative solutions chosen by different (groups of) decision-makers. It has two features with respect to the learning process: NetWeaver models provide detailed, permanent documentation about how we conceive that a system should be evaluated, and the modular architecture of logic models is well suited to the incremental, evolutionary adaptation of the models as knowledge changes [
24]. The overall analysis was performed in three steps:
Generation of the criteria space and selection of six alternative points on the Pareto frontier for subsequent analysis.
Aspatial multiple-criteria decision analysis of the aggregate performance of points on the Pareto frontier for the entire study region.
Spatial logic-based assessment of points on the Pareto frontier for individual stands of the study area.
Details of these three analytical steps are presented in the next three subsections and a flow chart of this methodology is provided (
Figure 2).
2.3. The Criteria Space and the Pareto Frontier
A multicriteria mathematical programming method within the web-based SADfLOR was used to define the forest management planning decision space. This was translated into the corresponding criteria space and Pareto frontiers by the interactive decision map technique. In this application, we emulated a negotiation process for the selection of six bundles of ecosystem services (six points in Pareto frontiers) for further analysis. The tool read this information and triggered the solution by the multicriteria mathematical programming method to find the corresponding six management plans.
2.4. Aspatial Multiple-Criteria Decision Analysis
A multiple-criteria decision analysis model within the CDP component of EMDS was used to compare the performance of multiple points on the Pareto frontier with respect to criteria and weights determined by a decision-maker or group of decision-makers. In this aspatial model, we considered the same criteria used in the first analytical step, but with solution variables aggregated into two primary sets of criteria: Wood Harvested and Other Ecosystem Services. In practice, the criteria weights will vary individually by decision-maker or by consensus in the case of a group. However, for purposes of demonstration in this study, we assumed a consensus across all decision-makers, in which the primary criteria are equally weighted. The key analysis metric for CDP is a measure of utility that expresses the preference among Pareto points, given the weight assigned to the primary criteria in the CDP model. Note that this a classical strategic MCDA analysis in which the alternatives are the aggregate results summarized by alternative solution points on the Pareto frontier.
2.5. Spatial Logic-Based Assessment of Points on the Pareto Frontier for Individual Stands
In the final step of the analysis, we used the weights developed in the second analytical step in a logic model that was executed by the NetWeaver component of EMDS to assess performance of the delivery of ecosystem services at the scale of individual stands. Each of the six nondominated points from the Pareto solution was associated with a map that was evaluated in EMDS, because the input data vary with the Pareto point selected. The NetWeaver model is structurally and computationally equivalent to the CDP model, but the key NetWeaver metric is a measure of strength of evidence. In this application, the topmost logic topic in the model was evaluating evidence that a particular Pareto solution delivers ecosystem services in a stand. One map of NetWeaver outcomes was generated in EMDS for each Pareto point. Maps are symbolized in terms of degree of evidence that the solution delivers ecosystem services. This analysis disaggregated the strategic MCDA analysis to individual stands.
The data inputs to NetWeaver included the ecosystem services target values in the six bundles selected by stakeholders and the corresponding management plans (e.g., one prescription per stand). Each data input to NetWeaver was evaluated by a fuzzy membership function that returned the strength of evidence that the observed stand-level value contributed to delivery of ecosystem services. Fuzzy membership functions in NetWeaver were scaled to the interval [–1, 1], in which –1 indicates no evidence, and 1 indicates full evidence. For each data input, the function was normalized to the observed data by setting the X-value corresponding to an evidence value of –1 as the minimum of the data distribution, considering the data range in the variable across all six Pareto solutions and all land stands. In general, data values were significantly right-skewed, so the X-value corresponding to an evidence value of 1 was set to the 90th percentile of the distribution for each input in order to provide a more useful distribution of outcomes on each input variable.
4. Discussion
The present study demonstrates how knowledge-based decision support approaches encapsulated in the EMDS can be used to complement the multicriteria decision methods within the Pareto frontiers tool of SADfLOR, when the two DSS are employed in a simple collaborative way. Importantly, no additional programming was required because the two EMDS components, CDP and NetWeaver, both received outputs more or less directly from SADfLOR, with only minimal data processing needed for data exchange between the two systems. To briefly summarize the analytical steps used in this study, we first obtained the six nondominated solutions from the multiple-criteria Pareto frontiers approach in SADfLOR (
Figure 3); next, an aspatial MCDA analysis with CDP was conducted to assess the aggregate performance of Pareto solutions A–F over all stands for the 90-year planning horizon with respect to delivery of ecosystem services (
Figure 4), and; finally, a NetWeaver logical model version of the CDP model (
Figure 6) was run in the EMDS system to assess the performance of the Pareto solutions at the scale of individual stands (
Figure 6 and
Figure 7).
Pareto frontier applications have been used to study various forestry and natural resources phenomena related to several ecosystem services and forest management and planning. Authors in Reference [
26] developed a Pareto frontier approach to assess trade-offs between timber-related criteria and spatial conditions correlated to other ecosystem services. Authors in Reference [
27] used a production possibility frontier-based approach to evaluate the effects of bioenergy production on timber production potential and biodiversity. SADfLOR was built from forest ecosystem decision support system architectures [
15,
16,
17] to evolve a Pareto frontier module that encapsulates the Feasible Goals Method/Interactive Decision Maps technique for interactive visualization of the web-based Pareto frontier [
23,
28]. Former developments of this multiple-criteria approach included the analysis of trade-offs between cork- and timber-related products at the regional level [
22] and the analysis of trade-offs between timber-related products and a wide range of nonwood forest products [
18] and regulatory services [
29] at the landscape level.
The multiple-criteria method was developed further by this research to include three wood products and a subset of other ecosystem services (e.g., biodiversity, average carbon stored). Moreover, its collaborative planning functionality [
18] was enhanced by the combination with other decision support approaches encapsulated in EMDS.
When there is only one decision-maker, or a group of decision-makers with similar preferences, the choice of only a single point on the Pareto frontier may be relatively easy. Nevertheless, in contexts characterized by multiple decision-makers, this may no longer be the case. Further negotiation may be needed to achieve a consensual solution. In this research, we have emulated such a situation where various (groups of) decision-makers with different preferences identify several optimal (nondominated) solutions in the criteria space (SADfLOR;
Figure 3). So, the follow-on application of the classical aspatial strategic decision model in CDP (
Figure 4) was used to elicit weights assigned by decision-makers to the ecosystem services and extend the analysis to compare the performance of six alternative Pareto solutions with respect to these preferences and the corresponding provisioning of ecosystem services (
Figure 5). Although the CDP model in our study was relatively small and simple, being limited to consideration of a few provisioning services, the approach is readily extendible to more complex models that also consider supporting, cultural, and regulating services.
In this study, for demonstration purposes, the CDP and NetWeaver analyses were each conducted with a single model. More realistically, the present work could be extended to account for different weighting preferences among stakeholder groups. For example, Delphi methods could be used in conjunction with the CDP model structure to distinguish weighting preferences among a few to several stakeholder groups. With the groups being identified, the landscape of the whole study area could be partitioned according to the stakeholder groups and customized analyses run for each stakeholder group. It might also be of interest to the full group of all landowners in the association to reassemble the separate maps of outcomes from NetWeaver into a single composite map. This is easy to do in EMDS.
Finally, this study demonstrated the collaborative application of the SADfLOR and EMDS systems, but this collaboration required a modest amount of human intervention in terms of moving data between the systems, so an interesting question to pose is, “Is it feasible to automate the interoperability of the two systems?” The answer is a qualified yes. The EMDS system was rearchitected between 2012 and 2015, and is now built on the Microsoft Windows Workflow Foundation [
30]. As a result, it is now relatively easy to integrate new analytical components into EMDS, albeit with a moderate amount of coding still required. For example, two new analytical engines were added into the EMDS architecture in the past year to support Bayesian network processing and Prolog-based decision trees. However, in the next year or two, EMDS will also integrate an open-source workflow editor, such as Trident [
31], which will greatly simplify the task of integrating third-party analytical components into its architecture [
32]. It is probably fair to say that open architecture DSS supported by workflow engines that simplify extensibility of DSS are the way of the future in continuing DSS development for environmental analysis and planning.