Extreme Longevity: Analysis of the Direct or Indirect Influence of Environmental Factors on Old, Nonagenarians, and Centenarians in Cilento, Italy

Longevity is rightly considered one of the greatest achievements of modern society, an achievement understood as the possibility of increasing the healthy part of life and not only its full duration. This study investigated the phenomenon of regional longevity in Cilento by analyzing the association between longevity indicators and some environmental factors, in order to understand if factors such as altitude, climate, UNESCO protected areas, and hinterland can directly or indirectly influence the measure of healthy living and lead to longevity. Demographic and environmental data were collected through the Archives of the Cilento municipalities, the National Institute of Statistics, the Italian Institute for Environmental Protection and Research, and the Italian National Commission for UNESCO. The Geographically Weighted Regression were used to determine the association between longevity indicators and environmental factors. Correlation analysis between the longevity indicators was investigated in order to have a complete picture of longevity in Cilento. It was discovered that Cilento longevity is mainly found in the central area of the territory and from there, by age groups, it widens towards the south-eastern area. This designated area is part of the hilly area of the Cilento, at an altitude between 400 and 700 m above sea level. The towns of this area are part of the UNESCO heritage and are characterized by a transitional climate between Mediterranean and temperate climate. Moreover, the correlation analysis between the six indicators of longevity has confirmed a linear relationship between the six variables and this indicates that in the years to come there will be the probabilities of a generational turn over between the old, great old, nonagenarians and current centenarians, provided that the SARS-CoV-2 pandemic subsides.


Introduction
Due to the decline in birth rate and the increase in life expectancy, the population worldwide is aging at an unprecedented speed. According to data from World Population Prospects [1], there are 727 million older people (aged 65 and older) and half a million centenarians (aged 100 and older) in the world and this trend will tend to increase in the future, provided that the SARS-CoV-2 pandemic subsides [2,3]. According to some studies, in fact, there is a greater vulnerability of the elderly but also a resilience of centenarians [4,5].
Longevity is rightly considered one of the greatest achievements of modern society. On the topic of longevity and aging, it is important to make a distinction between old, great old, nonagenarians, and centenarians. According to the current literature, the population understand if factors such as altitude, climate, UNESCO protected areas, and hinterland can directly or indirectly influence the measure of healthy living and lead to local longevity.

Description of the Case Study: The Cilento Region
Cilento is an area located in the Campania region, in southern Italy. The territorial area of Cilento includes 102 towns (Figure 1), covering about 490,000 hectares. From the environmental point of view, it is a homogeneous territory, characterized by the integration of different landscape types including coastal, hilly, and mountainous areas, the mouths of important watercourses, and the geomorphological and climatic variability creates a potential not easily found in other areas of the Italian Peninsula. In this area, man has been able to integrate harmoniously with the forms of the landscape [45]. [14,15,44]. Our hypothesis is that some environmental indicators, such as altitude, climate, and rural hinterland may be protective factors for longevity. Therefore, the aim of the present study was to investigate the phenomenon of regional longevity in Cilento by analyzing the association between longevity indicators and some environmental factors, in order to understand if factors such as altitude, climate, UNESCO protected areas, and hinterland can directly or indirectly influence the measure of healthy living and lead to local longevity.

Description of the Case Study: The Cilento Region
Cilento is an area located in the Campania region, in southern Italy. The territorial area of Cilento includes 102 towns (Figure 1), covering about 490,000 hectares. From the environmental point of view, it is a homogeneous territory, characterized by the integration of different landscape types including coastal, hilly, and mountainous areas, the mouths of important watercourses, and the geomorphological and climatic variability creates a potential not easily found in other areas of the Italian Peninsula. In this area, man has been able to integrate harmoniously with the forms of the landscape [45].
With a substantial old population base, Cilento is an ideal region to examine how regional longevity is associated with several natural factors.

Defining the Old Population
Cilento has a population of about 288,185 inhabitants and maintains relatively high proportions of long-lived people (old, nonagenarians, and centenarians) with a ratio of 43.05/1000 inhabitants [46] In this study, six longevity indicators (Centenarian ratio, Aging tendency, 85+ ratio, 90+ ratio, Longevity index, and Centenarity index) were chosen to fully reflect the extreme regional longevity, the longevity of nonagenarians and centenarians within the popula- With a substantial old population base, Cilento is an ideal region to examine how regional longevity is associated with several natural factors.

Defining the Old Population
Cilento has a population of about 288,185 inhabitants and maintains relatively high proportions of long-lived people (old, nonagenarians, and centenarians) with a ratio of 43.05/1000 inhabitants [46] In this study, six longevity indicators (Centenarian ratio, Aging tendency, 85+ ratio, 90+ ratio, Longevity index, and Centenarity index) were chosen to fully reflect the extreme regional longevity, the longevity of nonagenarians and centenarians within the population and the overall old population in Cilento. Two indices were used to reduce the effect of important phenomena such as declining birth rates and migration observed in regional territories: Centenarity index (CI%), given by the ratio of centenarians to those over 90, and Longevity index (LI%), given by the ratio of those over 90 to those over 65 [9]. The longevity indicators used are summarized in Table 1. Note: + = more than.

Environmental Variables Influencing the Structure of Old Age
Environmental factors were selected to understand their influence on regional longevity. The selection of these factors was based on review of relevant literature and data availability. The selected natural indicators include altitude, climate, UNESCO heritage area, and hinterland. Altitude was divided into plain, hill, and mountain based on meters of elevation, according to the altitude bands established by ISTAT [47]. For the climate of Cilento the mean annual temperature in degrees Celsius and the annual precipitation in millimeters (period of analysis year 2020) were analyzed [48], according to the climatic regions identified in Cilento by the directives of the Campania region [49]. In addition, the municipalities of Cilento were distinguished in number of municipalities present in UNESCO protected areas and municipalities located in areas defined as buffer zones because they are not included in the protected area [50]. Finally, the municipalities of Cilento were distinguished into municipalities present in the coastal area and municipalities in the hinterland [47]. These main factors and their subgroup sections are summarized in Table 2. For this study, population data were collected through the Archives of the Cilento municipalities, in the years 2019-2020 and compared with ISTAT [51] (2020) population statistics. Longevity indicators were calculated from the collected population data and used as dependent variables in the study. Environmental indicators, used in the research as independent variables, were collected through: the database of the Italian National Institute of Statistics [47] for altitude and coastal/no coastal areas; data from the Italian Institute for Environmental Protection and Research [48] for climate; and the data from the Italian National Commission for UNESCO [50] for the towns of Cilento belonging to UNESCO.

Statistical Analysis
This study used the STATA software [52] to calculate longevity indicators from the National Institute of Statistics 2020 and Cilento Municipal Archives 2019-2020. Correlations among the six longevity indicators were performed using STATA. QGIS 3.14.15 [53] was then used to produce spatial distribution maps of longevity and environmental variables. Multiple approaches were applied to address the research questions. Global Moran's I index and local Moran's I index (spatial autocorrelation analysis) were used to describe the characteristics of the global and local spatial cluster value of the longevity index, through ArcGIS [54]. An Ordinary Least Squares (OLS) regression and a Geographically Weighted Regression (GWR) model, with ArcGIS software, were used to study the spatial relationship between the explanatory variables and the longevity outcome variables in the given dataset. Finally, we compared the performance between the OLS regression and the GWR model. In the models presented, we eliminated variables that had a VIF value greater than 7.5 because GWR is sensitive to multicollinearity.

Spatial Autocorrelation Analysis
Moran's index (Moran's I) was used to measure the spatial autocorrelation of the indicators and was calculated using the following formula [55]: where n is the number of spatial units indexed by i and j; x is the variable of interest; x i and x j are the values of the observed variable at sites i and j; x is the mean of x; and the weights w ij are written in a (n × n) weight matrix. The weight matrix depicts the relationship between an element and its surrounding elements. Weight can be based on contiguity relationship or distance. The value of Moran's I usually ranges from −1 to +1. A positive Moran's I value denotes positive spatial autocorrelation, whereas a negative value denotes negative spatial autocorrelation. For statistical hypothesis testing, values of Moran's I were evaluated based on Z-score. A |Z| value higher than 1.96 indicates significant spatial autocorrelation at the 0.05 significance level, and a |Z| value higher than 2.58 indicates significant spatial autocorrelation at the 0.01 significance level.
In this study, spatial autocorrelation analysis (Moran's I) between the dependent variables, showed a value ranging from 0.08 to 0.35 and these variables had relatively positive Z and p values.
However, the global spatial autocorrelation index measures the overall spatial autocorrelation of our dataset and does not indicate where clusters of high or low longevity might occur. Therefore, local indicators of spatial autocorrelation are needed.
The local Moran's I [56] equation is as follows: where z i is the value of the variable z at location i; z is the average value of z with the sample number of n; z j is the value of the variable z at all the other locations (where j = i); σ 2 is the variance of variable z; and w ij is a weight which can be defined as the inverse of the distance d ij among locations i and j. The weight w ij can also be determined using a distance band: samples within a distance band are given the same weight, while those outside the distance band are given the weight of 0. A high positive local Moran's I value implies that the location under study has similar high or low values to its neighbors, so the locations are spatial clusters. Spatial clusters include high-high clusters (high values in a high-value neighborhood) and low-low clusters (low values in a neighborhood. A high negative local Moran's I value means that the location under study is a spatial outlier. Spatial outliers are those values that are obviously different from the values of their surrounding localities. Spatial outliers include high-low (a high value in a low-value neighborhood) and low-high (a low value in a high-value neighborhood) [57].

OLS Regression and GWR Model
To quantify the impacts of environmental factors on longevity indicators, we used OLS regression, which is the most commonly used method for analyzing relationships between two or more variables. OLS is a traditional regression method that estimates a global regression coefficient, which is constant over space. The equation for OLS is as follows: where y is the dependent variable, β 0 is the intercept, x i is the ith independent variable, β i is the ith regression coefficient, ε is the error term, and p is the number of independent variables. OLS cannot test spatial effects between observations and does not take into account spatial autocorrelation, so we applied GWR to measure altitude, climate, UNESCO listed area or buffer zones, coastal and inland areas in terms of spatial impact. GWR is an extension of OLS regression that allows locally varying parameters to consider a spatial nonstationary in a sample [58]. The GWR formula is as follows: where y i refers to the dependent variables such as Centenarian ratio, Aging tendency, 80+ ratio, 90+ ratio, Longevity index, and Centenarity index, at location i. (u i , v i ) means the coordinates of the centroid at location i. β 0 (u i , v i ) refers to the intercept for location i. βk(u i , v i ) means the local parameter for independent variable k at location i. x ik is the value of independent variable k at location i. ε i is the error term for location i.
where β(u i , v i ) means the local regression coefficient at location i. X refers to the matrix of the independent variables of the environmental factors. Y is the vector of the dependent variable of longevity indicators.

Comparison of OLS Regression and GWR Models
We compared the model performance between the OLS regression and GWR models based on the coefficient of determination (R 2 ) and corrected Akaike Information Criterion (AICc) values. The R 2 presents the prediction ability of a regression model to fit the measured values of the dependent variable, and the AICc is an indicator of the relative information lost by the model during the estimation process [58].

Evaluation of the Spatial Distribution of Longevity Indicators in Cilento, Italy
The longevity area in Cilento is clearly delineated in the central municipalities of the territory ( Figure 2) and from there, by age groups, it widens towards the south-eastern area. The towns with the highest Centenarian ratio (Figure 2a)   Influenced by migration and birth rates, Aging tendency, 85+ ratio, 90+ ratio, and Centenarian ratio cannot accurately explain the regional distribution of longevity in Cilento. To exclude the influence of migration and birth rates, Centenarity index and Longevity index were introduced [9].
The highest Longevity index (Figure 2e) is held by the central and north-western areas of Cilento, and the reference towns are: Sant'Angelo a Fasanella (17.5), Sessa Cilento Influenced by migration and birth rates, Aging tendency, 85+ ratio, 90+ ratio, and Centenarian ratio cannot accurately explain the regional distribution of longevity in Cilento. To exclude the influence of migration and birth rates, Centenarity index and Longevity index were introduced [9].
The six indicators analyzed all show high levels of longevity in the central and southeastern area while they have lower levels towards the borders with Basilicata and the coastal zone.
Based on the results of the global autocorrelation values of the longevity indicators, we can state that the spatial distribution of the longevity indices in the dataset is spatially clustered more than would be expected if the underlying spatial processes were random, so we can reject the null hypothesis. Indeed, Moran's I Z-score values indicate spatial autocorrelation at the 0.01 significance level for the Aging tendency, the 85+ ratio, and the 90+ ratio, and a 0.05 significance level for the Longevity index. No significance level occurs for the values of the Centenarian ratio and the Centenarity index. The results of the global spatial autocorrelation analysis are reported in Table 3.  Figure 3 shows the local spatial patterns (Anselin local Moran's I) of longevity indicators in Cilento, Italy. The Centenarian ratio clusters were highlighted mainly in the designed area High-High and Low-Low, which are located in the central area of Cilento. Furthermore, due to the Aging tendency, the 85+ ratio, the 90+ ratio, and the Longevity index, the aggregated areas with a high longevity value defined High-High, were distributed mainly in the central and southeastern part of Cilento. The Low-Low areas were grouped into the Northwest and Southwest areas. In addition, the high and low areas were randomly placed in the areas surrounding the High-High and Low-Low regions. These results indicate that longevity is concentrated in the central and southeastern areas of Cilento and may be influenced by surrounding environmental factors.

The Spatial Pattern of Environmental Factors That May Influence the Distribution of Longevity
The spatial distribution of indicators relating to altitude, climate, and UNESCO heritage areas and coastal/non coastal areas, which may have an impact on the distribution of longevity in Cilento, is shown in Figure 4a-d, respectively.
Furthermore, due to the Aging tendency, the 85+ ratio, the 90+ ratio, and the Longevity index, the aggregated areas with a high longevity value defined High-High, were distributed mainly in the central and southeastern part of Cilento. The Low-Low areas were grouped into the Northwest and Southwest areas. In addition, the high and low areas were randomly placed in the areas surrounding the High-High and Low-Low regions. These results indicate that longevity is concentrated in the central and southeastern areas of Cilento and may be influenced by surrounding environmental factors.

The Spatial Pattern of Environmental Factors That May Influence the Distribution of Longevity
The spatial distribution of indicators relating to altitude, climate, and UNESCO heritage areas and coastal/non coastal areas, which may have an impact on the distribution  In Figure 4a we can observe three different altitudes [47]: 1. the plain, located in the coastal area, between the towns: Pollica, Ascea, Pisciotta, San Mauro Cilento, Agropoli, Serre; 2. the hill, which covers most of the Cilento area; 3. the mountain, which is located in the northern and central-eastern area with Teggiano, Piaggine, Monte San Giacomo, Petrucci, Auletta, Caggiano.
One of the most interesting features of Cilento is certainly the geomorphologic and climatic variability hardly found in other areas of the peninsula. Figure 4b highlights the three climatic zones of Cilento: the Mediterranean climate, located in the coastal zone, with hot and dry summers and warm and humid winters, and an average temperature of 22-28 °C; the transitional climate, located in the central, centralwestern, and south-eastern hilly zone, with a temperate climate of transition to the Mediterranean and average temperatures around 20 °C; finally, the temperate climate, located towards the northern and central-eastern mountain zones, with cool summers and mild winters and average temperatures of 15-20 °C. The climatic variability of Cilento, which extends from the hottest and driest sectors to the coolest and most humid ones, is a source of exceptional value for the vegetation [48]. Figure 4c shows that Cilento also falls within the area of the Vallo di Diano National Park, which has been recognized as a UNESCO World Heritage Site for both the natural environmental and cultural value of the area. Therefore, the municipalities of Cilento are divided between the areas included in the UNESCO list and the buffer zones [50].
Finally, Figure 4d shows that the Cilento area is divided into 3/4 for the non-coastal zone and 1/4 for the coastal zone [47]. In Figure 4a we can observe three different altitudes [47]: 1. the plain, located in the coastal area, between the towns: Pollica, Ascea, Pisciotta, San Mauro Cilento, Agropoli, Serre; 2. the hill, which covers most of the Cilento area; 3. the mountain, which is located in the northern and central-eastern area with Teggiano, Piaggine, Monte San Giacomo, Petrucci, Auletta, Caggiano.
One of the most interesting features of Cilento is certainly the geomorphologic and climatic variability hardly found in other areas of the peninsula. Figure 4b highlights the three climatic zones of Cilento: the Mediterranean climate, located in the coastal zone, with hot and dry summers and warm and humid winters, and an average temperature of 22-28 • C; the transitional climate, located in the central, central-western, and south-eastern hilly zone, with a temperate climate of transition to the Mediterranean and average temperatures around 20 • C; finally, the temperate climate, located towards the northern and central-eastern mountain zones, with cool summers and mild winters and average temperatures of 15-20 • C. The climatic variability of Cilento, which extends from the hottest and driest sectors to the coolest and most humid ones, is a source of exceptional value for the vegetation [48]. Figure 4c shows that Cilento also falls within the area of the Vallo di Diano National Park, which has been recognized as a UNESCO World Heritage Site for both the natural environmental and cultural value of the area. Therefore, the municipalities of Cilento are divided between the areas included in the UNESCO list and the buffer zones [50].
Finally, Figure 4d shows that the Cilento area is divided into 3/4 for the non-coastal zone and 1/4 for the coastal zone [47].
Analyzing the spatial pattern of environmental factors that may influence the distribution of longevity, we find that the Cilento areas with higher longevity tend to be located in the hinterland areas, at an altitude between 400 and 700 m.

OLS Regression and GWR Models
OLS regression data and GWR models are reported in Table 4 and Figure 5, respectively. The OLS analysis revealed that longevity indicators in the studied area were dissimilarly influenced by different independent variables. In particular, a negative relationship was found between UNESCO buffer zones and Aging tendency, whereas a positive relationship was found between nonagenarians and the Cilento hinterland. Moreover, as show in Table 4, UNESCO buffer zones also show a negative relationship with the great old and nonagenarians, while the relationship between the longevity index with the Cilento hinterland and altitude emerged. In addition, to test for spatial autocorrelation, Moran's I index was used to identify whether case densities were randomly distributed across the study area. We found that longevity densities were clustered across environmental factors, with altitude (Moran's I = 0.26; z-score = 4.46; p = 0.0018), climate (Moran's I = 0.21; z-score = 3.73; p = 0.0008), UNESCO heritage (Moran's I = 0.27; z-score = 4.60; p = 0.0004), and littoral/non-littoral areas (Moran's I = 0.65; z-score = 10.87; p < 0.0001). For the GWR models, shown in Figure 5, only variables that were significant in the OLS regression were used.

Comparisons of the Performance between the OLS and GWR Models
Comparison of OLS regression and GWR models for all independent variables used in the analysis are shown in Table 5. (c) (d) Figure 5. Results of the GWR models for Aging tendency, 85+ ratio, 90+ ratio and Longevity index showing the local R 2 . Note: (a) local R 2 for Aging tendency; (b) local R 2 for 85+ ratio; (c) local R 2 for 90+ ratio; (d) local R 2 for Longevity index.
In addition, to test for spatial autocorrelation, Moran's I index was used to identify whether case densities were randomly distributed across the study area. We found that longevity densities were clustered across environmental factors, with altitude (Moran's I = 0.26; z-score = 4.46; p = 0.0018), climate (Moran's I = 0.21; z-score = 3.73; p = 0.0008), UNESCO heritage (Moran's I = 0.27; z-score = 4.60; p = 0.0004), and littoral/non-littoral areas (Moran's I = 0.65; z-score = 10.87; p < 0.0001). For the GWR models, shown in Figure  5, only variables that were significant in the OLS regression were used.

Comparisons of the Performance between the OLS and GWR Models
Comparison of OLS regression and GWR models for all independent variables used in the analysis are shown in Table 5.  Table 4; R 2 indicates R-Squared, coefficient of determination; AICc-Akaike's Information Criterion.

Correlation between the Six Longevity Indicators
Analysis of the correlation between the six longevity indicators in Cilento, showed Figure 5. Results of the GWR models for Aging tendency, 85+ ratio, 90+ ratio and Longevity index showing the local R 2 . Note: (a) local R 2 for Aging tendency; (b) local R 2 for 85+ ratio; (c) local R 2 for 90+ ratio; (d) local R 2 for Longevity index.

Correlation between the Six Longevity Indicators
Analysis of the correlation between the six longevity indicators in Cilento, showed that nonagenarians (90+ ratio) have a strong correlation with the Longevity index and the Aging tendency. The great old (85+ ratio) have a correlation with the Longevity index, the Aging tendency and the 90+ ratio. The Centenarian ratio (C/10,000) correlates with the Aging tendency, the 85+ ratio, and the 90+ ratio, the Centenarity index, and the Longevity index. The Centenarity index correlates only with the Centenarian ratio (Table 6).

Discussion
By analyzing the association between six longevity indicators and five environmental factors, this study aimed to investigate the phenomenon of regional longevity in Cilento. Based on our results, we are able to support the hypothesis that healthy living and longevity in Cilento may be influenced directly or indirectly by several environmental factors, such as altitude, climate, UNESCO protected areas, and hinterland. In our opinion, there are at least three findings that support our hypothesis, which can be summarized as follows: First, it has been observed that longevity in Cilento is clearly manifested in central municipal areas of the territory and from there, by age groups, it spreads towards the southeastern area. This result, in agreement with the study of Deng et al., [59] confirms that there are only limited regions that show a high longevity index, and that the phenomenon of longevity is not universal. Herein, we found that centenarians are located in the heart of the clustered area at an altitude between 440 and 600 m above sea level (hilly area). Cilento area is characterized by a transitional climate between Mediterranean and temperate climate, with mild temperatures around 20 • C, fairly wet winters and moderate summer drought. The mild climate seems to be an important factor for the longevity phenomenon in Cilento. The results of our study agree with other studies [60,61] in confirming that the old, great old, nonagenarians and centenarians could benefit from the climate. From a clinical point of view, this could have important implications. Indeed, Mathieu et al., and Wyder [60,61], who mainly traced the genesis of medical interest in altitude climate, showed how mountain climate could play a role as therapeutic agent to reduce the risk of cardiovascular diseases.
The second finding of our study is that the old, great old, and nonagenarians seem to be associated with the municipalities included in the UNESCO "World Heritage List". A potential explanation could be, at least in part, found considering the diet [62] of those populations. In particular, the UNESCO natural area includes a Mediterranean park par excellence where grows a species of olive tree, Olea europaea. This food, with protected designation of origin, contains substances capable of playing an important antioxidant role [63], preventing neoplastic, inflammatory, cardiovascular, and metabolic disorders [64]. Pes et al. [65] demonstrated that a higher intake of olive oil has a beneficial effect on selfperceived health, physical performance, and sense organ function. Regarding Cilento area, Scelzo et al. [14] who studied nonagenarians and centenarians in Cilento, showed that those who lived in the hilly rural area of the territory adopted specific lifestyles, such as hard work, love of the land, family, and religion, which allowed them to maintain mental well-being and made them particularly resilient and optimistic. Moreover, the study by Pizza et al. [15] on Cilento nonagenarians also emphasized the importance of lifestyle, diet, but also individual personality factors. However, it is interesting to note that long-lived populations living in the buffer zone seem to be less protected. Therefore, more studies are necessary to explain this phenomenon.
As third finding, using a correlation analysis between the six longevity indicators, we found a linear continuity of the long-lived population in Cilento, which could ensure generational turnover among the old, the great old, the nonagenarians, and the current centenarians. Based on the results of the global autocorrelation values of the longevity indicators, we found that the spatial distribution of the longevity indices in the dataset resulted spatially clustered more than would be expected, when the underlying spatial processes were random. However, it should be underlined that no single longevity indicator was sufficient to assess regional longevity. In the current study, a complete picture of longevity in Cilento was obtained by using a set of indicators to delineate the different dimensions of longevity.
Finally, all independent variables used in the analysis were compared using OLS regression and GWR models. Herein we found that longevity in Cilento is clearly manifested in the central municipalities of the territory and from there, by age group, it extends to the southeastern area. This result confirms that there are only limited regions that show a high longevity index, and that the longevity phenomenon is not universal [59]. Specifically, we found that longevous are negatively affected by UNESCO buffer zones, but positively affected by the hinterland areas of Cilento where UNESCO protected area municipalities are located. In addition, longevous benefited from the positive influence of altitude between 400 and 700 m. Climate seems to have no effect on longevity. Furthermore, in accordance with other studies, herein we found that GWR models achieved higher R 2 and AICc values than the OLS regression model, indicating better performance of GWR models than OLS regression models [66][67][68].
Herein, longevity indicators in the studied area were dissimilarly influenced by several independent variables. In particular, we found a negative relationship between UNESCO buffer zones and Aging tendency, while a positive relationship between nonagenarians and the Cilento hinterland. Furthermore, UNESCO buffer zones also show a negative relationship with the great old and nonagenarians, while a relationship between the longevity index with the Cilento hinterland and altitude emerged. This result supports our hypothesis that countries included in the UNESCO heritage list represent a possible natural protective factor for longevity.
From a political and economic point of view, considering that an aging society implies an increase in health care spending, it should be mandatory to consider more and more interventions to support healthy aging, selecting those factors that are able to ward off disease and make long-term care policies economically sustainable.
This study has some limitations. First, it analyzed longevity in Cilento by considering only a specific year without examining a longer period of time. Second, it investigates the phenomenon of longevity in Cilento by analyzing only specific environmental factors, such as altitude, climate, UNESCO-protected areas, and hinterland, that could directly or indirectly influence the extent of healthy living and long life in Cilento, but did not consider factors such as drinking water that might have traces of elements beneficial to human health. Third, despite that some environmental factors could be better studied by techniques, such as remote sensing or land-use data, in the present study the input layers for environmental factors were not studied.

Conclusions
To the best of our knowledge, this is the first study to define the area of longevity in Cilento through six longevity indicators and four environmental factors. The results obtained from the study showed that environmental factors, such as hill altitude, transitional climate between Mediterranean and temperate climate, UNESCO heritage areas, and hinterland areas, can contribute directly or indirectly to local longevity in Cilento and centenarians. Further studies are needed to better understand the relationship between health, diet, and living area. In addition, studies related to the association between environmental and genetic factors on local longevity in Cilento should be encouraged and investigated to understand what factors are able to modify the epigenome to establish healthy aging.

Funding:
The authors received no financial support for the research, authorship, and/or publication of this article.
Institutional Review Board Statement: Ethical approval from the Human Research Ethics Committee was not required for this study. Our study exploited the demographic database of Official Archives. This represents a retrospective study and is not directly associated with patients.
Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.

Data Availability Statement:
The data included in this manuscript were provided by the National Institute of Statistics and the Cilento Municipal Archives. Therefore, we are not authorized to share the data with third party organizations. However, the corresponding author is available to provide any explanation to the Editor if requested.