GIS-Based Aesthetic Appraisal of Short-Range Viewsheds of Coastal Dune and Forest Landscapes

: This paper demonstrates the possibilities of a Geographical Information System (GIS) for investigating and explicating the spatial variation of the short-range viewshed aesthetic appeal in a World Heritage coastal dune and forest area. The study pursues the following objectives: (1) develop and trial a GIS-based algorithm for computing the Aesthetic Appeal Index for a Short-Range Viewshed ( ¯æ ); (2) deliver an output map showing the spatial variation of the computed ¯æ values in the target territory and distribution of the zones with high scenic quality and potential aesthetic ecosystem services (PAES); and (3) assess management alternatives in zones with high PAES and high conservation value. This study combines two key innovative aspects. First, it integrates an objective digital map of habitats with subjective scenic preferences of coastal forest and dune landscapes based on psychophysical and cognitive perceptions of scenic beauty. Second, it applies a GIS-based algorithm to translate subjective scenic preferences to an output map of ¯æ . The study’s main conclusion is that the combined aesthetic appraisal of the immediate and foreground viewshed of coastal forests and dunes, by applying a specially created GIS algorithm, allows an assessment of the scenic quality of this landscape reliably in statistical terms.


Introduction
The aim of this study is to demonstrate the possibilities of a Geographical Information System for investigating and explicating the geographical variation of the aesthetic appeal of a UNESCO-listed coastal dune and forest landscape, interpreted in psychophysical and cognitive terms as Potential Aesthetic Ecosystem Services (PAES). Landscape perceptions result from the transactions between individuals and landscapes [1]. Experts have argued for many years that considering the public's aesthetic perceptions, judgments, needs, and demands could contribute to more effective and intelligent landscape planning and wiser resource use [2]. What is needed for that purpose is a visual aesthetic assessment of the geographic and temporal patterns of environmental conditions [3].
The role of aesthetic resources in the economic valuation of natural and semi-natural landscapes has increased in recent decades [4][5][6][7]. Currently, the research on the scenic quality and aesthetic value of landscapes, particularly those with a high recreational appeal, mainly focuses on aesthetic ecosystem services (AES) as part of a broader sphere of cultural ecosystem services (CES). Booth et al. [8] and Schirpke et al. [9] focus on modelling the scenic value of various regions in the USA and Europe as an essential component of AES (former) and CES (latter).

1.
Recreation and tourism: physical use of nature and landscape; 2.
Natural and cultural heritage: intellectual and scientific values.
However, Raumer et al. [13] argue that a CES approach separates aesthetics and recreation, impeding deeper insights into landscape research. Gosal and Ziv [14] and Kalinauskas et al. [15] concur, arguing that aesthetic values of landscapes interpreted as CES are context-specific and difficult to assess, generalize, and validate spatially. Hence, visual resources are a critical but often ignored type of CES [11]. On the other hand, Povilanskas et al. [16] and Urbis et al. [17,18] argue that, if properly assessed and interpreted, AES and CES could deliver technical tools for landscape conservation and management on the local and regional levels.
Traditionally, two contrasting approaches-'objectivist' and 'subjectivist'-have dominated methodology development in landscape aesthetics [19,20]. They delineate a crucial methodological divide concerning the answer to the critical question: "Is scenic quality an objective property inherent in the landscape or a subjective impression based on individual cognitive processes?" [19]. In other words, one approach presupposes that scenic quality is inherent in the actual setting, while the second approach assumes that scenic quality is "in the eyes of the beholder" [19][20][21].
Adherents of the objectivist (ecological/physical/formal aesthetic) paradigm characterize the landscape as an object. Consequently, descriptive inventories usually rely on the judgments and standards held by trained and skilled experts. This approach asserts that professionals can objectively determine the public's preferences with enough expertise without asking the public. Most scenic quality assessments conducted by forest and land management institutions worldwide maintain the objectivist paradigm. On the other hand, the "subjectivist" paradigm relies on scenic preferences resulting from the physical characteristics of the landscape in question, and the attitudes, expectations or perceptions of the observing individuals.
The subjective approach represents the perspective that the focus of a study is the individual perceiver, not the landscape itself. From the subjectivist viewpoint, the 'scenery' exists only as the individuals under scrutiny perceive it. Personal perception is the cornerstone of the "subjectivist" paradigm, with internalized attitudes and values comprising the units of study [22]. In other words, a perception-based approach emphasizes the human (subjective) view of the landscape and treats its biophysical features as stimuli that evoke relevant psychological responses. The subjectivist paradigm considers scenic quality as a human construct.
Such a human construct is based on the interpretation of what an observer perceives through memories, associations, imagination and any symbolism it evokes [19]. In perception-based assessments, the context implies that ecological objectives (e.g., diversity, sustainability) are not definable at any specific point but require consideration of a changing spatial mosaic of different environmental conditions. Determining relationships between the ecological quality and visual aesthetic quality of the landscape, then, must address a dynamic kaleidoscope of conditions that naturally and continually change in response to ecosystem processes [3].
In technical terms, the main difference between the 'objectivist' and 'subjectivist' approaches may be summarized as follows. The aesthetic valuation in the objectivist paradigm relies on the ecological, geographical and architectural criteria of the visual landscape units according to the structural and physiological principle of the scenic expression. Only suitably qualified experts can carry out such an assessment. The objectivist aesthetic valuation indicators determine the aesthetic appeal of the landscape based on the structure of individual spaces (open vs. closed, expressive vs. monotonous), and the ratio of various landscape elements (forested or urban spaces, water bodies, balmy microclimate). The obtained result is typically a comprehensive study describing many different parameters of landscape elements (forest cover, forest quality, visual transparence, the composition of visual elements and lines, forest age and structure).
The subjectivist aesthetic valuation relies on an emotional reaction to the scenic view here and now or its fixed photographic or IT-simulated image. In order to achieve reliable results, the research aims to address as many people as possible with different perspectives. Based on this paradigm, the methods of aesthetic-emotional assessment of the landscape reflect the dynamics of the variable aesthetic value of the society (community), depending on the user. The result obtained is presented in free emotional visual image descriptions or breakdowns according to the predicted emotional response scale. The scalability of the subjectivist aesthetic valuation is negligible; as the scale changes, the assessment changes little.
After scrutinizing 134 studies published between 2015 and 2021, we conclude that the objectivist paradigm primarily addresses the issues of landscape valorization on the national or regional scale [15,23]. In this way, it tends to overlook local landscape management units, which is a shortcoming of the objectivist approach. Meanwhile, the subjectivist paradigm emphasizes landscape perceptions by local people or visitors. In the age of social media, this may include a large group of virtual visitors, yet, the focus is still on local landscape units. Hence, the difficulty of scaling up results, generalizing and validating them is the main criticism of the subjectivist approach [15].
The psychophysical functions have proven to be robust to changes of landscapes and observers [26]. Therefore, the psychophysical approach is fruitful in scenic quality research, which is pertinent to forest recreation, since a perennial question on the academic tourism research agenda is the understanding of the tourist [33], and hence also her or his landscape perceptions and preferences [21].
Furthermore, Vining and Stevens [2] argued for relating the physical environment and the observer's psychophysical perception with her cognitive structure. They saw no reason to separate the psychophysical and cognitive approaches in studying the relationship between the physical traits of the landscape and the resulting perceptual response of the observer, as the latter's response is a function of the cognitive process. Experiencing involves the processes of seeing, knowing, believing and recognizing. These processes are cognitive and are based on individuals' background knowledge, learning, and reasoning capabilities [34].
The human factors affecting scenic preference include education, personality, profession, involvement, and socio-cultural groupings [20,35]. Hence, integration of the psychophysical and cognitive paradigms implies that the aesthetic quality of the landscape is a product of an interaction between the individual viewer and the landscape in which the human component comprises experience, future expectations and the sociocultural context, and the landscape component includes individual elements, scenes and habitats as entities [31,34]. In other words, the landscape relates the individual to their environment [36,37].
Such a complex view of landscape perception acknowledges underlying cognitive processes, such as environmental attitudes or acceptability of an impact, as well as the observer's particular frame of reference, which mediates the response to a visual stimulus [22]. Therefore, in landscape preference studies, we have to express quantitatively meanings and attitudes associated with landscapes bearing salient societal connotations and complex management implications, i.e., semi-natural forest and dune landscapes with high recreational and heritage value, like the Curonian Spit [38,39].
Such perspectives may provide a robust foundation for more cognitive perceptions [40]. In cognitively and emotionally charged cases of conservation and management of highquality scenery, quantitative methods of analysis and demonstration of the landscape structure and image are needed [41]. An integrated approach, which focuses on the psychophysical and cognitive nature of human-landscape interactions and their outcomes (including preference, values, satisfaction, well-being, predictions, and physical landscape elements or habitats) prevails in the scrutinized literature.
Furthermore, considering the effects of the human-centered integrated approach towards interpretation of the aesthetic quality of landscapes, the biophilic effect should be noted as technically resulting from the psychophysical and cognitive approaches. It is a technical toolbox used within landscape architecture to enhance the beneficial relationship of both residents and visitors with the environment through nature, professionally designed landscapes and habitats, and also, place and space conditions. It is argued that the biophilic effect has evidence-based health benefits for residents and visitors in various aesthetically appealing recreational environments [42][43][44].
The interest in spatial data and GIS for assessing the scenic quality of natural and seminatural landscapes started in the 1980s. Early attempts to apply GIS for the scenic quality assessment already incorporated visual values in sustainable landscape planning [45]. According to an observation by Bishop and Hulse [46] (p. 60): 'If visual values can be predicted using mapped data and the computational capabilities of a GIS, then there exists the potential for development of more objective and cost-effective procedures for assessment of visual qualities and the impact of change'.
Applying GIS for scenic quality valuation of landscapes with a recreational potential became popular in the 2000s with the rapid advance of the computational capabilities of GIS software, including statistical analysis, visualization and computer modelling [47][48][49]. Xu et al. [50] analyzed 92 studies of cultural and ecological landscape corridors in Europe, published between 1990 and 2018, focusing on their planning and management, the problems addressed, approaches and tools used, spatial scales and the stakeholders involved. They found that GIS analysis-often together with cost-distance analysis, graph-based analysis and spatial network analysis-was the most widely employed approach (82% of the surveyed studies) for planning and managing landscape corridors.
GIS-based procedures to value AES and, more generally, CES result from linking process-based models with non-market economic valuation methods [51][52][53]. For example, many case studies worldwide explore a framework that combines GIS with the Hedonic Pricing Method in analyzing the relationship of the value of green and open spaces with housing prices [54][55][56][57][58][59][60]. Therefore, the main aim of our work is to explore the relationship between the scenic dimensions, such as naturalness, heterogeneity and succession, and public perception of the scenic quality of semi-natural forest and dune landscapes with high recreational and heritage value in a new and innovative way.
Our study comprises the full range of scenes representing all forest and dune habitats and landscapes-from forest plantations to the natural dune habitats of interest for the European Union NATURA 2000 series necessary for the protection of biodiversity.
The study pursues the following concrete objectives: 1.
Create and trial a GIS-based algorithm for computing values of the Aesthetic Appeal Index for a Short-Range Viewshed (ǣ); 2.
Deliver a GIS map showing the spatial variation ofǣ values in the target territory and distribution of the zones with high scenic quality and, hence, Potential Aesthetic Ecosystem Services (PAES); 3.
Assess management alternatives at the sites with extraordinary PAES and high conservation value, suggesting the best hiking routes in strict nature reserves; 4.
Discuss the advantages and limitations of the methodology as a decision support process in coastal forest and dune management and tourism planning.
This study features two innovative elements: (1) Comprehensive analysis of scenic preferences for forest and dune landscapes in the immediate and foreground viewshed using a qualitative-quantitative combination of methods; and (2) creating and applying a valuation spreadsheet of the immediate and foreground habitat combinations of coastal forests and dunes and thus reducing the required hardware capacity and computation time. The study's main conclusion is that a combined valuation of the short-range, i.e., immediate and foreground viewshed of coastal forests and dunes on the Curonian Spit, applying an innovative GIS algorithm, allows assessment of the scenic quality of this series of habitats and landscapes reliably in statistical terms.

Study Area
The Curonian barrier spit separates the Curonian Lagoon from the open Baltic Sea ( Figure 1). The length of the Curonian Spit is 94 km and the width varies from 380 m to 4 km. It is the largest accumulative barrier sand spit in the Baltic Sea Region [61]. The varied and dynamic dune landscape with high biological diversity and a strong visual appeal is a distinctive feature of the Curonian Spit on a regional scale [62]. The 32.6-km long Great Curonian Dune Ridge of 40-60 m high mobile dunes is the longest coastal mobile dune ridge in Europe [63]. Politically, it is shared by Lithuania and Russia. The Curonian Spit extends from Zelenogradsk (Russia) in the southwest (54 •

Materials and Data Used
The three photos for each viewshed were obtained in the following manner. First, a full list of the dune, forest and successive transitional habitats of the Curonian Spit was compiled using a pre-existing national database of NATURA 2000 habitats. The photographer (Arvydas Urbis) then took a series of high resolution (1536 × 1024 pixels) monochromatic photos of each viewshed using a digital camera Nikon N5005 wide zoom with a 28-200 mm focal length Tamron lens, with a 75° degree field of view. To eliminate the The landscape of the Curonian Spit is symbolically significant for Lithuanians-for locals and tourists alike [64]. The specific semi-natural landscape with mobile dunes, untouched coastline and period cabins communicates the image archetypal to the Lithuanian Baltic coast. At the same time, the Curonian Spit, with its distinct nature and the sense of the tamed wilderness, is a favorite international tourist destination [39]. Due to its unique blend of natural, cultural and symbolic values, the entire Curonian Spit was included in the UNESCO World Heritage List as a transboundary cultural landscape of mobile and forested dunes of outstanding international importance [65].
The Great Curonian Dune Ridge of mobile dunes is also protected as a strict nature reserve within the Kurshskaya Kosa national park on the Russian part of the spit (est. 1987), and Kuršių nerija national park on the Lithuanian part (est. 1991). We have investigated the whole series of the dune habitats on the Curonian Spit, regarding their comparative aesthetic appeal and scenic quality, including dry sand heaths, mobile bare dunes ('white dunes'), fixed dunes with herbaceous vegetation ('grey dunes'), wet dune slacks, as well as all the transitional ecotones between these habitats and the forested ones.
The forest habitats included in our investigation comprise most natural forest and forest plantation types that are abundant on the Curonian Spit, due to a protracted history of dune afforestation [38]. Natural forests and forest plantations cover 62% of the total area of the spit. Referring to the results of our previous studies [17,18], the most visually appealing forest habitats for domestic Lithuanian visitors are: (i) mature Scots pine forest plantations with dry sand heaths in the background; (ii) mature Scots pine forest plantations; (iii) young Scots pine plantations with grey dunes in the background; (iv) middle-aged Silver birch stands; and (v) mature Norway spruce forest. The remaining forest habitats possess a much lesser visual appeal.

Materials and Data Used
The three photos for each viewshed were obtained in the following manner. First, a full list of the dune, forest and successive transitional habitats of the Curonian Spit was compiled using a pre-existing national database of NATURA 2000 habitats. The photographer (Arvydas Urbis) then took a series of high resolution (1536 × 1024 pixels) monochromatic photos of each viewshed using a digital camera Nikon N5005 wide zoom with a 28-200 mm focal length Tamron lens, with a 75 • degree field of view. To eliminate the potential bias of photographic representations showing a limited field of view, the photographer produced photographs presenting a viewshed within an angle of 150 • (double width wide zoom photos).
In the screening phase, a team of experts then unanimously selected the three photographs for each viewshed, taking care to ensure that all possible combinations of the surveyed landscape attributes (stages of succession, the naturalness and heterogeneity of scenery) were indeed represented by the images selected [16]. The selection process resulted in a list of 45 short-range viewsheds of various landscapes and habitats of the Curonian Spit. The remaining short-range viewsheds are so rare that visitors rarely see them during their visits to the spit.
We have conducted a pre-test of the distinguishability of the landscape images using a Chi squared-test and have proven (p < 0.00001) that lay visitors of the Curonian Spit successfully distinguish different landscapes and their combinations portrayed in the photographs [16]. Hence, the 135 images representing 45 forest and dune landscapes of the Curonian Spit, and their succession and management stages were used as visual stimuli for the paired-comparison survey. The respondents could see any ten of all possible pairs of the objects (N × (N−1)/2 = 990) if N = 45 sorted in random order.
The viewsheds were repeated 15 times, with each of the three photos representing each viewshed having equal opportunities to be paired with any of the photos representing the other viewsheds. This yielded a final sample of 14,850 pairs. Each respondent had equal opportunities to receive any of the 14,850 pairs of photographs for the paired choice provided in ten randomly sorted pairs in a flipper album. A single evaluator was asked to assess ten pairs of photographs, considering this number to be the biggest possible for an attentive and dedicated judgment by an individual [17,18], thus preventing the visual and mental evaluator's fatigue. We completed the paired comparison survey in July and August 2017. The population in the survey comprised all domestic summer visitors to the Lithuanian part of Curonian Spit, currently amounting to 387,000, according to the data from the toll-point. The sample (n = 1485) included randomly selected adult domestic visitors to the Curonian Spit (0.38% of the total visitors' population). The survey resulted in a ranking list of normalized values representing the relative aesthetic appeal values of 45 landscapes and habitats of the Curonian Spit. The aesthetic appeal values were computed using the method of paired comparisons for the normalized relative ranking of different landscapes of the Curonian Spit [5,[16][17][18][19]. We have ranked the photos of the examined landscapes and habitats relative to each other regarding their attractiveness. This was achieved by using a simple ranking technique and involved the computation of the proportion of times each photo was judged by the 1485 interviewees as more attractive than any other photo. The computation results were then transformed to normalized values representing the relative aesthetic appeal values of various landscapes. More technical details are given below in Section 2.3. Materials and data used. STAGE III (LANDSCAPE EVALUATION AND MANAGEMENT DECISION SUPPORT) is the final stage of the whole research project and comprises seven steps.

Step 11. Creating a Spreadsheet of the Normalized Aesthetic Appeal Values of the Immediate and Foreground Habitat Combinations of Coastal Forests and Dunes
First, we have interpreted all 45 ranked images in terms of their short-range viewshed information. Second, we have adapted the approach from [63][64][65] for a raster ArcGIS map. We have considered the central 20 m × 20 m cell as an observer's viewpoint, and the surrounding eight 20 m × 20 m cells as an immediate viewshed (i.e., the shortest-range viewshed, which is around the observer). Meanwhile, the following sixteen 20 m × 20 m cells were considered a foreground (i.e., the short-range viewshed next to the observer). The principal scheme of the full short-range viewshed is given in Figure 3. viewshed, which is around the observer). Meanwhile, the following sixteen 20 m × 20 m cells were considered a foreground (i.e., the short-range viewshed next to the observer). The principal scheme of the full short-range viewshed is given in Figure 3. We then rendered all of the data concerning the aesthetic appeal values of the immediate and foreground habitat combinations of the coastal forests and dunes (normalized from 0 to 1, where 1 is the highest score) into 45 columns and 45 lines of an Excel spreadsheet. The resulting cells of the spreadsheet were segregated into four types ( Figure 4): A. Aesthetic appeal values of the homogeneous habitats stretching over 100 m × 100 m area deduced directly from the field study results are given in the diagonal cells as the intersections of the same scenery views (Type A cells in Figure 4, 18 cases). The most common habitats include mobile (white) dunes, fixed (grey) dunes and mature Scots pine forest plantations. We then rendered all of the data concerning the aesthetic appeal values of the immediate and foreground habitat combinations of the coastal forests and dunes (normalized from 0 to 1, where 1 is the highest score) into 45 columns and 45 lines of an Excel spreadsheet. The resulting cells of the spreadsheet were segregated into four types ( Figure 4): A. Aesthetic appeal values of the homogeneous habitats stretching over 100 m × 100 m area deduced directly from the field study results are given in the diagonal cells as the intersections of the same scenery views (Type A cells in Figure 4, 18 cases). The most common habitats include mobile (white) dunes, fixed (grey) dunes and mature Scots pine forest plantations. B. Aesthetic appeal values of the ecotone habitats, bordering other different habitats at less than 30 m in any direction from a viewpoint cell when the immediate habitat around the viewpoint cell comprises impermeable thickets obscuring the foreground behind it (Type B cells in Figure 4, 5 cases). In that case we interpret the viewpoint cell as located within a homogeneous habitat, ignoring the foreground of the habitats behind it (Type B cells in Figure 4, five cases). Such cases include mature Mugo pine plantations, young Scots pine plantations, young stands of Silver birch, and Black alder and brushwood of various species. Finally, the aesthetic appeal values of the ecotone habitats, which border other habitats at less than 30 m in any direction from a viewpoint cell when the immediate habitat around the viewpoint cell provides a permeable viewshed and the combination of the immediate and foreground habitats and landscapes is not included into the 45 sceneries ranked by the participants in the field survey. In that case, an aesthetic appeal value of a habitat is given in a cell of the spreadsheet where the line represents an immediate habitat and the column represents a foreground habitat. In that case, the aesthetic appeal value of an immediate and foreground habitat combination is calculated using the formula (1) based on the aesthetic appeal value of both homogenous habitats comprising the whole short-range viewshed and taking the visual distance decay effect into account [66][67][68][69][70] (Type D cells in Figure 4, 331 cases): where V n is the resulting aesthetic appeal value of an ecotone habitat, V 1 is the measured aesthetic appeal value of the immediate habitat, and V 2 is the measured aesthetic appeal value of a foreground habitat behind the measured immediate habitat. The most common examples of Type D ecotone habitats include: mature Scots pine plantations in the immediate habitat with semi-mature Black alder stands in the foreground; mature Silver birch stands in the immediate habitat with young Scots pine plantations in the foreground; semi-mature Scots pine plantations in the immediate habitat with mature Mugo pine plantations in the foreground; mobile (white) dunes in the immediate habitat with semi-mature Silver birch stands in the foreground, as well as dry sand heaths in the immediate habitat with young Scots pine plantations in the foreground.  Finally, the aesthetic appeal values of the ecotone habitats, which border other habitats at less than 30 m in any direction from a viewpoint cell when the immediate habitat around the viewpoint cell provides a permeable viewshed and the combination of the immediate and foreground habitats and landscapes is not included into the 45 sceneries ranked by the participants in the field survey. In that case, an aesthetic appeal value of a habitat is given in a cell of the spreadsheet where the line represents an immediate habitat The most common examples of Type D ecotone habitats include: mature Scots pine plantations in the immediate habitat with semi-mature Black alder stands in the foreground; mature Silver birch stands in the immediate habitat with young Scots pine plantations in the foreground; semi-mature Scots pine plantations in the immediate habitat with mature Mugo pine plantations in the foreground; mobile (white) dunes in the immediate habitat with semi-mature Silver birch stands in the foreground, as well as dry sand heaths in the immediate habitat with young Scots pine plantations in the foreground.  [71]. The proposed algorithm is given in Figure  5 and explained below. First, an immediate and a foreground habitat and their types (A, B, C or D) have to be identified for each viewpoint in each of the surrounding 24 cells within the 360° field of view in a GIS-viewshed [70] in 16 directions ( Figure 6). If the viewpoint is in a homogeneous habitat stretching over an area of at least one hectare, the ǣ value of the viewpoint cell is considered equal to the scenic value for that habitat type quarried from the spreadsheet (Type A or B cells).  [71]. The proposed algorithm is given in Figure 5 and explained below. First, an immediate and a foreground habitat and their types (A, B, C or D) have to be identified for each viewpoint in each of the surrounding 24 cells within the 360 • field of view in a GIS-viewshed [70] in 16 directions ( Figure 6). If the viewpoint is in a homogeneous habitat stretching over an area of at least one hectare, theǣ value of the viewpoint cell is considered equal to the scenic value for that habitat type quarried from the spreadsheet (Type A or B cells).
If the viewpoint is located in an ecotone habitat, the computing ofǣ includes several additional conditions. In the directions where the immediate habitats are visually impermeable, only the scenic value of the immediate habitat, quarried from the spreadsheet (Type B cells), is used forǣ computation. In the directions where the immediate habitats are visually permeable, the scenic values of the foreground habitats are also included in the computation ofǣ. If this combination falls among the 45 images examined in the field study, we use the deduced values quarried from the spreadsheet (Type C cells); if not, we use the calculated ones (Type D cells).
Finally, for each viewpoint of the ecotone habitats, theǣ value is computed as the average of all deduced values aggregated over 100 m square grid cells and resulting from the scenic values of the whole short-range viewshed as an average of 16 directions (2): whereǣ n is a resulting aesthetic appeal value of n viewpoint of an ecotone habitat, and from V 1dir to V 16dir are aesthetic appeal values at each of the 16 directions representing the 360 • field of view of the short-range viewshed and taking into account all cases of the immediate and foreground habitat combinations ( Figure 6).  If the viewpoint is located in an ecotone habitat, the computing of ǣ includes several additional conditions. In the directions where the immediate habitats are visually impermeable, only the scenic value of the immediate habitat, quarried from the spreadsheet (Type B cells), is used for ǣ computation. In the directions where the immediate habitats are visually permeable, the scenic values of the foreground habitats are also included in the computation of ǣ. If this combination falls among the 45 images examined in the field study, we use the deduced values quarried from the spreadsheet (Type C cells); if not, we use the calculated ones (Type D cells).
Finally, for each viewpoint of the ecotone habitats, the ǣ value is computed as the average of all deduced values aggregated over 100 m square grid cells and resulting from the scenic values of the whole short-range viewshed as an average of 16 directions (2): ǣn = (V1dir + V2dir + V3dir + V4dir + V5dir + V6dir + V7dir + V8dir + V9dir + V10dir + V11dir + V12dir + V13dir + V14dir + V15dir + V16dir)/16, (2) where ǣn is a resulting aesthetic appeal value of n viewpoint of an ecotone habitat, and from V1dir to V16dir are aesthetic appeal values at each of the 16 directions representing the 360° field of view of the short-range viewshed and taking into account all cases of the immediate and foreground habitat combinations ( Figure 6).
To compute the ǣ value at each viewpoint of the ecotone habitats, a dynamic Excel matrix with 25 computational grid cells arranged in 16 directions is applied (Figure 7). The dynamic matrix allows computing of the ǣ value in all relevant short-range viewshed cases using the IF function of the Excel program for each grid cell array (3): where Vp(n) is the primary scenic value of a habitat, and Ci,j is the particular case of an immediate and foreground habitat combination. To compute theǣ value at each viewpoint of the ecotone habitats, a dynamic Excel matrix with 25 computational grid cells arranged in 16 directions is applied (Figure 7). The dynamic matrix allows computing of theǣ value in all relevant short-range viewshed cases using the IF function of the Excel program for each grid cell array (3): where V p(n) is the primary scenic value of a habitat, and C i,j is the particular case of an immediate and foreground habitat combination. To develop a GIS-based algorithm for computing the ǣ values, ArcGIS software (Esri, V10.1; Redlands, CA, USA) was applied. For saving the computation capacity and time of the available hardware, the target territory comprised only Nagliai Strict Nature Reserve of the Kuršių nerija national park and its buffer zone (ca. 395 hectares). The geographical extent of the study area in terms of latitude and longitude is as follows: from 55°24′55.23″ N, 21°2′50.65″ E in the northeast to 55°32′9.46″ N, 21°7′45.70″ E in the southwest.
In the study area, the topography varies from 1.0 m in dune slacks to 62.0 m at the highest peaks of the Great Dune Ridge. The average height of the study area is 4.7 m above sea level. The average height of the Great Dune Ridge in the analyzed area is 26.4 m above sea level. The target territory is the most representative of habitat and landscape diversity since it contains a whole range of dune and forest habitats and their combinations (Figure 8).
We have created the background map of the target territory from the spatial ArcGIS vector database of the State Forest Cadaster of the Republic of Lithuania for each forest  The target territory is the most representative of habitat and landscape diversity since it tains a whole range of dune and forest habitats and their combinations (Figure 8).
We have created the background map of the target territory from the spatial Ar vector database of the State Forest Cadaster of the Republic of Lithuania for each f plot given in the cadaster. According to the principles of national forest inventory in uania, dunes and dry sand heaths are also included in the forest cadaster. We hav ported the vector data of Nagliai Strict Nature Reserve and its buffer zone from the Forest Cadaster to a raster layer (cell size 20 × 20 m). As mentioned above, each cell o background map was interpreted in all three functions: as a viewpoint, an immediate itat of the viewpoint and a foreground habitat of the viewpoint. We have created the background map of the target territory from the spatial ArcGIS vector database of the State Forest Cadaster of the Republic of Lithuania for each forest plot given in the cadaster. According to the principles of national forest inventory in Lithuania, dunes and dry sand heaths are also included in the forest cadaster. We have exported the vector data of Nagliai Strict Nature Reserve and its buffer zone from the State Forest Cadaster to a raster layer (cell size 20 × 20 m). As mentioned above, each cell of the background map was interpreted in all three functions: as a viewpoint, an immediate habitat of the viewpoint and a foreground habitat of the viewpoint.

Step 14. Delivering an Output GIS Map of the Spatial Variation of theǣ Values in the Target Territory
The output GIS map of the spatial variation of theǣ values in Nagliai Strict Nature Reserve and its buffer zone was produced as an additional ArcGIS raster layer based on the computation of theǣ values for each raster cell of the territory done in Step 13 ( Figure 8). This helps to gain an overview of the spatial variation of theǣ values compared to the national park management regime layer, the extent of conservation policy restrictions in this protected coastal dune and forest area, and the location and availability of potential aesthetic ecosystem services (PAES) for tourists across the whole range of dune and forest landscapes. The output map resolution is 20 m × 20 m of each raster cell as for the background GIS map.

Step 15. Validating the Output GIS Map of the Spatial Variation of theǣ Values in the Target Territory
Suppose the statistical analysis shows a validated and significant correlation between GIS-based modelled aesthetic appraisal values of the short-range viewsheds and the visitors' scenic quality ratings. In that case, we can apply this output for the area-wide simulation of scenic quality using the GIS [72]. Therefore, to validate the output GIS map statistically, independently acquired data from a field study on aesthetic appraisal of short-range viewsheds of coastal dune and forest landscapes were used. The predictedǣ values for all the grid cells within the target territory of Nagliai Strict Nature Reserve were correlated with the elicited ones during a specially conducted on-site route survey.
The field data were collected via a questionnaire that was distributed to university students following the methodology proposed by Chhetri and Arrowsmith [73]. Field measurements of aesthetic appraisal of the short-range viewsheds were recorded by a convenience sample of 44 Klaipeda University students hiking along the walking trail or without any track in the open dunes and mature forest plantations (Figure 9). The students collected the field data by completing questionnaires. They had to evaluate the aesthetic appeal of the entire 360 • field of view on a ten-point Likert scale at the 100-stop series of observation points along the track (Figure 10).   The students were divided into four groups with 11 persons each. The first gr hiked in one direction on the first day in May of 2019 and the second group walked in reverse direction on the same day. The purpose was to see whether preceding events pacted upon responses given. The third and fourth groups repeated the survey, follow the same route and stopping at the same points on the second day, in May of 2019. B days were sunny with few clouds and air temperature hovering around 16° to 17°C, ical weather at the beginning of the tourist season on the Curonian Spit.
The length of the survey trail was 15.9 km (Figure 9). The stops for the 360° fiel view assessment of the aesthetic appraisal of the short-range viewshed were at each to 200 m aiming to choose the most representative dune and forest habitats sites. A gether, there were 100 valid points of the survey. The locations have been identified u the GPS receivers Garmin Oregon 600 with positioning information based on x and y thogonal coordinates in a WGS84 system. The position accuracy was 3 m. In addition two professional researchers guiding the groups recorded the biophysical details of locations in terms of forest succession, age and composition, and landscape features.
Hull and Stewart [74] and Chhetri et al. [34] observe that the mood of an aesth appraisal varies throughout a hiking survey. They propose the notion of 'experien landscape' containing three key elements: (1) encountered landscape, (2) sequence and feelings. Encountered landscape comprises the seen views, the people and the objects landscape. A sequence is an order in which these scenes or objects are encountered. F ings are the subjective emotions felt while encountering these views [34,74]. For valida the GIS-based model of the aesthetic appraisal of short-range viewsheds, at each po the students were asked to look around and judge the aesthetic appeal of the entire field of view, answering the question: "How attractive do you find the landscape rounding you?".
The last step in the validation of the GIS-based model of the aesthetic appraisa short-range viewsheds was the calculation of the correlation between the predicted thetic appeal ratings and the averaged aesthetic appeal ratings actually assessed in field at each of the 100 hiking event points using a bivariate correlation (Pearson's r). is a measure of linear association between the simulated and the actual, independent of ratings [75]. The calculation of the Pearson's correlation coefficient was completed standard Excel function PEARSON, which calculates the r values for the dataset autom ically. The students were divided into four groups with 11 persons each. The first group hiked in one direction on the first day in May of 2019 and the second group walked in the reverse direction on the same day. The purpose was to see whether preceding events impacted upon responses given. The third and fourth groups repeated the survey, following the same route and stopping at the same points on the second day, in May of 2019. Both days were sunny with few clouds and air temperature hovering around 16 • to 17 • C, typical weather at the beginning of the tourist season on the Curonian Spit.
The length of the survey trail was 15.9 km (Figure 9). The stops for the 360 • field of view assessment of the aesthetic appraisal of the short-range viewshed were at each 100 to 200 m aiming to choose the most representative dune and forest habitats sites. Altogether, there were 100 valid points of the survey. The locations have been identified using the GPS receivers Garmin Oregon 600 with positioning information based on x and y orthogonal coordinates in a WGS84 system. The position accuracy was 3 m. In addition, the two professional researchers guiding the groups recorded the biophysical details of the locations in terms of forest succession, age and composition, and landscape features.
Hull and Stewart [74] and Chhetri et al. [34] observe that the mood of an aesthetic appraisal varies throughout a hiking survey. They propose the notion of 'experienced landscape' containing three key elements: (1) encountered landscape, (2) sequence and (3) feelings. Encountered landscape comprises the seen views, the people and the objects in a landscape. A sequence is an order in which these scenes or objects are encountered. Feelings are the subjective emotions felt while encountering these views [34,74]. For validating the GIS-based model of the aesthetic appraisal of short-range viewsheds, at each point, the students were asked to look around and judge the aesthetic appeal of the entire 360 • field of view, answering the question: "How attractive do you find the landscape surrounding you?".
The last step in the validation of the GIS-based model of the aesthetic appraisal of short-range viewsheds was the calculation of the correlation between the predicted aesthetic appeal ratings and the averaged aesthetic appeal ratings actually assessed in the field at each of the 100 hiking event points using a bivariate correlation (Pearson's r). This is a measure of linear association between the simulated and the actual, independent sets of ratings [75]. The calculation of the Pearson's correlation coefficient was completed by a standard Excel function PEARSON, which calculates the r values for the dataset automatically.

Step 16. Interpreting the Spatial Variation of theǣ Values in the Target Area in Terms of Scenic Quality and PAES
To interpret the spatial variation of theǣ values in the target area in terms of scenic quality and PAES, a focus group was conducted in August 2021. Quoting a recent study by Gundumogula, who have summarized the main methodological principles of successful conducting of focus groups [76] (p. 301): 'Focus groups help to discover new aspects and information of one's research, as the participants own and contribute together much more and more diverse perspectives on the selected topic than the researcher could imagine alone. Focus groups have a high chance to catch and consider the peoples' feeling, ex-pressions, views, believes and responses while collecting the data.' The necessary information was obtained from an informal discussion with six practitioners and experts in a focus group comprised of people with first-hand knowledge of the sector, including an economist, a psychologist, a geographer, an environmentalist, and other professionals, in which this small number of experts were asked to informally discuss various aspects of scenic quality distribution in the target territory to include the practitioner's perspective. The aim of the focus group was to establish a conceptual link between the scenic quality and the PAES since the former notion is based on simulated and validated tourist offer-related aspects of the aesthetic appeal of coastal dunes and forests whilst the latter notion is tourist demand-related.
There are differing opinions on the optimal size of the focus groups ranging from four to twelve participants [76]. The moderator facilitated the intercourse within the focus group with a set of lightly structured open-ended questions. These questions, or rather topics, did not restrict the direction taken by each discussion stage. Participants were encouraged to discuss their experiences in applying the aesthetic appeal as an incentive to visit the coastal dunes and forests of the Curonian Spit. To facilitate the free flow and openness of the discussion, the note-keeper took written notes without audio or video recording of the session.

Step 17. Practical Decision-Support in Coastal Dune and Forest Management and Tourism Planning in the Target Territory
In order to determine the most effective landscaping measures by improving the aesthetic appeal of the landscape in the target territory (increasing the scenic value): 1.
Areas that are important to visitors for their aesthetic appeal (e.g., near tourist trails) have been identified. For this purpose, the existing or planned trail routes or the existing or planned beauty spots were marked on the forest management plot plans.

2.
In the target territory, with the help of the developed methodology, the areas of the lowest aesthetic value were determined. In cases where these areas are extensive, special localized scenic quality improvement measures have been proposed. The primary visual assessment threshold matched the area's clearly defined landscape elements (e.g., dune foot, dune slacks, hummocks or hollows). Furthermore, in selected areas of low aesthetic value, the GIS model helped determine which taxonomic parameters of the habitat or landscape reduce the short-range scenic value and which correction would allow for forming a more aesthetically appealing scenery.

3.
Specific landscaping measures (clear-cutting, thinning, planting) were selected in the identified plots. We have assessed how the essential parameters and scenery type may change. On this basis, landscaping measures were selected and implemented. After performing these, the plots were inventoried by determining the changed landscape scenic quality parameters, and the type of the newly formed scenery was specified. Finally, the calculations of newǣ values were performed for various forest and dune plots, and the resulting change in scenic quality value was estimated.

GIS Output Map of the Short-Range Viewshed Scenic Quality and Its Validation
As a result of computations using the developed GIS-based algorithm for computing the Aesthetic Appeal Index for a Short-Range Viewshed (ǣ), anArcGIS raster layer with a grid resolution of 20 m × 20 m cells based on the computation of theǣ values for each raster cell was created, showing the spatial variation of the computedǣ values in Nagliai Strict Nature Reserve and its buffer zone of Kuršių nerija national park and distribution of the zones with low, medium and high scenic quality ( Figure 11). This integrates an objective digital map of habitats with subjective scenic preferences of coastal forest and dune landscapes based on psychophysical and cognitive perceptions of scenic beauty and translates these subjective scenic preferences toǣ values.

GIS Output Map of the Short-Range Viewshed Scenic Quality and its Validation
As a result of computations using the developed GIS-based algorithm for computing the Aesthetic Appeal Index for a Short-Range Viewshed (ǣ), anArcGIS raster layer with a grid resolution of 20 m × 20 m cells based on the computation of the ǣ values for each raster cell was created, showing the spatial variation of the computed ǣ values in Nagliai Strict Nature Reserve and its buffer zone of Kuršių nerija national park and distribution of the zones with low, medium and high scenic quality ( Figure 11). This integrates an objective digital map of habitats with subjective scenic preferences of coastal forest and dune landscapes based on psychophysical and cognitive perceptions of scenic beauty and translates these subjective scenic preferences to ǣ values. We note a distinct positive Pearson's correlation between the modelled and observed short-range aesthetic appeal values at 100 viewpoints (r = 0.77, p < 0.0005). This result indicates high validity for a modelled short-range aesthetic appeal values using the created and applied GIS-based computation algorithm ( Figure 12). Roth et al. [23] indicate that, from a social sciences and psychology perspective, a value of r > 0.5 shows a strong effect, whereas from a geographical perspective, a value of 0.6 < r < 0.8 shows a distinct effect. For instance, Palmer [68] considers Pearson's r = 0.587 between the modelled scenic appeal ratings and the ratings made in the field at the 73 viewpoints used for visual simulations as sufficient [68]. We note a distinct positive Pearson's correlation between the modelled and observed short-range aesthetic appeal values at 100 viewpoints (r = 0.77, p < 0.0005). This result indicates high validity for a modelled short-range aesthetic appeal values using the created and applied GIS-based computation algorithm ( Figure 12). Roth et al. [23] indicate that, from a social sciences and psychology perspective, a value of r > 0.5 shows a strong effect, whereas from a geographical perspective, a value of 0.6 < r < 0.8 shows a distinct effect. For instance, Palmer [68] considers Pearson's r = 0.587 between the modelled scenic appeal ratings and the ratings made in the field at the 73 viewpoints used for visual simulations as sufficient [68].

GIS-Based Decision Support for Scenic Quality Management
The comparison of GIS data before and after the intervention measures into the dune and forest management, or changes that happened due to natural succession processes allowed us to create another ArcGIS raster layer indicating zones with the resulting scenic value changes of the near-range viewshed and their direction (positive-neutral-negative).
The obtained results allowed the analysis of the distribution of the scenic value of the landscapes in the territory based on the new ArcGIS raster layer reliably and to evaluate the suitability of the intervention measures that have taken place in terms of scenic quality ( Figure 13).

GIS-Based Decision Support for Scenic Quality Management
The comparison of GIS data before and after the intervention measures into the dune and forest management, or changes that happened due to natural succession processes allowed us to create another ArcGIS raster layer indicating zones with the resulting scenic value changes of the near-range viewshed and their direction (positive-neutral-negative). The obtained results allowed the analysis of the distribution of the scenic value of the landscapes in the territory based on the new ArcGIS raster layer reliably and to evaluate the suitability of the intervention measures that have taken place in terms of scenic quality ( Figure 13).

GIS-Based Decision Support for Scenic Quality Management
The comparison of GIS data before and after the intervention measures into the dune and forest management, or changes that happened due to natural succession processes allowed us to create another ArcGIS raster layer indicating zones with the resulting scenic value changes of the near-range viewshed and their direction (positive-neutral-negative). The obtained results allowed the analysis of the distribution of the scenic value of the landscapes in the territory based on the new ArcGIS raster layer reliably and to evaluate the suitability of the intervention measures that have taken place in terms of scenic quality ( Figure 13).  As mentioned, in the Nagliai Strict Nature Reserve, with the help of the developed methodology, the areas of low aesthetic value were determined. In cases where these areas are extensive, special localized scenic quality improvement measures have been proposed. The primary visual assessment threshold matched the area's clearly defined landscape elements (e.g., dune foot, dune slacks, hummocks or hollows). Furthermore, in selected areas of low aesthetic value, the GIS model helped to determine which taxonomic parameters of the habitat reduced the short-range scenic value, and correction measures were proposed to form a more aesthetically appealing scenery.
Last but not least, based on the research results, the methodology for planning, selecting and optimizing linear landscape elements was created and successfully tested. The scenic quality of the surrounding landscape of a self-guided nature trail in Nagliai Strict Nature Reserve (providing the only access for the visitors to the strict nature reserve) was assessed for three planned alternatives. Based on the results of computing the scenic quality, the route with the highest aesthetic value was proposed and successfully implemented as a pilot project ( Figure 14). It enabled balancing the psychophysical and cognitive aspects of the aesthetic appraisal of coastal dunes and forests with the conservation prescriptions of Kuršių nerija national park [39]. As mentioned, in the Nagliai Strict Nature Reserve, with the help of the developed methodology, the areas of low aesthetic value were determined. In cases where these areas are extensive, special localized scenic quality improvement measures have been proposed. The primary visual assessment threshold matched the area's clearly defined landscape elements (e.g., dune foot, dune slacks, hummocks or hollows). Furthermore, in selected areas of low aesthetic value, the GIS model helped to determine which taxonomic parameters of the habitat reduced the short-range scenic value, and correction measures were proposed to form a more aesthetically appealing scenery.
Last but not least, based on the research results, the methodology for planning, selecting and optimizing linear landscape elements was created and successfully tested. The scenic quality of the surrounding landscape of a self-guided nature trail in Nagliai Strict Nature Reserve (providing the only access for the visitors to the strict nature reserve) was assessed for three planned alternatives. Based on the results of computing the scenic quality, the route with the highest aesthetic value was proposed and successfully implemented as a pilot project ( Figure 14). It enabled balancing the psychophysical and cognitive aspects of the aesthetic appraisal of coastal dunes and forests with the conservation prescriptions of Kuršių nerija national park [39].

Discussion
The focus group helped to identify aesthetic criteria and attributes as one of the main reasons for choosing the Curonian Spit as a tourist destination at the seaside. Furthermore, it exposed why specific attributes were more important for tourists' satisfaction levels. As observed by Pérez-Hernández et al. [77], the scenic quality of mobile dunes with a natural and open background is usually much higher than that of any forested landscape. The aesthetic demand concepts were elicited through a focus group in which six practitioners with 20+ years of experience in tourism-related services and a background in marketing,

Discussion
The focus group helped to identify aesthetic criteria and attributes as one of the main reasons for choosing the Curonian Spit as a tourist destination at the seaside. Furthermore, it exposed why specific attributes were more important for tourists' satisfaction levels. As observed by Pérez-Hernández et al. [77], the scenic quality of mobile dunes with a natural and open background is usually much higher than that of any forested landscape. The aesthetic demand concepts were elicited through a focus group in which six practitioners with 20+ years of experience in tourism-related services and a background in marketing, psychology, landscape management, nature conservation, entertainment industry and geography were asked to informally discuss their involvement with dune aesthetics, their scenic preferences, and their connection to nature tourism.
As mentioned above, the interpretation of the spatial variation ofǣ values in the target area focused on establishing a conceptual link between the scenic quality and the PAES since the former notion relies on tourist offer-driven aspects of the aesthetic appeal of coastal dunes and forests whilst the latter notion is demand-driven. Participants discussed their experiences in applying the aesthetic appeal as an incentive to enjoy the coastal dunes and forests of the Curonian Spit. It became clear that, for some participants, the Curonian Spit meant primarily travelling to the mobile dune areas and nothing more, whereas for others, dune tourism follows naturally from enjoying the 'peace and quiet' of the Curonian Spit as part of their lifestyle.
Eight relishing aspects [78] have been elicited as the source of PAES of the dune and forest landscapes of the Curonian Spit from the focus group discussions (without rating any priority): (1) staginess; (2) visuality; (3) affectedness; (4) liminality; (5) apprehensiveness; (6) sense of togetherness; (7) uniqueness; and (8) knowledgeability. The focus group concurred that the perceived scenic quality leads to perceived PAES, which leads to visit intentions [79,80]. The concept of the image of the Curonian Spit as a visually unique destination did emerge in the focus group discussions. Still, the term brand did not. All focus group participants believed that the Curonian Spit image was an essential issue, including those who did not have any business or research interests there.

Conclusions
The most significant finding of our research was discovering the role of the double short-range viewshed in comprising the immediate viewshed and the foreground viewshed for the psychophysical investigations in landscape aesthetics. This discovery resulted in the practical development of an innovative GIS-based aesthetic appraisal methodology for coastal dunes and forests with high recreational and conservation value. Therefore, we may conclude that the combined aesthetic appraisal of the short-range viewshed, comprising the immediate and foreground viewshed of coastal forests and dunes, applying a specially created GIS algorithm, allows for assessing the scenic quality of this landscape reliably (Pearson's r = 0.77, p < 0.0005).
Relying on the results of our comprehensive research project on GIS-based aesthetic appraisal of the short-range viewsheds of coastal dune and forest landscapes, we created a decision-support system for the coastal dune and forest management in Kuršių nerija national park. This system converted the objective forest and dune habitat data into the scenic quality and PAES values. Therefore, the recommendations for further work focus on two interrelated avenues of research: (1) to investigate the role of the near background in aesthetic appraisal of short-range viewsheds; and (2) to incorporate the aesthetic valuation of water bodies in the foreground and background into the developed GIS-based algorithm of aesthetic appraisal of short-range watersheds.

Conflicts of Interest:
The authors declare no conflict of interest.