Effects of Terrain Parameters and Spatial Resolution of a Digital Elevation Model on the Calculation of Potential Solar Radiation in the Mountain Environment: A Case Study of the Tatra Mountains

: Solar radiation signiﬁcantly affects many processes on Earth. In situ measurements are demanding and require a dense network of sensors. A suitable alternative solution could be the modelling of potential solar radiation based on a digital elevation model (DEM) in geographic information systems. The key issue of this study is to determine the inﬂuence of the terrain parameters and the spatial resolution of a DEM on the calculation of potential solar radiation. The area of study is the Tatra Mountains (the highest mountains of the Carpathians). The DEM determined from light detection and ranging (LiDAR) was used. To determine the inﬂuence of the terrain, the following terrain parameters were applied: slope; aspect, represented by northness and eastness; elevation; and topographical position index using six different circular neighbourhoods (10 m, 30 m, 50 m, 100 m, 500 m, and 1000 m). The results indicate a moderate correlation (0.32–0.46) between the solar radiation calculation errors and the absolute values of the topographic position indices with small neighbourhoods (10 m–100 m). To show the impact of the spatial resolution, the calculation was performed based on four different DEM resolutions, namely 5 m, 10 m, 30 m, and 90 m. Mutual differences in potential solar radiation were quantiﬁed concerning the topographic position index. The result is also a model of potential annual solar radiation in the Tatra Mountains, calculated at a resolution of 5 m or 2 m.


Introduction
Solar radiation is a general term for the electromagnetic radiation emitted by the Sun [1]. It is a source of energy and significantly affects the environment, as well as geological and natural processes and the health and vitality of humans, animals, and plants [2][3][4][5]. It is also an important renewable energy source that helps sustain life on Earth. Solar energy can be actively used to generate electricity in power plants but also passively to generate electricity using photovoltaic systems on the roofs and facades of buildings. To optimize the distribution of these devices, it is necessary to know the distribution of solar radiation in the country. Because solar radiation is essential for plant growth, knowledge of its distribution also helps, for example, in planning the cultivation or protection of sun-demanding plants. In the mountain environment, it is specifically a matter of forest protection. Another application in the mountain environment, where it is necessary to know the distribution of solar radiation, is the prediction of avalanche hazards. Solar radiation that hits the georelief is the result of a complex interaction between the atmosphere and the Earth's surface. Its amount depends on many factors, for example, geographic location, season, time of day, local landscape, cloud cover, and other optical ISPRS Int. J. Geo-Inf. 2022, 11, 389 2 of 16 properties of the atmosphere [6]. In the case of considerable vertical fragmentation of the area, direct measurements are insufficient to express the spatial and temporal variability of solar radiation incidents on the surface of the georelief [7,8]. Although there are thousands of solar radiation monitoring stations around the world, for most geographic areas, there is no accurate data available on solar radiation [9]. Management of measurements in a mountain environment is even more problematic due to difficult access to rugged terrain. In situ measurements are also demanding and require a dense network of sensors. Therefore, in addition to the available measurements, it is necessary to use existing knowledge and create suitable models of solar radiation. Because modelling the influence of local weather on solar radiation modelling is very demanding and inaccurate, potential solar radiation is often calculated based on digital elevation models (DEMs) [10,11]. This is also due to the fact that many analyses require modelling of potential solar radiation with sufficient positional accuracy, albeit with lower temporal accuracy. In addition, for many spatial analyses, the value of solar radiation at a certain point in relation to its surroundings is often more important than absolutely accurate values in the whole area. A suitable parameter that meets these requirements is potential annual solar radiation calculated at high spatial resolution in geographical information systems (GIS).
The main spatial distribution factor of potential solar radiation is georelief, especially elevation, slope, and aspect. This paper is therefore devoted to determining the effect of the used DEM resolution, as well as the above-mentioned terrain parameters, on the calculation of potential solar radiation in GIS.
Currently, the most used software tools for the calculation of potential solar radiation are the Solar Radiation toolset implemented in ArcGIS [12,13] and r.sun [14], which was designed for the GRASS GIS [15]. Selected studies that describe the implementation of calculation of solar radiation based on 3D data are also presented in [9]. The article [9] states that one of the main models created in a GIS tool is Solar Analyst, developed by Fu and Rich [16][17][18][19][20] based on a theory called SolarFlux. The model was created as a module in ArcView 3.0 software created by ESRI (Environmental Systems Research Institute) and later implemented in the ArcGIS 10 extension known as the Solar Radiation toolset [12,13].
Many articles have been devoted to the calculation of potential solar radiation. Some of them focus on fundamental research of the calculation of potential solar radiation [5,[21][22][23][24]; others focus on its use in various applications [8,[25][26][27][28][29][30][31][32][33]. For example, [8] is devoted to the calibration and validation of the ArcGIS Solar Radiation tools for photovoltaic potential determination, focusing mainly on the input atmospheric parameters, diffusivity, and transmissivity but also pointing to the importance of spatial resolution in the calculation of potential solar radiation. According to [8], most GIS-based methods for calculating solar radiation are based on some form of geographic data, such as DEMs [10,11,21,22] or light detection and ranging (LiDAR) data [5,[24][25][26][27][28]. These methods also use different assumptions and therefore differ in terms of accuracy and performance [8]. For example, refs [5,24,25,28] describe tools for calculating potential solar radiation from 3D LiDAR data.
In the field of applications and uses of calculated potential solar radiation in mountain areas, the effect of solar radiation on forest damage by bark beetle attacks was described in some papers [34][35][36][37][38][39][40][41]. Specifically, in [30,39,40], potential solar radiation was applied as one of the site characteristics in the determination of the probability of bark beetle attack in the Tatra Mountains. In [30], the authors stated that potential solar radiation calculated based on the DEM better identified beetle infestations than commonly used meteorological variables. However, it should be emphasized that DEM at a resolution of 30 m (e.g., based on shuttle radar topography mission (SRTM) [41][42][43]) was mostly used in those analyses.
Although the influence of the DEM resolution on the calculation of potential solar radiation is known to be significant [8,10,11], it has not yet been quantified how accurate the input DEM should be or which error in the calculation of potential solar radiation would result from using a lower-quality DEM. The authors mostly point to different DEMs, which means that the results also show a difference resulting from another data source. In this study, we used only one data source and aggregated it by the average function to different ISPRS Int. J. Geo-Inf. 2022, 11, 389 3 of 16 resolutions. This leads to a real indication of the influence of the DEM spatial resolution on the calculation of potential solar radiation. We used an available DEM derived from LiDAR data.
Another issue is how terrain fragmentation affects the accuracy of the calculation, although we can assume that differences in potential solar radiation are more significant in mountainous environments. Therefore, the aim is also to determine the influence of terrain parameters on errors in the calculated potential radiation caused by the lower resolution of the input data. Because the values of solar radiation change most in rugged terrain, we focused mainly on the mountain environment, specifically on the Slovak side of the Tatra Mountains (the highest mountains of the Carpathians).
We assume that based on the results, we will be able to better determine in which areas it is important to use a DEM with higher resolution and where a lower resolution is sufficient to model potential solar radiation. In the experiment, we tested DEM resolutions of 5 m, 10 m, 30 m, and 90 m.
This paper is further organized as follows. Section 2 describes the data and methods used, including a description of the study area, available data sources, the design of the experiment, and a description of the software tools used. Section 3 presents the results of quantifying the effects of DEM resolution on the calculation of potential solar radiation, as well as correlations between terrain parameters and errors in the calculation of potential solar radiation caused by lower spatial resolution. Section 4 provides a discussion of the results, and Section 5 concludes the document.

Study Area
The study area is located in the Tatra Mountains (Central Europe), the highest range within the Carpathians, located on the border between Poland and Slovakia ( Figure 1). In 1949, the Tatra Mountains were declared the first national park in Slovakia (TANAP), and in 1993, together with the Polish part of the Tatras, they were declared the "Biosphere Reserve of the Tatras" by UNESCO [44]. The High Tatras include as many as 30 peaks more than 2500 m above sea level (a.s.l.), with the elevation reaching 2655 m a.s.l. (Gerlach). The Tatra Mountains are characterized by deep valleys and mountain ridges, so they are very suitable for performing an analysis of the effects of georelief on the calculation of potential solar radiation. The study area has a west-east elongated shape (Figure 1), and its area is 959 km 2 , corresponding to zone 26 from LiDAR data provided by the Geodesy Cartography and Cadastre Authority of the Slovak Republic (GCCA SR; the abbreviation in Slovak used to denote data is ÚGKK SR) [45] (Figure 2).

Data
For the Slovak side of the Tatra Mountains, solar radiation data are partly provided by the Slovak Hydrometeorological Institute [48]. Unfortunately, the network of ground stations is not dense enough for precise spatial analysis, and measurements are relatively expensive. The data on solar radiation is also provided, for example, by Global Solar Atlas [49]. These data are very good in terms of temporal accuracy and also consider meteorological conditions. However, their spatial resolution (pixel size: 250 m) is not sufficient for detailed spatial analysis. Therefore, it is more appropriate to use calculated potential solar radiation for many spatial analyses with higher spatial resolution. The calculation of potential solar radiation is based on DEM. In spatial analysis in the Tatra Mountains, SRTM [43] with a resolution of 10 m or more has been an available and often used DEM. However, since 2017, the GCCA SR has been continuously creating and providing a digital terrain model, DTM 5.0, and digital surface model, DSM 1.0, of the entire territory of the Slovak Republic [45]. The DTM 5.0 and DSM 1.0 models were created from aerial laser scanning (ALS) data, and the declared density of the data is at least five points per m 2 . Both models are provided in raster format at a spatial resolution of 1 m. The whole territory of the Slovak Republic is divided into 42 localities ( Figure 2). Scanning takes place gradually in individual locations from the west of Slovakia to the east. Available areas for users are highlighted in green in Figure 2. Number 26 indicates the study area. The parameters of ALS data collection from area 26 (Tatras) are presented in Table 1. We assume that these data provide a suitable basis for more accurate modelling of potential solar radiation.

Data
For the Slovak side of the Tatra Mountains, solar radiation data are partly provided by the Slovak Hydrometeorological Institute [48]. Unfortunately, the network of ground stations is not dense enough for precise spatial analysis, and measurements are relatively expensive. The data on solar radiation is also provided, for example, by Global Solar Atlas [49]. These data are very good in terms of temporal accuracy and also consider meteorological conditions. However, their spatial resolution (pixel size: 250 m) is not sufficient for detailed spatial analysis. Therefore, it is more appropriate to use calculated potential solar radiation for many spatial analyses with higher spatial resolution. The calculation of potential solar radiation is based on DEM. In spatial analysis in the Tatra Mountains, SRTM [43] with a resolution of 10 m or more has been an available and often used DEM. However, since 2017, the GCCA SR has been continuously creating and providing a digital terrain model, DTM 5.0, and digital surface model, DSM 1.0, of the entire territory of the Slovak Republic [45]. The DTM 5.0 and DSM 1.0 models were created from aerial laser scanning (ALS) data, and the declared density of the data is at least five points per m 2 . Both models are provided in raster format at a spatial resolution of 1 m. The whole territory of the Slovak Republic is divided into 42 localities ( Figure 2). Scanning takes place gradually in individual locations from the west of Slovakia to the east. Available areas for users are highlighted in green in Figure 2. Number 26 indicates the study area. The parameters of ALS data collection from area 26 (Tatras) are presented in Table 1. We assume that these data provide a suitable basis for more accurate modelling of potential solar radiation.  The data are freely available, but the user is obliged to state the source of ALS products when creating his own work and publishing it as follows: "ÚGKK SR". Data from smaller areas (up to 400 km 2 ) can also be downloaded via the ZBGIS Map Client application [50] ( Figure 2). In this case study, we used DTM 5.0 aggregated to a resolution of 5 m. Because DTM 5.0 is the name of a specific model that could be replaced by another DEM in our study, we prefer to use the general abbreviation DEM hereafter. We use the abbreviation DTM 5.0 only if we want to emphasize this model specifically. The data are freely available, but the user is obliged to state the source of ALS products when creating his own work and publishing it as follows: "Ú GKK SR". Data from smaller areas (up to 400 km 2 ) can also be downloaded via the ZBGIS Map Client application [50] (Figure 2). In this case study, we used DTM 5.0 aggregated to a resolution of 5 m. Because DTM 5.0 is the name of a specific model that could be replaced by another DEM in our study, we prefer to use the general abbreviation DEM hereafter. We use the abbreviation DTM 5.0 only if we want to emphasize this model specifically.

Experimental Design
In an experiment, we compared potential solar radiation calculated at the resolution corresponding to SRTM (10 m or 30 m) or lower (90 m) with values calculated at a finer resolution (5 m). We chose resolutions 3 times higher and 3 lower than the resolution of the commonly used SRTM model in the study area. The highest resolution used (5 m) was one with which it was possible to perform each of the designed and implemented calculations in the ArcGIS environment in a relatively short time (in this case, within two days). We used DTM 5.0 data aggregated from the original resolution of 1 m to a resolution of 5 m, 10 m, 30 m, and 90 m because we wanted to emphasize the effect of the input resolution used, not the DEM differences. Raster aggregation generates a reduced-resolution version of a raster. In this case, we used the aggregation function Mean. Therefore, each output cell contains the mean of the input cells that are encompassed by the extent of that cell. Therefore, we included all values from the original raster in the calculation, which is an advantage compared to using one of a resampling method, such as Nearest, Majority, Bilinear, or Cubic [52]. Using aggregation, we later also created models of potential solar radiation with resolutions of 10 m, 30 m, and 90 m. We aggregated them from a model calculated with a resolution of 5 m. For the above reason, we assume that it is more appropriate than calculating them directly at a lower resolution. We emphasize that our aim is to point out the effect of terrain parameters and spatial resolution of DEM on the calculation of potential solar radiation, not to calculate potential solar radiation as accurately as possible. Four different resolutions are sufficient for this purpose. We did not use the raster with the original resolution of 1 m due to the computational complexity of determining the potential solar radiation in such a large area. A model with a resolution of 2 m can also be calculated, but the calculation can take more than 10 days.
To assess the impact of terrain, we calculated the terrain parameters (aspect and slope) and the derived terrain parameters, such as northness, eastness [53], and topographic position indices (TPIs) [54,55]. TIPs compares the elevation of each cell in a DEM to the mean elevation of a specified neighbourhood around that cell [52]. Although aspect is an important terrain parameter, it is difficult to include it directly in the analysis due to its angular scale (0 • to 360 • , with values 0 • and 360 • being close to each other). Therefore, we suggest using the parameters of northness (northness = cos(aspect)) and eastness (eastness = sin(aspect) [53]) instead of aspect. The slope is expressed in degrees but only in the range of 0 • to 90 • , so we included it in the analysis in this form. To better express terrain edges and extremes, we suggest applying absolute values of TPI (i.e., we do not distinguish whether it is a higher or lower position of points, and only mutual differences are important). All parameters were calculated based on DTM 5.0 [45] at a spatial resolution of 5 m.
We subsequently calculated the correlations of the terrain parameters with the differences in potential solar radiation calculated at various resolutions. To exclude values from marginal areas, which may be inaccurate due to incomplete marginal input data, we excluded from the comparison the zone up to 200 m. Finally, we compared the calculated and aggregated values of potential annual solar radiation at multiple resolutions. The simplified experimental procedure is shown in Figure 3.

Software Tools and Data Processing
We calculated potential annual global solar radiation in the ArcGIS 10.8 software environment using the Area Solar Radiation tool (implemented in the Solar Radiation toolset) [12,13]. This tool calculates potential solar radiation over a geographic area based on the hemispherical viewshed algorithm explained in [18,19]. The total amount of radiation is calculated for a given location in units of Wh/m 2 [12,13]. The Area Solar Radiation tool takes location (latitude-longitude), elevation, slope, aspect, and atmospheric transmission as the most relevant inputs [8,13]. These parameters usually result from the input data. Other input settings are listed in Table 2. We used default values, except for latitude, which was automatically calculated for the site area.

Software Tools and Data Processing
We calculated potential annual global solar radiation in the ArcGIS 10.8 software environment using the Area Solar Radiation tool (implemented in the Solar Radiation toolset) [12,13]. This tool calculates potential solar radiation over a geographic area based on the hemispherical viewshed algorithm explained in [18,19]. The total amount of radiation is calculated for a given location in units of Wh/m 2 [12,13]. The Area Solar Radiation tool takes location (latitude-longitude), elevation, slope, aspect, and atmospheric transmission as the most relevant inputs [8,13]. These parameters usually result from the input data. Other input settings are listed in Table 2. We used default values, except for latitude, which was automatically calculated for the site area. For calculation of the terrain parameters and parameters derived from them, we used standard tools, such as Aspect, Slope, Focal Statistics, and Raster Calculator (Spatial Analyst) in ArcGIS 10.8 [56]. The input data representing variables were used in raster form at a spatial resolution of 5 m.

Results
The results of the analysis can be summarized in several points: (1) Differences in potential annual solar radiation calculated at resolutions of 10 m, 30 m, and 90 m compared to potential annual solar radiation calculated at a resolution of 5 m, supplemented by their statistical characteristics; (2) Correlations between absolute values of potential solar radiation differences and the terrain parameters to assess the influence of terrain fragmentation on the accuracy of the calculation of potential solar radiation; (3) Comparison of calculated and aggregated models of potential annual solar radiation; (4) The model of annual potential solar radiation in the Tatra Mountains with a cell size of 5 m (or 2 m), which is applicable in various spatial analyses. Table 3 contains the statistical characteristics (minimum, maximum, mean, and root mean square error (RMSE)) of the calculated differences in solar radiation at a resolution of 5 m compared to the potential solar radiation calculated at a resolution of 10 m, 30 m, and 90 m. The differences in potential annual solar radiation calculated at a resolution of 30 m compared to potential annual solar radiation calculated at a resolution of 5 m are shown in Figure 4a,b. Figure 4c shows TPI calculated with a circular neighbourhood of 100 m (TPI100) and illustrates, in detail, the relationship between the calculated differences shown in Figure 4a,b and TPI100. It is also important to note that the processing time increases significantly with increasing raster resolution. In our study, the processing times were 0 m:

Correlations between Absolute Values of Potential Solar Radiation Differences and the Terrain Parameters
Potential solar radiation, slope, aspect, and TPI100, all calculated at a resolution of 5 m, are shown in Figure 5. Other TPI values are shown as their absolute values to better demonstrate their relationship with differences in potential solar radiation calculated at different resolutions ( Figure 6). demonstrate their relationship with differences in potential solar radiation calculated at different resolutions ( Figure 6).    The results indicate independence of the potential solar radiation differences on eastness; very low dependence on |TPI1000|, |TPI500|, and northness; and low dependence on elevation. Of particular interest is the low dependence on aspect (northness and eastness), on which the potential solar radiation directly depends. Conversely, the dependence on absolute values of TPI with input neighbourhoods of 10 m to 100 m is relatively high ( Table 4). The correlation coefficient of these parameters reaches values of 0.33-0.46, suggesting a moderate correlation, although still sufficiently indicating the influence of TPI on errors in the calculation of potential solar radiation. The inaccuracy of its calculation on ridges and valleys also affects their surroundings, which may not be characterized by significant TPI values. Moderate correlations were also shown for slope (0.31-0.38) ( Table 4). Table 5 shows the intervals of absolute values of TPI50 and the associated RMSE of the potential solar radiation models calculated at a resolution lower than 5 m. On mountain ridges and valleys, RMSE reaches up to 326 321 Wh/m 2 (at resolutions of 5 m and 90 m). The most significant differences in the calculation of potential solar radiation resulting from different raster resolutions are in areas with high absolute values of TPI (for resolutions of 10 m, 30 m, or 90 m, values greater than 10 m lead to 4.7-6.2 times greater errors compared to errors in the flat area (|TPI10| < 5 m)) ( Table 5). A graphical representation of absolute values of TPI, together with differences in potential solar radiation (5 m vs. 30 m), is shown in Figure 6. This figure also confirms that the most considerable differences in potential solar radiation values are at and around peaks, ridges, and valleys.

Comparison of Calculated and Aggregated Models of Potential Annual Solar Radiation
The next result obtained in this study is a comparison of potential solar radiation calculated based on DEM at a resolution of 10 m, 30 m, and 90 m and values aggregated from a solar radiation model with a resolution of 5 m to a resolution of 10 m, 30 m, and 90 m, respectively. The statistical characteristics (RMSE) of the differences are shown in Table 6. Based on the results, we can confirm that it is more appropriate to aggregate solar radiation from a higher-resolution model. Table 6. Comparison of the differences in potential solar radiation calculated and aggregated at a spatial resolution of 10 m, 30 m, and 90 m. The model with which these models are compared is calculated with a resolution of 5 m.

Model of Annual Potential Solar Radiation in the Tatra Mountains with a Cell Size of 5 m
The last result obtained in our study is the model of annual potential global solar radiation in the Tatra Mountains with a spatial resolution of 5 m (Figure 5b) or 2 m. The reference coordinate system used is S-JTSK (EPSG: 5514) with a file size of 464 MB (or 2.9 GB for a spatial resolution of 2 m). The above results (Section 3.3) indicate the suitability of this model for analyses where lower spatial resolution is sufficient.

Discussion
Based on the obtained results, we can state that the resolution of input data in the calculation of potential solar radiation in the GIS environment has a significant impact on the results, especially in areas with high absolute values of TPI calculated with small neighbourhoods (10 m-100 m). This is especially evident in high mountain areas, which contain rugged terrain. Consequently, in areas with a |TPI50| or |TPI100| greater than 10 m, we recommend applying a high-resolution DEM to calculate potential solar radiation. In this study, we used a model of potential annual solar radiation with a spatial resolution of 5 m. In addition, we calculated a model at a resolution of 2 m, which we did not include in the analyses in this study due to computational complexity and time-consuming calculations. However, it can be used for other follow-up analyses, especially in rugged terrain. The second reason is the fact that four different resolutions are sufficient to show the effect of resolution on the result. The authors of [13] also stated that calculating insolation can be very time-consuming. Calculations for a large DEM can take several hours, and for a very large DEM, even days [13]. Similarly, the authors of [8] points to the high computational complexity of the calculation of potential solar radiation at a high raster resolution. In comparison, for a study area of about 1 km 2 , the processing time was recorded with a standard model run of 01 m: 12 s, 06 m: 22 s, and 10 m: 7 s for 30 m, 5 m, and 0.50 m, respectively [8]. The calculations were performed on a Windows computer with an Intel i5 processor with four cores and eight GB of RAM. It should also be emphasized that processing time increases when using shorter time intervals.
However, it should be noted that even a file size of more than 2 GB, in this case, is a significant limit. In this context, it is also important to note that the Area Solar Radiation tool is designed only for local landscape scales, where it is generally acceptable to use one latitude value for the whole area [13]. With larger datasets, such as for states, countries, or continents, the results will differ significantly at different latitudes (greater than one degree). For this reason, it is not sufficient to use only a higher raster resolution for a calculation of a larger area. It is also necessary to divide broader geographical regions into smaller zones with different latitudes [13]. Our analysis in the Tatra Mountains was related to the elongated shape of the area with a difference in latitude of 0.22 degrees, although with a range of longitude greater than 1 degree [51]. In this case, the error caused by the latitude setting leads to a potential solar radiation error of as much as 4000 Wh/m 2 .
The authors of [8] argue that the default model values of diffusivity and transmissivity lead to a substantial underestimation or overestimation of solar radiation. The authors further claim that model validation is necessary because actual values cannot be defined from atmospheric data prior to model implementation [8,57]. The values of solar radiation from meteorological stations are certainly useful for a comparison of the values of the calculated potential solar radiation. However, even the values of the modelled potential solar radiation are sufficient to perform many spatial analyses, especially if it is not necessary to know the exact values of solar radiation but rather their local changes in the area. The values calculated based on the DEM in our study are suitable for this purpose.
Therefore, we recommend using the potential annual solar radiation model calculated in this study with a resolution of 5 m (Figure 5a, Supplementary Material File S1) or higher for analysis in the Tatra Mountains (see Data Availability Statement). In the case of analyses in which a lower resolution is sufficient, the values of the calculated model can be aggregated to the required lower resolution (rather than calculating potential solar radiation with a lower resolution). We would like to add that model values need to be aggregated through the average function, not resampled, as resampling leads to less reliable values.

Conclusions
We found and confirmed that the most significant differences in the calculation of potential solar radiation resulting from different raster resolutions are in areas with high absolute values of TPI, i.e., mountainous areas (for resolutions of 10 m, 30 m, or 90 m, values greater than 10 m lead to 4-6 times greater errors compared to errors in a flat area).
According to our analysis, more accurate values of potential solar radiation at the required resolution (e.g., 30 m) are obtained by calculation at higher resolution (e.g., 5 m) and subsequent aggregation of values by the average function.
Based on the results, we recommend: • In mountainous areas, calculate potential solar radiation at as high a resolution as possible (in a flat area, this is not necessary); • For spatial analysis in the Tatra Mountains, we recommend using the solar radiation model proposed in this study (Figure 5a, Supplementary Materials File S1) or its aggregation at a lower resolution.
In future research, it could be useful to compare our results with the potential solar radiation calculated on the basis of SRTM at a resolution of 30 m or other DEMs used in spatial analyses in the Tatra Mountains. Another open issue is the calculation of potential solar radiation with even higher resolution, mainly due to analyses in smaller localities requiring more detailed knowledge of solar radiation.
Although our analysis was performed in a specific area, the results may be applicable to other areas. A high-spatial-resolution model of potential solar radiation can help improve the accuracy of various geological and ecological spatial analyses.

Data Availability Statement:
The data presented in this study are available from the corresponding author upon request, including the potential annual global solar radiation model calculated at a resolution of 2 m (2.9 GB). The potential annual global solar radiation model calculated at a resolution of 5 m is available in Supplementary Materials (131 MB). The input data are available at: https://www.geoportal.sk/en/zbgis/als_dmr/ (accessed on 21 February 2022). The data are freely available, but the user is obliged to state the source of aerial laser scanning products when creating his own work and publishing it as follows: "ÚGKK SR".