Environmental and Social Dynamics of Urban Rooftop Agriculture (URTA) and Their Impacts on Microclimate Change

Urban cities are facing the challenges of microclimatic changes with substantially warmer environments and much less access to fresh vegetables for a healthier food supply than in adjacent rural areas. In this respect, urban rooftop agriculture is considered as a green technology for city dwellers and the community to attain environmental and socioeconomic benefits in a city. For this purpose, a roof top of 216 square meters was selected as an experimental plot where 70% of the area was covered with the selected crops (Tomato, Brinjal, Chili, Bottle Gourd and Leafy vegetables such as Spinach, Red Spinach and Water Spinach; they were cultivated under fencing panels of Bottle Gourd). The microclimatic parameters such as air temperature, near roof surface temperature, indoor temperature and relative humidity and carbon dioxide concentration from different locations of the agricultural roof and from nearby bare roofs were observed during the whole experimental period (November 2018–May 2019). Five existing rooftop gardens with green area coverages of 40, 50, 60, 80, and 85% were selected, and 5 bare nearby roofs were also selected through field visits and questionnaire surveys of 200 existing rooftop gardens. The air and ambient temperature, cooling degree day and energy saving trends were assessed for the selected roofs. The economic assessment was carried out through the net present value and internal rate of return approach of urban rooftop agriculutre. The results showed that the temperature was reduced from 1.2 to 5.5% in different area coverages of agricultural roofs with plants compared to the nearest bare roofs. For the time being, the cooling load was decreased from 3.62 to 23.73%, and energy saving was increased significantly from 5.87 to 55.63% for agricultural roofs compared to bare roofs. The study suggested that the value of urban rooftop agriculture was high environmentally and economically compared to the traditional bare roof, which would be an added amenity by the city dweller’s individual motivations and state interests, and it could be aligned to achieve a more sustainable city.


Introduction
In a phase of global warming, the urban warming effect is likely to be amplified, especially increasing human discomfort during summer. The local warming, caused by the urban heat island (UHI), significantly increases temperatures as well as economic losses in addition to global warming [1]. The rapid urbanization process plays a key role in the formation of UHI as well as global warming, which impacts the urban quality of life [2,3]. Due to high economic growth and improved living standards, energy demand in urban cities is rising for increasing electricity consumption, mostly for using air conditioning systems in urban buildings [4][5][6][7]. Therefore, finding a way to reduce energy consumption with respect to the cooling load can significantly reduce the heat burden directly and greenhouse emissions indirectly. For colder days (25 • C or 77 • F), a 1 • C increase in daily temperature leads to a 14.5% increase in electricity consumption [8][9][10]. Urban warming could lead to double the economic losses expected from human-caused climate change, and it would be probably comparable to about half of the warming caused by climate change by the year 2050 [11][12][13][14]. Daily minimum temperature readings at related urban and rural sites frequently show that the urban sites are 6 • C to 11 • C (10 • F to 20 • F) warmer than the agricultural sites [15,16]. However, the overheating of both urban building roofs as well as wall surfaces negatively affect the indoor air temperature, which is a crucial factor impacting urban warming, building energy consumption and occupant welfare [17][18][19][20]. In summer, urban masonry and asphalt capture, store and reradiate more solar energy per unit area than the vegetation and soils of rural areas [21]. Furthermore, less of this energy can be used for evaporation in urban areas, which characteristically exhibit greater precipitation runoff from streets and buildings [22,23]. At night, radiative losses from urban buildings and street materials keep the city air warmer than that of rural areas. In addition, human cultural and economic activities have distinctive effects on urban climate warming [24].
Based on the observations of different studies, the air temperatures under green roofs are cooler than normal roofs at least by 3 • C to 4 • C. A study in the city of Toronto also found that an entire typical residential building with a green roof experienced a 25% cooling effect, while the floor below the green roof had a 60% cooling effect [25]. Hence, roof gardens have the potential to act as insulation for the roof because of the heat exchange with the outside environment [26]. Green roofs reduce a huge amount of heat absorption in summer. Thus, less energy is required to keep the indoor air cool [27,28].
Greenery vegetation activities such as urban rooftop agriculture (URTA) are turning out to be a very effective adaptation tool to cool down UHI [19,29]. URTA is a complex interaction between natural systems and human activities [30]. URTA helps to improve the environmental quality and the economic conditions of city dwellers and has the potential to provide fresh vegetables to people for a healthier food supply [31,32]. URTA not only produces and distributes food, but also reduces carbon footprints for urban activities through the energy-saving, efficient management of building resources [33]. Conversely, most of the existing research is based on green version approaches such as parks, gardens, green roofs, green façades/walls, porous the diversion pavements and green and blue belts in the context of urban development and of the potential impacts of global climate change from the city center [34].
Dhaka is a mega city that is about to become the sixth largest city in the world with a estimated population of about 27 million by 2030 [35]. Global warming will increase the severe UHI conditions in Dhaka by creating heat stress, and that will affect Dhaka's quality of life [36]. As a result, Dhaka is anticipating the challenges of having fresh food, water security, healthy lives and a cleaner city to achieve Sustainable Development Goals (SDGs) 2, 6 and 11 [37,38]. On the contrary, rooftop gardens/agriculture have very high social recognition (85%) and is commonly practiced in Dhaka, and it has good economic prospects [39]. However, there is a research gap about the cooling effect of different crop types in URTA, the crop area coverage percentage of URTA and the fresh food production potential of URTA. No previous study has testified based on either roof area coverages by plants or microclimatic parameters dynamics of cooling load.
Addressing the above research background, the objectives of this study were set to: (1) assess the dynamics of microclimatic parameters of ARs and BRs during the cooling period in Dhaka to mitigate the UHI effect; (2) quantify the role of URTA on energy saving through different area coverages of roofs; and (3) summarize the benefits and social impacts of URTA. This study aims to provide a widespread prospect assessment of URTA in the Dhaka Metropolitan Area (DMA).

Study Area and Research Design
Dhaka is a mega city and the capital of Bangladesh, having an area of 1463.60 sq km, boundaried by the Gazipur, Tangail, Munshiganj, Narayanganj, Manikganj and Faridpur districts (Figure 1a,b). The foregoing areal expansion of Dhaka with its dense population has triggered the processes of land transformation because of the growth of urbanization. It is also responsible for the physical and environmental instability of that area. Due to fast urbanization, it is facing the loss of natural vegetation, loss of open spaces and a general decline within the spatial extent and connectivity of wetlands and wildlife habitat. The city of Dhaka is rapidly growing in terms of both population and extent. It is becoming the center of the country's industrial, commercial, cultural, educational and political activities. That is why Dhaka is becoming a warmer city compared to the rural areas. In this regard, the present study was conducted through URTA within the DMA from November 2018 to May 2019. The experimental plot of AR was designed on the roof of the Department of Agriculture, Sher-e-Bangla Agricultural University (SAU), Dhaka, Bangladesh, to cultivate the selected crops (Tomato, Brinjal, Chili, and Bottle Gourd and leafy vegetables) (Figure 1c). The height of the institutional building roof was about 50 ft (about 15 m). The area of the experimental plot was 216 m 2 , 18-m long and 12-m wide, which is physically significant for such experiments ( Figure 2). Thus, a completely randomized design (CRD) was used for designing experimental plots of AR in such a way that at least 70% of the roof area was covered with the selected crops (Tomatoes-16%, Brinjall-15% Chilli-13%, Bottle Gourd-25% and Recreational room-2.3%) (Figure 3). Similarly, 216 m 2 BR was also selected near the experimental AR. However, two rooms were covered with the experimental AR below it and two rooms below the BR, which were also considered in this study. On the other hand, 5 existing rooftop gardens or ARs with an area coverage of 40, 50, 60, 80 and 85% by agriculture and 5 nearby BRs were selected through a survey of 200 existing private rooftop gardens within the DMA. However, the experimental AR was divided into 10 rows to organize URTA with those crops in 150 plastic drums. The drip irrigation method was used for cultivation, and Bottle Gourd was cultivated around the experimental plot using fencing panels. Below the fencing panel, leafy vegetables such as Spinach, Red Spinach and Water Spinach were cultivated in 8 wooden frames where the area of each wooden frame was 2.23 m 2 . Soil was specially prepared with cocoa dust and vermicompost (2:1:1) and plant center to center spacing of those crops was maintained as per the Bangladesh Agricultural Research Institute's (BARI) recommended guidelines.   This research was divided into the following two sections: Sections 1 and 2. Section 1 was devoted to the assessment of the thermal environment changing due to URTA according to the different area coverage of the roofs (temperature, relative humidity and carbon dioxide) and effects on energy saving using cooling degree days (CDD) technique. Section 2 was intended to carry out an assessment of the socioeconomic impacts of URTA with experimental study and also from open-end questionnaire surveys of 200 existing rooftop gardens. The survey was conducted within the DMA, and rooftop gardens were selected randomly in the specific area. It was done by survey through oral interviews of the building/rooftop garden owner to gather more in-depth information, opinions and preferences. The role of food production, priorities of involvement of women, choices of multiple crops, promotion or popularization, growth and trends and economic strengths or values of URTA, etc., were considered in the study to clarify the social dynamics of URTA.
On the other hand, a suitable and easy handling drip irrigation system was set in the entire plot area to deliver the right amount of water from the source of portable water storage with free pumping energy at the root zone of each plant at a regular interval [ Figure 2]. The drip irrigation system prevented the plants from suffering stress or strain of under and overwatering. The research study was implemented in four main phases: (i) plot and soil preparation, (ii) fieldwork, (iii) data collection and (iv) analysis. The study roof was selected as it was most suitable for carrying out the daily research and for collecting the regular data. Then, the growing medium, i.e., soil, was prepared according to the suggested ingredients and then put in the container. The selected existing five greenery roofs showed many similarities in plant types and materials used except for some characteristics as orientation and disguise phase. All the height of building roofs were also similar (50 ft from the ground) and were made of concrete bricks with a plaster coating and a waterproofing system. The microclimatic environmental variations in terms of temperature and relative humidity were measured for the selected greenery roofs. In addition, three main actors involved in this study as a means to understand the in-depth socio-spatial characteristics and specific social practices of URTA. Moreover, the social practices were also studied through continuous, direct observation by visiting the five selected existing gardens two times per week for three months. The primary data such as air temperature (AT), near roof surface temperature (RST), relative humidity (RH) and carbon dioxide (CO 2 ) data were collected from the middle part of the experimental AR plot and from four middle locations of the east, west, south and north edges of the experimental AR plot. Temperature-RH-CO 2 measuring data loggers (HUATO, S653) were used to measure temperature, RH and CO 2 concentration. It was done 5 days per week from November 2018 to May 2019 from the experimental AR and nearby BR at 9:00 a.m., 1:30 p.m. and 5:30 p.m. A digital compact infrared thermometer with a 4-h interval was also used for the measurement of roof surface temperature. The temperature was measured near the roof surface 1 m, 1.5 m and 2 m above from roof surface for both AR and BR for the same height. However, all data (except at 2 m) were collected from below the shaded of the canopy layer of plants so that direct solar radiation could be avoided. The temperature was also measured from the selected five existing ARs and top floor room at 1.5-m height above the roof surface. At the same time, temperatures were measured from the middle of the top floor room of the selected 6 BRs near the experimental AR. RH and CO 2 data were collected only at 1.5-m height above the roof surface for the URTA roof and comparatively to the adjacent BRs in this study. All data were collected at 1.5 m above the roof surface because the average minimum height and branch density of those selected crops varied from 0.9 m to 1.5 m due to human comfort breathing at this height in Bangladesh context [40]. The ambient temperature data from the top floor room under the experimental AR and the BR were also collected for comparing the thermal variation among them.

Temperature Trend Measurement of Soil in the Container and Air under URTA
The temperature trend in the container was measured in the URTA at the lower, middle and upper portion of the container. At the same time, the temperature in the land was measured for the same height of the lower middle, and upper portion of the container for comparing the trend of temperature and find out the cause of roof cooling by URTA. Temperatures were measured at 0.8 m, 1 m, 1.5 m and 2 m from the roof surface to find out the temperature trend in URTA during the whole experiment period for 9:30 a.m., 1:30 p.m. and 5:30 p.m. daily.

Cooling Degree Day (CDD) Calculation
The degree-day approach was directly proportional to a difference between the mean daily temperature of ambient air and indoor temperature. The higher the CDD, the higher was the energy requirement for cooling. Considering the average outside temperature of Dhaka during the summer season, the base temperature for cooling comfort (T base ) was 20 • C [41]. In this study, the principles of CDD were used to study the energy consumption of URTA-containing buildings in Dhaka. About 40-85% area coverage scenarios were explained and calculated for the trend of energy consumption. We determined the corresponding energy and cost savings in those selected URTA locations across Dhaka. The CDDs were calculated for six URTA roofs based on roof area coverage by plants. However, the optimal area coverage and plant density worked as an insulation thickness. A function of CDD and the pay back period of the URTA, as well as insulation costs and other costs, were analyzed, which addressed the comparison between URTA roofs and roofs without URTA. The following equation was used for the calculation of the CDD. If T max < T base , CDD = 0; If the average value of the minimum and maximum temperature below the base temperature, then the corresponding values of the daily and monthly CDD was calculated by the following formulas (Equations (1) and (2)): If, (T max + T min )/2 < T base , then CDD = (T max + T base )/4 (1) If, T min > T base , then CDD = (T max + T min )/2 − T base (2) where T max , T min , T base and CDD are the maximum temperature, minimum temperature, base temperature and cooling degree day, respectively.

Overall Heat Transfer Coefficient Calculation
The overall heat transfer coefficient, U, may change because of variations in inflow conditions and fluid properties. For steady-state conditions, the rate of heat flow per unit area through a compound element, such as in the AR, was estimated by the following Equation (3) and heat flow through the BR was calculated by the following Equation (4): where ∑R (m 2 K/W) is the total resistance (the sum of individual resistances), R 0 is the thermal resistance of BR, R s is the thermal resistance of soil with 40% moisture content, R cp is the thermal resistance of cement plaster, R sl is the thermal resistance of slab and R p is the thermal resistance of small plants. R-value of different layers of an AR is given in Table 1.

Cooling Load and Energy Saving Calculation
The total cooling load on a room or building consists of internal loads. The external loads contain heat transfer by conduction through the building walls, roofs, floors, doors, etc., heat transfer by radiation through fenestration such as windows and skylights. The load due to heat transfer through the envelope is named as the external load, while all other loads are called indoor loads. In the case of an internal load of a building, the cooling load is required, especially for internal heat-generating sources such as occupants, lights or appliances. The proportion of external versus internal load varies with building type, site climate and building design. Since the surrounding conditions are highly variable on any given day, the cooling load of an outside-loaded building varies extensively. Apparently, from the energy production and economics points of view, the system design approach for an externally loaded building is a very important issue. Peak load calculations evaluate the utmost load to size and choose the refrigeration equipment. The energy analysis program compares the entire energy use during a certain period with various alternatives so as to work out the optimum one. In this study, Cooling Load Temperature Differential (CLTD) through the roof (URTA roof and BR) was derived and used tabulated data to simplify the calculation process. The basic conduction equation for warmth gain is: where Q (W/m 2 ) is the rate of heat flow per unit area through a compound element and ∆T (K) is the temperature difference. For steady-state conditions, the rate of heat per unit area between each surface is the same. The heat gain is converted to cooling load using the space transfer functions (sol-air temperature) for the rooms with light, medium and heavy thermal characteristics. The equation is modified as [46]: where Q = cooling load (W/m 2 ) is the rate of heat flow per unit area through a compound element; U = Coefficient of heat transfer of roof or wall or glass, W/m 2 K; A = area of the roof in m 2 (% of the area covered by plants, soil, plants with soil and bare are shown in Table 2); and CDD = cooling degree day temperature difference (k). In this study, in the case of overall heat transfer coefficient (U) calculation, the area covered by plants and soil was considered from the survey of existing roofs and experimental roofs. Here, the area covered by soil with plants represent the total area of roof covered by soil with the container or growing medium of plants, and only plants represent the area covered by a leaf of plants that were free from the container or growing medium. Surface temperatures were also measured on the AR under the plants (t 2 ) and over the plant (t 1 ) cover and room below the roof's temperature (t 3 ). On the other hand, at the BR surface, temperatures measured on the roof surface (t 1 ) and the room below the roof's temperature (t 2 ) were used in the calculation of U. Total amount of energy consumption for air conditioning is calculated by the following equation: where Q (W/m 2 ) is the cooling load or heat transmission through the roof (W/m 2 ), and C op is the co-efficient of performance of air conditioning system and is the ratio of useful heating or cooling provided to work required. In this study, C op is calculated by the following Equations [47,48]: where T base is the preferred human comfort temperature as 20 • C, and T mean is the daily average value of room temperature.

Analysis of Yield and Commercial Value of Selected Crops in the URTA
The Yield cash flow quantifies the economic value of the outputs of rooftops' productive use, i.e., of food harvested or energy generated. In the food production scenarios for the selected crops, two seasons (summer and winter) were considered in this analysis. The values of food supply chains (long and short) in the URTA system, where crops were distributed from gardeners to vendors or consumers, were also calculated at the average local market selling price. Planting pots or growing mediums of existing URTA (typically consisting of shallow, free-standing blue plastic drums, wooden boxes, bottles, tins, drums, jutes and plastic bags) were also categorized. Bottle gourd, Tomato and Brinjal were cultivated in two seasons (winter and summer, from November 2018 to May 2019). After completing the winter season, the maximum root zone depth and area were measured, and the container size was selected in the next summer season and winter season based on the first year winter season findings. Similarly, the vegetative area and yield of those crops were also calculated and compared between these two setups.
The weight of each Tomato, Brinjal, Chili and Bottle gourd was measured by a digital weighing machine, and then the total weight of those selected crops was calculated later. For the second year experiment, the Bean was grown in those experiments using the same size of the plastic drum that was used during the first year experiment of Bottle gourd. Several leafy vegetables were considered for this study, including the Spinach, Red Spinach, Water Spinach and Green Spinach, for vertical cultivation as well as for bed cultivation of those vegetables since they were under the fencing panel of Bottle gourd and Bean. The potential crop yields of selected crops in the experimental AR were calculated, and an average sell price in BDT/m 2 from the yearly production was estimated for two seasons to assess the economic benefit and recognize the commercial value of URTA.

Economic Analysis of URTA
Environmental and social benefits of ARs require quantitative appraisals to estimate the financial benefits. However, information on such benefits is not yet commercialized and are mostly individually based on regular conventions. The assessment presented here took into account several accessible sources that justify the economic rewards of URTA. In this study, URTA was implemented with 70% area coverage (on experimental Ars) with the selected crops, as mentioned above. The cost and benefit included two levels: (a) direct effects incurred by the operators of the systems, i.e., investment costs, operation costs and profits generated from yields; and (b) societal effects on the local community, such as market impacts (household savings in food expenses, local jobs creation) and environmental impacts (AR-enhanced air quality, AR habitat creation and mitigation of UHI effect and energy savings) were considered for economic analysis ( Table 3). The potential costs and benefits of a proposed strategy or initiative and, ultimately, its feasibility was the output of an economic evaluation. A systematic process for decision making and trade-offs was also the alternative comparison of economic analysis of URTA. The tests of net present value (NPV), Profitability Index (PI) or Cost-Benefit Ratio (BCR), Internal Rate of Return (IRR) and the payback period, i.e., the period of time required before total revenues equal or surpass total costs for the first time, are the standard methods for economic analysis of any kind of project/firm/scheme. In this study, NPV and IRR were considered for evaluating the feasibility of URTA. NPV defined as the change between the present value (PV) of cash outflows and the PV of cash inflows over a period of time. IRR is a calculation used to estimate the profitability of potential investments or the discount rate at which NPV is zero. The preferred profitable and viable condition of a project is sustained when NPV > 0. In this study, the roof option with the highest NPV and highest IRR indicated the preferred option of URTA. However, Life Cycle Cost (LCC) was considered during the analysis, which was the cost that was associated with the rooftop agricultural firm from the beginning of the project/firm to the end of its useful life, which was considered as 30 years. In this study, NPV is calculated by the following formula [49]: where, F t = net cash inflow-outflows during a single period; r = discount rate or return that be earned in alternative investments; C i = initial investment cost of all setups of URTA; t = number of years within the time periods of first instalment of URTA, generally computed yearly for which the economic evaluation is desired (15 years for this study).
However, if the present value of future cash flows from a likely project using the internal rate as the discount rate, which is subtracted out from the original investment, the net present value would be zero. IRR will be bigger than the discount rate of return (r) for the accepted project. In this research, the discount rate was considered as 12% for calculating NPV based on the Bangladesh government project plan implementation guideline. The payback period means the period of time that a project requires for recovering the money invested in it as well as the life span of the URTA system including all installed material. The payback period is calculated by the following formula: The total irrigation cost has been calculated from the total water requirement (mm) of the selected vegetables and figured the amount of cost based on the tariff of Dhaka Water Supply and Sewerage Authority (DWASA), Bangladesh, from the total amount of water. The cost of yearly energy savings was computed through the multiplication of the simulated energy savings in kWh/m 2 , the total area of the roof, and the energy consumption tariff in Bangladesh. The net cash flows were computed yearly and were assumed to be constant over the investment lifetime. Labour requirements were also considered in the study on the economic evaluation of URTA. Table 3. Sources for net present value (NPV) computation parameters.

Variable Value Source
The installation cost of URTA system including (i) installation of irrigation system and fencing panels, Changing air temperature aspects were observed from different locations of the experimental AR (70% of roof area covered by agriculture with cultivation of Tomato, Chili, Brinjal, Bottle groud, Spinach, Red Spinach, Green Spinach) and from the nearby BR during the whole experimental period. The average temperatures of different times of a single day are shown in Figure 4. However, it has been detected that the temperature reached its maximum range during the month of May at 1.30 p.m., for both AR and BR (Table 4) plots, in comparison to 9.30 p.m. and 5.30 p.m. According to the descriptive statistics, the trend of the air temperature of the BR was always higher than the AR throughout the day.  According to    According to Table 4, the maximum temperatures of the BRs were 2 • C (9:30 a.m.), 5 • C (1:30 p.m.) and 4 • C (5:30 p.m.) higher compared to the AR (experimental roof). Consequently, the minimum temperatures of the BRs were 1 • C (9:30 a.m.), 2 • C (1:30 a.m.) and 2 • C (5:30 p.m.) higher compared to the ARs. The comparison of the mean temperature between the AR and the BR showed that the AR, with 70% of its area covered by plants, was 2.61 • C, 3.41 • C and 2.66 • C cooler than the BR at 9:30 a.m., 1:30 p.m. and 5:30 p.m., respectively. In the case of rainy days, the temperature differences of both roofs became very minimal. So, rain periods were avoided for temperature data analysis.
In this study, percentiles are used to understand the values of thermal dynamics in ARs and BRs, as well as to clearly recognize the advantage of URTA due to reduction in temperatures. From Table 5, it is clearly seen that the different percentile ranges of ARs were always higher than that of BRs. This means that of the 5, 10, 25, 50, 75, 90 and 95 percent temperature values, ARs had a range that was always less than BRs at the same temperature range. Similarly, the air temperature difference histogram, recorded from the five selected roofs (85, 80, 60, 50 and 40% roof area covered by agriculture) and the nearby BRs at 1:30 p.m. and 1.52 m above the roof surface, is shown in Figure 5. The histogram represents the mean value, standard deviation and normal distribution of the temperature difference frequencies of the different area covered roofs during the month of March. From the temperature differences analysis, it was exposed that during the month of March, the maximum frequencies of temperature differences were 5.5 • C, 4.5 • C, 3.5 • C, 2.3 • C, 1.2 • C and 0.45 • C in 85, 80, 60, 50 and 40% roof area covered by agriculture, respectively ( Figure 5).  The mean temperature differences in April and May of those selected ARs and BRs were found to be 4.76 • C, 4.29 • C, 3.37 • C, 2.19 • C, 1.18 • C and 0.41 • C, and 4.41 • C, 3.51 • C, 3.42 • C, 1.85 • C, 1.00 • C and 0 • C in 85, 80, 60, 50 and 40% covered roof area by agriculture, respectively. The minimum and maximum temperature differences in March were recoded as 3 • C and 6 • C, 2.10 • C and 7.20 • C, 2 • C and 6.5 • C, 1 • C and 3 • C, 0.1 • C and 1.8 • C and −1 • C and 1 • C of ARs and BRs, respectively. A 95% confidence interval for the mean was also calculated from the observed temperature differences between ARs and BRs through SPSS. It was revealed that the lower and upper bound of the 95% confidence interval for a mean temperature difference of those ARs and BRs were 4.55 • C and 5.08 • C, 3.90 • C and 5.32 • C, 2.83 • C and 3.91 • C, 1.79 • C and 2.26 • C, 0.85 • C and 1.18 • C and −0.36 • C and 0.27 • C, respectively. So, it is experiential that 50 and 40% covered roof area by URTA obtained a lower temperature reduction, recording a maximum of 1.8 • C and 1 • C and a minimum of 0.1 • C and −0.1 • C at 1:30 p.m. during the month of March to May, respectively. On the contrary, the 85, 80, 70 and 60% URTA were equally and highly effective on air temperature reduction compared to 50 and 40% roof area covered URTA roofs and BRs at the hottest time over the day during the summer season. However, the most noticeable difference is shown by the 85% roof area covered AR, which maintained its temperature variances and standard deviation as 0.34 • C and 0.58 • C compared to the other selected roofs. Thus, the temperature difference should vary on the percentage of area covered by rooftop agriculture persists during the day at 1:30 p.m.

Near Roof Surface Thermal Dynamics in AR and BR
The potential ranges of exterior roof surface temperature reduction were calculated for the experimental AR and BR. The winter and summer season variation of the near roof surface for the AR and the BR is presented as a Box-and-Whisker plot in Figure 6. This boxand-whisker plot shows the lowest value, highest value, median of surface temperature and performance on the roofs. From Figure 6, it is shown that the variability of the roof surface temperature was higher in the BRs for the month of December 2018 to May 2019 (both winter and summer seasons) than ARs. It should be noted that roof surface temperature in both roofs significantly varied at 1:30 p.m. and 5:30 p.m. comparatively with 9:30 a.m. due to the shadow-shading effect or the higher leaf density of plants on the building roof surface. The maximum difference was 12 • C during the month of May, and the average difference was found to be 6 • C on a specific day. It was also observed from one sample t-test that there was no significant temperature difference during the months from December to February at 9:30 a.m. However, from March to May, the temperature significantly differed by 3-12 • C at 9.30 a.m. Similarly, at 1.30 p.m. and 5.30 p.m., the temperature reduction was 4 to 12 • C. It was also reviewed that on semi-intensive green roofs, the roof surface temperature reduction was found to be 7 to 14 • C [53]. The results showed that the URTA roof was very effective in peripheral surface temperature reduction and thereby provided thermal shading to the building. The degree of surface temperature reduction by the URTA increased with the increased solar intensity, as a higher reduction was observed during the daytime at 1:30 p.m. Table 6 represents the temperature gradient at the north-south, middle and east-west side of the AR at 1.25 m above the roof surface from January to April 2019 daily at 9:30 a.m. and 1:30 p.m. The mean temperature incline, measured below the canopy shading at near roof surface 1 m, 2 m, 1.25 m and 2 m above the roof surface, varied within a confined average range of 0.74 • C to 2.32 • C and 1.15 • C to 3.37 • C on the AR during March at 9.30 a.m. and 1.30 p.m., respectively. On the BR, the air temperature incline at different heights varied within the range of 0.50 • C to 1.1 • C and 1 • C to 1.7 • C at 9:30 a.m. and 1:30 p.m., respectively, in March. It was also clearly observed that the temperature gradient was lower below the fencing panel at the north and south sides of the experimental roof compared to other locations, such as the middle, east and west side of the experimental AR where the fencing panel was 1.5 m above the roof surface. Leafy vegetables cultivated under the fencing panel worked as an additional input for lowering the heating effect. The canopy density and height of plants resulted in the temperature change at different heights during the day. It was an upward trend both at 9:30 a.m. and at 1:30 p.m. in the case of the BR. However, for the experimental AR, the temperature variance trend was upward at 9:30 a.m., whereas, at 1:30 p.m., the temperature variance moved downward. Irrigation was given in the experimental AR each morning and in the afternoon by drip irrigation system according to crops' water requirements. So, this impact was limited to affect the temperature changes at a different height. On the other hand, solar intensity was low in the morning compared to noon (1:30 p.m.), and air warming in the morning differred at different heights with little effect. It is notable that the maximum variances occurred on the near-surface at 1.2 m height compared to at 2 m height from the roof surface. Thus, the thermal variation at different heights was more effective in the AR than BRs during the peak solar intensity. It was shown that the highest reduction in the temperature gradient was relative to the density of the canopy and plant distance. This was caused by the difference between the dense greenery and the increased evapotranspiration on the roof. However, the temperature gradient at different sides and heights through URTA created a much cooler microclimate than that of the adjacent BR surface. This could be due to the heat absorbance by plants. Solar radiation and evapotranspiration caused an average temperature gradient of 7.12 • C on the roof surface of 4.58 • C, 2.79 • C and 2.01 • C at 1 m, 1.25 m and 2 m from the roof surface, respectively, compared to the BR. It could be described that there was a relatively cooler air layer on the AR up to 1 m and warmer air above 2 m from the surface around noon. However, at the BR, the temperature suppression was more pronounced at near roof surface and gradually decreased up to 1.25 m, which was also higher than the agricultural roof.

Relative Humidity Dynamics of AR and BRs
In order to assess the effect of rooftop agriculture upon microclimate changes, the average mean, minimum, maximum, standard deviation and variance of the relative humidity (RH) were analyzed from December 2018 to May 2020, collected daily at 9:30 a.m., 1:30 p.m. and 5:30 p.m. The trends of RH in the AR and BR are shown in Figure 7a-c. The maximum, minimum, mean and standard deviation (SD) of RH were observed as 88.83, 30.67, 56.92 and 11.44 for AR and 25.33, 78.67, 49.29 and 10.64% for BR, respectively. Minimum 5% and maximum 10% variations were found between the AR and BR, which indicates that rooftop agriculture is proficient in increasing the RH in the air layer compared to the BR, and it influenced the microclimate of the surrounding air by its evapotranspiration. So, Figure 7a-c demonstrates that the relative humidity changes in the AR were always higher then the BR, which played a significant role in the thermal behaviour of the roof in the daytime at 1:30 p.m. As humidity itself was a climatic variable, it also influenced other climatic variables. Thus, URTA would have positive impacts on the thermal comfort of the people living in urban cities through the reduction of air temperature. So, it is highly recommended to include URTA in the building code of Bangladesh to mitigate the UHI effect. It is applicable for all regions of the globe to reduce global microclimatic change during warmer seasons.

CO 2 Dynamics for the AR and BR
An average 1.63% reduction of CO 2 concentration was observed at 1 m above the roof surface in the experimental AR compared to the adjacent BR within the period from December 2018 to May 2019. The significances of the concentration were analyzed through regression analysis and shown in Figure 8. However, it was observed that the mean concentration of CO 2 (ppm) was 400 ppm and 406 ppm in AR near the plants and in BR, respectively. The maximum concentrations were found to be 431.00 ppm and 440 ppm in AR near the plant and in BR, respectively. Different percentiles of CO 2 concentration were also analyzed. It has been found that the 75, 90 and 95 percentiles of CO 2 concentrations were 404 ppm and 408 ppm, 413 ppm and 410 ppm and 414 ppm and 419 ppm in AR and in BR, respectively. So, the study found that the AR had a higher respiration rate from plants that cause the differences in CO 2 concentration compared to the nearby BRs same as green roof [54]. Hence, AR is able to mitigate the microclimatic changes in urban cities and the UHI effect by reducing heat-trapping gas concentration leading to thermal comfort at a local scale. The regression model was fitted to CO 2 concentration values of both roofs from December to May. A total of 141 days, which were represented by the sequence in the Figure 8, and the deviations from the fitted line to the observed values were noted. The linear line denotes the validity of all values with the dates, which had a negative linear correlation with CO 2 concentration. It indicated that there was no positive relation with date or month, but in the AR, CO 2 concentration was comparatively lower than BR. From February to April, the maximum values were close to the regression line, i.e., the maximum values were found to be close to 401 ppm and 410 ppm for the AR and BR, respectively. So, the URTA plays a very crucial role in the microclimatic changes and controls the temperature and CO 2 rises (1 to 10 ppm) with and around the roofs. This study showed that, for the lower percentage of area coverage by ARs, CDD was nearly same to the adjacent bare roofs and differs by only 0.18 • C. However, for the higher percentage of area coverage by ARs, the CDD difference was found to be 4.25 • C. AR was more suitable for decreasing ambient temperatures and for reducing the cooling load between 32% and 100% [55,56]. From Figure 9, it is found that the possible mean CDD difference was subordinate between the lower percentage of green area coverage and higher percentage of green area coverage of roofs. At least 60% of green area coverage of roofs could be chosen where CDD difference is increasing compare to BR. A paired sample t-test with a 95% confidence interval was used to compare the means of cooling load potential (KW/m 2 ) of selected ARs and nearby BRs with six pairs (Table 7). There was a significant difference in cooling load requirement for the different area coverage of ARs (M = 0.788 to 1.30, SD = 0.0062 to 0.01332) and nearby BRs (M = 1.369 to 1.387, SD = 0.0108 to 0.0280). These results suggested that AR had a substantial cooling effect and depended on the roof area coverage by agriculture. Our research results suggested that when the agricultural roof was covered more than seventy percent, cooling load requirement decreased. The maximum cooling load prerequisite was 1337.74 W/m 2 for the 40% area coverage AR, and the minimum cooling load prerequisite was 772.31 W/m 2 at the 85% area coverage AR compared to other ARs. However, it was detected that among these roofs, indispensable cooling load varied from 12.15 to 20.34%. Therefore, due to the increases in area coverage of URTA, the daily peak cooling load value would be decreased and cooling load saving increased significantly. The buildings with intensive, semi-intensive and extensive green roofs could save about 20-60, 10-45 and 20% energy consumption, respectively [57]. On the other hand, ARs could save 1 to 34% of the amount of total annual energy consumption, 10 to 33.33% of the space cooling load and 20 to 50% of the peak space load [3,58]. According to Tables 8 and 9, it was observed that energy consumption decreased in the high area covered ARs, and different percentile levels of energy savings were observed in all roofs with the increase of green areas. It was clearly observed that 85, 80, 70, 60, 50 and 40% of roof area covered ARs saved energy on top floor of a building by 59.45, 55.63, 39.81, 25.94, 18.88 and 5.87%, respectively. Sixty percent or below 60% area covered ARs saved energy up to 30%. In the case of 60, 50 and 40% of the roof area covered ARs, the 95th percentile of energy saving was as much as one-fourth to one-fifth compared to the 70-85% roof area covered ARs. The results showed that energy consumption differed in all roof options, which was closely related to the area covered by ARs, soil and density of leaves and plants. On the other hand, in the summer season, the energy consumption reduction was barely significant, which is shown in the quantile-quantile (Q-Q) plot normal distribution in Figure 10a-c. The Q-Q plots represent the probability distributions of all values of energy savings by plotting their quantiles against each other and creating a perfectly straight line for the 40, 70 and 85% roof area covered ARs. From Figure 10a-c, it was confirmed that there was a significant potential application of ARs as an energy conservation approach in buildings in hot and moist climatic conditions. However, it can be seen that the energy saving fluctuations of ARs were always less than the thermal fluctuations of BRs, especially in the warm months of the year. However, the equipment operation efficiency was not considered in this study to calculate the energy saving.

Gender-Sensitive Socialization of URTA
Out of the total of 200 rooftop agriculture owners, it could be observed that 51.94% of owners of ARs were female and 48.06% were male, and the ages of the male and female owners were different. However, it was noticed that the URTA was mostly female and elderly male sensitive, where women clearly had an important role to play in increasing the productivity of rooftop agriculture. Therefore, the sensitivity of gender-oriented social dynamics of URTA was multifarious, with individuals expressing degrees of perception towards the four different age factors ( Figure 11). However, personal socialization activities of URTA were identified most strongly among these four groups of age. So, age and sex were the most imperative social factors to put URTA into practice. Figure 11 represents the different year groups (1 = less than 40 years, 2 = 40-50 years, 3 = 50-60 years and 4 = more than 60 years) of males and females and their contribution to the implementation of URTA. Figure 11 characterizes that the 40-50-year-old group of females and the over 60-year-old male group were most perceived by the respondents group of URTA. On the other hand, from the questionnaire survey, it was observed that the maximum everyday jobs of URTA were done mostly by women, including soil preparation, fertilizer application and water management. Some responsibilities were shared with labour such as loading the soil and heavy material transferring such as bamboo, rod sit, containers, soil and organic fertilizer, caring of the roof garden, etc. From the data analysis, the overall skillfulness of women had been increased by 68.78% through rooftop agriculture. So, gender contribution was highly related with URTA and their understanding of agriculture was enhanced through regular involvement in the cultivation of different fruits, flowers, vegetables and other plants in rooftop agriculture. It was found that personal capabilities about the commercialization of URTA products came out as strong factor among the three dynamic parameters. Mandatory in the building code and proper monitoring (36.92%) and subsidies, incentives and bank loans from the government (50%) and training on the agricultural system (13.08%) were the most perceived by respondents. It has also been shown from the previous studies that health (53%) and education (62%), planning social welfare (40%), social group integration (40%), community recreation (35%) and social empowerment (25%) were professed by respondents [59].  Figure 12 represents the results of an economic assessment employing the NPV approach of URTA. NPV was close to zero at the end of the fifth year at a 12% discounted rate. However, NPV became positive, which led to a greater cash inflow compared to cash outflow at a 12% discounted rate at the end of the fifth year within the life period of 15 years of URTA, and at the end of 14th year, NPV was close to zero when the internal rate of return (IRR) is 21.59%. In this study, the 12% discount rate was considered according to Bangaldesh government development project proposal (DPP) appraisal. Due to the very highly sensitive productivity of URTA, the experiment led to a positive NPV after 5 years with proper carrying, including efficient water management techniques both in crops and leafy vegetables. Thus, it can also be concluded that the benefits depends on crop type, production and area covered of the roof by crops and would only be achieved towards the end of the life cycle of the first investment materials of URTA [60]. NPV results of the food production from the URTA scenarios revealed that first-year production was comparatively less than second-year production due to lack of technical knowledge and experience of organic food production on the roof.  However, organic and soil-less cultivation on the roof top led to a positive NPV for its growing capability for around the year and provided a fresh supply of agricultural products to the consumer in a sustainable way. It may contribute to the whole year short duration of food supply chains by as much as 30.07%. The results conluded that 98 ha vegetable gardens and 2539 ha arable land could satisfy the demand of about 63,700 and 321,000 consumers through vegetables and cereal products, respectively [61]. Table 10 represents the yearly benefit of URTA. This study observed that annual job creation advantage @BDT 138.90/m 2 (29.84%) of the total benefits of URTA [ Table 10]. Figure 13 represents the other benefits of the experimental URTA. Yearly energy savings were: BDT 6.04/KWH and BDT 184.44/m 2 /year, considering the fourth step of DPDC tariff, Bangladesh from 301 to 400 units, 6.04/KWh (19.81%). Annual AR-enhanced air quality advantage and annual mitigation of UHI effect are @BDT 2/m 2 (5.16%) and @BDT 67.17/m 2 (14.43%), respectively.

Conclusions
This study aimed at observing the impacts of URTA on microclimate change and socioeconomic dimensions in the urban areas. The output of this study represents the stage for the holistic assessment of alternative solutions, integrating environmental and socioeconomic dimensions and putting URTA into perspective by comparing it to alternative uses of roofs as vacant urban space. The findings of this work reveal that the maximum temperature differences between the ARs and the nearest BRs were 0.45 • C to 5.5 • C during the summer season. It was found that 60-85% roof area covered by URTA were equally and highly effective for air temperature reduction compared to 50% or below roof area covered by URTA (maximum 1 • C to 1.8 • C). The relative humidity was increased by a minimum of 5% and a maximum of 10% in the ARs compared to the BRs. The 70% covered ARs could decrease the Temperature Humidity Index (THI) by minimum of 8%, while the THI for the AR is 26.60 (comfortable), and for the BR, 29 (uncomfortable) in the dry season.
The results of this study also revealed that ARs were effective in reducing heat flow through the roof. Thus, the energy demand for cooling load in the top floor of the building was lowered. The URTA could achieve a saving of 3.62 to 32.28% the peak cooling load. It resulted in 5.87 to 59.45% energy saving with financial benefits compare to the adjacent BR. The increases in area coverage of URTA led to the decrease of the daily peak cooling load. It enhanced energy-saving significantly. The energy-saving fluctuated with ARs with the vegetated area, soil layer coverage and leaf area indices of plants.
However, it has also been found that URTA is mostly female friendly with the age group of 40-50 year. URTA becomes elderly male sensitive with the age group of over 60 years. It indicates that retired males are mostly involved with URTA. Economic sustainability of URTA depends on yields and prices. In this study, at a 12% discount rate, NPV becomes positive at the end of the fifth year, resulting in more cash inflow. URTA is an economically accountable process with financial benefits of yearly energy savings of 19.81%. Annual job creation is 29.84%, enhanced air quality advantage is 5.16% and the annual mitigation of heat island effect is 14.43%. So, the commercial dynamics of URTA refer to achieving financial success according to the demand of the local population. Investments and proper carrying of URTA can raise incomes and produce overall economic growths for longer-term food security and improved well-being. URTA also can provide sustainable, interactive community spaces for flat members or relatives and can enjoy health benefits through recreation and relaxation. URTA brings the unusable space into productive spaces and increases the property value of the building. It plays an important role in addressing a different range of micro-environmental challenges through adaptations.
URTA can play significant roles in producing fresh and affordable vegetables, enhancing cooling load and saving energy, improving urban micro-climatology through reduction of roof temperature and increasing the relative humidity and creating empowerment. The findings of this study may inspire urban planners and decision makers to recognize that URTA can provide measurable benefits both to the city dwellers and to the community to attain environmental and socioeconomic benefits in comparison to traditional urban roof uses.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy or ethical restrictions.