A New Method to Determine How Compaction A ﬀ ects Water and Heat Transport in Green Roof Substrates

Featured Application: The method presented in this paper aims to estimate the hydrodynamic and thermal properties of green roofs after settling has occurred. Abstract: Although compaction a ﬀ ects water and heat transport processes in porous media, few studies have dealt with this problem. This is particularly true for substrates, which are artiﬁcial porous media used for engineering and technological solutions, such as in vegetated or green roofs. We propose a methodology to study the e ﬀ ect of substrate compaction on the characterization of physical, hydrodynamic and thermal properties of ﬁve green roof substrates. The methodology consists in a parametric analysis that uses the properties of a substrate with known bulk density, and then modiﬁes the substrate properties to consider how compaction a ﬀ ects water and heat ﬂuxes. Coupled heat and water transport numerical simulations were performed to assess the impact of the changes in the previous properties on the hydraulic and thermal performance of a hypothetical roof system. Our results showed that compaction reduced the amplitude of the ﬂuctuations in the volumetric water content daily cycles, increasing the average water content and reducing the breakthrough time of the green roof substrates. Compaction changes the thermal behavior of the green roof substrates in di ﬀ erent ways for each substrate due to the dependence of the air, water and soil fraction of each substrate.


Introduction
Green roofs have been extensively adopted over the last few decades as a technological solution for sustainable development [1][2][3], as they integrate vegetation into buildings to minimize some negative impacts of urbanization [2][3][4][5]. A properly designed green roof can reduce pollution, noise levels, building energy consumption, and stormwater runoff [2,6,7]. At the same time, green roofs increase biodiversity and provide ecological services in urban environments [8,9].
Green roof overall performance strongly relies on the interaction between environmental conditions and the green roof's biophysical components [2,5]. A systematic comprehension of the processes occurring in every component of a green roof, under site-specific conditions, is required to understand the benefits described above. This systematic comprehension is extremely important in arid and semi-arid areas, where construction materials may be different and meteorological conditions significantly differ from humid climates, where majority of green roof investigations have been carried out to date [1,2,[10][11][12].
The vegetation and substrate layers are the central components that determine water and heat fluxes on green roofs [2,3,12,13]. This work focuses on green roof substrates, even when it is known that vegetation is relevant on green roof hydraulic and thermal performance. Substrates are engineered media with an improved performance compared to natural soils [14]. Substrates have been typically used in many civil engineering applications. The substrate provides resources that are critical for the vegetation, such as a support media, nutrients and water [2,15]. The thermal and water retention properties of a green roof greatly rest on the substrate's physical properties [16]. The heat flux through a green roof and the substrate's moisture content affect the buildings' heat gains and losses, which impacts the consumption of energy in the building. Also, substrate's moisture content is key for the survival of the vegetation, especially in arid and semi-arid areas [1,2]. Therefore, it is important to comprehend the physical, thermal, and hydrodynamic substrate properties and how they change through time.
Even when many investigations addressing green roofs thermal performance consider the impacts of moisture on the substrate thermal properties, they commonly neglect how water flow affects the energy balance [17,18]. Sailor [17] developed a widely used model that represents heat fluxes across the substrate and includes the atmosphere-vegetation-substrate energy exchange. Sailor's [17] model assumes that moisture in the substrate diffuses at a constant rate. However, it is widely known that there are highly non-linear processes occurring in this system [2,19]. Sailor and Hagos [18] and Zhao et al. [20] incorporated how compaction affects the substrate thermal properties. However, their model did not account for the effect of water flow on the energy balance.
Soil compaction plays a key role in solute transport, nutrient availability, crop productivity, and water movement [21]. Soil compaction has been typically investigated without considering its impact on the coupled water and heat transport in porous media [22,23]. There are many investigations that deal with the impact of compaction on water fluxes [21,[24][25][26], as well as many studies that investigate how compaction affects heat fluxes in soils [18,20]. Nonetheless, there is a lack of studies that have investigated the effect of substrate compaction on the combined water and heat transport in green roofs. Hence, substrate compaction in green roofs has been typically neglected. However, studies demonstrate that through time, the substrate settles and changes the thermal and hydraulic behavior of green roofs [7,12,27].
The objective of this study is to develop a methodology that can be used to represent the effect of compaction on the hydrodynamic and thermal properties of a porous medium. The proposed methodology combines the Richards' [28] and the energy equations to predict water and heat fluxes in a compacted medium. Using numerical simulations populated with experimental data obtained from green roof substrates, the proposed methodology is used to predict water and heat fluxes in a hypothetical roof system and to determine the heat needed to maintain a constant reference temperature in a room located below the roof.

Materials and Methods
In this section, we briefly introduce the water flow and heat transport theory in porous media. Then, we describe the theory used to incorporate the effect of compaction on the hydrodynamic and thermal properties of porous media, which can then be used to determine the impact of substrate compaction on water and heat transport. Later, we present the substrates used to apply this methodology. Finally, the numerical simulations used to investigate the impacts of substrate compaction on water and heat transport through a hypothetical roof system are described.

Water Flow and Heat Transport Theory in Porous Media
Water transport in porous media usually is represented using the two-dimensional version of the Richards' [28] where θ is the volumetric water content; K(θ) is the hydraulic conductivity curve; h is the pressure head; x and z are the horizontal and vertical directions, respectively; and t is the time. The water retention curve, h(θ), and the hydraulic conductivity curve, K(θ), are two important properties that control water flow in a partially saturated media [15], and are needed to solve Equation (1).
There are many expressions that can be used to represent the h(θ) [29][30][31]. The most used equation to represent the h(θ) is that proposed by van Genuchten [30]: where S e is the effective saturation; θ r , and θ s are the residual and saturated volumetric water contents, respectively; α is the inverse of the air-entry pressure; and n and m are empirical parameters. Assouline et al. [31] also proposed a water retention model that has been used to investigate the impact of compaction on both h(θ) and K(θ); therefore, it is a relevant model for this investigation. The h(θ) model proposed by Assouline et al. [31] is given in Equation (3): where ξ and µ are two fitting parameters, and h L = 15 bar is the pressure head corresponding to a very low water content. It is also interesting to note that the parameters of the h(θ) models of Assouline et al. [31] and van Genuchten [30] can be related using the following expressions: µ= 0.51 ln 1.1 , There are also many expressions to represent the K(θ). The most used expression is the van Genuchten-Mualem model with m = 1 − 1/n [32]: where K S is the saturated hydraulic conductivity, l is the pore-connectivity parameter (=0.5 [32]), and the S e is determined using the van Genuchten [30] model. Assouline and Tartakovsky [33], on the other hand, describes the K(θ) using the following expression: where γ and Γ are the incomplete and complete Gamma functions, respectively, and The heat flow in a porous medium has been typically described using the energy equation [34]: where C p (θ) and C w are the volumetric heat capacities of the moist media and water, respectively; T is the temperature; q lx and q lz are the liquid flux density in the horizontal and vertical directions, respectively; and λ(θ) is the porous media apparent thermal conductivity, defined as [35]: where β t is the thermal dispersivity; q l is the magnitude of the liquid flux density; and λ 0 (θ) is the thermal conductivity in absence of fluid flow, which in this work has been modeled using the Chung and Horton [36] method. The volumetric heat capacity of the moist soil, C p (θ), is described using the following equation [35]: where C is the volumetric heat capacity; a v is the air content; the subscripts n, a, and w refer to the solid, air, and liquid phases, respectively; and the subscript o refers to organic matter. To represent the water and heat flux through porous media, as described in Equations (1) and (8), both the hydrodynamic and thermal properties of the media must be known.

Impact of Compaction on the Hydrodynamic Properties
To consider the effect of compaction on the water retention curve, the method proposed by Assouline [37] was employed. In this method, a homogeneous porous medium with an initial bulk density, ρ, is compacted to a density ρ c . After the porous medium is compacted, the compacted saturated and residual water contents, θ sc and θ rc , respectively, are calculated as follows: where ρ s is the solid particle density of the substrate. When the porous media is compacted, Assouline [37] proposed to estimate the change in the remainder parameters of his water retention curve as follows: where the subscript c refers to the compacted state, and w represents the effect of the silt content (SC) and the clay content (CC) on the way that the porous medium is compacted. Equations (13) and (14) can be also written using the α and n parameters of the van Genuchten [30] model: Appl. Sci. 2019, 9, 4697

of 19
The effect of compaction on the saturated hydraulic conductivity was estimated using the approach presented by Assouline [38]. He proposed obtaining the saturated hydraulic conductivity for a compacted bulk density using the following expression: with ξ = ξ µ , and where η is a parameter that depends on the coefficient of variation of the water retention model of Assouline [38]: where σ 2 and r G are the second and first moments of Equation (3): The calculation of K Sc /K s requires the utilization of Equations (19)- (21) to find the parameters for the uncompacted and compacted state, and the values of µ and µ c from Equation (14). Using Equations (16) and (17), the parameters of the Assouline [38] model in the compacted state can be related to those of the van Genuchten [30] water retention curve.

Impact of Compaction on the Thermal Properties
After knowing the thermal conductivity curve as a function of S e in an uncompacted porous media, e.g., a relationship such as that presented in Equation (9), it is of interest to characterize the impact of the compaction in the thermal conductivity curve. We proposed to use different theoretical approaches combined with empirical information to include the impact of compaction on the thermal conductivity curve. The Johansen [39] model can be used to represent λ 0 . This relationship depends on S e and on the bulk density (ρ) of the media [40]: where λ 0s and λ 0d are the thermal conductivities in the saturated and the dry states, respectively; and λ e is the Kersten number, defined as [39,41]: where κ is an empirical parameter that depends on the soil type and does not change with ρ [42]. The λ 0s is estimated using the following expression [39]: where λ w is the water thermal conductivity and λ sol is the solid particles thermal conductivity. Equation (24) assumes that the porosity of the porous media is similar to its saturated volumetric water content. The λ 0d is determined using the Ball and Arp model [40], which is suitable for any type of soil: where a is an empirical parameter, ρ is the bulk density, ρ s is the particle density of the porous media, and λ air is the thermal conductivity of air. The thermal conductivity curve in absence of fluid flow for the compacted porous media, λ c (θ), needs some inputs from the hydrodynamic properties: first, to determine λ 0s , the θ s is required, as shown in Equation (24); and second, the lower and upper bounds of the thermal conductivity curve are related to the θ r and θ s that also change as compaction occurs. Also, note that the impact of compaction on the λ c (θ) is represented through the bulk density that is used to determine λ 0d , but also through the change in the θ s described in Equation (11), which modifies the λ 0s .
The compaction methodology proposed in this investigation required the use of empirical data to determine how the thermal conductivity changed as compaction occurred. First, the thermal conductivity curve must be known for a specific ρ, i.e., the values of λ 0 (S e ), λ 0s , and λ 0d were known. Combining this information with Equation (22) allowed the Kersten number (λ e ) to be found for different saturation degrees. Then, Equation (23) was used to find the value of κ that represents the best the λ e for each porous media. Recall that κ does not change with ρ [42]. According to Coté and Konrad [41] the value of κ depends on soil type; well-graded gravels and coarse sands have κ ≈ 4.6, silts and clays have κ ≈ 1.9, and peat has a κ ≈ 0.6. In this study, we used the κ value as a fitting parameter to obtain a λ e (S e ) for each substrate that agreed with the experimental information. Later, Equations (24) and (25) were used to find the value of λ sol that matched better the thermal conductivity curve of each porous medium, and to find the value of a that represented the behavior of all the porous media.
The change in the volumetric heat capacity of the porous media as compaction occurs could be readily considered using Equation (10) by modifying the saturated water content as described in Equation (11).

Substrates
To evaluate the proposed methodology that investigates the impact of compaction on the coupled water and heat flux, empirical information from five substrates was used. The substrates used in this study are those employed by Sandoval et al. [2]. Figure 1 presents photographs of each substrate and Table 1 provides a detailed characterization of their physical properties. These five substrates were chosen because they are commonly used in green roofs (S1-S3), and because they are easily found in the local market (S4,S5). A detailed characterization of the physical, hydrodynamic, and thermal properties of these substrates, for their uncompacted state, is presented elsewhere [2]. For the compacted state, the substrate's hydrodynamic and thermal properties were estimated using the theory presented in Sections 2.2 and 2.3. Table 1. Properties of the substrates: texture determined using the United States Department of Agriculture (USDA) soil classification system, organic matter (OM), soil particle density (ρ s ), maximum density (ρ max ), uniformity (C u ) and curvature (C c ) coefficients for the substrate's uncompacted state. Modified from Reference [2].   Table 1 provides a detailed characterization of their physical properties. These five substrates were chosen because they are commonly used in green roofs (S1-S3), and because they are easily found in the local market (S4,S5). A detailed characterization of the physical, hydrodynamic, and thermal properties of these substrates, for their uncompacted state, is presented elsewhere [2]. For the compacted state, the substrate's hydrodynamic and thermal properties were estimated using the theory presented in Sections 2.2,2.3.

Influence of Substrate Compaction on the Performance of a Hypothetical Roof System
To assess the effects of substrate compaction on water and heat transport through green roof substrates, we carried out numerical simulations in the same hypothetical roof system investigated by Sandoval et al. [2]. As the simulations reported in this study used the same domain and boundary conditions used by Sandoval et al. [2], here we only provide a brief description of the numerical model. Figure 2 presents the conceptual model utilized to investigate the effect of substrate compaction on roof functioning. The system was comprised by a substrate on top of a concrete layer that was above a room. As we focused on substrate compaction, and without diminishing the relevance that vegetation has on the overall green roof functioning, the vegetation layer of this system was not included. Instead, we put our attention in analyzing the effect of compaction on green roof functioning when the five substrates investigated in this study were employed. We also put our attention into a green roof located in a semiarid zone during summer (see meteorological conditions in Appendix A), so our results were focused on the response of the hypothetical roof system subject to irrigation for a period with no precipitation, in which it is relevant to know when the green roof begins to drain (the breakthrough time). For this reason, we did not analyze water retention efficiency, peak reduction, and delay of a precipitation event.
. 2019, 9, x 7 o luence of Substrate Compaction on the Performance of a Hypothetical Roof System assess the effects of substrate compaction on water and heat transport through green r tes, we carried out numerical simulations in the same hypothetical roof system investiga doval et al. [2]. As the simulations reported in this study used the same domain and bound ions used by Sandoval et al. [2], here we only provide a brief description of the numer . gure 2 presents the conceptual model utilized to investigate the effect of substrate compac f functioning. The system was comprised by a substrate on top of a concrete layer that w a room. As we focused on substrate compaction, and without diminishing the relevance tion has on the overall green roof functioning, the vegetation layer of this system was ed. Instead, we put our attention in analyzing the effect of compaction on green r ning when the five substrates investigated in this study were employed. We also put on into a green roof located in a semiarid zone during summer (see meteorological conditi endix A), so our results were focused on the response of the hypothetical roof system sub ation for a period with no precipitation, in which it is relevant to know when the green r to drain (the breakthrough time). For this reason, we did not analyze water reten cy, peak reduction, and delay of a precipitation event. e combined water and heat fluxes (Equations (1,8)) were solved using the HYDRUS re [43], and considering the uncompacted and compacted hydrodynamic and ther The combined water and heat fluxes (Equations (1) and (8)) were solved using the HYDRUS-2D software [43], and considering the uncompacted and compacted hydrodynamic and thermal properties for each simulation. We selected the HYDRUS-2D software to solve the governing equations, instead of more sophisticated simulation software such as EnergyPlus (http://energyplus.net/) or WUFI (https://wufi.de/), because it allows a detailed representation of the hydraulic and thermal processes occurring in the substrates. The uncompacted hydrodynamic properties of the substrates were defined according to the experimental values reported by Sandoval et al. [2], while the compacted properties were determined according to the proposed model. In all the simulations, a uniform and arbitrary initial moisture content of 0.3 cm 3 cm −3 was selected. The flow boundary conditions at the top of the domain were defined using the precipitation, irrigation, and evaporation measured at the Laboratory of Vegetative Infrastructure of Buildings at Pontificia Universidad Católica de Chile [1]. The irrigation schedule used 10 mm/day allocated in three pulses per day at 7:00, 14:00, and 22:00 hours (see Appendix A). The lateral boundary conditions were defined as no-flux boundaries, except where the drainage was located. The seepage boundary condition [44] was used to represent the drain, and the interface between the concrete layer (ceiling) and the room was represented with a boundary condition of no flow. The cumulative water flow through the drain was employed to compare the hydraulic functioning of the substrates.
The thermal initial condition was a uniform temperature (=23 • C). At the domain's top, the surface energy balance was utilized as a boundary condition for heat [44]. The lateral heat transport boundary conditions were adiabatic. The thermal functioning of the system was evaluated by considering that the roof is above a single room with adiabatic walls and floor, and using the cumulative thermal energy (cooling/heating) required to keep the room at 23 • C. This thermal energy balance was estimated as described by Causone et al. [45]. The reader is referred to the work of Sandoval et al. [2] for a detailed description of the simulations under the uncompacted state.  [30] model for different bulk densities, and Table 2 presents the parameters of the water retention models for the uncompacted and compacted states of the substrates. In general, as compaction occurs, the size of the pores decreased, and capillary forces increased.    Table 2 presents the parameters of the water retention models for the uncompacted and compacted states of the substrates. In general, as compaction occurs, the size of the pores decreased, and capillary forces increased.

Impact of Compaction on the Substrate's Thermal Properties
The parameters of the Chung and Horton model [36]-which was used to describe the thermal conductivity of the different porous media, λ 0 (θ)-are presented in Table 3.
To find how the thermal conductivity changes as compaction occurs, the first step in our methodology was to determine the λ e , λ sol , and κ for each substrate, and the a parameter for the set of substrates. Figure 4 presents the λ e as a function of the S e for each of the substrates investigated, estimated by fitting Equation (23) to the experimental thermal conductivity data to find the κ parameter for each substrate; and by fixing the a parameter to 0.1614, which yielded a good representation of the experimental values of λ 0d (Equation (25)) for all the substrates. In this process, λ sol was also fitted to represent correctly the thermal conductivity of each substrate and, in particular, the experimental values of λ 0s and λ 0d . Table 3 also presents the values of κ and λ sol obtained from this fitting process. Table 3. Thermal properties of the substrates studied in this investigation. The first row of each substrate refers to the uncompacted (UC) state, while the second row corresponds to the compacted (C) values obtained with the parametric analysis.

Substrate
State    Figure 5a-e presents the thermal conductivity of the five substrates for both the uncompacted and the compacted state. In general, for a specific moisture level, as compaction occurred the thermal conductivity of the substrates increased.

Influence of Substrate Compaction on the Performance of a Hypothetical Roof System
For conciseness, we solely present the water and heat fluxes exiting the hypothetical roof system. The temporal transition of the volumetric water contents and of the temperatures in the system are presented in Appendices B and C. Figure 6a-d shows the cumulative water flux that exits from the bottom of each substrate during the simulations. Light colors correspond to the uncompacted substrate, and dark compacted substrate. The S2 substrate did not reach its total volume storage in either state, so it is not shown. Figure 7 shows the cumulative thermal energy (heating and cooling) consumed to maintain 23 • C in the room for the entire simulation (30 days). The increment of the thermal daily amplitude, which is the result of substrate compaction, increased the heating and cooling energy due to natural convection. The biggest changes are related to substrates S1 and S2; in both cases the energy for heating was reduced and the cooling energy increased. This behavior occurs due to the mixed effects of the general increment in thermal conductivity with moisture and the high porosity that reduces the vertical diffusion.

Influence of Substrate Compaction on the Performance of a Hypothetical Roof System
For conciseness, we solely present the water and heat fluxes exiting the hypothetical roof system. The temporal transition of the volumetric water contents and of the temperatures in the system are presented in Appendices B,C. Figure 6a-d shows the cumulative water flux that exits from the bottom of each substrate during the simulations. Light colors correspond to the uncompacted substrate, and dark compacted substrate. The S2 substrate did not reach its total volume storage in either state, so it is not shown.  Figure 7 shows the cumulative thermal energy (heating and cooling) consumed to maintain 23 °C in the room for the entire simulation (30 days). The increment of the thermal daily amplitude, which is the result of substrate compaction, increased the heating and cooling energy due to natural convection. The biggest changes are related to substrates S1 and S2; in both cases the energy for heating was reduced and the cooling energy increased. This behavior occurs due to the mixed effects

Discussion
To truly quantify the urbanization negative impacts that a green roof mitigates, a quantitative assessment of each of its components must be performed, and the relevant physical processes must be included in this assessment. One aspect that has been typically neglected is the impact that substrate compaction has on the coupled water and heat fluxes across the roof. In this work, we developed a parametric methodology that used empirical substrate data to consider how compaction changed the substrate's hydrodynamic and thermal properties.
By analyzing the changes in the water retention curve when the substrate was compacted (Figure 3), it was clear that the pores decreased their size, diminishing the θ s and increasing the air-entry pressure (h b = α −1 ) due to an increase in the capillary forces [37]. These changes in the porous media water retention capacity are consistent with the observations published in the scientific literature [2,15,24,25,46,47]. For instance, in the S2 compacted water retention curve, the h b increased from 17 to 526 cm, while the h b of the S1 substrate increased from 9 to 53 cm. The increases in the h b values were on the same order of magnitude as those reported by Assouline [37]. However, the impact of compaction on the θ r was not the same for all the substrates. S1, S2 and S5 showed an increase between 0.05 and 0.07 cm 3 cm −3 in θ r as compaction occurred while S3 and S4 showed an insignificant change. The reason of this change is that the original density values of θ r of S3 and S4 were near to zero and the Equation (12) used the original density to determine the compacted values of θ r . In addition, the changes in n and small changes in θ r drove unexpected retention curves as shown in Figure 3c: the compacted water retention curve of S3 (ρ = 1 g cm −3 ) intercepted the uncompacted curve (ρ = 1.2 g cm −3 ) at~0.15 cm 3 cm −3 , meaning that below this value the water retention decreased. However, as the pore size decreases, an increase in the capillary forces is expected and thus one would expect a general increase in water retention [21]. This phenomenon most likely occurs because the properties of these artificial substrates are very different than those of natural soils, which are the media used to develop the algorithms that predict the effect of compaction on the hydrodynamic properties [37]. Therefore, for substrates S3 and S4 the proposed model may not be the best representation of how compaction influences the hydrodynamic properties. From the results obtained for substrates S1-S5 (Figure 3), we concluded that the proposed method resulted in physically possible values for volumetric water contents that are greater than 0.2 cm 3 cm −3 .
Regarding the thermal properties, the values of κ were within the range of values reported for natural soils, and were typically in the range reported for silt, clay, and peat [41]. The difference between the κ values obtained for the substrates used in this study and those reported for natural soils are most likely explained by the artificial nature of the substrates. The λ sol values depend on the fraction of gravel, sand, silt, clay, and percentage of quartz of each substrate. As a reference, soils with high organic matter have λ sol~0 .25 W m −1 K −1 , while silty loam soils have λ sol~2 .05 W m −1 K −1 , and crushed quartz has λ sol~7 .2 W m −1 K −1 [42]. The λ sol presented in Table 3 ranges between 1.103 and 2.124 W m −1 K −1 . The lower values obtained for the λ sol of each substrate, compared to natural soils, could be justified by the presence of organic matter and inorganic materials. The λ sol values were fitted so the simulated λ 0s and λ 0d were similar to the experimental measures. This agreement between the simulated and experimental values was achieved in all the substrates except for substrate S2, where the fitted λ 0d (0.047 W m −1 K −1 ) was lower than the measured λ 0d (0.101 W m −1 K −1 ). However, the fitted and measured λ 0s of this substrate agreed with a difference of~1%.
In general, for a specific moisture level, as compaction occurs the thermal conductivity of the substrates increase. This behavior is consistent with the findings of other investigations [18,20] and occurs because as the substrate is compacted, there is more physical contact between the solid particles, which increase thermal conduction [20,24,25,35]. Also, the magnitude of the thermal conductivity is similar to that presented in the literature [15,20]. Figure 5a-e also shows that substrates S2, S3, and S5 present a linear increase of thermal conductivity with moisture, contrary to the non-linear behavior observed in substrates S1 and S4. In general, the impact of compaction on the thermal conductivity curve is more relevant in the saturated zone of the curve, i.e., when S e approaches 1. For example, in substrate S2 the λ 0s increased by 0.17 W m −1 K −1 (~30% increase). In the dry zone of the curve, i.e., when S e approaches 0, the substrates' λ 0d presented slight increases, except for substrate S1, which exhibited an increment of 0.15 W m −1 K −1 . Similar effects are reported in other investigations, where the thermal conductivity curves of different soils shows increases of 0.15 and of 0.4 W m −1 K −1 for λ 0d and λ 0s when the soil is compacted [18]. However, previous investigations do not account for the hydrodynamic effects of compaction on the conductivity curves. The addition of the hydrodynamic properties on the compaction method could be clearly seen in substrates S1 and S2, as shown in Figure 5a,b. The highest change between the uncompacted and compacted thermal conductivity curves occurred in substrates S1 and S2. In S1, the increase in thermal conductivity was similar in all the range of the curve; contrary to what was observed for substrate S2, where the increment near the saturated part was too large while the increment near the dry part of the curve was too small. The small increment of the dry thermal conductivity in substrate S2 was due to its very low bulk density (Equation (25)). In addition, the large change of density caused a large increment in θ r (0.23 to 0.40 cm 3 cm −3 ) (Equation (12)). As a result, the compacted and uncompacted curves intersected.
Regarding the hydrodynamic conditions through the hypothetical green roof, our results showed that compaction reduced the amplitude of the fluctuations in the volumetric water content daily cycles (see Appendix B). Also, the cumulative water that exits the roof ( Figure 6) suggested that the breakthrough time occurred earlier in the compacted substrates, most likely due to a reduction in the pore spaces. After compaction, the reduction in the breakthrough time of the substrates ranged between one hour (S3) and three days (S4). This result highlights the nonlinearity effect that compaction has on water fluxes. In all the substrates the amount of water that drains the roof increased, especially in the S3 substrate where the cumulative water flux increased by more than three times. This significant increase in the cumulative water flux was the result of a decrease in the pore spaces and the changes in the response of each substrate to the atmospheric conditions.
The comparison among the thermal conductivity curves of the five substrates in the uncompacted and compacted state-near to their maximum densities as shown in Table 1-can help to understand how the heat fluxes change through time in a green roof (Figure 7). Before compaction, the substrates with higher thermal conductivity were S3 and S4 ( Figure 5). After compaction, the same pattern was observed, except when all the substrates were saturated. At this point, all the substrates reached thermal conductivity values of~0.9 W m −1 K −1 . In the case of the S2 substrate, the saturated thermal conductivity increased from 0.65 (uncompacted state) to 0.93 W m −1 K −1 (compacted state). In addition, the S2 substrate had the lowest thermal conductivity values for the entire range of moisture for both the uncompacted and compacted state. On the other hand, after compaction, the thermal conductivity of substrates S1 and S5 presented similar values for the entire range of volumetric water content. The numerical simulations allow the comparison of the resulting thermal behavior of the substrates (see also Appendix C). It can be seen in Figure 7 that the energy required for heating and cooling for substrate S2 was smaller than that for substrates S1 and S5. Therefore, substrate S2 was more effective in reducing the vertical heat diffusion when compaction occurred. As explained by Sandoval et al. [2], the effect of the volumetric air content on the reduction of the heat diffusion is very important. In the simulations presented in this work, the S2 substrate presented a low average water content (i.e., the pore spaces have more air than water) since the substrate never reached saturation.

Conclusions
In this study, a new methodology to understand how compaction affects the hydrodynamic and thermal properties of a porous medium was developed. This methodology was used to investigate the effects of compaction on the coupled water and heat transport in green roof substrates through numerical simulations. The results showed that as the substrate was compacted, the pore-size of the substrates was reduced. As a result, the capillary forces increased and the water storage volume decreased. These results agree with other investigations performed in natural soils. However, because the foundations of the new methodology are based on compaction methods developed for natural soils that have higher residual water contents and bulk densities compared to the artificial substrates, the new approach failed to produce reasonable results for some substrates. For instance, for substrates S3 and S4 the methodology predicted compacted residual water contents that were lower than those of the uncompacted state. For this reason, the proposed method could be used for green roof substrates that present residual volumetric water contents larger than 0.05 cm 3 cm −3 .
Regarding the substrate's thermal properties, the thermal conductivity curves increased their values as the bulk density increased. This increase in the thermal conductivity occurred because there was more physical contact between the solid particles of the media. The effect of compaction on the thermal conductivity curve was more relevant in the saturated part of the curves for all the substrates.
Also, when considering that compaction also affects the hydrodynamic properties of the substrates, the variations of the thermal conductivity curve are now restricted to the new limits of moisture, which reduces the possible fluctuations of the thermal conductivity. When the substrates have a very low bulk density (<0.5 g cm −3 ), the proposed method can present problems, so care must be taken when analyzing the results.
The numerical simulations showed that the water and thermal behaviors of the roof system changed with compaction, reducing the amplitude of the fluctuations in the volumetric water content daily cycles, increasing the average water content and reducing the breakthrough time of the green roof substrates.
When comparing the overall thermal behavior within the substrates studied, the S2 substrate reduced the vertical heat diffusion through the roof. Compaction changes the thermal behavior of the green roof substrates in different ways for each substrate due to the dependence of the air, water, and soil fraction of each substrate. The results showed that compaction affected the coupled water and heat fluxes in a hypothetical green roof. However, the compaction methodology needs an experimental validation for artificial substrates that have low bulk densities and residual water contents. Future research should focus on performing an appropriate experimental validation of this methodology to further confirm the results from the numerical simulations.    Figure A2a-e shows the temporal evolution of the volumetric water content during the first five days of simulation, at the atmosphere-substrate interface. Light colors correspond to the uncompacted substrate while dark colors are related to the compacted substrate. The water content evolution showed a daily cycle in all the substrates where the water content increased after the roof was irrigated, and then decreased due to evaporation. The effect of compaction on the moisture level is clear: after the substrates were compacted, the fluctuations in the volumetric water content were smaller. The reduction in the moisture fluctuations was clear for substrates S1 and S4. On the other hand, for substrates S2, S3, and S5 the mean water content increased most likely due to the reduction of water storage volume. Figure A2 demonstrates that there was a change in the hydrodynamic response of each substrate to the same meteorological conditions as compaction occurred.   Figure A3a-e shows the modeled temperature evolution of the five substrates in 30 days of simulation and at 34.5 cm below the surface, i.e., near the ceiling of the room. The changes at the roof ceiling when compaction occurs showed a general increment in the daily amplitude of the temperatures, which could be explained by the larger average water contents that can be seen in Figure A2 that provoked higher thermal conductivities and vertical heat diffusion. The change in temperatures at the ceiling due to compaction were mostly noticed in substrates S1 and S2. The temperature daily cycles of the uncompacted and compacted state were incremented by~1 • C in the case of S1 and S2. These changes in the thermal amplitude agreed with the changes in the thermal conductivity curves shown in Figure 5. Appendix C Figure A3a-e shows the modeled temperature evolution of the five substrates in 30 days of simulation and at 34.5 cm below the surface, i.e., near the ceiling of the room. The changes at the roof ceiling when compaction occurs showed a general increment in the daily amplitude of the temperatures, which could be explained by the larger average water contents that can be seen in Figure A2 that provoked higher thermal conductivities and vertical heat diffusion. The change in temperatures at the ceiling due to compaction were mostly noticed in substrates S1 and S2. The temperature daily cycles of the uncompacted and compacted state were incremented by ~1 °C in the case of S1 and S2. These changes in the thermal amplitude agreed with the changes in the thermal conductivity curves shown in Figure 5.