Regionalizing Aquatic Ecosystems Based on the River Subbasin Taxonomy Concept and Spatial Clustering Techniques
Abstract
:1. Introduction
2. Material and Methods
2.1. Study Area
2.2. Technical Flow
2.3. Regionalizing Objectives and Principles
- Principle of including land area. Terrestrial ecosystems, climate, geology, soil and other natural conditions as well as human activities in the watershed were the most important factors influencing or determining the composition, structure, pattern, process and function of aquatic ecosystems. In the watershed hydrological processes, a variety of nutrients and pollutants were transported into water and then affected the structure and function of aquatic ecosystems. That is to say, the watershed or subwatershed characteristics could control or influence the aquatic life in rivers, streams and other types of water. Therefore, the land area in subwatershed should be included as part of aquatic ecoregion.
- Principle of comprehensiveness and dominance. An aquatic ecoregion should not be delineated based only on some of the aquatic ecological components or their driving factors, but rather the comprehensive characteristics of the aquatic ecosystem or their driving factors. In the process of ecoregion delineation, the pattern of various aquatic ecosystem components as well as their similarities and differences of comprehensive features of aquatic ecosystem or their driving factors must be taken fully into account. Based on the comprehensive analysis, the dominant factors influencing the spatial differentiation of aquatic ecosystem should also be considered.
- Principle of conjugation. The boundaries of the same level aquatic ecoregions did not intersect each other, and the boundaries of adjacent aquatic ecoregions had no space left, and their relationships were seamless and continuous in space. That is to say, each aquatic ecoregion was a complete unit and there was no separation and overlap between each other.
- Principle of hierarchy. Aquatic ecoregions should have multiple levels such as from level I to level n, and be organized in a hierarchical framework and operated at different spatial scales. The watershed aquatic ecosystem should be delineated into multi-level aquatic ecoregions. High-level aquatic ecoregions should contain low-level ecoregions and low-level aquatic ecoregions should be embodied in high-level ecoregions.
2.4. Regionalizing Indicators
2.5. Regionalizing Using Spatial Clustering Technique
- Delineating Subwatersheds. Subwatersheds were used as basic clustering units for delineating level I and II aquatic ecoregions in the Taihu Lake watershed. The SRTM DEM data with pixel spatial resolution of 90 m was used to delineate the subwatersheds in the Taihu Lake watershed. A watershed analysis on the terrain model for the Taihu Lake watershed was performed to generate data on flow direction, flow accumulation, stream definition, stream segmentation and watershed delineation using hydrology analysis tool. After the above several processing steps, sketch maps of subwatersheds were obtained. However, the subwatersheds delineated by the hydrology module were inconsistent with the actual stream networks, so further manually-assisted modification was carried out based on the distribution of DEM and the actual river networks. Thus 1,107 small subwatersheds were obtained based on surface drainage patterns.
- Mapping the Spatial Distribution of Regionalizing Indicators. The 1,107 subwatersheds in the Taihu Lake watershed were used as basic calculation units for each indicator. Based on the raster data of regionalizing indicators, the spatial distribution maps of the six key regionalizing indicators, including elevation, drainage density, percent of construction land area, percent of cultivated land area, soil type and slope, were produced using the Zonal Statistics Tool, which calculated statistics on values of raster data within each subwatershed. Figure 2 represents the spatial distribution maps of average elevation, drainage density, percent of construction land area, percent of cultivated land area, soil type and slope of 90 × 90 m raster units and 1,107 subwatersheds in the Taihu Lake watershed.
- Spatial Clustering and Manual-Assisted Optimization. Spatial clustering, which groups similar spatial objects into classes, is an important component of spatial data mining, and it can be used in the identification of areas of similar land usage in an Earth observation database or in merging regions with similar weather patterns, etc. [27]. Spatial clustering exceeds the ability of traditional multivariate cluster analysis technique. The two-step cluster method is a scalable cluster analysis algorithm designed to handle very large data sets, and it can handle both continuous and categorical variables. The two-step cluster analysis approach has the advantage of automatically determining the optimal number of clusters according to the clustering criterion with a rapid computation speed and less subjectivity and randomness. The two-step cluster analysis algorithm contains two stages: (1) preclustering and (2) hierarchical clustering. The precluster stage groups the respondents into several small clusters. The cluster stage uses the small clusters as input and groups them into larger clusters.
2.6. Aquatic Life Survey and Evaluation
3. Results
3.1. Aquatic Ecoregions Scheme
3.2. Aquatic Ecoregions Characteristics
- Aquatic Life Characteristics. To further validate and assess the scientific, reliability and validity of regionalizing aquatic ecosystems, various indicators of aquatic life, including the dominant species of fish, benthic density, biomass, dominant species, Shannon-Wiener diversity index, Margalef species richness index, Pielou evenness index and ecological dominance, were used to quantitatively compare the spatial differentiations across aquatic ecosystems in different level I aquatic ecoregions and level II aquatic subecoregions.
4. Discussion and Conclusions
- The level I aquatic ecregions and level II aquatic subecregions in the Taihu Lake watershed were delineated using the 1,107 subwatersheds as the basic clustering units based on the principle of including land area from the view of effects of driving factors. It reflected the land surface physical processes influencing the water quality and aquatic life under the action of the overland flow. The regionalizing method based on spatial clustering technique was feasible in the delineation of level I aquatic ecregions and level II aquatic subecregions in the Taihu Lake watershed.
- Using the river subbasin taxonomy concept and spatial clustering approach to delineate the level I aquatic ecregions and level II aquatic subecregions in the Taihu Lake watershed was operable and acceptable. In the delineation process, the impacts of people’s subjective experiences or ideas were minimized to the greatest extent. Compared with the traditional regionalizing method based on indicators’ GIS overlay analysis, the proposed regionalizing method using spatial clustering technique based subwatersheds had the advantages of convenience and automation.
- The Taihu Lake watershed could be delineated into two level I aquatic ecoregions, including Ecoregion I1 and Ecoregion I2, and five level II aquatic subecoregions, including Subecoregion II11, Subecoregion II12, Subecoregion II21, Subecoregion II22 and Subecoregion II23.
Acknowledgments
- Conflict of InterestThe authors declare no conflict of interest.
References
- Crowly, J.M. Biogeography. Canadian Geographer 1967, 11, 312–326. [Google Scholar]
- Bryce, S.A.; Clarke, S.E. Landscape-level ecological regions: Linking state-level ecoregion frameworks with stream habitat classifications. Environ. Manag. 1996, 20, 297–311. [Google Scholar]
- Olson, D.M.; Dinerstein, E.; Wikramanayake, E.D.; Burgess, N.D.; Powell, G.V.N.; Underwood, E.C.; D’amico, J.A.; Itoua, I.; Strand, H.E.; Morrison, J.C.; Loucks, C.J.; Allnutt, T.F.; Ricketts, T.H.; Kura, Y.; Lamoreux, J.F.; Wettengel, W.W.; Hedao, P.; Kassem, K.R. Terrestrial ecoregions of the world: A new map of life on earth. BioScience 2001, 51, 933–938. [Google Scholar]
- Wikramanayake, E.; Dinerstein, E.; Loucks, C.; Olson, D.; Morrison, J.; Lamoreux, J.; McKnight, M.; Hedao, P. Ecoregions in ascendance: Reply to Jepson and Whittaker. Conservat. Biol. 2002, 16, 238–243. [Google Scholar]
- Hessburg, P.F.; Salter, R.B.; Richmond, M.B.; Smith, B.G. Ecological subregions of the interior Columbia Basin, USA. Appl. Veg. Sci. 2000, 3, 163–180. [Google Scholar]
- McMahon, G.; Gregonis, S.M.; Waltman, S.W.; Omernik, J.M.; Thorson, T.D.; Freeouf, J.A.; Rorick, A.H.; Keys, J.E. Developing a spatial framework of common ecological regions for the conterminous United States. Environ. Manag. 2001, 28, 293–316. [Google Scholar]
- Clarke, S.E.; White, D.; Schaedel, A.L. Oregon, USA, ecological regions and subregions for water quality management. Environ. Manag 1991, 15, 847–856. [Google Scholar]
- Jenerette, D.; Lee, J.; Waller, D.; Carlson, R.E. Multivariate analysis of the ecoregion delineation for aquatic systems. Environ. Manag 2002, 29, 67–75. [Google Scholar]
- Abell, R.; Thieme, M.L.; Revenga, C.; Bryer, M.; Kottelat, M.; Bogutskaya, N.; Coad, B.; Mandrak, N.; Balderas, S.C.; Bussing, W.; Stiassny, M.L.J.; Skelton, P.; Alleb, G.R.; Unmack, P.; Naseka, A.; Ng, R.; Sindorf, N.; Robertson, J.; Armijo, E.; Higgins, J.V.; Heibel, T.J.; Wikramanayake, E.; Olson, D.; López, H.L.; Reis, R.E.; Lundberg, J.G.; Pérez, M.H.S.; Petry, P. Freshwater ecoregions of the world: A new map of biogeographic units for freshwater biodiversity conservation. BioScience 2008, 58, 403–414. [Google Scholar]
- Loveland, T.R.; Merchant, J.M. Ecoregions and ecoregionalization geographical and ecological perspectives. Environ. Manag. 2004, 34, s1–s13. [Google Scholar]
- Spalding, M.D.; Fox, H.E.; Allen, G.R.; Davidson, N.; Ferdana, Z.A.; Finlayson, M.; Halpern, B.S.; Jorge, M.A.; Lombana, A.; Lourie, S.A.; Martin, K.D.; McManus, E.; Molnar, J.; Recchia, C.A.; Robertson, J. Marine ecoregions of the world: A bioregionalization of coast and shelf areas. BioScience 2007, 57, 573–583. [Google Scholar]
- Xie, G.D.; Lu, C.X.; Zeng, L.; Cao, S.Y.; Zhang, Y.S.; Leng, Y.F. Objective, progress and methodology of spatial function zoning. Geogr. Res 2009, 28, 561–570, in Chinese. [Google Scholar]
- Fu, B.J.; Liu, G.H.; Chen, L.D.; Ma, K.M.; Li, J.R. Scheme of ecological regionalization in China. Acta Ecol. Sinica 2001, 21, 1–6, in Chinese. [Google Scholar]
- Omernik, J.M. Ecoregions of the conterminous United States. Ann. Assoc. Am. Geogr. 1987, 77, 118–125. [Google Scholar]
- US Geological Survey (USGS). Aquatic Ecoregions in the United States GIS Data for the EPA Clean Air Mapping and Analysis Program (C-MAP); USGS: Golden, CO, USA, 1987. Available online: http://www.data.gov/geodata/g601930 accessed on 15 October 2011.
- Australian and New Zealand Environment and Conservation Councill (ANZECC). An Examination of an Aquatic Ecoregion Protocol for Australia; ANZECC: Canberra, Australia, 1997. Available online: http://www.environment.gov.au/water/publications/environmental/rivers/nrhp/pubs/ecoregion1.pdf accessed on 15 October 2011.
- Whittier, T.R.; Larsen, D.P.; Hughes, R.M.; Rohm, C.M.; Gallant, A.L. Ohio Stream Regionalization Project: A Compendium of Results; United States Environmental Protection Agency: Washington, DC, USA, 1988. Available online: http://nepis.epa.gov/Exe/ZyPURL.cgi?Dockey=2000HUQW.txt accessed on 15 October 2011.
- Cohen, P.; Andriamahefa, H.; Wasson, J.G. Towards a regionalization of aquatic habitat: Distribution of mesohabitats at the scale of a large basin. Regul. Rivers: Res. Manag. 1998, 14, 391–404. [Google Scholar]
- Bryce, S.A.; CLARKE, S.E. Landscape-level ecological regions: Linking state-level ecoregion frameworks with stream habitat classifications. Environ. Manag. 1996, 20, 297–311. [Google Scholar]
- Meng, W.; Zhang, Y.; Zheng, B.H. Aquatic ecological region approach and its application in China. Adv. Water Sci. 2007, 18, 293–300, in Chinese. [Google Scholar]
- Meng, W.; Zhang, Y.; Zheng, B.H. Study of aquatic ecoregion in Liao River Basin. Acta Sci. Circum. 2007, 27, 911–918, in Chinese. [Google Scholar]
- Gao, Y.N.; Gao, J.F. Delineation of aquatic ecoregions in Taihu Lake Basin. Geogr. Res 2010, 29, 111–117, in Chinese. [Google Scholar]
- Zhou, B.; Zheng, B. Research on aquatic ecoregions for lakes and reservoirs in China. Environ. Monit. Assess 2008, 147, 339–350. [Google Scholar]
- Huang, Y.; Cai, J.L.; Zheng, W.S.; Zhou, F.; Guo, H.C. Research progress in aquatic ecological function regionalization and its approach at watershed scale. Chin. J. Ecol. 2009, 28, 542–548, in Chinese. [Google Scholar]
- Li, Y.M.; Zeng, W.L.; Zhou, Q.X. Research progress in water eco-functional regionalization. Chin. J. Appl. Ecol 2009, 20, 3101–3108, in Chinese. [Google Scholar]
- Gao, Y.N.; Gao, J.F. Delineation of aquatic ecoregions in Taihu Lake Basin. Geogr. Res 2010, 29, 111–117, in Chinese. [Google Scholar]
- Han, J.; Kamber, M. Data Mining: Concepts and Techniques, 1st ed; Morgan Kaufmann: San Francisco, CA, USA, 2000; Volume 1, p. 550. [Google Scholar]
- Jin, X.C.; Tu, Q.Y. The Criterion of Eutrophic Lake Survey; China Environmental Science Press: Beijing, China, 1990; pp. 1–302, in Chinese. [Google Scholar]
- European Committee for Standardization. Water Quality—Sampling of Fish with Multi-Mesh Gillnets; European Committee for Standardization: Brussels, Belgium, 2005. Available online: http://surveygillnet.com/Z-files/multimesh%20standard.pdf accessed on 15 July 2011.
- Gao, F.; Deng, J.C.; Xu, Z.B.; Ning, Y.; Yin, H.B.; Gao, J.F. Ecological characteristics of macrobenthic communities in the Chaohu Lake Basin and their relationship with environmental factors. J. Anim. Vet. Adv. 2011, 10, 627–634. [Google Scholar]
- Huang, S.G.; Lin, S.Z.; Cao, G.Q.; Wu, S.F.; Chen, J.Y. Primary study on the characteristics of species diversity in Chinese-fir and broad-leaved mixed forests in different habitats. Biodivers. Sci. 2001, 9, 162–167, in Chinese. [Google Scholar]
Regionalizing levels | Regionalizing indicators | The role of indicators |
---|---|---|
Level I | Drainage density | To reflect the potential effects of spatial heterogeneous distribution of surface water resources on aquatic ecosystems. |
Elevation | To reflect the potential effects of regional terrain, which decides the spatial distribution of various factors such as precipitation, temperature, surface runoff, and other macro-scale factors in the Taihu Lake watershed, on the spatial variations of watershed aquatic ecosystems. | |
Level II | Percent of construction land area | To reflect the potential effects of point source and life diffused pollution load intensity on aquatic ecosystems. |
Percent of cultivated land area | To reflect the potential effects of agricultural non-point source pollution load intensity on aquatic ecosystems. | |
Soil type | To reflect the potential effects of spatial distribution of soil types heterogeneity on aquatic ecosystems. | |
Slope | To reflect the potential scour strength effects of undulating terrain, influencing the transport of nutrients and pollutants caused by land use activities and soil conditions, on aquatic ecosystems. |
Item | Aquatic ecoregion I1 | Aquatic ecoregion I2 |
---|---|---|
Total area | The total area is 11.0 thousand km2, accounting for 29.81% the total area of the Taihu Lake watershed | The total area is 25.9 thousand km2, accounting for 70.19% the total area of the Taihu Lake watershed |
Involved administrative regions | Zhenjiang, Changzhou, Nanjing, Wuxi, Xuancheng, Huzhou, Hangzhou | Zhenjiang, Changzhou, Wuxi, Suzhou, Shanghai, Jiaxing, Huzhou, Hangzhou |
Topography | Undulating surface, elevation rangs from 0 to 1,575 m with an average of 103.22 m, mean slope is 6.8° with a range difference of 66.89°. The area that elevation equal to and height above 10 m accounts for 73.97% the total area of aquatic ecoregion I1. | Flat surface with scattered hills, elevation rangs from 0 to 351 m with an average of 5.56 m, mean slope is 0.85° with a range difference of 39.79°. The area that elevation is less than 10 m accounts for 95.15% the total area of aquatic ecoregion I2 |
socioeconomic condition | The relatively smaller population density, slower industrial and agricultural development, and lower economic level, the main land use types are forest land and cultivated land | The bigger population density, the faster industrial and agricultural development and the higher economic level, the main land use types are construction land and irrigated cultivated land |
Water quality | The surface water quality is relatively better, mainly belonging to Grade II and Grade III according to the Chinese environmental quality standards for surface water (GB 3838-2002) | The surface water was seriously polluted with a poor water quality, mainly belonging to Grade IV, Grade V and Worse than Grade V according to the Chinese environmental quality standards for surface water (GB 3838-2002). |
Aquatic life | The best living conditions for aquatic life, with the nutritional status of mild eutrophication in some large lakes and reservoirs in the region, dominated by Carassius auratus and Hemiculter leucisculus fish species, tubificida and chironomidae on benthos. There are Limnodrilus hoffmeisteri, Glyptotendipes lobiferus, Bellamya aeruginosa, Branchiura sowerbyi, Limnodrilus grandisetosus, etc. | The worst living conditions for aquatic life, with the nutritional status of moderate eutrophication in some large lakes and reservoirs in the region. There are more than 60 fish species and 60 benthic macroinvertebrates species, dominated by Corbicula fluminea, Limnodrilus hoffmeisteri and Stenothyra glabra, etc. |
Item | Aquatic subecoregion II11 | Aquatic subecoregion II12 | Aquatic subecoregion II21 | Aquatic subecoregion II22 | Aquatic subecoregion II23 |
---|---|---|---|---|---|
Total area | The total area is 4.3 thousand km2, accounting for 11.65% the total area of the Taihu Lake watershed | The total area is 6.9 thousand km2, accounting for 18.70% the total area of the Taihu Lake watershed | The total area is 8.7 thousand km2, accounting for 23.58% the total area of the Taihu Lake watershed | The total area is 2.4 thousand km2, accounting for 6.50% the total area of the Taihu Lake watershed | The total area is 14.6 thousand km2, accounting for 39.57% the total area of the Taihu Lake watershed |
Involved administrative region | Zhenjiang, Changzhou, Nanjing, Wuxi, Xuancheng, | Huzhou, Hangzhou, Xuancheng | Zhenjiang, Changzhou, Wuxi, Suzhou | Changzhou, Wuxi, Suzhou | Wuxi, Suzhou, Shanghai, Jiaxing, Huzhou, Hangzhou |
Topography | Undulating surface, elevation ranges from 0 to 586 m with an average of 37.62 m, mean slope is 3.04° with a range difference of 41.30° | Undulating surface, elevation ranges from 0 to 1,575 m with an average of 140.92 m, mean slope is 8.55° with a range difference of 66.89° | Flat surface with scattered hills, elevation ranges from 0 to 353 m with an average of 5.96 m, mean slope is 1.00° with a range difference of 35.27° | Water surface with a small hill, elevation ranges from 0 to 324 m with an average of 1.57 m, mean slope is 0.33° with a range difference of 31.91° | Relative flat surface with scattered hills, elevation ranges from 0 to 331 m with an average of 5.55 m, mean slope is 0.79° with a range difference of 39.79° |
socioeconomic condition | The fourth biggest population density (more than 500 persons per sq km), the slowest industrial and agricultural development and the lowest economic level, the main land use types are cultivated land an forest land | The third biggest population density (more than 600 persons per sq km), the third fastest industrial and agricultural development and the third highest economic level, the main land use types are forest land and cultivated land | The second biggest population density (more than 1,000 persons per sq km), the faster industrial and agricultural development and the higher economic level, the main land use type is construction land and irrigated cultivated land | the main land use type is water body | The biggest population density (more than 2,600 persons per sq km), the fastest industrial development and the highest economic level, the main land use ypes are construction land and irrigated cultivated land |
Water quality | Relatively good water quality, most belonging to Grade III according to the Chinese environmental quality standards for surface water (GB 3838-2002) with a higher proportion meeting water quality standards. Grade III water environment functional area accounts for 82%, and Grade IV 18% | Relatively good water quality, the Tiaoxi river system and the Sianxi river system have the best surface water quality, and they belong to Grade III except the southern part of the region The best surface water quality, Water quality compliance rate of Dongtiaoxi is 60, and Xitiaoxi 100%, the monitoring sections into Taihu Lake have good water quality, most belonging to Grade III | The surface water was seriously polluted with a poor water quality, mainly belonging to Worse than Grade V and Grade V according to the Chinese environmental quality standards for surface water (GB 3838-2002). The major pollution indicators were NH3-N and TN. | The surface water most belongs to the Worse than Grade V with a nutritional status of moderate eutrophication, and high frequency of algal bloom, seriously affecting the safety of drinking water | The surface water was seriously polluted with a poor water quality, mainly belonging to Worse than Grade V and Grade V. Nitrogen and phosphorus pollution load is much higher than the standard of water environment capacity, and the living conditions of aquatic organisms is very poor |
Aquatic life | The relatively good living conditions for aquatic life. There are Limnodrilus hoffmeisteri, Glyptotendipes lobiferus, Bellamya aeruginosa, Branchiura sowerbyi, Limnodrilus grandisetosus, etc. | The best living conditions for aquatic life, dominated by Carassius auratus and Hemiculter leucisculus among fish species, tubificida and chironomidae in the benthos | The relatively poor living conditions for aquatic life, dominated by Limnodrilus hoffmeisteri and Bellamya aeruginosa in benthos | The relatively poor living conditions for aquatic life, with the nutritional status of moderate eutrophication in the local water body. There are more than 60 fish species and 59 benthic macroinvertebrates species in Taihu Lake, dominated by Corbicula fluminea, Limnodrilus hoffmeisteri and Stenothyra glabra, etc. | The relatively poor living conditions for aquatic life, dominated by tubificida and chironomidae and accompanied by planorbidae; trumpet snail, gammaridae and erpobdellidae, etc. |
Season | Ecoregion | Species1 | Species2 | Species3 | Species4 |
---|---|---|---|---|---|
Normal | I1 | C. auratus | H. leucisculus | ||
season | I2 | C. auratus | H. leucisculus | Pseudorasbora | C. carpio |
Wet | I1 | C. auratus | H. leucisculus | C. carpio | |
season | I2 | C. auratus | H. leucisculus | S. macrops | C. carpio |
Ecoregions and subecoregions | Benthic density (ind/m2) | Benthic biomass (g/m2) | Dominant species |
---|---|---|---|
I1 | 5,544 | 204.82 | B. aeruginosa, C. plumosus |
I2 | 3,777 | 117.61 | L. hoffmeisteri, B. aeruginosa |
II11 | 12,435 | 475.33 | B. eruginosa, L. hoffmeisteri |
II12 | 376 | 1.94 | L. hoffmeisteri, C. plumosus |
II21 | 8,339 | 56.19 | L. hoffmeisteri, B. aeruginosa |
II22 | 45 | 0.13 | L. hoffmeisteri, Nais sp. |
II23 | 2,117 | 167.90 | L. hoffmeisteri, B. aeruginosa |
© 2011 by the authors; licensee MDPI, Basel, Switzerland This article is an open-access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/3.0/).
Share and Cite
Gao, Y.; Gao, J.; Chen, J.; Xu, Y.; Zhao, J. Regionalizing Aquatic Ecosystems Based on the River Subbasin Taxonomy Concept and Spatial Clustering Techniques. Int. J. Environ. Res. Public Health 2011, 8, 4367-4385. https://doi.org/10.3390/ijerph8114367
Gao Y, Gao J, Chen J, Xu Y, Zhao J. Regionalizing Aquatic Ecosystems Based on the River Subbasin Taxonomy Concept and Spatial Clustering Techniques. International Journal of Environmental Research and Public Health. 2011; 8(11):4367-4385. https://doi.org/10.3390/ijerph8114367
Chicago/Turabian StyleGao, Yongnian, Junfeng Gao, Jiongfeng Chen, Yan Xu, and Jiahu Zhao. 2011. "Regionalizing Aquatic Ecosystems Based on the River Subbasin Taxonomy Concept and Spatial Clustering Techniques" International Journal of Environmental Research and Public Health 8, no. 11: 4367-4385. https://doi.org/10.3390/ijerph8114367
APA StyleGao, Y., Gao, J., Chen, J., Xu, Y., & Zhao, J. (2011). Regionalizing Aquatic Ecosystems Based on the River Subbasin Taxonomy Concept and Spatial Clustering Techniques. International Journal of Environmental Research and Public Health, 8(11), 4367-4385. https://doi.org/10.3390/ijerph8114367