Models and heuristics for forest management with environmental restrictions

Tese de doutoramento, Estatística e Investigação Operacional (Otimização), Universidade de Lisboa, Faculdade de Ciências, 2018The main focus of this thesis was to develop mathematical models and methods in integer programming for solving harvest scheduling problems with environmental restrictions. Constraints on maximum clearcut area, minimum total habitat area, minimum total core area and inter-habitat connectivity were addressed for this purpose. The research was structured in a collection of three papers, each one describing the study of a different forest harvest scheduling problem with respect to the environmental constraints. Problems of papers 1 and 2 aim at maximizing the net present value. A bi objective problem is considered in paper 3. The objectives are the maximization of the net present value and the maximization of the inter-habitat connectivity. The tree search methods branch-and-bound and multiobjective Monte Carlo tree search were designed specifically to solve the problems. The methods could be used as heuristics, as a time limit of 2 hours was imposed. All harvest scheduling problems were based on the socalled cluster formulation. The proposed models and methods were tested with sixteen real and hypothetical instances ranging from small to large. The results obtained for branch-and-bound and Monte Carlo tree search show that these methods were able to find solutions for all instances. The results suggest that it is possible to address the environmental restrictions with small reductions of the net present value. With respect to the forestry fragmentation caused by harvestings, the results suggest that, although clearcut size constraints tend to disperse clearcuts across the forest, compromising the development of large habitats, close to each other, the proposed models, with the other environmental constraints, attempt to mitigate this effect

Forest planning utilizing high spatial resolution data

This thesis presents planning approaches adapted for high spatial resolution data from remote sensing and evaluate whether such approaches can enhance the provision of ecosystem services from forests. The presented methods are compared with conventional, stand-level methods. The main focus lies on the planning concept of dynamic treatment units (DTU), where treatments in small units for modelling ecosystem processes and forest management are clustered spatiotemporally to form treatment units realistic in practical forestry. The methodological foundation of the thesis is mainly airborne laser scanning data (raster cells 12.5x12.5 m2), different optimization methods and the forest decision support system Heureka. Paper I demonstrates a mixed-integer programming model for DTU planning, and the results highlight the economic advances of clustering harvests. Paper II and III presents an addition to a DTU heuristic from the literature and further evaluates its performance. Results show that direct modelling of fixed costs for harvest operations can improve plans and that DTU planning enhances the economic outcome of forestry. The higher spatial resolution of data in the DTU approach enables the planning model to assign management with higher precision than if stand-based planning is applied. Paper IV evaluates whether this phenomenon is also valid for ecological values. Here, an approach adapted for cell-level data is compared to a schematic approach, dealing with stand-level data, for the purpose of allocating retention patches. The evaluation of economic and ecological values indicate that high spatial resolution data and an adapted planning approach increased the ecological values, while differences in economy were small. In conclusion, the studies in this thesis demonstrate how forest planning can utilize high spatial resolution data from remote sensing, and the results suggest that there is a potential to increase the overall provision of ecosystem services if such methods are applied

Optimization Model of Integrated Sustainable Forest Management Planning for Hydropower Power Plant

Forest management planning is placed not only on the management of timber production but also on the maintenance of forest ecosystem values, such as water conservation. Thus, the meaning of sustainable forest management has expanded from sustainable wood management to ecosystem-based multipurpose forest management. This paper discusses the need to conserve water by forests in such a way that water will flow into rivers to be used as a power source for hydroelectric power generation. The problem is then modelled as a multi-objective integer program. In this research, the model is referred to as the planning model. The planning model is filled with data and completed by appropriate procedures. The solution method is often a numerical optimization algorithm. The algorithm used to solve the model also affects the formulation of the problem, namely the type of planning model used. Based on the data analysis, the results revealed that the objectives have met in planning sustainable and integrated forest management for hydropower such as: 1) maximize the Net Present Value (NPV) of the forest within the planning horizon, 2) maximize the equity of volume harvested in each harvesting period (subject to policy constraints). In addition to 3) maximize area control criteria to determine regulated forest or age forest, 4) maximize the final inventory criteria ensuring a solution where the forest inventory will be greater than or equal to the initial inventory, 5) maximize the total carbon balance in the planning period, 6) maximize protected areas for water protection

Designing decision support tools for Mediterranean forest ecosystems management: a case study in Portugal

Original ArticleThe effectiveness of Mediterranean forest ecosystem management calls for the conceptualization and implementation of adequate decision support tools. The proposed decision support system encompasses a management information system, a prescription simulator, a constraint generator and a set of management models designed to solve decision problems. Emphasis is on the architecture of the prescription simulator and its linkage to the three other modules, as well as on methods for reporting and visualizing solutions. Results are discussed for a real world test case – Serra de Grândola, a management area with about 18 600 ha comprising 860 cork oak (Quercus suber L.) land units. Cork oak silviculture adds complexity to the traditional forest management problem. Results show that the devised system is able to address effectively the integration of ecosystem data, silviculture, growth-and-yield and management models. They further suggest that the proposed system architecture may help address the complexity of Mediterranean ecosystem management problemsinfo:eu-repo/semantics/publishedVersio

Optimal Forest Strategies for Addressing Tradeoffs and Uncertainty in Economic Development under Old-Growth Constraints

In Canada, governments have historically promoted economic development in rural regions by promoting exploitation of natural resources, particularly forests. Forest resources are an economic development driver in many of the more than 80% of native communities located in forest regions. But forests also provide aboriginal people with cultural and spiritual values, and non-timber forest amenities (e.g., biodiversity, wildlife harvests for meat and fur, etc.), that are incompatible with timber exploitation. Some cultural and other amenities can only be satisfied by maintaining a certain amount of timber in an old-growth state. In that case, resource constraints might be too onerous to satisfy development needs. We employ compromise programming and fuzzy programming to identify forest management strategies that best compromise between development and other objectives, applying our models to an aboriginal community in northern Alberta. In addition to describing how mathematical programming techniques can be applied to regional development and forest management, we conclude from the analysis that no management strategy is able to satisfy all of the technical, environmental and social/cultural constraints and, at the same time, offer aboriginal peoples forest-based economic development. Nonetheless, we demonstrate that extant forest management policies can be improved upon.forest-dependent aboriginal communities, boreal forest, compromise and fuzzy programming, sustainability and uncertainty, International Development, Resource /Energy Economics and Policy, R11, Q23, Q01, C61,

Source-tracking cadmium in New Zealand agricultural soils: a stable isotope approach

Cadmium (Cd) is a toxic heavy metal, which is accumulated by plants and animals and therefore enters the human food chain. In New Zealand (NZ), where Cd mainly originates from the application of phosphate fertilisers, stable isotopes can be used to trace the fate of Cd in soils and potentially the wider environment due to the limited number of sources in this setting. Prior to 1997, extraneous Cd added to soils in P fertilisers was essentially limited to a single source, the small pacific island of Nauru. Analysis of Cd isotope ratios (ɛ114/110Cd) in Nauru rock phosphate, pre-1997 superphosphate fertilisers, and Canterbury (Lismore Stony Silt Loam) topsoils (Winchmore Research Farm) has demonstrated their close similarity with respect to ɛ114/110Cd. We report a consistent ɛ114/110Cd signature in fertiliser-derived Cd throughout the latter twentieth century. This finding is useful because it allows the application of mixing models to determine the proportions of fertiliser-derived Cd in the wider environment. We believe this approach has good potential because we also found the ɛ114/110Cd in fertilisers to be distinct from unfertilised Canterbury subsoils. In our analysis of the Winchmore topsoil series (1949-2015), the ɛ114/110Cd remained quite constant following the change from Nauru to other rock phosphate sources in 1997, despite a corresponding shift in fertiliser ɛ114/110Cd at this time. We can conclude that to the present day, the Cd in topsoil at Winchmore still mainly originates from historical phosphate fertilisers. One implication of this finding is that the current applications of P fertiliser are not resulting in further Cd accumulation. We aim to continue our research into Cd fate, mobility and transformations in the NZ environment by applying Cd isotopes in soils and aquatic environments across the country

On the Use of Imaging Spectroscopy from Unmanned Aerial Systems (UAS) to Model Yield and Assess Growth Stages of a Broadacre Crop

Snap bean production was valued at $363 million in 2018. Moreover, the increasing need in food production, caused by the exponential increase in population, makes this crop vitally important to study. Traditionally, harvest time determination and yield prediction are performed by collecting limited number of samples. While this approach could work, it is inaccurate, labor-intensive, and based on a small sample size. The ambiguous nature of this approach furthermore leaves the grower with under-ripe and over-mature plants, decreasing the final net profit and the overall quality of the product. A more cost-effective method would be a site-specific approach that would save time and labor for farmers and growers, while providing them with exact detail to when and where to harvest and how much is to be harvested (while forecasting yield). In this study we used hyperspectral (i.e., point-based and image-based), as well as biophysical data, to identify spectral signatures and biophysical attributes that could schedule harvest and forecast yield prior to harvest. Over the past two decades, there have been immense advances in the field of yield and harvest modeling using remote sensing data. Nevertheless, there still exists a wide gap in the literature covering yield and harvest assessment as a function of time using both ground-based and unmanned aerial systems. There is a need for a study focusing on crop-specific yield and harvest assessment using a rapid, affordable system. We hypothesize that a down-sampled multispectral system, tuned with spectral features identified from hyperspectral data, could address the mentioned gaps. Moreover, we hypothesize that the airborne data will contain noise that could negatively impact the performance and the reliability of the utilized models. Thus, We address these knowledge gaps with three objectives as below: 1. Assess yield prediction of snap bean crop using spectral and biophysical data and identify discriminating spectral features via statistical and machine learning approaches. 2. Evaluate snap bean harvest maturity at both the plant growth stage and pod maturity level, by means of spectral and biophysical indicators, and identify the corresponding discriminating spectral features. 3. Assess the feasibility of using a deep learning architecture for reducing noise in the hyperspectral data. In the light of the mentioned objectives, we carried out a greenhouse study in the winter and spring of 2019, where we studied temporal change in spectra and physical attributes of snap-bean crop, from Huntington cultivar, using a handheld spectrometer in the visible- to shortwave-infrared domain (400-2500 nm). Chapter 3 of this dissertation focuses on yield assessment of the greenhouse study. Findings from this best-case scenario yield study showed that the best time to study yield is approximately 20-25 days prior to harvest that would give out the most accurate yield predictions. The proposed approach was able to explain variability as high as R2 = 0.72, with spectral features residing in absorption regions for chlorophyll, protein, lignin, and nitrogen, among others. The captured data from this study contained minimal noise, even in the detector fall-off regions. Moving the focus to harvest maturity assessment, Chapter 4 presents findings from this objective in the greenhouse environment. Our findings showed that four stages of maturity, namely vegetative growth, budding, flowering, and pod formation, are distinguishable with 79% and 78% accuracy, respectively, via the two introduced vegetation indices, as snap-bean growth index (SGI) and normalized difference snap-bean growth index (NDSI), respectively. Moreover, pod-level maturity classification showed that ready-to-harvest and not-ready-to-harvest pods can be separated with 78% accuracy with identified wavelengths residing in green, red edge, and shortwave-infrared regions. Moreover, Chapters 5 and 6 focus on transitioning the learned concepts from the mentioned greenhouse scenario to UAS domain. We transitioned from a handheld spectrometer in the visible to short-wave infrared domain (400-2500 nm) to a UAS-mounted hyperspectral imager in the visible-to-near-infrared region (400-1000 nm). Two years worth of data, at two different geographical locations, were collected in upstate New York and examined for yield modeling and harvest scheduling objectives. For analysis of the collected data, we introduced a feature selection library in Python, named “Jostar”, to identify the most discriminating wavelengths. The findings from the yield modeling UAS study show that pod weight and seed length, as two different yield indicators, can be explained with R2 as high as 0.93 and 0.98, respectively. Identified wavelengths resided in blue, green, red, and red edge regions, and 44-55 days after planting (DAP) showed to be the optimal time for yield assessment. Chapter 6, on the other hand, evaluates maturity assessment, in terms of pod classification, from the UAS perspective. Results from this study showed that the identified features resided in blue, green, red, and red-edge regions, contributing to F1 score as high as 0.91 for differentiating between ready-to-harvest vs. not ready-to-harvest. The identified features from this study is in line with those detected from the UAS yield assessment study. In order to have a parallel comparison of the greenhouse study against the UAS study, we adopted the methodology employed for UAS studies and applied it to the greenhouse studies, in Chapter 7. Since the greenhouse data were captured in the visible-to-shortwave-infrared (400-2500 nm) domain, and the UAS study data were captured in the VNIR (400-1000 nm) domain, we truncated the spectral range of the collected data from the greenhouse study to the VNIR domain. The comparison experiment between the greenhouse study and the UAS studies for yield assessment, at two harvest stages early and late, showed that spectral features in 450-470, 500-520, 650, 700-730 nm regions were repeated on days with highest coefficient of determination. Moreover, 46-48 DAP with high coefficient of determination for yield prediction were repeated in five out of six data sets (two early stages, each three data sets). On the other hand, the harvest maturity comparison between the greenhouse study and the UAS data sets showed that similar identified wavelengths reside in ∼450, ∼530, ∼715, and ∼760 nm regions, with performance metric (F1 score) of 0.78, 0.84, and 0.9 for greenhouse, 2019 UAS, and 2020 UAS data, respectively. However, the incorporated noise in the captured data from the UAS study, along with the high computational cost of the classical mathematical approach employed for denoising hyperspectral data, have inspired us to leverage the computational performance of hyperspectral denoising by assessing the feasibility of transferring the learned concepts to deep learning models. In Chapter 8, we approached hyperspectral denoising in spectral domain (1D fashion) for two types of noise, integrated noise and non-independent and non-identically distributed (non-i.i.d.) noise. We utilized Memory Networks due to their power in image denoising for hyperspectral denoising, introduced a new loss and benchmarked it against several data sets and models. The proposed model, HypeMemNet, ranked first - up to 40% in terms of signal-to-noise ratio (SNR) for resolving integrated noise, and first or second, by a small margin for resolving non-i.i.d. noise. Our findings showed that a proper receptive field and a suitable number of filters are crucial for denoising integrated noise, while parameter size was shown to be of the highest importance for non-i.i.d. noise. Results from the conducted studies provide a comprehensive understanding encompassing yield modeling, harvest scheduling, and hyperspectral denoising. Our findings bode well for transitioning from an expensive hyperspectral imager to a multispectral imager, tuned with the identified bands, as well as employing a rapid deep learning model for hyperspectral denoising

Computer-based tools for supporting forest management. The experience and the expertise world-wide

Report of Cost Action FP 0804 Forest Management Decision Support Systems (FORSYS)Computer-based tools for supporting forest management. The experience and the expertise world-wide answers a call from both the research and the professional communities for a synthesis of current knowledge about the use of computerized tools in forest management planning. According to the aims of the Forest Management Decision Support Systems (FORSYS) (http://fp0804.emu.ee/) this synthesis is a critical success factor to develop a comprehensive quality reference for forest management decision support systems. The emphasis of the book is on identifying and assessing the support provided by computerized tools to enhance forest management planning in real-world contexts. The book thus identifies the management planning problems that prevail world-wide to discuss the architecture and the components of the tools used to address them. Of importance is the report of architecture approaches, models and methods, knowledge management and participatory planning techniques used to address specific management planning problems. We think that this synthesis may provide effective support to research and outreach activities that focus on the development of forest management decision support systems. It may contribute further to support forest managers when defining the requirements for a tool that best meets their needs. The first chapter of the book provides an introduction to the use of decision support systems in the forest sector and lays out the FORSYS framework for reporting the experience and expertise acquired in each country. Emphasis is on the FORSYS ontology to facilitate the sharing of experiences needed to characterize and evaluate the use of computerized tools when addressing forest management planning problems. The twenty six country reports share a structure designed to underline a problem-centric focus. Specifically, they all start with the identification of the management planning problems that are prevalent in the country and they move on to the characterization and assessment of the computerized tools used to address them. The reports were led by researchers with background and expertise in areas that range from ecological modeling to forest modeling, management planning and information and communication technology development. They benefited from the input provided by forest practitioners and by organizations that are responsible for developing and implementing forest management plans. A conclusions chapter highlights the success of bringing together such a wide range of disciplines and perspectives. This book benefited from voluntary contributions by 94 authors and from the involvement of several forest stakeholders from twenty six countries in Europe, North and South America, Africa and Asia over a three-year period. We, the chair of FORSYS and the editorial committee of the publication, acknowledge and thank for the valuable contributions from all authors, editors, stakeholders and FORSYS actors involved in this project