Exploring Switzerland’s Land Cover Change Dynamics Using a National Statistical Survey

: Timely and reliable Land Use and Cover change information is crucial to efﬁciently mitigate the negative impact of environmental changes. Switzerland has the ambitious objective of being a sustainable country while remaining an attractive business location with a high level of well-being. However, this aspiration is hampered by increasing pressures that are signiﬁcantly impacting the environment and putting serious demands on land. In the present study, we used the national Land Cover (LC) dataset, named ArealStatistik , produced by the Federal Statistical Ofﬁce, to explore the spatiotemporal patterns of Land Cover in Switzerland, providing a comprehensive assessment of land cover change at the national scale. Results indicate that, in general, Switzerland has undergone small, spatially dispersed, dynamic, and gradual change trends, with high rates of transition between low growing Brush Vegetation and forest LC classes in recent years. These pixel-level trends are more important in the lower altitude plateau and Jura regions, while greater changes in the spatial conﬁguration of LC are observed in the alpine regions. However, ﬁndings also suggest that identifying drivers and understanding the rate of change are limited by the spatial resolution and temporal update frequency of the ArealStatistik . The ability to understand these drivers would beneﬁt from a high-resolution annual LC dataset. Such a data product can be produced using the ArealStatistik together with dense satellite data time-series and Machine/Deep Learning techniques.


Introduction
The world is significantly impacted by human-induced environmental changes from local to global scales. Human well-being is critically dependent on Earth's natural systems. Despite the agreed goal of halting land degradation, it is projected to increase in the twenty-first century under all development scenarios [1]. Consequently, land degradation is recognized as a major global environmental issue with important linkages to climate change or biodiversity loss [2]. The drivers of land degradation are predominantly local, so actions to address them should be based on the understanding of the local interplay of a range of factors.
To understand these drivers, it is crucial to apprehend how people utilize the land for socio-economic activities and to precisely map the changes of observed (bio)physical cover and use of the Earth's surface (e.g., Land Use and Cover-LUC) [3]. Land cover (LC) is included in the list of Essential Climate Variables (ECV), on account of its relevance for assessment of changes in the climate system and its availability in the form of multiple coherent global datasets [4]. LC has been further identified as a possible Essential Biodiversity Variable [5], recognizing the potential for its detrimental impact on biodiversity of land conversion processes [6], and contributes to reporting on 14 of the 17 Sustainable Development Goals (SDGs) [7]. Moreover, LUC has been increasingly recognized as a Land 2023, 12, 1386 3 of 20 essential for human well-being [23,24] but also indicate that most ecosystems are being degraded and fragmented, and therefore not providing efficient and effective services as healthy ecosystems.
To track these changes in Switzerland, the official LUC data source is the Land Use Statistics (Arealstatistik) of the Federal Office for Statistics (FSO) obtained by visually interpreting aerial images and assigning a LC as well as a LU category of the lower-left corner of each sample point from a regular 100 m grid cell, corresponding to more than 4 million points over the country following three nomenclatures: standard (72 categories); land cover (27), and land use (46) over four-time periods (1979-85, 1992-97, 2004-09, 2013-18) [25]. Using this dataset, various studies have well documented Land Use change, such as [26,27]. Most of the national studies concentrate on Land Use change and how it can affect carbon stocks [28], species habitat [29] or urbanization and land abandonment [30]. However, less attention has been provided to Land Cover change and, currently, a comprehensive assessment of Land Cover change across all categories of the Arealstatistik at the national scale is still missing. Indeed, most research concentrates on a specific LC class, such as forest cover change [31] or cropland and grassland [32], producing information from LC data like a nationwide habitat map [33], LC degradation [34] or studying the effects on Gross Primary Production [35] Based on these previous considerations, this study aimed to explore the spatiotemporal patterns of Land Cover change in Switzerland using the official LC data for Switzerland, therefore filling the identified gap by providing a first extensive assessment. To achieve this objective, the following questions were investigated: (1) what are the spatiotemporal patterns of major LC change in the country between 1985 and 2018; (2) what are the key LC transitions; (3) what is the intensity of change?

Materials and Methods
Hereafter, we present the study area, the data, and the different analytical steps to determine Land Cover change spatiotemporal patterns (Figure 1). The overall objective is to produce a comprehensive assessment of Land Cover Change (LCC) at the national and regional scales over the entirety of Switzerland, and distinguish between random and systematic changes, as well as quantifying change in spatial patterns.
Land 2023, 12, x FOR PEER REVIEW 3 of 22 regional scale (European Green Deal), and through international policies (Sustainable Development Goals). They all agree that land can offer multiple valuable ecosystem services that are essential for human well-being [23,24] but also indicate that most ecosystems are being degraded and fragmented, and therefore not providing efficient and effective services as healthy ecosystems.
To track these changes in Switzerland, the official LUC data source is the Land Use Statistics (Arealstatistik) of the Federal Office for Statistics (FSO) obtained by visually interpreting aerial images and assigning a LC as well as a LU category of the lower-left corner of each sample point from a regular 100 m grid cell, corresponding to more than 4 million points over the country following three nomenclatures: standard (72 categories); land cover (27), and land use (46) over four-time periods (1979-85, 1992-97, 2004-09, 2013-18) [25]. Using this dataset, various studies have well documented Land Use change, such as [26,27]. Most of the national studies concentrate on Land Use change and how it can affect carbon stocks [28], species habitat [29] or urbanization and land abandonment [30]. However, less attention has been provided to Land Cover change and, currently, a comprehensive assessment of Land Cover change across all categories of the Arealstatistik at the national scale is still missing. Indeed, most research concentrates on a specific LC class, such as forest cover change [31] or cropland and grassland [32], producing information from LC data like a nationwide habitat map [33], LC degradation [34] or studying the effects on Gross Primary Production [35] Based on these previous considerations, this study aimed to explore the spatiotemporal patterns of Land Cover change in Switzerland using the official LC data for Switzerland, therefore filling the identified gap by providing a first extensive assessment. To achieve this objective, the following questions were investigated: (1) what are the spatiotemporal patterns of major LC change in the country between 1985 and 2018; (2) what are the key LC transitions; (3) what is the intensity of change?

Materials and Methods
Hereafter, we present the study area, the data, and the different analytical steps to determine Land Cover change spatiotemporal patterns (Figure 1). The overall objective is to produce a comprehensive assessment of Land Cover Change (LCC) at the national and regional scales over the entirety of Switzerland, and distinguish between random and systematic changes, as well as quantifying change in spatial patterns.    (Figure 2), which are based on zones of homogeneity of flora and fauna, as well as hydrologic and topographic features present [36,37].

Study Area
The study area covers all of Switzerland between the latitudes of 45.82° and 47.81° and the longitudes of 5.96° and 10.49°. The total land area is 41,285 km 2 . Switzerland can be divided into 6 biogeographical regions (Figure 2), which are based on zones of homogeneity of flora and fauna, as well as hydrologic and topographic features present [36,37]. These regions represent a diversity of landscapes and ecological features characterized by different climatic, geological, and vegetation conditions. Based on the altitude, these regions can be grouped into three major areas: the Alps (Northern Prealps, Eastern and Western Central, Southern), which is characterized by high elevations and accounting for 60% of the total surface of the country, and the Plateau (30%) and the Jura (10%), which are at a low elevation.
The climate of Switzerland is strongly affected by the Alps and its proximity to the Atlantic Ocean. The south side of the Alps is mostly influenced by the Mediterranean Sea. It can be defined as moderately continental in the plateau, alpine in the mountain regions, and more temperate in the Southern Alps. Switzerland has four distinct seasons with varying temperatures and precipitation patterns following these different climate conditions [38].

Data
In this study, we used the official source of LUC information, produced and made freely available by the Federal Office of Statistics (FSO). The ArealStatistik dataset is obtained by the visual interpretation of aerial images. Operators are assigned a LC and LU category for each sample point, following three nomenclatures: standard (72 categories); land cover (27), and land use (46), over four survey periods, each covering a 6-year timeframe (1979-1985, 1992-1997, 2004-2009, 2013-2018) [25,26]. For simplicity, the datasets are henceforth referred to by the last survey year (1985, 1997, 2009, and 2018). Gridded LC maps at a resolution of 100 m were obtained by rasterizing the sampling points of the ArealStatistik dataset. Scheduling of the aerial survey ensures that each point is sampled at the same time during the survey period, so that each data set can be considered a single LC map. Any change in LC is considered to have taken place at the same position These regions represent a diversity of landscapes and ecological features characterized by different climatic, geological, and vegetation conditions. Based on the altitude, these regions can be grouped into three major areas: the Alps (Northern Prealps, Eastern and Western Central, Southern), which is characterized by high elevations and accounting for 60% of the total surface of the country, and the Plateau (30%) and the Jura (10%), which are at a low elevation.
The climate of Switzerland is strongly affected by the Alps and its proximity to the Atlantic Ocean. The south side of the Alps is mostly influenced by the Mediterranean Sea. It can be defined as moderately continental in the plateau, alpine in the mountain regions, and more temperate in the Southern Alps. Switzerland has four distinct seasons with varying temperatures and precipitation patterns following these different climate conditions [38].

Data
In this study, we used the official source of LUC information, produced and made freely available by the Federal Office of Statistics (FSO). The ArealStatistik dataset is obtained by the visual interpretation of aerial images. Operators are assigned a LC and LU category for each sample point, following three nomenclatures: standard (72 categories); land cover (27), and land use (46), over four survey periods, each covering a 6-year timeframe (1979-1985, 1992-1997, 2004-2009, 2013-2018) [25,26]. For simplicity, the datasets are henceforth referred to by the last survey year (1985, 1997, 2009, and 2018). Gridded LC maps at a resolution of 100 m were obtained by rasterizing the sampling points of the ArealStatistik dataset. Scheduling of the aerial survey ensures that each point is sampled at the same time during the survey period, so that each data set can be considered a single LC map. Any change in LC is considered to have taken place at the same position in the interim period between surveys, following the methodology of FOEN (2021). Within the nomenclature used, Land Cover is classified into 27 'Basic Categories'. For better statistical reliability, especially for the use of this dataset on a small scale, these categories are also aggregated into 6 Principal Domains ( Figure 3). Details of the classifications are included in Table 1.

Change Detection
LCC is commonly performed using either pre-classification change detection, in which changes in the original remote sensing images are analyzed, or post-classification change detection, which compares the resulting classified images [41]. Usually, post-classification leads to a more accurate assessment of change and provides greater information as to the processes leading to change [42]. We applied this approach for change detection. Initial analysis of LCC involved the creation of a cross-tabulation matrix that represents the overlay of two spatial datasets from different times. The sum of horizontal entries represents the size of classes at t 0 , and the sum of vertical entries represents the size of classes at time t 1 . Values on the diagonal show the size of persistence of a class between the two time periods, whereas all other values indicate gain or loss. Visual representations of these matrices were produced in the form of Sankey diagrams to communicate the broad-scale LC transitions across multiple time periods and many thematic classes [43]. Land Cover data were downloaded directly from the FSO website as a Comma-Separated Value (CSV) file and then rasterized at 100 m resolution in QGIS 3.22. In addition, national border and biogeographical regions were obtained from the Federal Office for Topography (swisstopo) in the swissBOUNDARIES3D landscape model as shapefiles. All subsequent analyses were performed using R Version 4.2.1 [39] and Python version 3.10.7 [40].

Change Detection
LCC is commonly performed using either pre-classification change detection, in which changes in the original remote sensing images are analyzed, or post-classification change detection, which compares the resulting classified images [41]. Usually, post-classification leads to a more accurate assessment of change and provides greater information as to the processes leading to change [42]. We applied this approach for change detection.
Initial analysis of LCC involved the creation of a cross-tabulation matrix that represents the overlay of two spatial datasets from different times. The sum of horizontal entries represents the size of classes at t0, and the sum of vertical entries represents the size of classes at time t1. Values on the diagonal show the size of persistence of a class between the two time periods, whereas all other values indicate gain or loss. Visual representations of these matrices were produced in the form of Sankey diagrams to communicate the broad-scale LC transitions across multiple time periods and many thematic classes [43].

Intensity of LC Change
Quantification of systematic LCC was performed using intensity analysis, as detailed in Aldwaik and Pontius (2012), using the R package intensity.analysis (https://cran.r-project.org/web/packages/intensity.analysis/index.html, accessed on 10 January 2023). LCC

Intensity of LC Change
Quantification of systematic LCC was performed using intensity analysis, as detailed in Aldwaik and Pontius (2012), using the R package intensity.analysis (https://cran.rproject.org/web/packages/intensity.analysis/index.html, accessed on 10 January 2023). LCC was assessed on three additive levels (between time intervals, within categories, and transitions between categories), against the null hypothesis of spatiotemporally uniform change. The use of this methodology means that the focus is on systematic and stable transitions and spurious or ephemeral changes within single time-intervals are given less importance [44]. To simplify visualization and analysis, category level changes under 1% are considered inconsequential and are not included in the results.
The annual rate of change of the study region in the interval between datasets at different times is calculated as: This is then compared to the uniform rate of change for the study region, which represents the expected value of change if the annual rate of change remained constant throughout the study period (33 years). At the category level, annual rate of gains and losses can be estimated.
The results of Equations (3) and (4) can then also be compared to the result of Equation (2) for the interval in question to distinguish classes that undergo significant change from those which change under random processes, as would be expected under the normal evolution of the landscape. This avoids an overemphasis of the salt-and-pepper effect often observed with LUC, which may distract from larger-scale processes in action. As such, categories where annual change intensity during a time interval is greater than the uniform intensity for that time interval undergo active change. Categories where annual change intensity in a time interval is less than the uniform intensity for that time interval is considered dormant.
Finally, the intensity of transitions between categories is calculated to provide insight into which categories are targeted and which are avoided when a category gains or loses total surface area. The intensity of transitions to category n from category i in each time interval is calculated as: The uniform intensity of the transition to category n, i.e., the expected rate of transition if category n gained from all other categories in a uniform manner, is calculated as: The intensity of transitions R tin is calculated independently for gain and loss and are compared to uniform intensity W tin to distinguish classes that are avoided or targeted by a certain class when it changes.

Assessment of Change in Spatial Pattern of Land Cover
Of additional interest to the degree of change and intensity of transitions between categories is the change in the spatial pattern of land cover. This can provide insight into processes such as landscape fragmentation, which influence both soil carbon stocks and aboveground biomass [20], as well as provide an idea of the key LCC processes at a local scale. Pattern-based spatial analysis was performed using the R package motif [45] and applying the Jensen-Shannon Divergence (JSD) as a quantitative measurement of pattern similarity between gridded data from two-time intervals. Measurement of pattern similarity incorporates both the change in land cover composition, the proportion of different classes represented, and in land cover configuration; the layout of these classes are within a certain area [46]. The JSD provides a normalized measure of dissimilarity between the co-occurrence histograms of two maps with a value of 0 indicating that the two maps have identical co-occurrences, and a value of 1 indicating an entirely different distribution of LC classes. JSD was calculated on tiles of the Land Cover data created for windows of 50 pixels by 50 pixels, or 25 km 2 , to provide an extent that permits local comparison while capturing the granularity of the LC classes. Based on the work of [46], the similarity threshold of 0.012 was applied, with values lower than this indicating negligible change, to exclude landscapes considered to be 'unchanged'.  (Table S1).

Change Intensity over Time
At all levels of LC nomenclature, and for all time intervals considered, the annual intensity of change is greater than the uniform intensity for the study area, indicating that the changes occurring can be considered as relatively fast [44]. The annual rate of LC change is higher in the period from 2009 to 2018, and this pattern can be observed across all regions (Figures 5 and 6). Basic Categories are more dynamic than the Principal Domains, with a greater percentage of the map area undergoing change due to subtle changes within the broader classifications that are more easily triggered. Land Cover is, in general, more dynamic than Land Use on both the national and regional scale. The Western Central Alps are an exception to this, where the annual percentage change in Land Use was greater than for Land Cover in the first two intervals, driven by gains in non-agricultural forests and losses in unproductive areas and alpine grazing areas.  Figure S1).

Change Intensity over Time
At all levels of LC nomenclature, and for all time intervals considered, the annual intensity of change is greater than the uniform intensity for the study area, indicating that the changes occurring can be considered as relatively fast [44]. The annual rate of LC change is higher in the period from 2009 to 2018, and this pattern can be observed across all regions (Figures 5 and 6). Basic Categories are more dynamic than the Principal Domains, with a greater percentage of the map area undergoing change due to subtle changes within the broader classifications that are more easily triggered. Land Cover is, in general, more dynamic than Land Use on both the national and regional scale. The Western Central Alps are an exception to this, where the annual percentage change in Land Use was greater than for Land Cover in the first two intervals, driven by gains in non-agricultural forests and losses in unproductive areas and alpine grazing areas.

Category Intensity
As a contrast to gains in buildings that are only above uniform intensity in the first interval, annual percentage losses in the minor categories of artificial areas are greatest in

Category Intensity
As a contrast to gains in buildings that are only above uniform intensity in the first interval, annual percentage losses in the minor categories of artificial areas are greatest in

Category Intensity
As a contrast to gains in buildings that are only above uniform intensity in the first interval, annual percentage losses in the minor categories of artificial areas are greatest in 2009-2018, in all regions ( Figure 6). The Jura and Plateau regions show much greater annual change intensities than the alpine regions, driven by changes within vegetated LC classes (Figures 7 and 8).
nual change intensities than the alpine regions, driven by changes within vegetated LC classes (Figures 7 and 8).
Changes in most tree vegetation categories are also greatest in this interval across all regions, but to a greater degree in the Jura and Plateau regions where vegetation changes are generally more dynamic (Figures 7 and 8). Watery Areas undergo more minor changes, with fewer regions showing changes above uniform intensity; however, loss of glaciers and perpetual snow can be observed in the Eastern Central Alps and the Southern Alps. Bare Land categories occupy a much larger total surface area in the alpine regions, and this is reflected in the low intensity of change that remains lower than uniform intensity in all intervals.
Greater changes in marginal forest areas, such as Forest Strips and Forest Edges, can be observed in the alpine regions compared to the Plateau and Jura, and a lower rate of change in open forest reflects the difference in management strategies [31].  Changes in most tree vegetation categories are also greatest in this interval across all regions, but to a greater degree in the Jura and Plateau regions where vegetation changes are generally more dynamic (Figures 7 and 8). Watery Areas undergo more minor changes, with fewer regions showing changes above uniform intensity; however, loss of glaciers and perpetual snow can be observed in the Eastern Central Alps and the Southern Alps. Bare Land categories occupy a much larger total surface area in the alpine regions, and this is reflected in the low intensity of change that remains lower than uniform intensity in all intervals.
Greater changes in marginal forest areas, such as Forest Strips and Forest Edges, can be observed in the alpine regions compared to the Plateau and Jura, and a lower rate of change in open forest reflects the difference in management strategies [31].

Key LC Transitions
A transition between two classes, A and B, can be viewed as a dominant signal of change if class A systematically loses surface to class B, while at the same time class B systematically gains surface from class A [18]. At the national scale, the most important change signals are the gain in Tree Vegetation from Brush Vegetation categories, for which there is a clear increase in the 2009-2018 period compared to other years ( Figure 9). The intensity of gain of Tree Vegetation from Brush Vegetation is not equivalent to the intensity of loss from Brush Vegetation to Tree Vegetation, as these values are calculated based on the initial size of the changing category in question.
The transition from Watery Areas to Bare Land reflects the loss of permanent snow and glacier in alpine areas within the Basic Categories. However, the largest transitions are all weighted in one direction. While the loss of Grass and Herb Vegetation to Artificial Areas and the gain of Artificial Areas from Grass and Herb Vegetation are both systematic transitions above the rate of uniform change, the loss of Grass and Herb Vegetation to Artificial Areas occurs at a much higher intensity. As a result, the signals of the change identified are weak.
The central trends identified for the Tree and Brush Vegetation categories are from Brush classes to Tree classes, with the strongest signals of change in the 2009-2018 period being systematic transitions from Shrubs and Brush Meadows to Closed Forest.
At the lower level of thematic precision, alongside the transition from tree to brush vegetation, another important change is a gain of bare land from watery areas. This is observed particularly in the alpine areas where the Basic Categories indicate this is the loss of glaciers and permanent snow.

Key LC Transitions
A transition between two classes, A and B, can be viewed as a dominant signal of change if class A systematically loses surface to class B, while at the same time class B systematically gains surface from class A [18]. At the national scale, the most important change signals are the gain in Tree Vegetation from Brush Vegetation categories, for which there is a clear increase in the 2009-2018 period compared to other years (Figure 9). The intensity of gain of Tree Vegetation from Brush Vegetation is not equivalent to the intensity of loss from Brush Vegetation to Tree Vegetation, as these values are calculated based on the initial size of the changing category in question.
The transition from Watery Areas to Bare Land reflects the loss of permanent snow At the lower level of thematic precision, alongside the transition from tree to brush vegetation, another important change is a gain of bare land from watery areas. This is observed particularly in the alpine areas where the Basic Categories indicate this is the loss of glaciers and permanent snow.

Change in Spatial Patterns
Over the whole study period, most pixels show a non-negligible change in the spatial pattern. The highest difference in the spatial pattern is observed in the Plateau and northern Prealps regions, notably in the 1997-2009 period ( Figure 10). When grouped by the Basic Categories, all intervals undergo a similar distribution of change, however, more pixels show a JSD value > 0.012, reflecting the higher dynamism of LC change at the higher level of category disaggregation. Amongst the areas of highest distance, the most important changes to the spatial configuration of the landscape can be observed during 1997-2009, primarily in the form of changes towards less dense vegetation, transitioning from closed forest to shrubs, open forest, and clusters of trees.
Between 1985 and 2018, the loss of permanent snow and glaciers produces the highest difference in spatial configuration of the landscape, with major differences at the intersection of the Western and Eastern Central Alps, and the Southern Prealps ( Figure 11).

Change in Spatial Patterns
Over the whole study period, most pixels show a non-negligible change in the spatial pattern. The highest difference in the spatial pattern is observed in the Plateau and northern Prealps regions, notably in the 1997-2009 period ( Figure 10). When grouped by the Basic Categories, all intervals undergo a similar distribution of change, however, more pixels show a JSD value > 0.012, reflecting the higher dynamism of LC change at the higher level of category disaggregation. Amongst the areas of highest distance, the most important changes to the spatial configuration of the landscape can be observed during 1997-2009, primarily in the form of changes towards less dense vegetation, transitioning from closed forest to shrubs, open forest, and clusters of trees.
Between 1985 and 2018, the loss of permanent snow and glaciers produces the highest difference in spatial configuration of the landscape, with major differences at the intersection of the Western and Eastern Central Alps, and the Southern Prealps ( Figure 11).

Land Cover Dynamics in Switzerland
At the national scale, the gross values of LC gain and loss over the study period would suggest that Switzerland's landscape is relatively static, with just over 5% of land cover changing in the Principal Domain between each dataset, and low levels of difference in configuration and composition of land cover at the scale of 25 km 2 . However, exploration of pixel-level statistics shows that the land cover has changed above what would be expected under the hypothesis of random change, and particularly at the levels of higher disaggregation, land cover has been highly dynamic. Furthermore, many of these changes have been equally dynamic in terms of both LC gain and loss, which, while this may be reported as a net persistence in LUCC statistics, indicates a more dynamic landscape. Changes are of a lower magnitude than those reported for recent LU changes [26].
The low values of difference at the landscape scale, however, highlight that LC changes have taken the form of mostly multiple small-scale transformations rather than large contiguous areas. This is further emphasized by the weak, and mostly one-sided results, achieved when exploring the dominant signals of LC change. Change at spatial scales above the pixel level is most important for areas of glacier retreat and loss of permanent snow and ice, reflecting the significant loss of low-altitude snow and ice with climate change [26]. As indicated by official Land Use reporting [26], the transition from glacier to newly vegetated areas, while reported in the literature [47], is not yet significant enough to be detected through LUCC intensity analysis. The low-absolute levels of differences reported, however, reflect that these areas remain relatively intact, with loss occurring at the edges of the glaciers, rather than in the form of fragmentation of the landscape.
While built surfaces have shown to continue to increase the surface area with each survey period, intensity analysis confirms the published results showing the rate at which these surfaces' gain area has decreased [26]. The continued expansion of forest areas over former agricultural surfaces in alpine areas [26] and observed 'greening' of the plateau region [8] is also reflected.
All regions undergo a higher annual rate of land cover change between the 2009 and 2018 datasets. At the category level, spatial disaggregation of change is more evident, with the Jura and Plateau regions showing highly dynamic rates of change for vegetated categories. Investigation at the level of category transitions show that the key transition occurring within the Plateau region is to and from the closed forest, a trend which is more important in this region compared to the alpine regions where a relative change to expansive forest areas is less dramatic. The similarity in the intensity of gain and loss in the forest areas of the plateau region reflects the findings of the National Forest Inventory, which reports that the total forest area in this area was stable over the study period. The results presented indicate that the reason for the increased rate of change in the final time interval is due, in part, to the transition from brush vegetation to forest classes. However, this increased rate of change in the Plateau region is at odds with reports of the forested land in this region remaining static in recent years [26]. This increased rate of transition could indicate the emergence of new pressures on LC, or a change of intensity of certain pressures [19] that are not yet occurring at a magnitude detectable in absolute LC change.

The Need for Higher Spatial & Temporal National Land Cover Data
Accurate, reliable, and timely land cover data are crucial for the understanding and modeling of environmental processes [48], as well as for informing socio-economic policy at the right pace for action [49]. As such, the low update frequency of the ArealStatistik is at odds with the data needs for quantifying variables that are known to be as dynamic as land cover and land use [50]. Analysis of annual land cover maps would facilitate a greater understanding of the processes and drivers that result in LUCC [19,51] and allow for consideration of both subtle and rapid changes that can both indicate changes to LUC pressures.
Understanding ecological processes, including disturbances and climate feedback, which are key to the comprehensive assessment of the drivers and pressures facing the environment, requires the utilization of multiple time-steps of annual or even sub-annual land cover data [51,52]. This need is particularly evident considering the increased rate of change identified in some categories during the 2009-2018 interval and the difficulty in establishing a cause of this.
The knowledge gap surrounding the dynamics of LUCC in Switzerland is, to an extent, perpetuated by the limited temporal resolution of the datasets used. At the current resolution, the question remains whether the observed increased intensity of change between the two most recent survey periods is indicative of real changes in Land Cover change drivers, or merely an artifact of the periodicity of sampling that may be less relevant upon production of the next LC dataset. The analysis of absolute trends over the 9-12 years intervals presented here provides an indication of changes or persistence, however, there may be intermediate growths and declines in categories that manifest as an overall non-trend and are not effectively captured using the ArealStatistik. This hinders the effective use of time-series change detection and trend algorithms such as Breaks For Additive Season and Trend (BFAST) [53][54][55], Continuous Change Detection and Classification (CCDC), or LandTrendr [56]. Moreover, the timing of the aerial image acquisition can have an influence on the consistency of classification and observed changes. Indeed, for the last survey covering the period 2013 to 2018, it took 6 years to cover the entire country, meaning that each year, approximately one-sixth of the territory is covered and processed according to the schema of acquisition presented on https://www.bfs.admin.ch/bfs/fr/home/statistiques/ espace-environnement/enquetes/area/bases-donnees/photographies-aeriennes.html, accessed on 29 May 2023. Once the entire country is covered, the dataset is then released for the given survey period. With such a sequential approach, lots of LUC change can occur from the time the first tiles are processed until the last one. Consequently, LUC transitions can be missed, lowering the proportion of observed changes. This reinforces the need for yearly, consistent, accurate, and high-resolution LC maps for Switzerland.
It can be further argued that analysis of LUCC benefits from a consideration of different spatial scales, as landscape patterns sometimes show a high degree of variability depending on the spatial scale used [19], and dynamic but small-sized landscape areas can contribute greatly to observed changes. The relatively coarse spatial resolution of the ArealStatistik may mean that this dynamism is not fully captured. Refining LC data to a smaller spatial resolution facilitates the capture of processes occurring at the scale of human activity [57].
Recent efforts to downscale LUC data for Switzerland have improved the spatial resolution to 25 m, however, there remains a gap in the temporal dimension of datasets produced [21]. Several initiatives exist to produce global land cover datasets at high-spatial and temporal resolutions [58][59][60], however, their use at a local-national scale can be limited by variability in accuracy across space and a relatively low number of thematic classes compared to national datasets [48,61].

Perspectives
Although the ArealStatistik is useful and, thematically, more precise than commonly used classifications, it suffers from various limitations (i.e., spatial (1-hectare) and temporal (every 6 years) resolutions for consistent environmental monitoring. Consequently, this survey-based method constrains researchers and policymakers to work with data that do not necessarily reflect on the ground realities impending the ability to correctly capture detailed landscape features, qualities, particularities, and configurations, as well as the entire dynamics of change (e.g., low temporal granularity allows for capturing slow onset changes and not the potential rapid change of the landscape).
The big and constantly increasing amount of freely available EO data (Big Earth Data [62,63]) with enhanced temporal, spatial, and spectral resolution, together with rapid development in remote sensing data storage, handling, and processing technologies (e.g., Data Cube, Machine Learning), offers new opportunities for LUC mapping at finer spatial, temporal, and thematic resolutions [57,64,65]. The increasing availability of openaccess, high-resolution remote sensing data has vastly increased the potential for the development of land cover datasets from local to global scales [50]. Satellite imagery provides a consistent dataset of earth observations that is spatially continuous and contains the temporal resolution necessary to identify classes with strong temporal dynamics [48,61].
However, among the major identified issues concerning the increase in spatial and temporal resolution when producing LC data are (1) the possibility to discriminate a higher number of features, which, at the same time, leads to a general increase of variability in LC classes and, consequently, a decrease in accuracy level when using traditional classification methods; and (2) the lack of ground-reference data for validating the results [66][67][68].
To overcome these issues, Machine/Deep Learning (ML/DL) algorithms (e.g., convolutional neural networks, gradient-boosted trees) have demonstrated their potential to generate results of higher accuracy than those produced with parametric classifiers, being able to deal with higher dimensionality (e.g., space and time) and map classes with complex characteristics. This allows for the development of supervised classification models that correlate different information (spectral, spatial, temporal, and ancillary) to an area of interest [69].
The "ArealStatistik" is a unique dataset set providing 4 million ground-referenced points for each period (4 in total). Part of this data set can be used for both the train, test, and validate stages of classification. For training, this can be used as reference data to identify a sample of pixels of known classes; then in the allocation phase, pixels will be assigned to a given class according to their statistical similarities; finally, the accuracy of the classification can be tested by comparing a sample of pixels in both classified and reference data and generate an error matrix [70,71]. Recent advances in machine learning techniques, such as multitask learning, multi-instance learning, and multiview learning, could be also explored to tackle issues such as landscape heterogeneity in space and time (e.g., vegetation), potential lack of training data, and rarity of land cover change while developing a suitable model for LC and LCC [72,73].
Therefore, to improve the spatial and temporal resolution of land cover data of Switzerland, while keeping the thematic richness of the "Arealstatistik", the availability of more than 39 years of satellite EO Analysis Ready Data over Switzerland, made available by the Swiss Data Cube (https://www.swissdatacube.ch, accessed on 29 May 2023) [74][75][76], together with ML/DL techniques [77,78], allow the envisioning of innovative approaches to produce a yearly 10/30 m consistent time-series of LC and its changes, informing on class stability and transitions [79,80], to support effective national environmental monitoring [81,82].

Conclusions
The results presented here highlight the continued dynamism of LC in Switzerland over a period of 33 years. Assessment of spatial configuration identified low degrees of change across most of Switzerland in all time periods assessed, reflecting the many pixel-level transitions which have occurred. The findings, however, reflect the difficulty in assessing systematic LC change at the national and regional level, as the large range of small, spatially dispersed and gradual trends required aggregation to be detected, but this means that any locally important drivers or trends may have been masked. A clear dynamic over the study period is the high rates of annual change in low-growing Brush Vegetation, in the form of meadows and shrubs in the low-altitude regions, for which an increased intensity of change was detected in recent years. This could be a sign that drivers of LC change are evolving, but with the currently available data, it is not possible to deduce the drivers of this. The increased rate of transition between Brush Vegetation and Tree Vegetation LC classes in the Jura and Plateau regions are more dynamic, as there is a greater human influence in LUC, but there is not a clear pathway at the national or regional scale as to why shrub vegetation would be changing more than it has previously. Trends such as an increased rate of change to and from forest LC classes can be identified, but an understanding of reasons behind this is limited by infrequent updates to LC data. This work highlights that although the ArealStatistik is very useful thanks to its thematic richness, neither its low spatial resolution or update frequency allows for providing accurate and timely information to depict and understand the detailed dynamics of LCC and the related impact across the country. Accurate LC change assessment and effective LC projections require higher spatial (e.g., 10-30 m) and temporal (e.g., yearly) data products to build consistent time series. This restricts the comprehensive understanding of the gains and losses of LUC types, including the magnitudes, locations, and timings of transitions. Understanding the drivers of observed changes, particularly a potential increase in the rate of change towards Tree Vegetation classes, is important to inform both reflections on past environmental policies and the development of future management strategies. The ability to understand these drivers would benefit from a high-resolution annual LC dataset. Advances in ML/DL techniques, combined with increased availability of satellite data time series available in the Swiss Data Cube, present the potential for the development of such a dataset, in combination with the ArealStatistik, which provides a geographically specific LC classification scheme.