The Necessity of Maintaining the Resilience of Peri-Urban Forests to Secure Environmental and Ecological Balance: A Case Study of Forest Stands Located on the Romanian Sector of the Pannonian Plain

: Climate change’s negative effects, such as rising global temperatures and the disruption of global ecological ecosystems as a direct effect of rising carbon emissions in the atmosphere, are a signiﬁcant concern for human health, communities, and ecosystems. The condition and presence of forest ecosystems, especially those in peri-urban areas, play an essential role in mitigating the negative effects of climate change on society. They provide direct beneﬁts to the residents of large cities and their surrounding areas, and they must be managed sustainably to protect all their component ecosystems. This research was carried out in the forests of Lunca Muresului Natural Park and Bazos Arboretum, located in the Romanian sector of the Pannonian Plain, near urban agglomerations. The results showed high variability in the stands. Using the height-to-diameter ratio indicator concerning dbh and species, a strong Pearson correlation was registered (between 0.45 and 0.82). These values indicate the high stability of these stands, providing positive human–nature interactions such as recreational or outdoor activities (and a complementary yet indirect use value through attractive landscape views). Protecting these ecosystems offers a so-called insurance policy for the next generations from a climate change standpoint.


Introduction
Nowadays, one of the most critical environmental concerns of society is the increasing impact of greenhouse gas emissions, with the most essential being CO 2 , which subsequently influences the greenhouse effect in the atmosphere [1].Many consider it climate change [2][3][4] or global warming [5].As Europe has been projected to be the third-largest carbon emitter in the world [6] and forests act as carbon sinks [7], the paradigm is seen to shift slowly but steadily toward a socio-climate focused on national and international policies with increasing pressure on forests in many cases [8][9][10][11].While climate change has been seen to influence drastic changes in current policy frameworks, policy-makers have also considered other factors, such as urban expansion, migration, and the increasing need for sustainable resource provision and agriculture [12,13].As migration and urban expansion increase so does the need for more food and resources; thus, forests are faced with an old yet controversial problem of conversion to the agricultural land [14,15] more commonly seen in plain areas.
Forests absorb and store CO 2 through photosynthesis, reducing the amount of greenhouse gases in the atmosphere and mitigating the effects of climate change.Forests also store carbon in the form of wood, leaves, and other organic matter, helping to keep carbon locked away for long periods.When forests are cut down or burned, the carbon stored in the trees and other vegetation is released into the atmosphere, contributing to the buildup of greenhouse gases.Thus, forests can also act as a source of greenhouse gas emissions when disturbed through deforestation, forest fires, and other land use changes.Therefore, it is essential to manage forests sustainably to maximize their ability to act as a sink for greenhouse gases and minimize their contribution as a source.In this context, a possible solution for the sustainable management of forest ecosystems and the fulfillment of global policy guidelines resides in more effective management of forest areas [16][17][18].Last but not least, forest management must also consider preserving the beneficial effects of forest ecosystems [19,20], such as clean water or recreation.
Forests with special protection functions, in some cases, predominantly social or water, land, and soil protection, have an essential role in protecting forest ecosystems [21,22] while also securing material resources and exercising the provision function of a forest ecosystem [23,24].According to the forest legislative framework in Romania, only conservation cuttings are allowed in these areas.To a certain extent, without influencing the main protective role of the forest, these types of cuttings are defined by certain interventions, which are applied in forest stands, that fulfill special protection functions and aim to allow the establishment of a new tree group using low-intensity actions with a maximum of 10-15% of standing volume per decade [25].Since 1954, Romania's area of forests with special protection functions has increased from 12.7% of national forest land to approximately 66% in 2021.
Natural resources suffer the most in the face of urban land expansion into the surroundings mainly due to being often affected or displaced.Moreover, urban expansion is believed to directly influence the transformation or loss of forest areas and their management.While urbanization near forested areas increases, it affects the lives of most people in the transition, as forests act as a home for social and ecological services that improve human health and well-being [26].At the same time, Henwood et al. [27] suggested that management can evolve into ensuring extended benefits regarding trees, woods, and forests which society might need at a personal level with deep implications toward communities and the cultures in which they live.
Addressing the protection of stable forests, measurable adequate evidence-based support is believed to enhance forest stability, which is similar to having a life insurance policy providing certain security for future generations [28].Forests can partially mitigate natural hazards, such as erosion, floods, or rock falls [29]; out of these, river and coastal floods are the most critical hazards in many floodplains and coastal stretches of Central and Eastern Europe, with an expected annual damage (EAD) in Romania of USD 289 million up to 2050 [4].These natural hazards can be predicted by analyzing tree and stand stability, and forest managers can take appropriate measures to mitigate or prevent such events [30][31][32].Therefore, stand stability is a critical characteristic [33,34] regarding forest resilience.The research in this area focused on a clear correlation and prediction of height-to-diameter ratio (h/d ratio), in which species mixtures were also considered to provide insights into the overall stability generated by human-forest interaction [35,36].Recent research [35] has shown that this relationship between diameter and height is an important factor in forest research.Through this relationship, stand volumes and biomass can be estimated [37][38][39], but it can also give information on tree and stand growth [40], as well as information on stand productivity [41].
The aim of this study was to highlight the importance of peri-urban forests and their role in society.The objectives of this research were (i) to analyze the stand structure of several peri-urban forests located in the eastern sector of the Pannonian Plain in Romania; (ii) to highlight resilience through the stability of these forest ecosystems; and (iii) to generally assess the urban land expansion in the region and its implications for current and future generations.

Materials and Methods
This research was conducted in the Lunca Muresului Natural Park and the Bazos Arboretum (Figure 1) located in the western part of Romania on the Western Plains, which is the Romanian sector of the Pannonian Plain.The study area consists of forests located in the peri-urban region, composed of common oak (Quercus robur sp.) with special protection interest and scientific interest in conserving the gene pool and other ecosystems characterized by natural elements with special value.
Diversity 2023, 15, x FOR PEER REVIEW 3 of 14 The aim of this study was to highlight the importance of peri-urban forests and their role in society.The objectives of this research were (i) to analyze the stand structure of several peri-urban forests located in the eastern sector of the Pannonian Plain in Romania; (ii) to highlight resilience through the stability of these forest ecosystems; and (iii) to generally assess the urban land expansion in the region and its implications for current and future generations.

Materials and Methods
This research was conducted in the Lunca Muresului Natural Park and the Bazos Arboretum (Figure 1) located in the western part of Romania on the Western Plains, which is the Romanian sector of the Pannonian Plain.The study area consists of forests located in the peri-urban region, composed of common oak (Quercus robur sp.) with special protection interest and scientific interest in conserving the gene pool and other ecosystems characterized by natural elements with special value.According to their geographical position, these forest stands are found in the area of the continental climate, with sub-Mediterranean influence, where winters are milder, shorter, and relatively cold, while the summers are humid and long.According to the Kὅppen Climatic Classification, they are located in the Cfax climate type and described by a mild, humid climate with rainfall all year, where temperatures in the hottest month of the year (June) exceed 22 °C, and maximum rainfall is recorded at the end of spring with the minimum rainfall falling at the end of winter.
The humidity index has an average value of 60-80%, with lower values recorded during fall and winter, while higher values were reported in spring and summer.The average temperature of the year ranged between 9 and 11 °C , with amplitudes rising steadily to 22-23 °C.The average temperature in January exceeds −2 °C , while the average in July reaches 20-21 °C ; these temperatures are characteristic for the forest-steppe climate.
The forest stands in the areas studied are included in zones with special protection functions in which only special conservation and sanitary cuttings are allowed to some extent.The management of these types of forests is focused on the conservation of the living and nonliving ecosystem components, which results, in general, to less human influence over time.According to their geographical position, these forest stands are found in the area of the continental climate, with sub-Mediterranean influence, where winters are milder, shorter, and relatively cold, while the summers are humid and long.According to the K Diversity 2023, 15, x FOR PEER REVIEW 3 of 14 The aim of this study was to highlight the importance of peri-urban forests and their role in society.The objectives of this research were (i) to analyze the stand structure of several peri-urban forests located in the eastern sector of the Pannonian Plain in Romania; (ii) to highlight resilience through the stability of these forest ecosystems; and (iii) to generally assess the urban land expansion in the region and its implications for current and future generations.

Materials and Methods
This research was conducted in the Lunca Muresului Natural Park and the Bazos Arboretum (Figure 1) located in the western part of Romania on the Western Plains, which is the Romanian sector of the Pannonian Plain.The study area consists of forests located in the peri-urban region, composed of common oak (Quercus robur sp.) with special protection interest and scientific interest in conserving the gene pool and other ecosystems characterized by natural elements with special value.According to their geographical position, these forest stands are found in the area of the continental climate, with sub-Mediterranean influence, where winters are milder, shorter, and relatively cold, while the summers are humid and long.According to the Kὅppen Climatic Classification, they are located in the Cfax climate type and described by a mild, humid climate with rainfall all year, where temperatures in the hottest month of the year (June) exceed 22 °C, and maximum rainfall is recorded at the end of spring with the minimum rainfall falling at the end of winter.
The humidity index has an average value of 60-80%, with lower values recorded during fall and winter, while higher values were reported in spring and summer.The average temperature of the year ranged between 9 and 11 °C, with amplitudes rising steadily to 22-23 °C.The average temperature in January exceeds −2 °C, while the average in July reaches 20-21 °C; these temperatures are characteristic for the forest-steppe climate.
The forest stands in the areas studied are included in zones with special protection functions in which only special conservation and sanitary cuttings are allowed to some extent.The management of these types of forests is focused on the conservation of the living and nonliving ecosystem components, which results, in general, to less human influence over time.
ppen Climatic Classification, they are located in the Cfax climate type and described by a mild, humid climate with rainfall all year, where temperatures in the hottest month of the year (June) exceed 22 • C, and maximum rainfall is recorded at the end of spring with the minimum rainfall falling at the end of winter.
The humidity index has an average value of 60-80%, with lower values recorded during fall and winter, while higher values were reported in spring and summer.The average temperature of the year ranged between 9 and 11 • C, with amplitudes rising steadily to 22-23 • C. The average temperature in January exceeds −2 • C, while the average in July reaches 20-21 • C; these temperatures are characteristic for the forest-steppe climate.
The forest stands in the areas studied are included in zones with special protection functions in which only special conservation and sanitary cuttings are allowed to some extent.The management of these types of forests is focused on the conservation of the living and nonliving ecosystem components, which results, in general, to less human influence over time.
In order to conduct this study, four sample plots of one hectare each (100 m × 100 m) were installed and used to measure forest stands.The research plots were installed between 2014 and 2018 and are located in areas with no slope and at an altitude between 107 and 111 m above sea level.The dominate soil type is preluvosols, and over time, it has not been affected by gelling processes following the drop in groundwater levels.We collected information about the tree species, as well as the position and shape of the trees; biometric characteristics, such as diameter at breast height (dbh) of all the trees that exceed 2 cm; and the total height and assessed wood quality class of the tree.In order to assess the wood quality, 4 quality classes were taken into consideration (I, II, III, and IV) in relation to the percentage of industrial wood from total height of the tree (I > 0.5%; II = 0.25-0.50%;III = 0.1-0.25%;IV < 0.1%).
The total above ground tree volume was determined (including branches above five centimeters) using Equation (1) [42][43][44].Moreover, the total stand volume was determined by adding the individual tree volumes: log v = a 0 + a 1 log d +a 2 log d 2 + a 3 log h + a 4 log h 2  (1) where h represents the tree's total height expressed in meters, d represents the dbh expressed in centimeters, v represents the volume of the tree expressed in cubic meters, and a 0 , a 1 , a 2 , a 3 , and a 4 are regression coefficients customized for each species from yield tables [42,45].
To summarize the samples and analyses to be carried out, the main statistical parameters of stand characterization were determined using the R PASTECS package from R software version 1.3.21[46].
In order to evaluate the forest stand stability and adapt current forest management to the concerns and needs of society, an analysis of the optimal structures of the stands in the areas studied is required to establish the most appropriate theoretical functions to describe these stands.This was undertaken using the theoretical functions Normal [47], the 2-parameter Gamma [47], and Weibull [48] because they are characterized by high flexibility and accuracy in describing the stands [49].The goodness of fit was tested using the chi-square criterion and the Kolmogorov-Smirnov [50] test.
The graphical analyses were performed using R Studio's plot function [51].

Descriptive Statistics
The statistical indicators reveal that in all four plots, there is a relatively large variation, evidenced by the values of standard deviation, variance, and coefficient of variation (Table 1).In all four research plots, the number of trees (per hectare) varies from 203 (Bazos) to 722 (Bezdin); most of them contain more than 607 trees, except for the Bazos research plot where the stand is older (160 years old according to the management plan), larger in size, and has fewer trees.The minimum diameter varies from 2.35 (Ceala plot) to 11.50 cm (Bazos plot), the maximum diameter varies from 55.00 (Popin plot) to 108.00 cm (Bazos plot), and the average diameter ranges from 18.20 (Ceala plot) to 47.07 cm (Bazos plot).The standard deviation ranged from 11.49 (Popin plot) to 22.69 (Bazos plot), and the minimum coefficient of variation was 55% (Popin plot), and the maximum was 98% (Bazos plot).The variance values ranged from 132.08 (Popin plot) to 514.97 (Bazos plot).

Fitting of Experimental Diameter Distribution
The experimental and theoretical diameter distributions of the four research surfaces are graphically represented.From Figure 2, it can be seen that the four stands possess different structures, according to this analysis.Using theoretical functions with a high degree of diversity, it was possible to capture as well as possible the shape of the experimental distribution of the stands, while it could be seen that the structure of none of the stands was described by the theoretical normal function, showing that the stands do not have the characteristics of an even-aged stand.

Fitting of Experimental Diameter Distribution
The experimental and theoretical diameter distributions of the four research surfaces are graphically represented.From Figure 2, it can be seen that the four stands possess different structures, according to this analysis.Using theoretical functions with a high degree of diversity, it was possible to capture as well as possible the shape of the experimental distribution of the stands, while it could be seen that the structure of none of the stands was described by the theoretical normal function, showing that the stands do not have the characteristics of an even-aged stand.Using goodness-of-fit tests, the results obtained indicate that following the application of the chi-square criterion within the Bazos and Popin plots, there were no significant differences between the experimental and theoretical values (Table 2).Applying the same significance test for the Bezdin and Ceala plots, significant differences were recorded only for the normal distribution.No significant differences were recorded for the other distributions (Weibull and Gamma).
Applying the Kolmogorov-Smirnov significance test to all the plots for standardized data, significant differences between the theoretical and experimental values were recorded.However, it could be observed that the differences between the experimental and theoretical values were lower for the Weibull and Gamma functions, but it should be mentioned that the Kolmogorov-Smirnov test is much more rigorous than the chi-square criterion.Using goodness-of-fit tests, the results obtained indicate that following the application of the chi-square criterion within the Bazos and Popin plots, there were no significant differences between the experimental and theoretical values (Table 2).Applying the same significance test for the Bezdin and Ceala plots, significant differences were recorded only for the normal distribution.No significant differences were recorded for the other distributions (Weibull and Gamma).Applying the Kolmogorov-Smirnov significance test to all the plots for standardized data, significant differences between the theoretical and experimental values were recorded.However, it could be observed that the differences between the experimental and theoretical values were lower for the Weibull and Gamma functions, but it should be mentioned that the Kolmogorov-Smirnov test is much more rigorous than the chi-square criterion.

Stand Stability Height-to-Diameter Ratio (h/d Ratio) in Relation to Diameter and Species
Our analysis revealed that the description of the curve between h/d ratio in relation to diameter and species for the four survey areas shows a decreasing trend (Figure 3).In conjunction with the value of the correlation coefficient with values between 0.45 (Bezdin) and 0.82 (Bazos), it can be affirmed, with some certainty, that there is a high to very high correlation between the mentioned characteristics for almost all stands.An exception is the Bezdin plot, where ash (Fraxinus excelsior sp.) and elm (Ulmus minor sp.) trees (mixture species) are more distributed in the lower diameter categories, and only oak trees are found in the upper diameter categories, which are older, and as they were the first to settle in the research area, they manage to overlap other species.

Stand Stability Height-to-Diameter Ratio (h/d Ratio) in Relation to Wood Quality
For all the research areas, the comparative analyses between the h/d ratio and wood quality were performed (Figure 4) and showed a very close relationship between the two stand characteristics, expressed by correlation coefficient values ranging between 0.975 and 0.977.In addition, the results revealed that for all research plots, a higher wood quality was recorded for lower values of the h/d ratio.A small exception was observed in the Ceala research plot, where the fourth (lower) quality class recorded slightly lower h/d ratio values compared to the third (middle to lower) quality class, which is largely due to the species in the stand composition, which has lower longevity (e.g., Populus sp.).
to diameter and species for the four survey areas shows a decreasing trend (Figure 3).In conjunction with the value of the correlation coefficient with values between 0.45 (Bezdin) and 0.82 (Bazos), it can be affirmed, with some certainty, that there is a high to very high correlation between the mentioned characteristics for almost all stands.An exception is the Bezdin plot, where ash (Fraxinus excelsior sp.) and elm (Ulmus minor sp.) trees (mixture species) are more distributed in the lower diameter categories, and only oak trees are found in the upper diameter categories, which are older, and as they were the first to settle in the research area, they manage to overlap other species.

Stand Stability Height-to-Diameter Ratio (h/d Ratio) in Relation to Wood Quality
For all the research areas, the comparative analyses between the h/d ratio and wood quality were performed (Figure 4) and showed a very close relationship between the two stand characteristics, expressed by correlation coefficient values ranging between 0.975 and 0.977.In addition, the results revealed that for all research plots, a higher wood quality was recorded for lower values of the h/d ratio.A small exception was observed in the Ceala research plot, where the fourth (lower) quality class recorded slightly lower h/d ratio values compared to the third (middle to lower) quality class, which is largely due to the species in the stand composition, which has lower longevity (e.g., Populus sp.).
It can also be added that the h/d ratio values recorded for wood quality class 1 (superior) ranged from 50 (Bazos plot) to 76 (Popin plot).For the fourth (lower) quality class,

Discussion
The analysis of the main statistical indicators revealed high variability in diameters in all four research plots.Hence, the values of the coefficient of variation of the diameters were recorded to be 20-30% higher than the specific values for the even-aged stands [52].At the same time, the amplitude between minimum and maximum dbh indicates an une- It can also be added that the h/d ratio values recorded for wood quality class 1 (superior) ranged from 50 (Bazos plot) to 76 (Popin plot).For the fourth (lower) quality class, the minimum value recorded for this indicator was 81 (Ceala plot), and the maximum was 101 (Popin plot).The optimal values recorded for this characteristic of the researched stands ranged from 70 (Bazos plot) up to 85 (Bezdin and Popin plots).

Discussion
The analysis of the main statistical indicators revealed high variability in diameters in all four research plots.Hence, the values of the coefficient of variation of the diameters were recorded to be 20-30% higher than the specific values for the even-aged stands [52].At the same time, the amplitude between minimum and maximum dbh indicates an unevenaged structure of these stands [53][54][55].It should be added that uneven-aged stands are generally known to confer a high degree of stability compared to even-aged stands [56], thus providing protection against natural hazards [4] and increasing the resilience of forests overall.
By fitting experimental distributions of diameters, a descending trend of diameters was observed, with most of them found in the small diameter categories.This pattern is characteristic of uneven-aged stands [57].Further confirmation of the uneven-aged structure is the fact that the experimental distribution of diameters differed significantly from the theoretical Normal distribution, which is considered an indicator of the unevenaged structure [58].
As Weibull distribution is known to provide a more flexible diameter distribution description [59], it has been possible to capture the characterization of the trees in the stands in more detail.Despite the fact that the results showed differences in the null hypothesis test using the Kolmogorov-Smirnov test, it should be pointed out that this test is more demanding than the chi-square criterion [60].At the same time, the differences between the theoretical and experimental values recorded were not very large.However, the tree diversity recorded is high, highlighting their importance, mostly because it is known that urban and peri-urban forests with low tree diversity have a high risk of potential ecological disturbances [61].
Even though the term stability has no generally accepted definition, authors such as Grimm and Wissel [62] suggested that constancy, resilience, and persistence are the main stability properties, whereas [63] stressed that the vital stability properties in protective forests are resistance and elasticity and looked to assess the forests' stability properties through such forest responses to disturbances using silvicultural means.
A tree shape, expressed by the h/d ratio, is a good indicator of tree susceptibility to windthrow and snowthrow [64,65] and has been successfully applied in the past for this description [66].As in other studies [67,68], a trunk shape in the surveyed plots was influenced by silvicultural treatment.
Stand stability plays an important role in the evolution of forest ecosystems, and each tree contributes to the stability of the whole stand.The elements that significantly influence stand stability include component species, growth stage, density, structure, heightto-diameter ratio (h/d ratio), wood quality, phytosanitary conditions, and the presence of wounds caused by biotic factors [34,69].Recent research has shown a close relationship between two of these elements, especially the height-to-diameter ratio (h/d ratio) and wood quality [67,70,71].
These elements are easy to obtain, especially the height-to-diameter ratio (h/d ratio), and can also provide valuable information on stand stability [67,70] and, for managed stands, indicate the appropriate application of silvicultural treatments [72].
However, perceiving stability as a general ecosystem function might raise some issues regarding forest management due to its diverse implications.In addition to the stability they offer, the peri-urban forests investigated offer direct use value (recreation) and indirect use value (scenic views) [73].Thus, stability becomes a factor of specific forest characteristics of high importance to society, acting toward protection against various drivers of disturbances, such as changes in the structure and quality of the environment and issues related to biotic factors [4].
The stability of trees and forests near urban areas comes as an informal request of society to access the forest for different recreational reasons safely.In Eastern Europe, they not only exercise their perceived rights but act as a consequence of changes in public access to forests since the dawn of communism and the modern influences of the EU policies.Romania, along with Croatia and Serbia, had no regulations during the 1990s regarding public access to private forests, yet after the communist era and since 2008, private forest owners in Romania gained the legal right to exclusion, prohibiting public access to private forests, even on forest roads [74].However, among the other rights of owners (access, withdrawal, management, exclusion, and alienation), exclusion rights are the only real freedom that the private owner has; this issue is stressed by Nichiforel et al. [75], suggesting that private owners are becoming creative through rent-seeking activities which imply innovation in the case of Romania and are indirectly producing "rights" for themselves.
Nevertheless, this study is located in a state-owned forested area under protection; the public access is open to a certain extent and is restricted only under extreme conditions (forest fires, floods, wind throws, etc.).Whilst private owners are seeking opportunities to exercise and extend their rights, the state provides ecosystem services, such as recreational, aesthetic, and other similar societal needs, free of cost in most cases, based on the services provided.Whether private owners would take the opportunity and set aside part of their forest due to increased urbanization and societal needs is problematic.Setting aside part of one's property requires incentives under the form of compensations, whereas Romania's compensation schemes are highly regulated, overly bureaucratic, and time-consuming, and forest owners are becoming reticent to such voluntary actions.
Interestingly enough, by forecasting global urban expansion, Seto et al. [76] generally predicted for Europe that the average city would have an annual urban expansion rate of 2.50%, with around 86% attributed to GDP (gross domestic product) per capita growth and 4% attributed to population growth.Globally, the above-mentioned "average" city in the recorded study exhibited an urban population growth rate of 2.18% and an urban land expansion growth rate of 4.84%.This indicates that each city in the study added almost 46,000 urban dwellers per year and approximately 13.5 km 2 of new urban land.According to the Romanian National Institute of Statistics database [77], the overall population growth in urban areas in the studied regions decreased by more than 1% on average, while in rural areas, it increased by 1.5% on average.Regarding the dwellings nearby cities, the overall population growth in rural areas in a 10-year span exceeds 8%.However, in Romania, the 39% growth in GDP per capita in the Western region since 2009 [78] might have influenced urban expansion into rural areas.
Taking into account that forests are constantly mentioned within the ecological development discourse, specifically for rural areas with urban characteristics, the role of forests becomes more diverse rather than just production or afforestation/reforestation.Near urban areas, forestry falls out of its traditionalism and purposes and reinvents itself as urban forestry, where the quality of life aspects shade out the forest impact aspects [79].
The preservation of these ecosystems and the role that their stability plays in providing a so-called life insurance policy for future generations must be considered, mainly due to the low predictability of natural hazards (e.g., windthrow and snowthrow); however, an increase in tree and stand stability can mitigate the impact of natural hazards.By acquiring this knowledge, forest managers can adapt their management to mitigate the impact of such natural hazards.

Conclusions
Seen as a significant challenge to countries all over the globe, greenhouse gas (GHG) emissions have gained increased attention over recent decades due to the social problems caused, such as climate change, higher sea levels, and certain fauna extinction.Researchers are trying to find practical solutions to limit their negative effects and impacts on society while, at the same time, addressing the challenges raised by recent greenhouse gas policies on a national and international level.
In the past few years, climate change has gained the notoriety of a "climate crisis"; it became a social issue due to its social impact, with the right to access forests being highly considered in policy formulation.The ongoing increase in large-scale droughts, numerous forest disturbances, and, more importantly, land use changes has led to severe human migration, creating or expanding communities.Therefore, the management of these ecosystems must be adaptive and, at the same time, long-term.
Urbanization and land expansion have led to and will intensify migration from urban to rural areas.As in other European countries, in Romania, urban areas are expected to become larger and overpopulated, and it is believed this will harm people's well-being and health.
Altering the dynamics of human-nature interactions, urbanization influences changes in spatial (overpopulation in cities), temporal (increasing negative interactions over time), and socio-economic inequalities (increased positive interactions among wealthier individuals).The population growth in rural areas comes as a reaction to increased urbanization, and based on the growth in GDP per capita, the environment is faced with much more pressure, especially in terms of ecological and social functions, as households from more socio-economically advantaged backgrounds are expected to make the transition to rural areas and are usually known to be involved in solving environmental and social issues in the local community and ask for their rights by law.
Due to the increase in human-nature interactions in the past decade, beneficial environmental effects are commonly sought in the form of resources, protection from predators, stability and resilience to disturbances, well-being, health, and safety.However, these benefits usually come with a cost, as in the case of recreational activities in protected areas, which can exert negative impacts on local biodiversity.The future management of these peri-urban forests must consider their conservation status and the expansion of these areas where possible concerning the communities.Thus, it implies the establishment of integrated forest management where all ecosystem services are balanced through a participatory approach and fair involvement, satisfying all user needs and mitigating conflicts.

Figure 3 .
Figure 3. Diameter distribution in relation to species and height-to-diameter ratio.

Figure 3 .
Figure 3. Diameter distribution in relation to species and height-to-diameter ratio.

Figure 4 .
Figure 4. Relationship between h/d ratio and wood quality.

Figure 4 .
Figure 4. Relationship between h/d ratio and wood quality.

Table 1 .
Descriptive statistics of diameters.

Table 2 .
Results of applying statistical tests to the experimental distribution of diameters.

Table 2 .
Results of applying statistical tests to the experimental distribution of diameters.