Research Papers

Influence of Green Roofs on Early Morning Mixing Layer Depths in Mexico City

[+] Author and Article Information
Williams Vázquez Morales

Centro de Ciencias de la Atmósfera,
Universidad Nacional Autónoma de México,
CDMX, 04510, México
e-mail: williams.vazquez@unicach.mx

Arón Jazcilevich

Centro de Ciencias de la Atmósfera,
Universidad Nacional Autónoma de México,
CDMX, 04510, México
e-mail: jazcilev@unam.mx

Agustín García Reynoso

Centro de Ciencias de la Atmósfera,
Universidad Nacional Autónoma de México,
CDMX, 04510, México
e-mail: agustin@atmosfera.unam.mx

Ernesto Caetano

Centro de Ciencias de la Atmósfera,
Universidad Nacional Autónoma de México,
CDMX, 04510, México
e-mail: caetano@unam.mx

Gabriela Gómez

Instituto de Geografía,
Universidad Nacional Autónoma de México,
CDMX, 04510, México
e-mail: gabyg@igg.unam.mx

Robert D. Bornstein

Department of Meteorology and Climatology,
San Jose State University,
One Washington Square,
San Jose, CA 95192-0104
e-mail: pblmodel@hotmail.com

1Corresponding author.

Contributed by the Solar Energy Division of ASME for publication in the JOURNAL OF SOLAR ENERGY ENGINEERING: INCLUDING WIND ENERGY AND BUILDING ENERGY CONSERVATION. Manuscript received January 12, 2016; final manuscript received September 15, 2016; published online October 13, 2016. Assoc. Editor: Jorge E. Gonzalez.

J. Sol. Energy Eng 138(6), 061011 (Oct 13, 2016) (10 pages) Paper No: SOL-16-1023; doi: 10.1115/1.4034807 History: Received January 12, 2016; Revised September 15, 2016

An urbanized version of MM5 (uMM5) was used at a 500 m horizontal grid-resolution to study effects on morning urban mixing depths and near roof-top stability from use of extensive green roofs in Mexico City, which is characterized by large Bowen ratios and high building storages. The model uses urban-morphology data, while building hydrothermal uMM5 input parameters were obtained from measurements over green and nearby conventional roofs. Evaluation of uMM5 predicted values against rooftop and planetary boundary layer (PBL) observations from extensive field measurement campaigns showed that the model performed reasonably well. Additional simulations were carried assuming that the roofs in entire urban neighborhoods were greened. Predicted mixing depths from these simulations, along with observed air pollution concentrations, were then used in a simple box model to evaluate potential green roof impacts on concentration. Results showed that green roofs produced an early morning (7–10 LST) cooling of up to 1.2 °C at rooftop levels, which reduced mixing depths during that period. Effects were greater on a day with weak synoptic forcing that on one 48 h later with strong synoptic forcing. The mixing-depth decreases produced increased box-model pollutant concentrations. While the green roofs did not elevate the observed concentrations of CO, SO2, and NO2 above World Health Organization (WHO) health standards, they did increase PM10, values (which were already above its standard) by as much as 8% from 7 to 9 LST, when local populations are normally exposed to peak concentrations. This study has applications in the analyses of building energy efficiency.

Copyright © 2016 by ASME
Your Session has timed out. Please sign back in to continue.


DeNardo, J. C. , Jarrett, A. R. , Manbeck, H. B. , Beattie, D. J. , and Berghage, R. D. , 2005, “ Stormwater Mitigation and Surface Temperature Reduction by Green Roofs,” Trans. ASAE, 48(4), pp. 1491–1496. [CrossRef]
Gartland, L. , 2008, Heat Islands: Understanding and Mitigating Heat in Urban Areas, Earthscan, London.
Castleton, H. F. , Stovin, V. , Beck, S. B. M. , and Davison, J. B. , 2010, “ Green Roofs; Building Energy Savings and the Potential for Retrofit,” Energy Build., 42(10), pp. 1582–1591. [CrossRef]
Sailor, D. J. , Elley, T. B. , and Gibson, M. , 2012, “ Exploring the Building Energy Impacts of Green Roof Design Decisions—A Modeling Study of Buildings in Four Distinct Climate,” J. Build. Phys., 35(4), pp. 372–391. [CrossRef]
Gutierrez, E. , González, J. E. , Arend, M. , Bornstein, R. , and Martilli, A. , 2015, “ On the Anthropogenic Heat Fluxes Using an Air Conditioning Evaporative Cooling Parameterization for Mesoscale Urban Canopy Models,” ASME. J. Sol. Energy Eng., 137(5), p. 051005. [CrossRef]
Redon, E. , Lemonsu, A. , Musy, M. , Masson, V. , Munck, C. D. , and Chancibault, K. , 2015, “ Modeling of Urban Vegetation as a Thermal Regulator and Management of Associated Water Resources for Neighbourhood to City-scale Applications,” 9th International Conference on Urban Climate With 12th Symposium on the Urban Environmental, p. 3.
Taha, H. , 2013, “ Meteorological, Emissions and Air-Quality Modeling of Heat-Island Mitigation: Recent Findings for California, USA,” Int. J. Low-Carbon Technol., 10(1), pp. 3–14. [CrossRef]
Dupont, S. , Otte, T. L. , and Ching, J. K. , 2004, “ Simulation of Meteorological Fields Within and Above Urban and Rural Canopies With a Mesoscale Model,” Boundary-Layer Meteorol., 113(1), pp. 111–158. [CrossRef]
Martilli, A. , Clappier, A. , and Rotach, M. W. , 2002, “ An Urban Surface Exchange Parameterization for Mesoscale Models,” Boundary-Layer Meteorol., 104(2), pp. 261–304. [CrossRef]
Hitchcock, D. , 2005, “ Evaluating Meteorological Impacts of Urban Forest and Albedo Changes in the Houston—Galveston Region: A Fine Resolution (UCP) Meso-Urban Modeling Study of the August-September 2000 Episode,” Report C, Houston Advanced Research Center.
Taha, H. , 2008, “ Meso-Urban Meteorological and Photochemical Modeling of Heat Island Mitigation,” Atmos. Environ., 42(38), pp. 8795–8809. [CrossRef]
Taha, H. , 2008, “ Episodic Performance and Sensitivity of the Urbanized MM5 (uMM5) to Perturbations in Surface Properties in Houston TX,” Boundary Layer Meteorol., 127(2), pp. 193–218. [CrossRef]
Taha, H. , 2008, “ Urban Surface Modification as a Potential Ozone Air-Quality Improvement Strategy in California: A Mesoscale Modeling Study. Boundary-Layer Meteorology,” Boundary Layer Meteorol., 127(2), pp. 219–339. [CrossRef]
Malakooti, H. , Musson-Genon, L. , and Sportisse, B. , 2009, “ Numerical Simulations of Boundary Layer Structure During the High Pollution Episode in Tehran Region,” EGU General Assembly Conference, Vienna, Austria, Apr. 19–24, Vol. 11, p. 7885.
Bass, B. , Krayenhoff, S. , Martilli, A. , and Stull, R. , 2002, “ Mitigating the Urban Heat Island With Green Roof Infrastructure,” Urban Heat Island Summit, Toronto, Canada.
Claverie, R. , Bouyer, J. , Sabre, M. , Schwager, J. , Mucig, C. , Ramier, D. , Berthier, E. , Sellami, E. , Versini, P.-A. , and de Gouvello, B. , 2015, “ The TERRACES Project-A Collaborative Work to Understand the Role of Vegetative Green Roof in Refreshing the Urban Ambiances,” 9th International Conference on Urban Climate With 12th Symposium on the Urban Environmental, p. 27.
Oke, T. R. , Spronken-Smith, R. A. , Jáuregui, E. , and Grimmond, C. S. , 1999, “ The Energy Balance of Central Mexico City During the Dry Season,” Atmos. Environ., 33(24), pp. 3919–3930. [CrossRef]
Matsumoto, Y. , Valdés, M. , Urbano, J. A. , Kobayashi, T. , López, G. , and Peña, R. , 2014, “ Global Solar Irradiation in North Mexico City and Some Comparisons With the South,” Energy Procedia, 57, pp. 1179–1188. [CrossRef]
Molina, L. T. , Madronich, S. , Gaffney, J. S. , Apel, E. , Foy, B. D. , Fast, J. , and Zavala, M. , 2010, “ An Overview of the MILAGRO 2006 Campaign: Mexico City Emissions and Their Transport and Transformation,” Atmos. Chem. Phys., 10(18), pp. 8697–8760. [CrossRef]
Velasco, E. , 2006, “ Measurements at the SIMAT Site,” Molina Center for Energy and the Enviroment(MCE2), personal communication.
Doran, J. C. , Abbott, S. , Archuleta, J. , Bian, S. , Chow, J. C. , Coulter, R. L. , de Wekker, S. F. , Edgerton, S. A. , Elliot, S. , Fernandez, A. , Fast, J. D. , Hubbe, J. M. , King, C. , Langley, D. , Leach, J. , Lee, J. T. , Martin, T. J. , Martinez, S. , Martinez, D. , Martinez, J. L. , Mercado, G. , Mora, V. , Mulhearn, M. , Pena, J. L. , Petty, R. , Porch, W. , Russel, C. , Salas, R. , Shannon, J. D. , Shaw, W. J. , Sosa, G. , Tellier, L. , Templeman, B. , Watson, J. G. , White, R. , Whiteman, C. D. , and Wolfe, D. ,1998, “ The IMADA-AVER Boundary Layer Experiment in the Mexico City Area,” Bull. Am. Meteorol. Soc., 79(11), pp. 2497–2508. [CrossRef]
Dudhia, J. , Gill, D. , Yong, R. G. , and Manning, K. , 2005, “ PSU/NCAR Mesoscale Modeling System Tutorial Class Notes and User's Guide: MM5 Modeling System Version 3,” NCAR, Boulder, CO.
Burian, S. J. , Brown, M. J. , and Augustus, N. , 2007, “ Development and Assessment of the Second Generation National Building Statistics Database,” Seventh Symposium on the Urban Environment, pp. 10–13.
Ching, J. , Brown, M. , McPherson, T. , Burian, S. , Chen, F. , Cionco, R. , and Williams, D. , 2009, “ National Urban Database and Access Portal Tool,” Bull. Am. Meteorol. Soc., 90(8), pp. 1157–1168. [CrossRef]
Oke, T. R. , 1987, Boundary Layer Climates, 2nd ed., Vol. 452, Routledge, London.
National Centers for Environmental Prediction/National Weather Service/NOAA/U.S. Department of Commerce, 2000, “ NCEP FNL Operational Model Global Tropospheric Analyses, Continuing From July 1999,” Research Data Archive at the National Center for Atmospheric Research, Computational and Information Systems Laboratory, Boulder, CO, accessed Nov. 19, 2009, http://rda.ucar.edu/datasets/ds083.2
Shaw, W. J. , Pekour, M. S. , Coulter, R. L. , Martin, T. J. , and Walters, J. T. , 2007, “ The Daytime Mixing Layer Observed by Radiosonde, Profiler, and Lidar During MILAGRO,” Atmos. Chem. Phys. Discuss., 7(5), pp. 15025–15065. [CrossRef]
De Haij, M. , Wauben, W. , and Baltink, H. K. , 2006, “ Determination of Mixing Layer Height From Ceilometer Backscatter Profiles,” Proc. SPIE, 6362, p. 63620R.
Seinfeld, J. H. , and Spyros, N. P. , 1998, Atmospheric Chemistry and Physics: From Air Pollution to Climate Change, Wiley, New York, pp. 1197–1200.
Sportisse, B. , 2001, “ Box Models Versus Eulerian Models in Air Pollution Modeling,” Atmos. Environ., 35(1), pp. 173–178. [CrossRef]
Willmott, C. J. , 1981, “ On the Validation of Models,” Phys. Geogr., 2(2), pp. 184–194.
SMA-GDF., 2006, “ Informe Climatológico Ambiental del Valle de México 2006,” SMA-GDF, México D.F., accessed Feb. 10, 2009 http://www.sma.df.gob.mx
Bornstein, R. D. , 1968, “ Observations of the Urban Heat Island Effect in New York City,” J. Appl. Meteorol., 7(4), pp. 575–582. [CrossRef]


Grahic Jump Location
Fig. 1

Urban land covering within domain uMM5 innermost D4, including SIMAT site (5, star) and urban garden areas: (1) Bosque de Chapultepec, (2) Parque España, (3) Parque México, and (4) Cementerio Francés; modeling results will be shown along vertical cut L–L′

Grahic Jump Location
Fig. 2

(a) Satellite photograph showing measurement locations on green and conventional roofs in botanical garden area on UNAM campus and (b) green roof with succulent cactacea plants

Grahic Jump Location
Fig. 3

Central Mexico MM5 modeling domains D1–D3 and uMM5 domain D4

Grahic Jump Location
Fig. 4

Box model schematic, where all symbols are defined in text

Grahic Jump Location
Fig. 5

Surface anticyclonic synoptic conditions over central Mexico (black circle) on Mar. 10, 2006, showing horizontal winds (1 barb = 10 m s−1) and omega vertical velocities (scale bar at right; downward motion is positive, Pa·s−1) at (a) 03, (b) 06, and (c) 09 LST

Grahic Jump Location
Fig. 6

Measured Mar. 22, 2009 (a) virtual potential temperatures and (b) RHs, with solid line for green and dashed line for conventional roof simulations, same for (c) and (d), but for March 20.

Grahic Jump Location
Fig. 7

Comparison of Mar. 8–16, 2006 average hourly (a) wind speed and (b) temperature for measurements at SIMAT station (circles), uMM5 (solid line), and MM5 (dashed line)

Grahic Jump Location
Fig. 8

Modeled and measured (SIMAT site) fluxes Mar. 10, 2006 for (a) sensible and (b) latent heats, both for conventional roof case (CR dots), green roof case (GR triangles), and measurements (dashes); also shown are modeled CR (c) ground and (d) storage heat fluxes

Grahic Jump Location
Fig. 9

Ceilometer backscatter and uMM5 PBL heights (km, solid line and then dots) for (a) March 12 and (b) March 14

Grahic Jump Location
Fig. 10

Comparison of 06 to 11 LST uMM5 PBL heights for conventional (CR, dashes) and green roof (GR solid) cases on (a) Mar. 12, 2006 and (b) Mar. 14, 2006; enlarged view in (a) shows period with most significative differences

Grahic Jump Location
Fig. 11

Simulated uMM5 potential temperature (°C) and PBL height (horizontal lines) for conventional (CR) and green (GR) roof cases for (a) March 12 at 0910 LST and (b) March 14 at 0840 LST; (c) and (d) show corresponding modeled TKE values (J/kg)

Grahic Jump Location
Fig. 12

Percent reduction of maximum green roof uMM5 PBL heights for March 7–16, with time (LST) of maximum reduction indicated

Grahic Jump Location
Fig. 13

Vertical cross section (along line L–L′ of Fig. 1) of modeled mixing heights for conventional (CR) and green (GR) roof cases for (a) March 12 and (b) March 14; location of site SIMAT also shown

Grahic Jump Location
Fig. 14

Trends of CO, SO2, NO2, and PM10 concentrations at 07–10 LST for green (solid line) and observed/conventional (dashed line) roof cases for (a) March 12 and (b) March 14; horizontal dotted lines are WHO health norms




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In