A Model for Accurate Determination of Environmental Parameters in Indoor Zoological and Botanical Gardens Supporting Efﬁcient Species Management

: The detailed evaluation of environmental parameters can be a great tool for the optimal selection and location of vegetable species, not only in vegetable production facilities and greenhouses but also in zoological and botanical gardens, which frequently maintain delicate and exotic plant species with strict environmental requirements in immersive exhibits where conditions can vary remarkably. This study, developed at an indoor zoological garden (Biodomo—Parque de las Ciencias de Granada, Spain), evaluates a sampling protocol for the determination of seven environmental parameters: daily light integral (DLI) was determined at nine different locations of the facility using a portable Light Quantum SQ-500 sensor; air temperature, atmospheric pressure, and air relative humidity were measured using a ﬁxed ATMOS14 sensor; and soil temperature, soil water content, and soil conductivity were determined using a ﬁxed TEROS12 sensor. Values recorded for DLI showed statistically signiﬁcant variations across the nine different sampling locations, as well as between the different months in all sampling spots. Signiﬁcant variations were also detected across the 12 months of study for the rest of environmental parameters evaluated, and correlations were found between the studied parameters, with the correlation between soil and air temperature the strongest (r s = 0.758) and soil temperature signiﬁcantly superior to air temperature. The methodology described in this study can be easily reproduced in similar indoor zoological and botanical facilities, increasing the knowledge of the environmental conditions, and allowing corrections that could improve species selection, location, and management.


Introduction
The curiosity and desire to build artificial spaces that recreate nature within a safe environment for its contemplation and enjoyment has existed in our civilization for thousands of years. There are records of zoological and botanical gardens as far back as ancient Egypt [1]. Since the onset of long-distance travel and trade, humans have been fascinated by plant and animal species from other parts of the globe, so different from those present in their traditional environment. This has been accompanied by the desire to bring those foreign species back home to enjoy their natural properties, colors, and exotic shapes [2]. Transporting species from tropical areas to regions with colder climates indirectly demands ornamental greenhouses to allow their survival and development [3].
When focusing on vegetable species, the maintenance in indoor facilities such as botanical gardens has evolved enormously, from the first greenhouses that consisted of a simple protective cover from inclement weather to the sophisticated climate control systems that currently exist, where practically all parameters can be controlled including lighting supplementation [4,5]. Plants within the zoological facilities play a fundamental role, since they act as vital support for the animals, enriching the space where they live in various ways and making it easier for these animals to enjoy a life that is as natural as possible across the 24 h period. Furthermore, they create an environment with natural vegetation specific to a certain region, and therefore help to "transport" the visitors to another habitat and provide a feeling of being immersed in a tropical destination [6].
Studies evaluating the effect of light on plant growth and development focus on photosynthetically active radiation (PAR), which corresponds to the spectral range of radiation comprised between 400 and 700 nanometers wavelength, and is the radiation used by photosynthetic organisms during photosynthesis [7]. The number of photosynthetically active particles of light (photons) of PAR, delivered to a specific area over a 24-h period, is studied as daily light integral (DLI). The variable DLI is frequently used to describe plant light requirements, and there are multiple studies determining DLI in different types of horticultural and ornamental crops [8][9][10][11]. However, there are currently no published studies determining DLI values in interior zoological facilities. Furthermore, the development of DLI determination studies has also proven of great use, not only in terrestrial species but also in aquatic animal and plant species [12]. The design of a habitat where plant and animal species coexist entails that a balance is reached. These habitats, in which animals successfully develop natural behaviors while plants thrive and resist animal erosion, require detailed planning and consideration of environmental parameters [13]. In interior facilities, an environmental parameters design is essential for adequate species choice and distribution; if optimal lighting, temperature, humidity, etc., are not provided, plants will not develop and grow successfully, and the desired aesthetic potential will not be achieved [14]. Considering this, the quantification of light input and environmental parameters such as relative humidity, air, and soil temperature determination acquire a great importance, not only in vegetable production units but also in zoological facilities, in order to evaluate which lighting and environmental installation will be necessary, as well as the effective selection of species based on their environmental requirements [12,14].
The use of environmental sensors for the measurement of environmental parameters is already carried out in indoor horticultural and livestock farms; this has demonstrated the importance of environmental parameter (PAR light, temperature, air relative humidity, etc.) determination for the accurate management of these facilities [15,16]. On the other hand, studies determining environmental parameters in zoological and botanical facilities are still very limited, with little attention paid to the selection, use and management of the vegetable species, and with few studies reporting the importance of certain environmental parameters for the optimal health and development of animals [6,[17][18][19]. In livestock and horticultural production facilities, environmental parameter control studies are mainly focused on increasing production and nutritional quality, while in zoological and botanical gardens the main objective of environmental studies should be to improve animal and plant welfare, as well as to generate better criteria for species selection based on these parameters [15,16,20].
Because of this, the objective of our study was to develop and describe a method for environmental parameters measurement for indoor zoological gardens, combining different sensors of scientific quality used in indoor horticultural and livestock farms as well as research centers, which could be used easily and safely in mixed indoor installations in which plants and animals coexist. Our study aims to evaluate the potential use of these sensors and the sampling model developed, to obtain baseline data on the main environmental parameters that could provide information to improve and help adjust the climate and lighting control systems of this type of facility, as well as the selection, location and management of both plant and animal species within it.

Environmental Conditions
This study was developed at Biodomo-Parque de las Ciencias de Granada (Av. de las Ciencias. s/n. Granada, Spain; https://www.parqueciencias.com/biodomo (accessed on 10 September 2022)). Biodomo is a zoological and botanical facility inaugurated in 2016, which recreates terrestrial and aquatic habitats of the equatorial and tropical regions of the globe (Figure 1). Terrestrial ecosystems correspond to the Amazon rainforest, Southeast Asia and Madagascar; aquatic ecosystems represent fluvial, mangrove, coastal, and marine ecosystems. The present study focuses on the terrestrial multispecies enclosures of Biodomo, which house different animal and plant species living together, creating an ecosystem where they coexist in balance. With a total volume of 18,000 m 3 , a total surface of 4252 m 3 and total exhibit surface of 2700 m 2 , Biodomo is an interior zoological and botanical facility where the contribution of natural light through the roof is through polycarbonates, which constitute 30% of the roof surface. The rest of the roof is constituted by metal and thermal isolation where light does not penetrate. J. Zool. Bot. Gard. 2022, 3, FOR PEER REVIEW 3 and lighting control systems of this type of facility, as well as the selection, location and management of both plant and animal species within it.

Environmental Conditions
This study was developed at Biodomo-Parque de las Ciencias de Granada (Av. de las Ciencias. s/n. Granada, Spain; https://www.parqueciencias.com/biodomo (accessed on 10 th September 2022)). Biodomo is a zoological and botanical facility inaugurated in 2016, which recreates terrestrial and aquatic habitats of the equatorial and tropical regions of the globe (Figure 1). Terrestrial ecosystems correspond to the Amazon rainforest, Southeast Asia and Madagascar; aquatic ecosystems represent fluvial, mangrove, coastal, and marine ecosystems. The present study focuses on the terrestrial multispecies enclosures of Biodomo, which house different animal and plant species living together, creating an ecosystem where they coexist in balance. With a total volume of 18,000 m 3 , a total surface of 4252 m 3 and total exhibit surface of 2700 m 2 , Biodomo is an interior zoological and botanical facility where the contribution of natural light through the roof is through polycarbonates, which constitute 30% of the roof surface. The rest of the roof is constituted by metal and thermal isolation where light does not penetrate. To enrich the natural light input, lighting is supplemented by light-emitting diode (LED) screens (Sequoia Cultiva Wall 50-60 Hz 224 W IGNIA GREEN led screens), which emit photosynthetic active radiation (PAR) light spectrum as well as white light LED screens (Konak 200 W 5000 K 150° IP65 SECOM) with 12 h light: 12 darkness programmed periods The location of the PAR and white light LED screens throughout Biodomo has been optimized during the last four years and is represented in Figure 2. Environmental parameters at Biodomo are regulated by an air conditioning system made up of two water-water chillers and two air treatment units. Air treatment units purify the air using fans and filters (consisting of G4 Class Air 10 micron filter, carbon filters, and F7 EU7 Class Air filter which are renewed periodically), as well as climatize by heating, or refrigerating the air to achieve the desired temperature depending on the season. Air humidity is To enrich the natural light input, lighting is supplemented by light-emitting diode (LED) screens (Sequoia Cultiva Wall 50-60 Hz 224 W IGNIA GREEN led screens), which emit photosynthetic active radiation (PAR) light spectrum as well as white light LED screens (Konak 200 W 5000 K 150 • IP65 SECOM) with 12 h light: 12 darkness programmed periods The location of the PAR and white light LED screens throughout Biodomo has been optimized during the last four years and is represented in Figure 2. Environmental parameters at Biodomo are regulated by an air conditioning system made up of two waterwater chillers and two air treatment units. Air treatment units purify the air using fans and filters (consisting of G4 Class Air 10 micron filter, carbon filters, and F7 EU7 Class Air filter which are renewed periodically), as well as climatize by heating, or refrigerating the air to achieve the desired temperature depending on the season. Air humidity is maintained using a fogging system controlled by a humidistat (SMD4500 controller; ELLIWELL Ibérica; 46980 Valencia, Spain) Furthermore, the facility also has several additional auxiliary systems for climate control. For instance, some areas of Biodomo are equipped with fabric awnings on the ceiling below the roof that can be deployed to create shadows in the summer months when solar radiation is very strong. The roof of the Biodomo also has two movable skylights, which can be opened to increase sunlight input and air exchange when exterior conditions are optimal for animals and plants.
J. Zool. Bot. Gard. 2022, 3, FOR PEER REVIEW 4 maintained using a fogging system controlled by a humidistat (SMD4500 controller; EL-LIWELL Ibérica; 46980 Valencia, Spain) Furthermore, the facility also has several additional auxiliary systems for climate control. For instance, some areas of Biodomo are equipped with fabric awnings on the ceiling below the roof that can be deployed to create shadows in the summer months when solar radiation is very strong. The roof of the Biodomo also has two movable skylights, which can be opened to increase sunlight input and air exchange when exterior conditions are optimal for animals and plants.

Light Sensor
For this study, one Light Quantum Sensor Model SQ-500 (Apogee Instruments. Inc., Logan, UT 84321, USA) was used. Sensor sensitivity was 0.01 mV/µmol·m −2 s −1 , measurement range was 0 to 4000 µmol m −2 s −1 , long-term drift (non-stability) was under 2% per year, non-linearity less than 1%, and response time less than 1 ms, a field of view 180 • , spectral range 389 to 692 nm ± 5 nm, spectral selectivity less than 10% from 412 to 682 ± 5 nm, azimuth error less than 0.5%, tilt error less than 0.5%, and temperature response −0.11 ± 0.04% per • C. This sensor registered one DLI reading every 3 min, and the daily values were computed as the total DLI per squared meter per day.
To increase the reliability of the data obtained and determine the DLI that the plants within the facility were being exposed to, the light sensor was located in the enclosure at the vegetation level. As the different animal species maintained at the enclosure could interact with the sensor throughout the day, and in order to protect it, a protective cage was designed using wood on the sides and base, and a 150 mm methacrylate lid was fixed to the cage using a nylon flange ( Figure 3). While most avian species seemed not affected by the introduction of the sensor in their enclosures, some mammal species showed an initial interest in the sensor protective cage, such as ring-tailed lemurs (Lemur catta), the white-faced saki (Pithecia pithecia), and the white-lipped tamarin (Saguinus labiatus). These animals are provided weekly with different enrichments and initially inspected the cage in detail for food, though animals did not show further interest after the first week and during further data collections throughout the 12-month study.
within the facility were being exposed to, the light sensor was located in the enclosure at the vegetation level. As the different animal species maintained at the enclosure could interact with the sensor throughout the day, and in order to protect it, a protective cage was designed using wood on the sides and base, and a 150 mm methacrylate lid was fixed to the cage using a nylon flange ( Figure 3). While most avian species seemed not affected by the introduction of the sensor in their enclosures, some mammal species showed an initial interest in the sensor protective cage, such as ring-tailed lemurs (Lemur catta), the white-faced saki (Pithecia pithecia), and the white-lipped tamarin (Saguinus labiatus). These animals are provided weekly with different enrichments and initially inspected the cage in detail for food, though animals did not show further interest after the first week and during further data collections throughout the 12-month study.
The photosynthetic active radiation (PAR) light (measured in photosynthetic photon flux density (PPFD)) loss percentage associated with the protective cage was calculated by taking ten repetitive readings, with and without the cage, to determine the light loss percentage caused by the protective cage, which resulted in a mean loss of 8.18% ± 0.53 SD (min 7.42%, max 9.28%). The PPFD results provided in this study have been expressed as corrected data.

Air and Soil Sensors
For temperature, ambient humidity, and atmospheric pressure measurement, this study used one ATMOS14 sensor (Meteor group, Inc., USA. 2365 NE Hopkins Court, Pullman, WA 99163). Temperature measurement range was −40 to 80 °C, resolution of ± 0.1 °C, accuracy ± 0.5 °C, equilibration time < 400 s, long-term drift < 0.004 °C/year. Relative Humidity (RH) measurement range was 0 to 100% RH (0.00-1.00), resolution was 0.1% The photosynthetic active radiation (PAR) light (measured in photosynthetic photon flux density (PPFD)) loss percentage associated with the protective cage was calculated by taking ten repetitive readings, with and without the cage, to determine the light loss percentage caused by the protective cage, which resulted in a mean loss of 8.18% ± 0.53 SD (min 7.42%, max 9.28%). The PPFD results provided in this study have been expressed as corrected data.

Light (DLI) Sampling Protocol
To carry out this study, nine sampling points were selected for DLI determination. These sampling points are represented in Figure 2 and were designated because they are representative of the different enclosures and vegetation environments of Biodomo. The sensor, which registered one reading every 3 min, determined a total DLI value per square meter per day; it was moved daily across the different sampling locations from point number 1 to 9, thus rotating through all the points from 1 January 2021 to 31 December 2021 with a 9-day cycle and for 365 days. This portable data logger and the designed rotating system allowed us to obtain three to four monthly DLI measurements at each sampling point and then averaged these measurements to obtain the monthly DLI.

Sampling Protocol for the Rest of the Environmental Parameters
As represented in Figure 2, ATMOS14 and TEROS12 probes were located in a fixed position at the Asian installation at Biodomo and were not rotated throughout the year like the light sensor. The ATMOS14 sensor was located hidden between the foliage of a Port Jackson fig (Ficus rubiginosa) at 3 m above the ground, and the TEROS12 sensor was buried at a depth of 50 cm, under the same tree and between the roots. While the light sensor was portable, the data logger shared between ATMOS14 and TEROS12 was non-

Light (DLI) Sampling Protocol
To carry out this study, nine sampling points were selected for DLI determination. These sampling points are represented in Figure 2 and were designated because they are representative of the different enclosures and vegetation environments of Biodomo. The sensor, which registered one reading every 3 min, determined a total DLI value per square meter per day; it was moved daily across the different sampling locations from point number 1 to 9, thus rotating through all the points from 1 January 2021 to 31 December 2021 with a 9-day cycle and for 365 days. This portable data logger and the designed rotating system allowed us to obtain three to four monthly DLI measurements at each sampling point and then averaged these measurements to obtain the monthly DLI.

Sampling Protocol for the Rest of the Environmental Parameters
As represented in Figure 2, ATMOS14 and TEROS12 probes were located in a fixed position at the Asian installation at Biodomo and were not rotated throughout the year like the light sensor. The ATMOS14 sensor was located hidden between the foliage of a Port Jackson fig (Ficus rubiginosa) at 3 m above the ground, and the TEROS12 sensor was buried at a depth of 50 cm, under the same tree and between the roots. While the light sensor was portable, the data logger shared between ATMOS14 and TEROS12 was non-portable, and the strategic point was determined for the placement of both sensors. While light depends on the structure of the building and varies across the different areas and enclosures, the air environmental parameters collected by the ATMOS14 sensor are representative of the entire indoor facility. The soil probe TEROS12 was not translocated either, and its fixed location was maintained throughout the study as it was buried in the ground and lacked a portable datalogger. Both sensors were configured to register measures every 10 min. Data recording started on 27 May 2021 and ended on 25 May 2022; these dates correspond to the availability of the sensor for sampling at Biodomo.

Data and Statistical Analysis
The distribution of the DLI data was evaluated using a Kolmogorov-Smirnov normality test with Lilliefors correction, which showed that the data followed a non-Gaussian distribution [21]. To evaluate the existence of significant differences between the different sampling locations and between the different months at each sampling point, a Kruskal-Wallis test and Dunn's post hoc tests with Bonnferroni corrections were used, and significance was set at a p-value of 0.05. Statistical comparison between air and soil temperature was performed using a Mann-Whitney U test [22]. Correlation between the different environmental parameters was evaluated by calculating the Spearman's rank correlation coefficient (r s ) with a confidence interval of 0.95. A principal component analysis (PCA) was performed for the six variables measured using the ATMOS14 and TEROS12 sensors; daily light integral analysis has not been included in the PCA model due to a different sampling methodology. All statistical analyses were performed using the statistical software package RStudio R (Version 1.2.504; RStudio Team. 2020 Boston, MA 02210, USA; www.rstudio.com (accessed on 20 September 2022)).

Variations in DLI
Monthly DLI values and descriptive statistics for the different sampling locations are reported in Table 1. Values recorded for DLI showed statistically significant variations across the different sampling locations, as well as between the different months. Statistically significant differences (p < 0.05; Dunn's post hoc with Bonnferroni corrections) were detected between the DLI values recorded in the months of June vs. December, June vs. January, June vs. November, July vs. November, May vs. November, May vs. December, and May vs. January. Furthermore, as it can be appreciated in Figure 5, important differences were detected for DLI annual measurements between the locations ( Figure 5). These differences were statistically significant (Dunn's post hoc with Bonnferroni corrections) between: Amazon rainforest spot 1 vs. Amazon rainforest spot 2 (p < 0.05), Amazon rainforest spot 1 vs. Amazon rainforest spot 2 (p < 0.05), Amazon rainforest spot 3 vs. Asiatic island (p < 0.01), Amazon rainforest spot 1 vs. Asia spot 1 (p < 0.05), Amazon rainforest spot 2 vs. Asiatic Island (p < 0.01), Asia spot 1 vs. Asiatic island (p < 0.01), Asia enclosure spot 2 vs. Asiatic island (p < 0.05), Madagascar enclosure spot 1 vs. Asiatic island (p < 0.01), Burmese ruins vs. Asiatic island (p < 0.01), Amazon rainforest spot 3 vs. Madagascar enclosure spot 2 (p < 0.01), Amazon rainforest spot 2 vs. Madagascar enclosure spot 2 (p < 0.01), Asia spot 1 vs. Madagascar enclosure spot 2 (p < 0.01), Amazon rainforest spot 1 vs. Burmeese ruins (p < 0.05), Madagascar enclosure spot 2 vs. Burmeese ruins (p < 0.05), and Madagascar enclosure spot 1 vs. Madagascar enclosure spot 2 (p < 0.05). Statistical comparison between the different months across the sampling spots was not possible due to a small number of samples per month and spot (3-4 samples) associated with the rotating system and limited by the availability of only one DLI sensor.

Variations in Air Temperature, Air Relative Humidity, Ambient Pressure, Soil Water Content, Soil Temperature and Conductivity
Descriptive statistics for air temperature (°C), relative humidity (RH), ambient pressure (kPa), soil water content (m 3 /m 3 ), soil temperature (°C), and conductivity (mS/cm) at the fixed sampling spot for the different months are reported in Table 2. When focusing on air temperature, statistically significant differences (p < 0.01; Dunn's post hoc with Bonnferroni corrections) were detected between all months with the exception of December vs. November (p = 0.095) and April vs. February (p = 0.205); no differences were detected between August and July (p = 1.000). For relative humidity, statistically significant differences (p < 0.01; Dunn's post hoc with Bonnferroni corrections) were detected between all sampled months except for: November vs. September (p = 0.051), January vs. September (p = 0.016), November vs. October (0.019), May vs. October (p = 0.515), and April vs. October (p = 0.166). No differences were detected between May and November,

Variations in Air Temperature, Air Relative Humidity, Ambient Pressure, Soil Water Content, Soil Temperature and Conductivity
Descriptive statistics for air temperature ( • C), relative humidity (RH), ambient pressure (kPa), soil water content (m 3 /m 3 ), soil temperature ( • C), and conductivity (mS/cm) at the fixed sampling spot for the different months are reported in Table 2. When focusing on air temperature, statistically significant differences (p < 0.01; Dunn's post hoc with Bonnferroni corrections) were detected between all months with the exception of December vs. November (p = 0.095) and April vs. February (p = 0.205); no differences were detected between August and July (p = 1.000). For relative humidity, statistically significant differences (p < 0.01; Dunn's post hoc with Bonnferroni corrections) were detected between all sampled months except for: November vs. September (p = 0.051), January vs. September (p = 0.016), November vs. October (0.019), May vs. October (p = 0.515), and April vs. October (p = 0.166). No differences were detected between May and November, February vs. November, February vs. March, August vs. March, and August vs. February (p = 1.000). Something similar happened with atmospheric pressure, and statistically significant differences (p < 0.01; Dunn's post hoc with Bonnferroni corrections) were detected between all months except for July vs. November (p = 0.193). Soil water content also showed statistically significant differences (p < 0.01; Dunn's post hoc with Bonnferroni corrections) between all studied months except for August vs. September (p = 0.382). The same happened with soil temperature, and all months showed statistically significant differences (p < 0.01; Dunn's post hoc with Bonnferroni corrections) except for August vs. July (p = 0.055) and December vs. February (p = 0.087). Conductivity also varied significantly (p < 0.01; Dunn's post hoc with Bonnferroni corrections) across the different months, except for July vs. November (p = 0.158) and June vs. November (p = 0.035). Values for the main environmental parameters recorded every 10 min by ATMOS14 and TEROS12 probes are represented in Figure 6. Spearman's rank correlation coefficient (r s ) varied greatly when comparing the different environmental parameters measured by the ATMOS14 and TEROS12 probes. The correlation between the different environmental parameters measured by ATMOS14 and TEROS12 probes is presented in Table 3. Principal component analysis (PCA) is represented in Figure 7 and shows the correlation between air temperature and soil temperature, as well as between soil conductivity and soil water content. Principal component analysis also showed that the variables air temperature, soil temperature, soil conductivity, and soil water content were responsible for most of the variance in the data, while ambient pressure and relative humidity had a smaller contribution to this variance. As it can be appreciated in Table 4, principal component 1 only represented 38.82% of the total variance, principal component 2 only represented 24.26% of the variance, and principal component 3 only represented 17.325% of the variance. If we use the correlations between the principal components and the original variables to interpret these principal components (Table 5), we found that principal component 1 was strongly correlated with the variables soil temperature, air temperature, and soil conductivity; principal component 2 was strongly correlated with atmospheric pressure and soil water content; and principal component 3 was strongly correlated with relative humidity.    As it can be appreciated in Figure 8, soil temperature was always superior to air temperature. Statistically significant differences (p < 0.01; Mann-Whitney U test) were detected between soil and air temperature for all studied months, except for May (p = 0.027).  As it can be appreciated in Figure 8, soil temperature was always superior to air temperature. Statistically significant differences (p < 0.01; Mann-Whitney U test) were detected between soil and air temperature for all studied months, except for May (p = 0.027).  Figure 8. Mean ± SD air and soil temperatures at Biodomo-Parque de las Ciencias de Granada. Statistically significant differences (p < 0.01; Mann-Whitney U test) were detected between air and soil temperature in all months except May (p = 0.027).

Discussion
This study shows the potential of a detailed and accurate monitoring of the environmental parameters such as light intensity, air temperature, relative humidity, atmospheric pressure, soil temperature, soil water content, and soil conductivity. Our manuscript presents an affordable methodology that can provide the institution with important data. It also shows that many of the environmental parameters are correlated, demonstrating the potential use of the obtained information for improving ambient conditions and the efficient selection of vegetable species to form part of a zoological and botanical garden.
In our study, statistically significant differences were detected for annual DLI measurements not only in different animal enclosures within the building, but also between sampling spots in the same enclosure (For instance, between the sampling spots "Madagascar enclosure spot 1" vs. "Madagascar enclosure spot 2" and "Amazon rainforest spot 1" vs. "Amazon rainforest spot 2", as well as "Amazon rainforest spot 1" vs. "Amazon rainforest spot 3" (p < 0.05). These significant variations reveal the importance of developing studies determining DLI in the main areas where vegetation is located at the facility, to allow the efficient selection and location of the vegetable species within the enclosure depending on their DLI requirements. The DLI quantification allows us to accurately determine the most suitable location within the enclosure for the selected species, facilitating its adaptation and development [9,12,14]. Furthermore, quantifying DLI allows the correction of light input if DLI values are not the ones desired, which can be increased by

Discussion
This study shows the potential of a detailed and accurate monitoring of the environmental parameters such as light intensity, air temperature, relative humidity, atmospheric pressure, soil temperature, soil water content, and soil conductivity. Our manuscript presents an affordable methodology that can provide the institution with important data. It also shows that many of the environmental parameters are correlated, demonstrating the potential use of the obtained information for improving ambient conditions and the efficient selection of vegetable species to form part of a zoological and botanical garden.
In our study, statistically significant differences were detected for annual DLI measurements not only in different animal enclosures within the building, but also between sampling spots in the same enclosure (For instance, between the sampling spots "Madagascar enclosure spot 1" vs. "Madagascar enclosure spot 2" and "Amazon rainforest spot 1" vs. "Amazon rainforest spot 2", as well as "Amazon rainforest spot 1" vs. "Amazon rainforest spot 3" (p < 0.05). These significant variations reveal the importance of developing studies determining DLI in the main areas where vegetation is located at the facility, to allow the efficient selection and location of the vegetable species within the enclosure depending on their DLI requirements. The DLI quantification allows us to accurately determine the most suitable location within the enclosure for the selected species, facilitating its adaptation and development [9,12,14]. Furthermore, quantifying DLI allows the correction of light input if DLI values are not the ones desired, which can be increased by adding supplemental light or reduced by providing shade. This is important, as a linear relationship has been found between plant growth and the cumulative DLI, being an accurate indicator for growth and development [23].
The differences between sampling spots and the possibility to determine DLI levels at the desired study locations within the enclosure facilities allows the detection of DLI variations which could be produced by the building and enclosure design, which may not be appreciable to the naked eye [8,24]. The use of the method described in this study can allow the detection of the anomalies caused by the aesthetics of the building to be detected and corrected, perhaps by placing an artificial light source for the months in which the design of the building and the incidence of the sun produce shadows. An example of the detection of these anomalies can be observed at sampling point 1 (Asiatic island): during the months of October, November, and December, the sun's incidence makes the sun rays directly hit the island through some windows, producing a significant increase in DLI measurements during those months, with higher DLI levels than those recorded during the summer months. This is because when the sun is higher, the building structure creates a shadow in this location, reducing DLI levels. This information is of great use, as it will allow the technical and gardening team at Biodomo to place a PAR LED panel that will supplement this location of the building during the months when the structure creates this shadow, to achieve the desired DLI.
The implementation of studies such as the one proposed by the authors will allow a more efficient use of artificial light, as artificial light panels will only be used in areas where DLI are under the desired levels, promoting energy saving. This allows the maximum use of sunlight, and the placement and use of artificial light panels only at the spot/season/hour required. The DLI data would allow a separation in the artificial light panels control and the design of a lighting system that would work "like clockwork", varying daily and between seasons, depending on the design of the building and the vegetable species' requirements. This would avoid having all the light panels turned on continuously, many of them probably even located in spots where DLI is already over the desired levels. This whole process would allow energy savings, as well as a more efficient, economic, and sustainable light management in zoological and botanical gardens.
The experimental design and the methodology proposed in this study, with the use of high precision sensors for measuring environmental variables, opens the possibility of future studies with animal species maintained in indoor facilities, which also have certain DLI and UV necessities [25,26]. The knowledge of environmental DLI at the indoor facility will allow the optimal selection of the most appropriate species for the environmental light levels, as well as determining if an increase or decrease in light intensity is necessary for the enclosure depending on the species housed. A detailed study of these parameters could help to make a better choice of animal species or a better adaptation of the facility for the existing animals, thus improving their health and well-being [25,27]. Since the DLI sensor is submersible, its usefulness for designing correct lighting in aquatic environments for corals and aquatic plants in aquariums is significant, based on the fact that each species has certain DLI needs [12,28].
An important limitation of the light evaluation in this study was that only one DLI sensor was available for sampling, and therefore a rotating system was developed to determine the DLI levels in the selected spots within Biodomo. Future studies could include more than one sensor, which could be placed at the pre-determined sampling locations and therefore provide continuous DLI data for every sampling spot without the need to rotate the sensor, reducing sampling time and avoiding errors produced by the daily cloud and DLI variations.
When focusing on relative humidity, despite statistically significant variations being detected for this parameter across the year in our study, no statistically significant differences were detected for air relative humidity between August (which is one of the warmest and driest months in Granada, Spain) and February (which is one of the coldest months and accumulates most of the precipitation in Granada) [29]. The absence of statistically significant differences between these months can be explained as the environmental humidity control system (described in the Methodology section and based on automated air conditioning and fog systems at Biodomo) maintaining a relatively consistent relative humidity at Biodomo. Previous studies revealed that the maintenance of constant and appropriate relative humidity levels is essential not only for optimal plant development (as it affects photosynthesis, leaf growth, and disease incidence) but also for guaranteeing animal health and welfare [30,31]. Because of this, the monitorization of air relative humidity across the year in indoor facilities maintaining animal and plant species is of great importance to provide the most appropriate environmental conditions and avoid the occurrence of diseases.
Despite the variations observed in the temperatures registered in our study across the year, the indoor air and soil temperature results, with a maximum of 29.10 and 26.20 • C, respectively, and a minimum of 17.15 and 22.20 • C, respectively, are relatively stable compared to the continental climate of Granada, which in 2021 had maximum temperatures of 45 • C during August and minimum temperatures of −5 • C during February, revealing the efficiency of the automated climate control system design at Biodomo [29]. The maintenance of a constant temperature with mild variations in this parameter is important, as zoological and botanical gardens frequently host tropical species highly sensitive to temperature variations. The presence of sensors continuously registering environmental temperature is important to detect possible variations in temperature, which may exceed the desired maximum and minimum limits, and together with variations in air humidity can lead to disease in the animal and plant collections [32,33]. One factor to consider when keeping plant and animal species at indoor facilities is the possible alteration of circadian rhythms, since we are creating an artificial environment where it would be difficult to replicate the cycles of light and temperature that occur in nature [34,35]. This problem is exacerbated in mixed facilities where we house species from such diverse locations under the same roof. This marks an important guideline in the species selection criteria, considering that alterations in these cycles can cause certain problems, not only in the animals but also in the plants, such as alterations in the immune system or changes in their morphology [36,37]. In the care and management of certain animal species where seasonal cycles are very important, being able to emulate these cycles would be essential for their well-being and behavioral development [38].
The maintenance of plants in interior facilities such as indoor zoological gardens usually has particularities, with the limitation of space for the development of plants and their roots as an important concern. In Biodomo, all the plants are cultivated in buried pots, to avoid uncontrolled root proliferation, to have more precise control on irrigation, and to facilitate plant translocation. This study shows that the soil sensor can provide relevant information for daily maintenance and decision making, such as when to apply irrigation correctly. This agrees with previous studies evaluating the potential use of soil sensors to define efficient irrigation schedules [39].
This study shows that soil temperature and air temperature, while being strongly correlated (r s = 0.758), differ significantly from each other throughout the year, being soil temperature significantly higher at the sampled spot when compared to air temperature in all months except May. Biodomo does not count with a soil heating and refrigeration system, and the soil, due to its solid condition, accumulates thermal energy [40]. As can be observed in Figure 8, during spring the air temperature increases rapidly, reaching soil temperature in May, and increasing both soil and air temperatures sharply during the summer months, being that soil temperature is always higher due to solar radiation and energy accumulation. At the end of August, both air and soil temperatures drop, maintaining soil always at a higher temperature during autumn and winter. The determination of soil temperature is of great interest, as previous studies already demonstrated that these parameters directly affect the growth of the root system, and therefore the optimal development of the plant [41].
When looking at conductivity, our study detected a marked relationship between conductivity and soil water content (with a correlation of r s = 0.572). Furthermore, Figure 6 allows the visual detection of different peaks in the graphs for both parameters, which is related to the irrigations carried out throughout the year, observing a higher irrigation frequency during the months of March to August and a lower frequency from August to February. The live determination of soil parameters such as temperature, water content, and conductivity will also allow the early diagnosis and prevention of fungal pathologies since there is a direct relationship between these parameters and the appearance of fungal diseases [42]. A limitation to this study is that only one TEROS12 soil sensor was available during sampling, and therefore we could only determine soil parameters for a fixed location at Biodomo. Future studies with soil sensors placed across different locations will allow a more precise evaluation of soil conditions and if plant requirements are met at the different enclosures. This will allow the design of a detailed and sustainable irrigation schedule. For instance, this efficient irrigation system based on independent soil sensors would provide water only in the locations with a water deficit, always guaranteeing the optimal soil conditions and promoting plant development while saving water.
By relating atmospheric pressure and relative humidity we can obtain another useful parameter for plant management: evapotranspiration. This parameter allows the calculation of the amount of water that is being lost both by the plants and the soil [43]. Evapotranspiration, together with the data obtained by a soil probe, allows us to accurately and reliably determine plant water needs and establish more efficient irrigation schedules [43,44]. The objective of this study was not to present in detail the differences in atmospheric pressure and therefore in evapotranspiration, but to show an effective methodology to determine these parameters, which vary significantly over time. Future studies determining the variations in atmospheric pressure in indoor and outdoor zoological facilities will allow a more adequate interpretation of these parameters and their usefulness beyond greenhouse crops.
Together with the previously described correlations between variables, the PCA analysis for the six environmental parameters recorded by the fixed sensors (ATMOS14 and TEROS12) showed that air relative humidity contributed very poorly to the data variance, atmospheric pressure contributed moderately, while soil water content, soil conductivity, soil temperature and air temperature contributed greatly to data variance ( Figure 7). As can be appreciated in Table 4, even when evaluating the principal components with a variance over 1.000 (principal components 1, 2, and 3), they only represent 80.408% of the total variance, being a limitation of the PCA analysis in our study. It should also be considered that data variance, and the contribution of each variable to the principal components, could vary greatly depending on the conditions of each facility. In our study, the results of the PCA analysis could have been influenced by the climate control tools available at Biodomo: air relative humidity was strictly regulated by a climate control console linked to a fogging system, reducing data variance. Atmospheric pressure was not regulated, and its variance could significantly differ between studies. Despite air temperature being regulated automatically by a climate control system, this variable, together with the rest of environmental parameters studied, were greatly dependent on the daily management of the facility and the non-automated environmental control elements previously described in the methodology section of this manuscript. Although this was not the objective of our study, further trials could be performed evaluating in detail how the different auxiliary environmental parameters control tools affect each of the measured variables, so that the use of these elements can be optimized.
Together with the evaluation of environmental parameters described in this manuscript, further environmental studies performed in zoological and botanical gardens could benefit from the determination and quantification of both natural and anthropogenic air pollu-tants, particularly particulate matter [45]. These complementary studies will assist in the evaluation of indoor air quality, as indoor pollution has been shown to have severe health impacts and should be considered, especially in those facilities built near large urban areas. In addition, modern indoor zoological and botanical gardens such as Biodomo should operate with air treatment units equipped with air filters, which help control these harmful particles.
Ideally, desired environmental parameters should be strongly considered before constructing the building that will host the animal and plant collections, as it will facilitate the posterior maintenance of the optimal ambient conditions [24]. Once the facility is built, the development of environmental studies, such as the one proposed in this article, will allow determination of the effectiveness of the design and the need to implement artificial light sources, or modify the rest of the environmental conditions. Further studies determining environmental parameters such as DLI in facilities housing delicate animal species which are highly dependent on their environment such as corals, fish, reptiles, amphibians, and many avian and mammal species will provide important information for their successful care and management.

Conclusions
This study provides a simple and efficient method for the evaluation of environmental parameters in zoological and botanical gardens housing plants and animal species in indoor facilities. Studies determining environmental parameters are frequently developed in livestock and horticultural production facilities, while this type of studies is still uncommon in zoological/botanical facilities. This study applies a combination of the environmental monitorization systems used in plant and animal production facilities, for the design and improvement of climate management in indoor zoological and botanical gardens. The detailed methodology description provided in this study can be useful in modern zoological facilities, in which the optimal development and welfare of the species maintained is increasingly important. Results provided in this study show the importance of each parameter determination for the optimization of the facility design, species selection and location. Data Availability Statement: Data provided in this study will be made available upon request to the authors.