scieee AI-readable full text Open interactive document viewer

Evaluation of environmental comfort in a social housing prototype with bioclimatic double-skin in a tropical climate

López Escamilla, Álvaro; Herrera-Limones, Rafael; León-Rodríguez, Ángel Luis

Abstract

The aim of this paper is to broaden understanding of the social housing prototype that students and researchers at the University of Seville developed in response to tropical climate conditions to achieve the maximum level of hygrothermal comfort through the use of passive retrofit strategies. This prototype, known as Aura 1.0, was presented at the Solar Decathlon competition where hygrothermal conditions inside the prototype were monitored over a ten-day period. Nevertheless, as the monitoring period of the environmental conditions was so reduced, the aim of this research is to assess the environmental behaviour of the prototype for a full year, under tropical climate conditions. To do so, energy and environmental simulation instruments will evaluate the performance of the architectural proposal and its passive construction systems, in particular the double skin designed with different objectives (architectural, construction, socio-economic). The simulation model has been verified by the data monitored during the exhibition phase and the performance of the building. Various hypotheses have been analysed which demonstrate that employing a double skin in the dwelling contributed to improving interior thermal comfort, especially in the roof area as opposed to the facades, given the sun’s path at a latitude close to the equator.

Full text

Building and Environment 218 (2022) 109119 Available online 26 April 2022 0360-1323/© 2022 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). Evaluation of environmental comfort in a social housing prototype with bioclimatic double-skin in a tropical climate ´ Alvaro L´ opez-Escamilla a , b , Rafael Herrera-Limones a , b , * , ´ Angel Luis Le´ on-Rodríguez a , c a Instituto Universitario de Arquitectura y Ciencias de la Construcci´ on, Escuela T´ ecnica Superior de Arquitectura, Universidad de Sevilla, Av. Reina Mercedes 2, Seville 41012, Spain b Research Group HUM 965, Universidad de Sevilla, Spain c Research Group TEP 130, Universidad de Sevilla, Spain ARTICLE INFO Keywords: Thermal comfort Solar Decathlon Simulation Envelope Tropical climate Adaptive comfort Double skin façade ABSTRACT The aim of this paper is to broaden understanding of the social housing prototype that students and researchers at the University of Seville developed in response to tropical climate conditions to achieve the maximum level of hygrothermal comfort through the use of passive retrofit strategies. This prototype, known as Aura 1.0, was presented at the Solar Decathlon competition where hygrothermal conditions inside the prototype were monitored over a ten-day period. Nevertheless, as the monitoring period of the environmental conditions was so reduced, the aim of this research is to assess the environmental behaviour of the prototype for a full year, under tropical climate conditions. To do so, energy and environmental simulation instruments will evaluate the performance of the architectural proposal and its passive construction systems, in particular the double skin designed with different objectives (architectural, construction, socio-economic). The simulation model has been verified by the data monitored during the exhibition phase and the performance of the building. Various hypotheses have been analysed which demonstrate that employing a double skin in the dwelling contributed to improving interior thermal comfort, especially in the roof area as opposed to the facades, given the sun’s path at a latitude close to the equator. 1. Introduction The AURA 1.0 prototype was designed and built by a team made up of teachers, researchers, and final year students and recent graduates from the University of Seville’s School of Architecture (Spain) [1], in collaboration with the University of Santiago de Cali, (Colombia). This prototype was developed to participate in the first Latin American edition of the Solar Decathlon competition which took place in Santiago de Cali (Colombia) in 2015. Various editions of this contest had already been held in the US (origin of the contest) [2], as well as in Europe [3] and Asia [4]. In this competition, dwelling prototypes are built and then exhibited [5]. They are presented and explained by the participating university teams who are previously chosen in a competition process. These prototypes undertake 10 tests whose scores are either awarded by a jury of experts (Architecture, Engineering and Construction, Sustainability, Communication and Marketing, Urban design, Innovation) or through objective measurements (Comfort conditions, Energy efficiency, Energy balance, House functioning). This foments interdisciplinarity as researchers and students from different branches of knowledge work together with a common objective: design, build and test a sustainable dwelling with optimal interior environmental conditions [6]. This means that the dwelling prototypes can be used to carry out research into the suitability of different prefabricated construction systems [7], broadening visual comfort through natural lighting [8], the behaviour of specific installation systems from different prototypes in the same edition of the competition [9], or testing whether a particular system of energy production for self-consumption is marketable from an investment-performance perspective [10]. Other research has aimed for an adequate integration of energy efficient systems into the developed architectural concept [11]. The competition is also a good context for the introduction of reflections and debates as to whether comfort in a sustainable dwelling * Corresponding author. Instituto Universitario de Arquitectura y Ciencias de la Construcci´ on, Escuela T´ ecnica Superior de Arquitectura, Universidad de Sevilla, Av. Reina Mercedes 2, Seville 41012, Spain. E-mail address: [email protected] (R. Herrera-Limones). Contents lists available at ScienceDirect Building and Environment journal homepage: www.elsevier.com/locate/buildenv https://doi.org/10.1016/j.buildenv.2022.109119 Received 14 February 2022; Received in revised form 11 April 2022; Accepted 18 April 2022 Building and Environment 218 (2022) 109119 2 should be entrusted solely to passive strategies or whether, on the other hand, it is worth employing technology for active retrofitting when environmental conditions are extreme [12,13]. This aspect has been analysed in the different climates where the competition has been held: both tropical, in Colombia [14], and continental Mediterranean climates, such as Madrid (Spain) [15], highlighting that a combination of active-passive retrofitting strategies gives better results. Nevertheless, Solar Decathlon has a few limitations related to the long term analysis of housing behaviour [16]. This is due to the fact that the prototypes built for the competition are assessed in two phases: a first phase which tests the prefabricated construction system designed by each team because the assembly of each prototype dwelling unit must be carried out in a maximum of fifteen days in the so-called “solar city”; the second, 10-day long phase when the prototype is exhibited to the public and when the different contests are assessed, among them the comfort conditions. This means that environmental behaviour is only monitored for ten days, which also, logically, coincides with a very specific climate period, depending on the time of year and the climate in which the competition is being held. . The AURA 1.0 prototype took part in the first edition of the Latin American and Caribbean Solar Decathlon competition which was held in Santiago de Cali (Colombia) in 2015. The monitored data in two of its interior spaces (living room and bedroom) were very appropriate in terms of hygrothermal comfort [17] which led to this prototype winning the “Comfort Conditions” contest. Nevertheless, this is a social housing project, the aim of which is to achieve long term hygrothermal comfort for inhabitants, thereby minimizing energy poverty while guaranteeing that people with fewer resources can reside in an adequate space from the perspective of health and hygiene [18]. Consequently, the main objective of this article is to analyze what the environmental behaviour of the AURA 1.0 prototype would be over a year-long period. To this end, a simulation model was developed, verified by the monitored data, which allowed us to carry out calculations for the prototype’s hygrothermal behaviour based on five different hypotheses. Firstly, these hypotheses consist of the evaluation of the prototype with simulation instruments in the same conditions of use, interior loads and natural ventilation as were given in the competition (exhibition phase open to the public) but extrapolating them to a full year. Next, its behaviour during normal residential use (dwelling for a family of three or four occupants) will be analysed. In this case natural ventilation is limited as openings are adjusted to the usual use of the dwelling. Finally, the effect that one of the prototype’s key bioclimatic design strategies produces: the double skin that acts as a solar protection system and was designed for both the roof areas and the south facade. The aim is to quantify the impact of this double skin if it is used solely on the roof or solely on the south facade. As this is a temporary structure which was disassembled after the competition, there is no monitored data available for a year-long period which would give greater knowledge of behaviours. Neither is it possible to monitor other hypotheses related to thermal load or design. Therefore, the aim of this research centres around the analysis of the comfort conditions in the AURA 1.0 prototype in the city of Cali, with a tropical climate, according to K¨ oppen climate classification [19], in the long term and using proven simulation instruments which make it possible to calculate the different proposed hypotheses. During the exhibition period in Cali, the organisation’s assessment of the interior environmental conditions of the habitable spaces was made in a simple manner. Interior environmental conditions (dry temperature and relative humidity) only had to be maintained in a specific range. Outside of this static range, scores for the contest were lower or void. For the temperature, the range for maximum scores was between 24 ◦C and 28 ◦C. Relative humidity had to be below 60%. For this research, assessment of environmental conditions for the three simulated hypotheses has been carried out with widely recognised, adaptive comfort criteria which are set out in the European Standard EN 16798-1:2019 “Energy performance of buildings. Ventilation for buildings. Part 1” [20]. 1.1. Description of AURA 1.0: from urban scale to dwelling unit The AURA 1.0 prototype is a dwelling unit which resulted from the urban planning of one of the most socially and environmentally excluded areas of the city of Santiago de Cali. It was here that the AURA strategy was developed, focusing on reconditioning deprived urban spaces, with the premise of progressive and sustainable growth, and participative, guided self-construction. AURA 1.0 is developed at four scales: urban scale, building scale, housing scale, and exhibition prototype. A network of basic services is woven between them (again, at different scales) which makes the two base strategies of the project possible: sustainable growth in phases and participatory construction. Thus, the prototype which was built in the Solar City is the representation of a dwelling unit extracted from a project whose aim is the urban and social regeneration of the urban fabric in one of the most deprived areas in Cali. This dwelling has a 50 m 2 surface area and an 18 m 2 private hallway which opens onto the common areas. Access to the dwelling is through the north facade, where there is a covered gallery. The south facade is protected from solar radiation thanks to a double skin which acts as ventilated façade made out of cane (Fig. 1). This double skin also extends over the roof. The inside of the dwelling can be configured in different ways thanks to its flexible design in which spaces can have different functions throughout the day, depending on domestic needs. Thanks to the multifunctional furniture (Fig. 2), a space could be a bedroom at night, but a study or games room during the day, or it could serve to increase the space in the living or dining room. However, during the competition, the dwelling was exhibited as a completely diaphanous space; the furniture was set out in a way that meant the rooms were not compartmentalised. The AURA 1.0 prototype was therefore designed to overcome the characteristic tropical climate conditions in Santiago de Cali [21,22], located in the Northern hemisphere at a latitude of 3◦27′′ N, a longitude of 76◦32′′ W and 995 m above sea level. Based on the Givoni diagram [23] (Fig. 3), it is possible to analyze graphically the conditioning strategies required [24,25] in the climatic conditions of Santiago de Cali: on the one hand, internal gains will be necessary in the hours where the temperature is lower and, on the other hand, solar protection mechanisms must be available in the hours where there is greater radiation. With that in mind, a series of passive retrofit strategies, consistent with the location, were included such as protecting the south facade from the sun or creating a covered but open space like the entranceway which acted as both a sound barrier between houses and encouraged social relationships among neighbours. The competition was held in December, for this reason two solar protection strategies were designed in the prototype: to the North (facade without solar radiation in this month) a flight that served as protection to the entrance gallery to the house was arranged and, to the South (facade that receives direct solar radiation in the period of the year), a double skin was designed that continues on the roof, due to the verticality of solar radiation at this latitude. 1.2. Double skin on the south facade and roof The design of a double skin facade (DSF) is a passive conditioning strategy successfully tested in tropical climates [26,27], not only as solar protection but also to promote natural ventilation [28], although its performance is directly related to the materials used [29]. However, in the case of Santiago de Cali, not only the tropical character of its climate should be taken into account, but also its latitude. ´ A. L´ opez-Escamilla et al. Building and Environment 218 (2022) 109119 3 The impact of the equatorial latitude on architecture [30], as in the case of Cali, makes it necessary to analyze, among other aspects, the repercussion of the double skin on the roof [31,32], due to the high inclination of solar radiation, in addition to assessing the best orientation where to place the double skin on the facade (DSF), since both the North and South facades will receive direct solar radiation in different months of the year (Fig. 4). In the AURA 1.0 building scale design, a skin was developed and handcrafted from local materials (Fig. 5) which gave uniformity to the south elevation of a building based on safe, guided self-construction, guaranteeing that the final elevation would not be personalised by the tenant of each dwelling. Furthermore, this skin fomented the local economy, using natural materials which are not usually being used to generate wealth among the population. As a consequence, value is added to the work of the local carpenter who made the pieces for the facade. These pieces consist of different sized frames which make up the design of the facade. They were made from carved natural wood to which pieces of cane were attached mechanically, having previously been treated with antibacterial varnish and linseed oil. Fig. 1. South elevation of the AURA 1.0 prototype and cane panels as solar protection. Fig. 2. Multifunctional furniture. ´ A. L´ opez-Escamilla et al. Building and Environment 218 (2022) 109119 4 2. Materials and methods 2.1. Data monitored by the competition organisers The measurement process was carried out by the competition organisers during the 10 days that the contests lasted (from 5th to December 14, 2015). To do so, two temperature sensors were installed inside the prototype and three humidity sensors were placed in strategic locations in the dwelling (Fig. 6). During the competition, the prototype was considered an “exhibition pavilion” as it was visited in particular timeslots by a variable number of people, both expert and non-expert in the field. During the period of public exhibition, it was a completely open, ventilated space. Outside of this period, the prototype was closed, and this is when its interior environmental conditions were assessed. As energy consumption was measured and was scorable in another of the contests (Energy Balance), the heating and cooling systems were activated, or not, depending on needs. The best buildings, from an environmental aspect, were those that maintained interior conditions in a range with the lowest possible use of the active systems. The constructive solutions used (Table 1) were designed with the objective of using local resources, carrying out a study of existing companies within a radius of less than 100 km in order to propose solutions that could be transferred to society. Outside of visiting times, the dwelling was usually occupied by members of the team, for example, during the celebration of one of the competition’s contests. On the other hand, these timetables and this activity changed between weekdays and weekends. The prototype timetable on exhibition days is set out in Fig. 7. The number of visitors in the dwelling was not the same at 10am as 3pm, so although this timetable established by the competition organisers helps us to interpret the data obtained and verify the simulation data, it should be considered that they were achieved under very variable conditions. On the other hand, to measure exterior conditions, the organisation had a weather station which collected the values of temperature, humidity and lighting in real time on the site where the competition was held. It should be highlighted that during the competition period, there Fig. 3. Givoni diagram for Santiago de Cali’s climate. Fig. 4. Recorrido solar en Santiago de Cali. ´ A. L´ opez-Escamilla et al. Building and Environment 218 (2022) 109119 5 was no precipitation. Fig. 8 shows the starting data for the development of the validated simulation model: outdoor temperature and indoor temperature measured by the sensors placed by the organization. An average was made with the interior temperature values taken from each of the sensors located in the dwelling. 2.2. Simulation model To address the proposed objective, a simulation model of the study’s prototype, AURA 1.0 was developed. This model has been verified according to the temperature data obtained from the measurements taken by the competition organisers via a system of monitoring implemented in all the participating prototypes. SG SAVE version 2.9.2.1 was used to generate the simulation model. This uses the EnergyPlus [33,34] calculation engine, considered to date, the most advanced software in existence. It was developed by the U.S. Department of Energy’s National Renewable Energy Laboratory [35, 36]. The calculation methodology used by EnergyPlus, as well as the governing equations for heat and mass transfer [37,38] are available online [39]. SG SAVE uses Sketchup 2017 modelling software which makes it possible to model the spaces to be simulated in a very intuitive manner. A single space has been generated for the simulation model because, as previously explained in the prototype description, the rooms were not compartmentalised. The corresponding openings and shade elements used on the south facade have been added to the single space. To verify the simulation model and calculate the different hypotheses, a climate record from the Alfonso Bonilla Arag´ on International Airport located 20 km to the northwest of Cali (Colombia) was used. This is available on the web http://climate.onebuilding.org managed by the researchers Dr. Drury B. Crawley (Bentley Fellow/Director Building Performance Research at Bentley Systems, Inc., ASHRAE Director-At- Large, Vice President at IBPSA) and Dr. Linda Lawrie (Principal Software Design Engineer, Vice President at IBPSA). ASHRAE: American Society of Heating, Refrigerating & Air- Conditioning Engineers. IBPSA: International Building Performance Simulation Association. 2.3. Verification of the simulation model To calibrate the model, the temperature values collected by the competition organisers from the weather station, on the days chosen to verify the simulation model were used. Working with a space whose conditions, in regard to interior loads, were very variable for the previously mentioned reasons, the days on which the building had free evolution without active heating and cooling systems have been taken as the reference to carry out the verification. The simulation instrument has recreated the prototype’s operational conditions on the chosen days, bearing in mind the occupation and loads that these provided, as well as the renewal of interior air and the equipment and lighting loads. The exterior conditions used in the climate record were the same as those registered by the competition organisers’ weather station on the days selected for the verification of the model. Fig. 9 shows the comparison of temperatures measured in situ and the simulated temperatures. Although the graph shows how the data extracted from the simulation model follows the trend of the data of temperatures measured in the prototype, an analytical verification is also carried out. This is to ensure the model’s verification complies with the demands of ASHRAE Guideline 14:202 (ASHRAE) [40], which indicates that the calibrated model must have a normalized mean error (NMBE) of 55 and a Coefficient of Variation of Root-Mean Squared Error (CV)(RMSE) of 15% when monthly data is used for calibration, or requirements of ±10% or 30% respectively when hourly data is used, as in this case [41]. To calculate both coefficients, the expressions stated in section 5.2.11.3 of the guideline. CVRMSE =100 ×[∑(yi−yi)2/(n−p)]1/2/y NMBE =∑ n (yi−yi) (n−p) × y×100 The results of these calculations for the values obtained in the calibrated simulation model are: NMBE =- 0,393% CV(RMSE) =2,357% Fig. 5. Construction section of the double skin façade. ´ A. L´ opez-Escamilla et al. Building and Environment 218 (2022) 109119 6 Therefore, and according to the values mentioned previously, we can confirm that the simulation model is verified by the said method. 2.4. Hypothesis for the case study Once the simulation model is verified with the results obtained from the monitoring carried out during the exhibition period in the “Solar Decathlon City” in Cali (Colombia), different operational conditions were defined in accordance with each hypothesis: •Hypothesis nº1 (Exhibition mode): Prototype AURA 1.0 in the same operational conditions as during the competition. Building with high occupancy and open to the exterior. Parameters: 2 ren/h constant infiltration - Maximum ventilation 60 ren/h but variable depending on time of day. •Hypothesis nº2 (Dwelling Mode): Operational conditions of residential use. Limited and controlled openings. Low occupancy. Parameters: Constant ventilation: 0.85 ren/h - Summer night-time ventilation: 3.15 ren/h •Hypothesis nº3 (Effect of the double skin solely on the south facade): This hypothesis is based on hypothesis nº 2 (closed mode dwelling), but solar protection is limited to the south facade. Parameters: Same as Hypothesis nº2 Fig. 6. Location of sensors in the AURA 1.0 prototype. Table 1 Materials and thermal properties of the construction elements. Element Layer U (W/m 2 K) Espesor (m) Enclosure Fiber cement panel 0.30 0.15 Insulation OSB Board Floor Corrugated sheet 0.29 0.159 Cellulose insulation Corrugated sheet Foamed polyethylene sheeting Fiber cement panel Wood Ceiling Asphalt blanket 0.29 0.165 Fiber cement panel Insulation OSB Board Textile false ceiling ´ A. L´ opez-Escamilla et al. Building and Environment 218 (2022) 109119 7 •Hypothesis nº4 (Effect of the double skin solely on the roof): This hypothesis is based on hypothesis nº 2 (closed mode dwelling), but solar protection is limited to the roof. Parameters: Same as Hypothesis nº2 •Hypothesis nº5 (Total elimination of the double skin): This hypothesis is a combination of nº 3 and nº 4 and therefore the double skin is eliminated from both the south facade and the roof. Parameters: Same as Hypothesis nº2 In the 5 hypotheses that were studied, the behaviour of the prototype is analysed on the basis of the following parameters: occupancy, ventilation and design of the solar protection on the roof and south facade. All hypotheses have been studied for a full one-year period to analyze the thermal condition behaviour in the long term and in different seasonal Fig. 7. Activity times in the prototype. Fig. 8. T.outdoor and T.indoor measured by the weather station and sensors inside the dwelling during the monitoring period. Fig. 9. T.outdoor and T.indoor in free evolution, in comparison to the simulation model T.indoor. ´ A. L´ opez-Escamilla et al. Building and Environment 218 (2022) 109119 8 weather conditions. However, it should be pointed out that given the characteristics of a tropical climate, there are no large hygrometric oscillations between seasons (Fig. 10). Hypothesis nº1 aims to assess the AURA prototype over a year long period under the same operational conditions as the exhibition prototype and with which the verification model was carried out. For hypothesis nº2, the dwelling’s ventilation stream was modified with respect to hypothesis nº1. A maximum occupancy of 20 m 2 /person (in this case that would correspond to 3 people), and a 70 W/person load were determined appropriate for this type of housing. The profile occupancy timetable was as shown in Fig. 11 for weekdays, while a maximum constant occupancy was considered at weekends. This would be the standard hypothesis for social housing use. Hypothesis nº3, nº4 and nº5 are based on the operational conditions of hypothesis nº2, but different parts of the solar protection have been eliminated to see the effectiveness of this on the prototype’s thermal envelope (facade and roof). Specifically, in Hypothesis nº 3, the double skin would only be on the south facade. In Hypothesis nº4, the double skin would only be on the roof. Finally, in Hypothesis nº5 the double skin would be eliminated from the façade and the roof (Fig. 12). The results of the interior temperature obtained in these hypotheses are evaluated according to the adaptive comfort methodology set out in European Standard EN 16798-1:2019 “Energy performance of buildings. Ventilation for buildings. Part 1” [20], and a range of comfort is established supposing a value of 10% dissatisfaction. Regarding the relative humidity results, the same standard establishes a maximum relative humidity of 60% and minimum relative humidity of 25%. 3. Results and discussion After the calculation of the 5 hypotheses, the data obtained for temperature and humidity are analysed. This analysis will be carried out on a monthly basis, since atmospheric conditions are very stable throughout the year and there are no significant differences between the different seasons. However, it is convenient to take into account the solar path, since depending on the seasons of the year, the facade that receives direct solar radiation will vary: South in the months closest to the winter solstice, and North in the months close to the summer solstice. 3.1. Hypothesis nº1 In this hypothesis, a simulation is carried out with an open use of the dwelling, as it was exhibited in the competition. The results obtained simulating the thermal behaviour of the AURA 1.0 prototype over a full year are as follows (Fig. 13): The percentage of hours of comfort per month are shown in Table 2, indicating that behaviour is usually homogenous. From the thermal perspective, October and November stand out while July and August are more relevant when looking at relative humidity data. In this hypothesis, the thermal behaviour of the dwelling guarantees an average number of thermal comfort hours per year of over 40%. Relative interior humidity should stay, where possible, under 65%, this is achieved in more than 55% of the hours of the year. It is a relevant fact that in the months close to the solstices (December–January and June–July), the percentage of hours of thermal comfort is lower than the annual average. It should be pointed out that the occupancy load is very high in this hypothesis. This is due to the fact that the aim is to evaluate the prototype with the same load conditions as the exhibition phase where the average flow of public was high. 3.2. Hypothesis nº2 In this case, a simulation is carried out for a conventional use of a Fig. 10. Graph showing hygrothermal condition statistics in the city of Cali over a one-year period, obtained from the climate record detailed in point 2.2. Fig. 11. Occupancy profile on weekdays. ´ A. L´ opez-Escamilla et al. Building and Environment 218 (2022) 109119 9 dwelling and a much lower occupancy load than in hypothesis nº1 (20 m 2 /persona). The result, graphically (Fig. 14), is as follows: According to the data shown in Table 3, for hypothesis nº2, where a more closed dwelling use is calculated, the thermal behaviour of the dwelling guarantees a percentage of thermal comfort hours of around 55%, where the dwelling behaves better in the summer months. Regarding relative interior humidity of the dwelling, there is also a percentage of almost 55% of hours per year in which the dwelling has a humidity of less than 65%. In this hypothesis, the months close to the summer solstice have better thermal performance, and in the months close to the winter solstice (when the south façade receives direct solar radiation) the prototype has a worse thermal performance. 3.3. Hypothesis nº3 In this case, a simulation is carried out with the same operational Fig. 12. Exterior geometry of the 5 hypotheses in this study. Fig. 13. Hypothesis nº1, thermal behaviour over a full year. Table 2 Monthly hygrothermal behaviour of hypothesis nº1. nº hours within the thermal comfort range % Hours within the thermal comfort range nº hours RH <65 % Hours RH <65 January 307 41.26 457 61.42 February 266 39.58 448 66.67 March 357 47.98 312 41.94 April 331 45.97 331 45.97 May 353 47.45 323 43.41 June 300 41.67 402 55.83 July 289 38.84 544 73.12 August 325 43.68 519 69.76 September 314 43.61 447 62.08 October 358 48.12 436 58.60 November 350 48.61 281 39.03 December 276 37.10 370 49.73 Total 3826 43.68 4870 55.59 ´ A. L´ opez-Escamilla et al.