Associations between forest characteristics and socio-economic development: a case study from Portugal

The integration of socio-economic and environmental objectives is a major challenge in developing strategies for sustainable landscapes. We investigated associations between socio-economic variables, landscape metrics and measures of forest condition in the context of Portugal. The main goals of the study were to 1) investigate relationships between forest conditions and measures of socio-economic development at national and regional scales, 2) test the hypothesis that a systematic variation in forest landscape metrics occurs according to the stage of socio-economic development and, 3) assess the extent to which landscape metrics can inform strategies to enhance forest sustainability. A ranking approach and statistical techniques such as Principal Component Analysis were used to achieve these objectives. Relationships between socio-economic characteristics, landscape metrics and measures of forest condition were only significant in the regional analysis of municipalities in Northern Portugal. Landscape metrics for different tree species displayed significant variations across socio-economic groups of municipalities and these differences were consistent with changes in characteristics suggested by the forest transition model. The use of metrics also helped inform place-specific strategies to improve forest management, though it was also apparent that further work was required to better incorporate differences in forest functions into sustainability planning.     

Regionalization of Landscape Pattern Indices Using Multivariate Cluster Analysis

Regionalization, or the grouping of objects in space, is a useful tool for organizing, visualizing, and synthesizing the information contained in multivariate spatial data. Landscape pattern indices can be used to quantify the spatial pattern (composition and configuration) of land cover features. Observable patterns can be linked to underlying processes affecting the generation of landscape patterns (e.g., forest harvesting). The objective of this research is to develop an approach for investigating the spatial distribution of forest pattern across a study area where forest harvesting, other anthropogenic activities, and topography, are all influencing forest pattern. We generate spatial pattern regions (SPR) that describe forest pattern with a regionalization approach. Analysis is performed using a 2006 land cover dataset covering the Prince George and Quesnel Forest Districts, 5.5 million ha of primarily forested land base situated within the interior plateau of British Columbia, Canada. Multivariate cluster analysis (with the CLARA algorithm) is used to group landscape objects containing forest pattern information into SPR. Of the six generated SPR, the second cluster (SPR2) is the most prevalent covering 22% of the study area. On average, landscapes in SPR2 are comprised of 55.5% forest cover, and contain the highest number of patches, and forest/non-forest joins, indicating highly fragmented landscapes. Regionalization of landscape pattern metrics provides a useful approach for examining the spatial distribution of forest pattern. Where forest patterns are associated with positive or negative environmental conditions, SPR can be used to identify similar regions for conservation or management activities.     

Examining Landscape Dynamics at a Watershed Scale Using Landsat TM Imagery for Detection of Wintering Hooded Crane Decline in Yashiro,Japan

It is usually inappropriate to define rectangular land areas or administrative units as the extent for quantifying landscapes that possess hierarchical structure. As a functional unit established by geophysical relationships, the watershed is one of many natural scales in the hierarchical landscape. We examined the dynamics of the Yashiro watershed of Japan at the landscape level using pattern metrics based on Landsat thematic mapper (TM) imagery from 1985 to 1998. This watershed provides important habitats for the hooded crane (Grus monachus), a vulnerable species. While its world population has remained stable, the number wintering at Yashiro has declined in recent years. Changes in landscape metrics reveal that the spatial pattern within the watershed underwent homogenization due to depopulation of local people and shifts in local energy requirements and forest management policy at Yashiro. Specific changes include: a decrease in bare land area from 6.2% to 1.0% of the landscape, increased forest cover from 69.2% to 76.1%, reduction in patch number from 1194 to 616 and enlarged mean patch size, and a decrease in total edge from 223,740 m to 158,040 m. The rate of change in landscape metrics indicates a rapid change towards homogeneity in the landscape since 1990. The temporal changes in hooded crane populations corresponded to the changes in landscape. An alternative explanation has been proposed that decline of the species is influenced by landscape dynamics affecting both habitat selection and food resources. Conservation at the watershed scale is suggested to be complementary to the current conservation measures of the species.     

Assessing Impacts of Typhoons and the Chi-Chi Earthquake on Chenyulan Watershed Landscape Pattern in Central Taiwan Using Landscape Metrics

The Chi-Chi earthquake (M_L = 7.3) occurred in the central part of Taiwan on September 21, 1999. After the earthquake, typhoons Xangsane and Toraji produced heavy rainfall that fell across the eastern and central parts of Taiwan on November 2000 and July 2001. This study uses remote sensing data, landscape metrics, multivariate statistical analysis, and spatial autocorrelation to assess how earthquake and typhoons affect landscape patterns. It addresses variations of the Chenyulan watershed in Nantou County, near the earthquake’s epicenter and crossed by Typhoon Toraji. The subsequent disturbances have gradually changed landscape of the Chenyulan watershed. Disturbances of various types, sizes, and intensities, following various tracks, have various effects on the landscape patterns and variations of the Chenyulan watershed. The landscape metrics that are obtained by multivariate statistical analyses showed that the disturbances produced variously fragmented patches, interspersed with other patches and isolated from patches of the same type across the entire Chenyulan watershed. The disturbances also affected the isolation, size, and shape-complexity of patches at the landscape and class levels. The disturbances at the class level more strongly affected spatial variations in the landscape as well as patterns of grasslands and bare land, than variations in the watershed farmland and forest. Moreover, the earthquake with high magnitude was a starter to create these landscape variations in space in the Chenyulan watershed. The cumulative impacts of the disturbances on the watershed landscape pattern had existed, especially landslides and grassland in the study area, but were not always evident in space and time in landscape and other class levels.     

URBAN IMPACTS ON PHYSICAL STREAM CONDITION: EFFECTS OF SPATIAL SCALE,CONNECTIVITY, AND LONGITUDINAL TRENDS1

ABSTRACT: An assessment of physical conditions in urban streams of the Puget Sound region, coupled with spatially explicit watershed characterizations, demonstrates the importance of spatial scale, drainage network connectivity, and longitudinal downstream trends when considering the effects of urbanization on streams. A rapid stream assessment technique and a multimetric index were used to describe the physical conditions of multiple reaches in four watersheds. Watersheds were characterized using geographic information system (GIS) derived landscape metrics that represent the magnitude of urbanization at three spatial scales and the connectivity of urban land. Physical conditions, as measured by the physical stream conditions index (PSCI), were best explained for the watersheds by two landscape metrics: quantity of intense and grassy urban land in the subwatershed and quantity of intense and grassy urban land within 500 m of the site (R²= 0.52, p > 0.0005). A multiple regression of PSCI with these metrics and an additional connectivity metric (proximity of a road crossing) provided the best model for the three urban watersheds (R²= 0.41, p > 0.0005). Analyses of longitudinal trends in PSCI within the three urban watersheds showed that conditions improved when a stream flowed through an intact riparian buffer with forest or wetland vegetation and without road crossings. Results demonstrate that information on spatial scale and patterns of urbanization is essential to understanding and successfully managing urban streams.     

Impacts of Transportation Routes on Landscape Diversity: A Comparison of Different Route Types and Their Combined Effects

A comparison of different transportation route types and their combined effects on landscape diversity was conducted within Tiaoxi watershed (China) between 1994 and 2005. Buffer analysis and Mann–Kendall’s test were used to quantify the relationships between distance from transportation routes (railway, highway, national, and provincial road) and a family of landscape diversity parameters (Simpson’s diversity index, Simpson’s evenness index, Shannon’s diversity index, and Shannon’s evenness index). One-way ANOVA was further applied to compare influences from different route types and their combined effects. Five other landscape metrics (patch density, edge density, area-weighted mean shape index, connectance index, and Euclidean nearest neighbor distance) were also calculated to analyze the associations between landscape diversity and landscape pattern characteristics. Results showed that transportation routes exerted significant impacts on landscape diversity. Impact from railway was comparable to that from highway and national road but was more significant than that from provincial road. The spatial influential range of railway and national road was wider than that of highway and provincial road. Combined effects of routes were nonlinear, and impacts from different route types were more complex than those from the same type. The four landscape diversity metrics were comparably effective at the buffer zone scale. In addition, landscape diversity can be alternatively used to indicate fragmentation, connectivity, and isolation at route buffer scale. This study demonstrates an applicable approach to quantitatively characterize the impacts from transportation routes on landscape patterns and has potential to facilitate route network planning.     

Effects of Impervious Cover at Multiple Spatial Scales on Coastal Watershed Streams1

Abstract: The spatial scale and location of land whose development has the strongest influence on aquatic ecosystems must be known to support land use decisions that protect water resources in urbanizing watersheds. We explored impacts of urbanization on streams in the West River watershed, New Haven, Connecticut, to identify the spatial scale of watershed imperviousness that was most strongly related to water chemistry, macroinvertebrates, and physical habitat. A multiparameter water quality index was used to characterize regional urban nonpoint source pollution levels. We identified a critical level of 5% impervious cover, above which stream health declined. Conditions declined with increasing imperviousness and leveled off in a constant state of impairment at 10%. Instream variables were most correlated (0.77 ≤ |r| ≤ 0.92, p < 0.0125) to total impervious area (TIA) in the 100‐m buffer of local contributing areas (～5‐km² drainage area immediately upstream of each study site). Water and habitat quality had a relatively consistent strong relationship with TIA across each of the spatial scales of investigation, whereas macroinvertebrate metrics produced noticeably weaker relationships at the larger scales. Our findings illustrate the need for multiscale watershed management of aquatic ecosystems in small streams flowing through the spatial hierarchies that comprise watersheds with forest‐urban land use gradients.     

Scale Effects on Spatially Varying Relationships Between Urban Landscape Patterns and Water Quality

Scientific interpretation of the relationships between urban landscape patterns and water quality is important for sustainable urban planning and watershed environmental protection. This study applied the ordinary least squares regression model and the geographically weighted regression model to examine the spatially varying relationships between 12 explanatory variables (including three topographical factors, four land use parameters, and five landscape metrics) and 15 water quality indicators in watersheds of Yundang Lake, Maluan Bay, and Xinglin Bay with varying levels of urbanization in Xiamen City, China. A local and global investigation was carried out at the watershed-level, with 50 and 200 m riparian buffer scales. This study found that topographical features and landscape metrics are the dominant factors of water quality, while land uses are too weak to be considered as a strong influential factor on water quality. Such statistical results may be related with the characteristics of land use compositions in our study area. Water quality variations in the 50 m buffer were dominated by topographical variables. The impact of landscape metrics on water quality gradually strengthen with expanding buffer zones. The strongest relationships are obtained in entire watersheds, rather than in 50 and 200 m buffer zones. Spatially varying relationships and effective buffer zones were verified in this study. Spatially varying relationships between explanatory variables and water quality parameters are more diversified and complex in less urbanized areas than in highly urbanized areas. This study hypothesizes that all these varying relationships may be attributed to the heterogeneity of landscape patterns in different urban regions. Adjustment of landscape patterns in an entire watershed should be the key measure to successfully improving urban lake water quality.     

Landscape Trends in Mid-Atlantic and Southeastern United States Ecoregions

Landscape pattern and composition metrics are potential indicators for broad-scale monitoring of change and for relating change to human and ecological processes. We used a probability sample of 20-km × 20-km sampling blocks to characterize landscape composition and pattern in five US ecoregions: the Middle Atlantic Coastal Plain, Southeastern Plains, Northern Piedmont, Piedmont, and Blue Ridge Mountains. Land use/land cover (LULC) data for five dates between 1972 and 2000 were obtained for each sample block. Analyses focused on quantifying trends in selected landscape pattern metrics by ecoregion and comparing trends in land cover proportions and pattern metrics among ecoregions. Repeated measures analysis of the landscape pattern documented a statistically significant trend in all five ecoregions towards a more fine-grained landscape from the early 1970s through 2000. The ecologically important forest cover class also became more fine-grained with time (i.e., more numerous and smaller forest patches). Trends in LULC, forest edge, and forest percent like adjacencies differed among ecoregions. These results suggest that ecoregions provide a geographically coherent way to regionalize the story of national land use and land cover change in the United States. This study provides new information on LULC change in the southeast United States. Previous studies of the region from the 1930s to the 1980s showed a decrease in landscape fragmentation and an increase in percent forest, while this study showed an increase in forest fragmentation and a loss of forest cover.     

A Multi-scale Spatial Analysis of Native and Exotic Plant Species Richness Within a Mixed-Disturbance Oak Savanna Landscape

Impacts of human land use pose an increasing threat to global biodiversity. Resource managers must respond rapidly to this threat by assessing existing natural areas and prioritizing conservation actions across multiple spatial scales. Plant species richness is a useful measure of biodiversity but typically can only be evaluated on small portions of a given landscape. Modeling relationships between spatial heterogeneity and species richness may allow conservation planners to make predictions of species richness patterns within unsampled areas. We utilized a combination of field data, remotely sensed data, and landscape pattern metrics to develop models of native and exotic plant species richness at two spatial extents (60- and 120-m windows) and at four ecological levels for northwestern Ohio’s Oak Openings region. Multiple regression models explained 37–77 % of the variation in plant species richness. These models consistently explained more variation in exotic richness than in native richness. Exotic richness was better explained at the 120-m extent while native richness was better explained at the 60-m extent. Land cover composition of the surrounding landscape was an important component of all models. We found that percentage of human-modified land cover (negatively correlated with native richness and positively correlated with exotic richness) was a particularly useful predictor of plant species richness and that human-caused disturbances exert a strong influence on species richness patterns within a mixed-disturbance oak savanna landscape. Our results emphasize the importance of using a multi-scale approach to examine the complex relationships between spatial heterogeneity and plant species richness.     

A GIS Approach to Prioritizing Habitat for Restoration Using Neotropical Migrant Songbird Criteria

Restoration efforts to increase wildlife habitat quality in agricultural landscapes have limited funding and are typically done on a first come, first serve basis. In order to increase the efficiency of these restoration efforts, a prioritized ranking system is needed to obtain the greatest increase in habitat quality possible for the fewest amount of hectares restored. This project examines the use of a GIS based multi-criteria approach to prioritize lands for reforestation along the Kaskaskia River in Illinois. Loss of forested area and corresponding increase in forest fragmentation has decreased songbird habitat quality across the Midwestern United States. We prioritized areas for reforestation based on nine landscape metrics: available agricultural land, forest cover gaps, edge density, proximity to river, 200 m corridor area, total forest core area, fringe core area, distance to primary core value, and primary core area. The multi-criteria analysis revealed that high priority areas for reforestation were most likely to be close to the riparian corridor and existing large blocks of forest. Analysis of simulated reforestation (0, 0.5, 1.0, 5.0 10.0, 25.0, and 50.0% of highest priority parcels reforested) revealed different responses for multiple landscape metrics used to quantify forest fragmentation following reforestation, but indicated that the study area would get the greatest rate of return on reforestation efforts by reforesting 10.0% of the highest priority areas. This project demonstrates how GIS and a multi-criteria analysis approach can be used to increase the efficiency of restoration projects. This approach should be considered by land managers when attempting to identify the location and quantity of area for restoration within a landscape.     

Common Core Themes in Geomorphic, Ecological, and Social Systems

Core themes of geomorphology include: open systems and connectivity; feedbacks and complexity; spatial differentiation of dominant physical processes within a landscape; and legacy effects of historical human use of resources. Core themes of ecology include: open systems and connectivity; hierarchical, heterogeneous, dynamic, and context-dependent characteristics of ecological patterns and processes; nonlinearity, thresholds, hysteresis, and resilience within ecosystems; and human effects. Core themes of environmental governance include: architecture of institutions and decision-making; agency, or ability of actors to prescribe behavior of people in relation to the environment; adaptiveness of social groups to environmental change; accountability and legitimacy of systems of governance; allocation of and access to resources; and thresholds and feedback loops within environmental policy. Core themes common to these disciplines include connectivity, feedbacks, tipping points or thresholds, and resiliency. Emphasizing these points of disciplinary overlap can facilitate interdisciplinary understanding of complex systems, as well as more effective management of landscapes and ecosystems by highlighting drivers of change within systems. We use a previously published conceptual framework to examine how these core themes can be integrated into interdisciplinary research for human–landscape systems via the example of a river.     

An Ecoregional Context for Forest Management on National Wildlife Refuges of the Upper Midwest,USA

To facilitate forest planning and management on National Wildlife Refuges, we synthesized multiple data sources to describe land ownership patterns, land cover, landscape pattern, and changes in forest composition for four ecoregions and their associated refuges of the Upper Midwest. We related observed patterns to ecological processes important for forest conservation and restoration, with specific attention to refuge patterns of importance for forest landbirds of conservation priority. The large amount of public land within the ecoregions (31–80%) suggests that opportunities exist for coarse and meso-scale approaches to conserving and restoring ecological processes affecting the refuges, particularly historical fire regimes. Forests dominate both ecoregions and refuges, but refuge forest patches are generally larger and more aggregated than in associated ecoregions. Broadleaf taxa have increased in dominance in the ecoregions and displaced fire-dependent taxa such as pine (Pinus spp.) and other coniferous species; these changes in forest composition have likely also affected refuge forests. Despite compositional changes, larger forest patches on refuges suggests that they may provide better habitat for area-sensitive forest landbirds of mature, compositionally diverse forests than surrounding lands if management continues to promote increased patch size. We reason that although fine-scale research and monitoring for species of conservation priority is important, broad scale (ecoregional) assessments provide crucial context for effective forest and wildlife management in protected areas.     

Landscape Characterization Integrating Expert and Local Spatial Knowledge of Land and Forest Resources

In many developing countries, political documentation acknowledges the crucial elements of participation and spatiality for effective land use planning. However, operative approaches to spatial data inclusion and representation in participatory land management are often lacking. In this paper, we apply and develop an integrated landscape characterization approach to enhance spatial knowledge generation about the complex human–nature interactions in landscapes in the context of Zanzibar, Tanzania. We apply an integrated landscape conceptualization as a theoretical framework where the expert and local knowledge can meet in spatial context. The characterization is based on combining multiple data sources in GIS, and involves local communities and their local spatial knowledge since the beginning into the process. Focusing on the expected information needs for community forest management, our characterization integrates physical landscape features and retrospective landscape change data with place-specific community knowledge collected through participatory GIS techniques. The characterization is established in a map form consisting of four themes and their synthesis. The characterization maps are designed to support intuitive interpretation, express the inherently uncertain nature of the data, and accompanied by photographs to enhance communication. Visual interpretation of the characterization mediates information about the character of areas and places in the studied local landscape, depicting the role of forest resources as part of the landscape entity. We conclude that landscape characterization applied in GIS is a highly potential tool for participatory land and resource management, where spatial argumentation, stakeholder communication, and empowerment are critical issues.     

Watershed vulnerability predictions for the Ozarks using landscape models

Forty-six broad-scale landscape metrics derived from commonly used landscape metrics were used to develop potential indicators of total phosphorus (TP) concentration, total ammonia (TA) concentration, and Escherichia coli bacteria count among 244 sub-watersheds of the Upper White River (Ozark Mountains, USA). Indicator models were developed by correlating field-based water quality measurements and contemporaneous remote-sensing-based ecological metrics using partial least squares (PLS) analyses. The TP PLS model resulted in one significant factor explaining 91% of the variability in surface water TP concentrations. Among the 18 contributing landscape model variables for the TP PLS model, the proportions of a sub-watershed that are barren and in human use were key indicators of water chemistry in the associated sub-watersheds. The increased presence and reduced fragmentation of forested areas are negatively correlated with TP concentrations in associated sub-watersheds, particularly within close proximity to rivers and streams. The TA PLS model resulted in one significant factor explaining 93% of the variability in surface water TA concentrations. The eight contributing landscape model variables for the TA PLS model were among the same forest and urban metrics for the TP model, with a similar spatial gradient trend in relationship to distance from streams and rivers within a sub-watershed. The E. coli PLS model resulted in two significant factors explaining 99.7% of the variability in E. coli cell count. The 17 contributing landscape model variables for the E. coli PLS model were similar to the TP and TA models. The integration of model results demonstrates that forest, riparian, and urban attributes of sub-watersheds affect all three models. The results provide watershed managers in the Ozark Mountains with a broad-scale vulnerability prediction tool, focusing on TP, TA, and E. coli, and are being used to prioritize and evaluate monitoring and restoration efforts in the vicinity of the White River, a major tributary to the Mississippi River and Gulf of Mexico.     

Planning methodology for predicting spatial patterns of risk potential from industrial land uses

The nature of environmental risk is often determined by the location patterns of industrial firms in a locale. We develop a methodology for analyzing the capacity of toxics-generating industries to leave toxic residuals on the landscape, in the context of long-term master planning. To understand an area's risk profile, we first characterize the location pattern of risk-generating firms and develop ways to represent the risk potentials of these shapes. To this end, we develop a geometric approach for characterizing the spatial patterns of these clusters of dirty industries, using new measures for analyzing spatial densities and compactness. We then introduce context (i.e. actual zoning patterns, transportation corridors) and explain how this relates to the spatial patterns found. We illustrate this analytical method with application to Orange County, California, USA and point out how it affords a deeper understanding of the connections between industry and environmental risk. We end the article with a discussion of how these analytical methods might be used for land use planning.     

Variability of Lotic Macroinvertebrate Assemblages and Stream Habitat Characteristics Across Hierarchical Landscape Classifications

Streams are naturally hierarchical systems, and their biota are affected by factors effective at regional to local scales. However, there have been only a few attempts to quantify variation in ecological attributes across multiple spatial scales. We examined the variation in several macroinvertebrate metrics and environmental variables at three hierarchical scales (ecoregions, drainage systems, streams) in boreal headwater streams. In nested analyses of variance, significant spatial variability was observed for most of the macroinvertebrate metrics and environmental variables examined. For most metrics, ecoregions explained more variation than did drainage systems. There was, however, much variation attributable to residuals, suggesting high among-stream variation in macroinvertebrate assemblage characteristics. Nonmetric multidimensional scaling (NMDS) and multiresponse permutation procedure (MRPP) showed that assemblage composition differed significantly among both drainage systems and ecoregions. The associated R-statistics were, however, very low, indicating wide variation among sites within the defined landscape classifications. Regional delineations explained most of the variation in stream water chemistry, ecoregions being clearly more influential than drainage systems. For physical habitat characteristics, by contrast, the among-stream component was the major source of variation. Distinct differences attributable to stream size were observed for several metrics, especially total number of taxa and abundance of algae-scraping invertebrates. Although ecoregions clearly account for a considerable amount of variation in macroinvertebrate assemblage characteristics, we suggest that a three-tiered classification system (stratification through ecoregion and habitat type, followed by assemblage prediction within these ecologically meaningful units) will be needed for effective bioassessment of boreal running waters.     

Spatial modeling of harvest constraints on wood supply versus wildlife habitat objectives

We studied the effects of spatial and temporal timber harvesting constraints on competing objectives of sustaining wildlife habitat supply and meeting timber harvest objectives in a boreal mixedwood forest. A hierarchical modeling approach was taken, where strategic and tactical level models were used to project blocking and scheduling of harvest blocks. Harvest block size and proximity, together with short- and long-term temporal constraints, were adjusted in a factorial manner to allow creation of response-surface models. A new measure of the habitat mosaic was defined to describe the emergent pattern of habitat across the landscape. These models, together with multiple linear regression, were used to provide insight on convergence or divergence between spatial objectives. For example, green-up delay (defined as time required before a harvest block adjacent to a previously logged block can be scheduled for harvest) had an adverse effect on the amount of annual harvest area that could be allocated and blocked spatially, and habitat supply responded in an opposite direction to that of wood supply, where caribou, moose wintering, and marten habitat supply increased when harvest blocks were further apart, maximum block size smaller, and both a green-up delay and mesoscale stratification were applied. Although there was no "solution space" free of conflicts, the analysis suggests that application of the mesoscale stratification, together with a diversity of harvest block sizes and a between-harvest block proximity of 250 m, will perform relatively well with respect to wood supply objectives, and at the same time create a less fragmented landscape that better reflects natural forest patterns.     

豫西黄河流域土地利用景观空间格局分析

运用景观生态学原理,借助GIS技术,选取景观指数对豫西黄河流域土地利用景观格局进行分析。结果表明,林地和以旱地为主的农业景观是豫西黄河流域的主体景观类型;各景观类型间斑块密度值差距较大,斑块形状指数差异较明显,总体上较简单;研究区域景观异质性和多样性较高,各景观类型所占比例差异较大,农田和林地占有较大优势。该区域景观破碎化程度较小,特别是西部地区地形起伏较大,人类干扰活动相对较小。     

Integrating Landscape Analysis and Planning: A Multi-Scale Approach for Oriented Management of Tourist Recreation

Tourism and landscape are interdependent concepts. Nature- and culture-based tourism are now quite well developed activities and can constitute an excellent way of exploiting the natural resources of certain areas, and should therefore be considered as key objectives in landscape planning and management in a growing number of countries. All of this calls for careful evaluation of the effects of tourism on the territory. This article focuses on an integrated spatial method for landscape analysis aimed at quantifying the relationship between preferences of visitors and landscape features. The spatial expression of the model relating types of leisure and recreational preferences to the potential capacity of the landscape to meet them involves a set of maps showing degrees of potential visitor satisfaction. The method constitutes a useful tool for the design of tourism planning and management strategies, with landscape conservation as a reference.