case study method in geography

• DOI: 10.1080/03098269708725423
• Corpus ID: 129098785

A Claim for the Case Method in the Teaching of Geography.

• Published 1997
• Geography, Education
• Journal of Geography in Higher Education

88 Citations

The comparison of a thematic versus regional approach to teaching a world geography course, analysis of cultural and political geography learning concept design based on case method and team-based project, using the internet to integrate thematic and regional approaches in geographic education, living cases: authentic learning in action, making the case for case study learning, the use of case study teaching to promote autonomous learning, yawning sixth formers: an action research project examining how we can move beyond passive learning in sixth form teaching of case studies in urban management.

• Highly Influenced

World Geography: Organizational and Teaching Strategies

What a difference a case makes: teaching economic evaluation of public policies to non-economists, developing, implementing and evaluating case studies in materials science, 20 references, interactive lectures: a case study in a geographical concepts course, a scheme for the effective use of role plays for an emancipatory geography, using inquiry to enhance the learning and appreciation of geography., teaching geography in higher education: a manual of good practice, geography 16–19: some implications for higher education, beyond the lecture: case teaching and the learning of economic theory, the case method as a strategy for teaching policy analysis to undergraduates., in sight of the tunnel: the renaissance of geography education, situated cognition and the culture of learning, geography and education: north america, related papers.

Showing 1 through 3 of 0 Related Papers

Key Methods in Geography

Student resources, welcome to the companion website.

Welcome to the companion website for Research Methods in Geography, Third edition,  by Nicholas Clifford, Meghan Cope, Thomas Gillespie and Shaun French . The resources on the site have been specifically designed to support your study.

For students

• SAGE journal articles
• Further reading
• Supplementary material

About the book:

Key Methods in Geography  is the perfect introductory companion, providing an overview of qualitative and quantitative methods for human and physical geography. The third edition of this essential and accessible primer features:

• 12 new chapters  representing emerging themes including online, virtual and digital geographical methods
• Real-life case study  examples
• Summaries and exercises  for each chapter
• Free online access  to full text of  Progress in Human Geography  and  Progress in Physical Geography  Progress Reports

The teaching of research methods is integral to all geography courses:  Key Methods in Geography, 3rd edition  explains all of the key methods with which geography undergraduates must be conversant.

Acknowledgement: Thank you to our Online Content Editor: James D.A. Millington.

Disclaimer:

This website may contain links to both internal and external websites. All links included were active at the time the website was launched. SAGE does not operate these external websites and does not necessarily endorse the views expressed within them. SAGE cannot take responsibility for the changing content or nature of linked sites, as these sites are outside of our control and subject to change without our knowledge. If you do find an inactive link to an external website, please try to locate that website by using a search engine. SAGE will endeavour to update inactive or broken links when possible. 

Methodological Approaches in Physical Geography

• © 2022
• Firuza Begham Mustafa   ORCID: https://orcid.org/0000-0002-8152-6623 0

Department of Geography, University Malaya, Kuala Lumpur, Malaysia

You can also search for this editor in PubMed   Google Scholar

• consists of a diverse and detailed research methodology focus on physical geography
• with contributions from prominent geographers, professors and other subject experts from South Asia and East Asia
• Covers a wide range of physical geography methods

Part of the book series: Geography of the Physical Environment (GEOPHY)

4350 Accesses

4 Citations

1 Altmetric

This is a preview of subscription content, log in via an institution to check access.

Access this book

Subscribe and save.

• Get 10 units per month
• Download Article/Chapter or eBook
• 1 Unit = 1 Article or 1 Chapter
• Cancel anytime
• Available as EPUB and PDF
• Read on any device
• Instant download
• Own it forever
• Compact, lightweight edition
• Dispatched in 3 to 5 business days
• Free shipping worldwide - see info
• Durable hardcover edition

Tax calculation will be finalised at checkout

Other ways to access

Licence this eBook for your library

Institutional subscriptions

About this book

Geography science aims to observe the dynamics in describing earth's surface as a place and space for humans to carry out their lives, starting from simple identification using recording and sketching models, then utilizing tools such as maps, satellite imagery, statistics and Geographic Information Systems (GIS). In the development of geography science, it is appropriate to explain phenomena of the earth in the present context along with the process of developing science and technology using suitable and effective methods. Physical geography is the branch of natural science that deals with the study of processes and patterns in the natural environment such as the atmosphere, hydrosphere, lithosphere and biosphere. This book covers the methodology of the study for all aspects of physical geography, biosphere, hydrosphere, lithosphere, and atmosphere. A comprehensive geography textbook consists of a detailed research methodology for physical geography research including a few selected case studies in Asia. The uniqueness of this book is due to the contribution of several professors and subject experts from South East and East Asia with special particular reference to cases studies from a particular region. This book covered selected methodological approaches for hydrology, climatology and geomorphology including the discovery of the best method for exploring and assessing mysterious physical phenomena using a diversity of methodologies. This book explains the principal concept, basic method, optional method, detailed description of each method, and the advantages and disadvantages of the various methods. The technique of data selection, data acquisition, method of analysis, data interpretation and data analysis techniques with a specific focus on deterministic modeling, geography techniques, geospatial modeling with Geographic Information System (GIS), Artificial Intelligence (AI), Analytic Hierarchy Process (AHP), and Automated machine techniques and combination of statistical analysis. This book attempts to explore different approaches, methodological possibilities and challenges in conducting geographical research in physical geography. New digital geographic data sources and GIS applications can help researchers to receive clearer concepts and obtain better measurements of the relevant attributes changes in the physical environment. Opportunity to critically examine the conceptualization and identification of the field in geographical research and how digital media has not only expanded the scope of what constitutes the field but has redefined the field in itself as well as the practices of observing, knowing, and analyzing the real world.

Similar content being viewed by others

Portray of Methodologies in Geography: An Inquiry

An Introduction to Anthropogeomorphology and Geospatial Technology

Geospatial Information Research: State of the Art, Case Studies and Future Perspectives

• methods in physical geography
• digital geography
• methods in geomorphology
• flood susceptibility mapping
• monitoring in geography
• spatial analysis
• site selection methods

Table of contents (10 chapters)

Front matter, epistemology of geography.

• Dedi Hermon

Site Selection Method Using the Geographic Information System (GIS) and Analytic Hierarchy Process (AHP)

• Benjamin Bwadi Ezekiel, Firuza Begham Mustafa, Gabriel Temitope Adelalu, Bakoji Mohammed Yusuf

The Application of a Data-Driven Method for Spatial Analysis and Prediction of Gully Erosion Susceptibility

• Didams Gideon, Firuza Begham Mustafa

Methods and Approaches of Flood Susceptibility Assessment and Mapping: A Review in Geographical Perspective

• Khadija Bibi, Fareeha Siddique, Shehla Gul, Atta-ur Rahman, Firuza Begham Mustafa

Digital Geography and Its Methods

• Aparajita De

Automated in Situ Water Quality Monitoring—Characterizing System Dynamics in Urban-Impacted and Natural Environments

• Kim N. Irvine, Lloyd H. C. Chua, Cameron A. Irvine

Research Methods and Techniques in Physical Geography

• Virendra Nagarale, Subhash Anand, Piyush Telang

The Methodological Approaches in Physical Geography

• K. W. G. Rekha Nianthi, A. K. Wickramasooriya, Lalitha Dissanayake, C. S. Hettiarachchi, R. M. G. N. Rajapaksha

Urban River Restoration: A Methodological Discourse with Examples from Kerala, India

• Srikumar Chattopadhyay

The Methodological Approach of Assessing Urban Vertical Expansion Using Satellite Remote Sensing Techniques

• L. Manawadu, V. P. I. S. Wijeratne

Editors and Affiliations

Firuza Begham Mustafa

About the editor

Bibliographic information.

Book Title : Methodological Approaches in Physical Geography

Editors : Firuza Begham Mustafa

Series Title : Geography of the Physical Environment

DOI : https://doi.org/10.1007/978-3-031-07113-3

Publisher : Springer Cham

eBook Packages : Earth and Environmental Science , Earth and Environmental Science (R0)

Copyright Information : The Editor(s) (if applicable) and The Author(s), under exclusive license to Springer Nature Switzerland AG 2022

Hardcover ISBN : 978-3-031-07112-6 Published: 02 September 2022

Softcover ISBN : 978-3-031-07115-7 Published: 03 September 2023

eBook ISBN : 978-3-031-07113-3 Published: 01 September 2022

Series ISSN : 2366-8865

Series E-ISSN : 2366-8873

Edition Number : 1

Number of Pages : XI, 178

Number of Illustrations : 14 b/w illustrations, 57 illustrations in colour

Topics : Physical Geography , Geography, general , Environment, general

• Publish with us

Policies and ethics

• Find a journal
• Track your research

• Human Geography
• Physical Geography
• Economic Geography
• Urban Geography
• Political Geography
• The United States and Canada
• Latin America
• Russia and the Republics
• Australia, New Zealand, and Antarctica
• Arctic Ocean
• Atlantic Ocean
• Indian Ocean
• Pacific Ocean
• Climate Classification

Issues and Challenges of the Case Study Approach

Use of the case study approach in both teaching and research has been criticized for more than five decades. Indeed, unlike this reference book, many research and pedagogy guidebooks fail to include the case study approach in their discussions. The reasons for this are twofold. First, case study research has been criticized for: (1) its unscientific nature (because findings cannot be replicated) and (2) reliance on overgeneralizable findings. Key to overcoming this first limitation is triangulating a rigorous set of mixed method approaches to data collection and analysis and maintaining a chain of evidence to argue a case. The second criticism is best mitigated by using the findings from the case study to address and contribute to larger questions, issues, and theories in  human geography .

Since almost all case studies involve the use of interview methodologies or ethnographic work, one of the strongest arguments for their validity is to emphasize that larger scale data sets often overlook or blur the significance of individual stories. Therefore, case study research has the potential to capture and analyze the lived experiences of people, and understand more about particular places on the ground. In sum, despite criticism related to studies that focus on specific places, groups, or issues, scaling up the findings from small scale projects to respond to larger research questions makes case study research and teaching critical in helping link local issues to larger global challenges.

• Using the Case Study Approach to Teach Human Geography
• Methods Useful in Case Study Research
• Types of Case Studies
• Case Study Approach
• Capital’s Consumption Spatiality
• Capital and Space: Capital’s Crisis-Spatiality
• Capital and Space: Capital’s ‘Normal’ Spatiality
• How Does Capital Work?: Mechanisms of Capital
• Capital and Space

• Find My Rep

You are here

Key Methods in Geography

• Nicholas Clifford - Loughborough University, UK
• Meghan Cope - University of Vermont, USA
• Thomas Gillespie - University of California, Los Angeles, USA
• Description

Key Methods in Geography is the perfect introductory companion, providing an overview of qualitative and quantitative methods for human and physical geography. The fourth edition of this essential and accessible primer covers the breadth of the discipline and offer critical and contextual perspectives on research methods. New coverage takes account of newer technologies and practice, and 9 new chapters bring greater diversity of positionality and perspective to the volume, including decolonial methods, predicting, visualizing and modelling climate and environmental change, and writing up research. Case study examples, summaries and exercises have been included in each chapter to enable learning.

This is vital reading for any student undertaking a Geography Methods module as well as a valuable resource for any student embarking on independent research as part of their degree.

An in-depth inquisition of the complex dynamic landscape associated with human and physical geography and their interactions, expertly articulated and packaged by an ecosystem of experts, key methods in geography is an ideal reference book for novice and experienced researchers.  Faith Njoki Karanja - Geoinformation Expert

This usefully revised book should be essential reading for any of our undergraduates (and even postgraduates), particularly those undergrads in their final year that are undertaking independent diss/project work. As I have done with previous editions, I will certainly recommend it to all my supervisees and will suggest to colleagues that they do the same. Many colleagues probably already do so, either through the tutorial programme, the dissertation/project modules or fieldwork modules. [One thought: in future editions, there will perhaps be a need for a dedicated chapter on use of AI in research. I notice that AI appears in the index of this volume but understandably doesn't cover the latest developments in the field (e.g. release of ChatGPT and the like).]

It is important for students to have a disciplinary book that confirm all of them how to make their research dissertation.

Preview this book

For instructors.

Please select a format:

Select a Purchasing Option

• Electronic Order Options VitalSource Amazon Kindle Google Play eBooks.com Kobo

How to revise geography case studies

Molly, one of my readers, wrote to me to ask:

I was just wondering if you had any tips for revising and remembering geography case studies?

When I've asked around some of the things that students find most difficult about revising geography case studies are:

• Condensing all the information
• Remembering all the statistics
• Knowing what you need to know, and what you can afford to forget

Having been a bit of an ace at this kind of thing myself (I got an A* at GCSE, A at A-Level and a degree in the subject) I thought I'd share some of my top tips on how to revise geography case studies today.

1. Make sure you understand the case study

The first step in remembering anything is understanding it. You need to have a clear model in your mind of how the case study works. This includes how it's laid out in space (a mental map), who the people were who were involved and the context of the case study (historical, political, social, economic and environmental. These tips will help you with this:

• Make sure you've seen a map of the place. In this day and age this is easy with google maps, google earth and google streetview. All of these things can help you understand both the 2-D and 3-D landscape of the case study.
• Find newspaper articles and pictures to give you some background and also help you to visual the place
• Watch videos if they exist. For some case studies there are amazing clips of films (Kibera, the Nairobi shanty town at the beginning of The Constant Gardner springs to mind). For others there will be great video clips on YouTube to help you.
• If you can, visit the place. Nothing is as powerful as this in fully understanding a place.

2. Condense your notes

Once you've thoroughly understood the case study it's time to condense your notes. There are various ways you can do this.

• Create an A3 annotated map of the area. Colour code things like causes and effects or social, economic, environmental and political factors. Have a key. You can even have flaps. Stick the map up on the wall and look at it frequently. The great thing about this is that the finite size of the page forces you to condense the information.
• Create a table. You could put things like the social, economic, political and environmental factors along one side and background, causes and effects along the other.  Inside The Extraordinaries Club I have some grids and guidelines for you to download and use. These are exclusively for members. Find out more about the club here .
• Create index cards. This was one of my favourites as it was quite a kinaesthetic way of separating the information into bite-sized chunks. It was a great format to give other people to test me so that I could learn all those facts and figures. You can even have different coloured index cards for different topics.
• Traditional revision notes . In my opinion this is a bit boring, and can also be a bit intimidating when you see reams of notes that you have to memorise. I'd go with one of the other options if I was you.

Download a free extract of my book, The Ten Step Guide to Acing Every Exam You Ever Take , to find out how you learn best. Enter your details in the box below to get your extract.

DOWNLOAD YOUR FREE CHAPTER

Of the ten step guide to acing every exam you ever take.

Identify how YOU learn best and make your independent study as successful as it possibly can be

Your personal data will be treated with respect. View my privacy policy here .

3. Memorise

Now you've condensed your class notes you need to memorise them. Good memorisation, in my experience comes down to two things:

• Using the information in different formats.

I'd advise you to do a combination of the following:

• Read index cards out loud, cover and test yourself.
• Get other people to test you.
• Act it out.
• Make up songs or rhymes
• Whatever else works for you…

4. Teach someone else about the case study

Teaching someone else is one of the best ways there is to a) check your understanding (because they'll never understand it if you don't) and b) practice putting what you know into words so that someone else understands it.

5. Do Past Papers

The final step is to do past papers. I strongly recommend that you do this in the format of Revision Power Hours.

If you do power hours, and make a point of marking your work, you'll not only do lots of repetition of the case study you've been learning, you'll start to learn to think like an examiner and also get a brilliant insight into exactly what they expect you to know in terms of facts and statistics.

I will say this. I used to remember literally hundreds of stats for my case studies. When I became a teacher it surprised me how few students actually needed to know in order to get good marks. However, this comes with a word of warning. You need to have a good insight into what your exam board expects you to know.

Over to you

That's pretty much a masterclass in how to revise geography case studies. Now it's up to you to put it into practice

In the comments below I'd love to know:

• What you're finding difficult about revising geography case studies
• What other subjects you need revision help with

Need more help with your revision?

It works because you develop a personalised study system tailored to the way you work, what motivates you and your ambitions.

Exam success with way less stress.

You might also like...

Why your revision isn’t working

Why you need to make past paper practise part of your revision

3 Exam Technique Mistakes GCSE and A-Level Students Commonly Make

How to revise GCSE English with 7 weeks to go

“Is it too late for revision to make a difference for this year’s GCSEs and A-Levels?”

Why finding your study style is like finding your fashion sense

Lucy Parsons

Click here to leave a comment below, leave a reply:.

Save my name, email, and website in this browser for the next time I comment.

Geographic Information Systems/Science

• The Basics of GIS
• Getting started with GIS
• Selected book titles for GIS
• ArcGIS software
• Open source software
• GIS manuals
• Computer programming with GIS
• Articles & journals
• ArcGIS for Office
• GIS for the Mac
• Case studies in GIS
• Archaeology & GIS
• Census & GIS
• History & GIS
• Urban applications
• Spatial analysis & modelling
• Public health
• Internet resources for GIS
• GIS Education Blog
• ArcGIS Online
• Web mapping This link opens in a new window
• Digital mapping
• Google Earth & Maps
• Geovisualization
• Managing citations
• Scholarly communication This link opens in a new window

Subject Librarian

Introduction

    Case studies are specific examples of how GIS was used to solve a problem or made information sharing easier in a particular industry. The Library has a collection of books which illustrate a number of case studies. These case studies can also provide you with inspiration for your own GIS projects. This is not the complete list of case studies available, but those listed should help you with ideas.

Looking for statistics or data to actually map? Look here .

Finding library resources

• GIS and case studies This is a keyword search in the online catalog.
• geographic information systems case studies
• geographic information systems

Selected case studies

In addition to the books listed below, ESRI has published a series of booklets online to show best practices in different fields using ArcGIS.

Internet resource(s)

Looking for GIS case studies on the Internet? Just do a search for "gis case studies." That will generate a long list. Then you need to decide what you want from the case study. That will help you narrow your search.

• Case studies from ESRI GIS Education Community These case studies are rich, real-world stories of GIS in education – across several educational categories. These stories represent best practices of GIS planning, implementation or evaluation in education.
• Case studies from Lake County, Florida

• GIS Case Studies from the County of Los Angeles, CA
• << Previous: GIS for the Mac
• Next: Archaeology & GIS >>
• Last Updated: Mar 22, 2024 2:00 PM
• URL: https://researchguides.dartmouth.edu/gis

An official website of the United States government

The .gov means it’s official. Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

The site is secure. The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

• Publications
• Account settings

Preview improvements coming to the PMC website in October 2024. Learn More or Try it out now .

• Advanced Search
• Journal List
• Sage Choice
• PMC11334375

Methodologic and Data-Analysis Triangulation in Case Studies: A Scoping Review

Margarithe charlotte schlunegger.

1 Department of Health Professions, Applied Research & Development in Nursing, Bern University of Applied Sciences, Bern, Switzerland

2 Faculty of Health, School of Nursing Science, Witten/Herdecke University, Witten, Germany

Maya Zumstein-Shaha

Rebecca palm.

3 Department of Health Care Research, Carl von Ossietzky University Oldenburg, Oldenburg, Germany

Associated Data

Supplemental material, sj-docx-1-wjn-10.1177_01939459241263011 for Methodologic and Data-Analysis Triangulation in Case Studies: A Scoping Review by Margarithe Charlotte Schlunegger, Maya Zumstein-Shaha and Rebecca Palm in Western Journal of Nursing Research

We sought to explore the processes of methodologic and data-analysis triangulation in case studies using the example of research on nurse practitioners in primary health care.

Design and methods:

We conducted a scoping review within Arksey and O’Malley’s methodological framework, considering studies that defined a case study design and used 2 or more data sources, published in English or German before August 2023.

Data sources:

The databases searched were MEDLINE and CINAHL, supplemented with hand searching of relevant nursing journals. We also examined the reference list of all the included studies.

In total, 63 reports were assessed for eligibility. Ultimately, we included 8 articles. Five studies described within-method triangulation, whereas 3 provided information on between/across-method triangulation. No study reported within-method triangulation of 2 or more quantitative data-collection procedures. The data-collection procedures were interviews, observation, documentation/documents, service records, and questionnaires/assessments. The data-analysis triangulation involved various qualitative and quantitative methods of analysis. Details about comparing or contrasting results from different qualitative and mixed-methods data were lacking.

Conclusions:

Various processes for methodologic and data-analysis triangulation are described in this scoping review but lack detail, thus hampering standardization in case study research, potentially affecting research traceability. Triangulation is complicated by terminological confusion. To advance case study research in nursing, authors should reflect critically on the processes of triangulation and employ existing tools, like a protocol or mixed-methods matrix, for transparent reporting. The only existing reporting guideline should be complemented with directions on methodologic and data-analysis triangulation.

Case study research is defined as “an empirical method that investigates a contemporary phenomenon (the ‘case’) in depth and within its real-world context, especially when the boundaries between phenomenon and context may not be clearly evident. A case study relies on multiple sources of evidence, with data needing to converge in a triangulating fashion.” 1 (p15) This design is described as a stand-alone research approach equivalent to grounded theory and can entail single and multiple cases. 1 , 2 However, case study research should not be confused with single clinical case reports. “Case reports are familiar ways of sharing events of intervening with single patients with previously unreported features.” 3 (p107) As a methodology, case study research encompasses substantially more complexity than a typical clinical case report. 1 , 3

A particular characteristic of case study research is the use of various data sources, such as quantitative data originating from questionnaires as well as qualitative data emerging from interviews, observations, or documents. Therefore, a case study always draws on multiple sources of evidence, and the data must converge in a triangulating manner. 1 When using multiple data sources, a case or cases can be examined more convincingly and accurately, compensating for the weaknesses of the respective data sources. 1 Another characteristic is the interaction of various perspectives. This involves comparing or contrasting perspectives of people with different points of view, eg, patients, staff, or leaders. 4 Through triangulation, case studies contribute to the completeness of the research on complex topics, such as role implementation in clinical practice. 1 , 5 Triangulation involves a combination of researchers from various disciplines, of theories, of methods, and/or of data sources. By creating connections between these sources (ie, investigator, theories, methods, data sources, and/or data analysis), a new understanding of the phenomenon under study can be obtained. 6 , 7

This scoping review focuses on methodologic and data-analysis triangulation because concrete procedures are missing, eg, in reporting guidelines. Methodologic triangulation has been called methods, mixed methods, or multimethods. 6 It can encompass within-method triangulation and between/across-method triangulation. 7 “Researchers using within-method triangulation use at least 2 data-collection procedures from the same design approach.” 6 (p254) Within-method triangulation is either qualitative or quantitative but not both. Therefore, within-method triangulation can also be considered data source triangulation. 8 In contrast, “researchers using between/across-method triangulation employ both qualitative and quantitative data-collection methods in the same study.” 6 (p254) Hence, methodologic approaches are combined as well as various data sources. For this scoping review, the term “methodologic triangulation” is maintained to denote between/across-method triangulation. “Data-analysis triangulation is the combination of 2 or more methods of analyzing data.” 6 (p254)

Although much has been published on case studies, there is little consensus on the quality of the various data sources, the most appropriate methods, or the procedures for conducting methodologic and data-analysis triangulation. 5 According to the EQUATOR (Enhancing the QUAlity and Transparency Of health Research) clearinghouse for reporting guidelines, one standard exists for organizational case studies. 9 Organizational case studies provide insights into organizational change in health care services. 9 Rodgers et al 9 pointed out that, although high-quality studies are being funded and published, they are sometimes poorly articulated and methodologically inadequate. In the reporting checklist by Rodgers et al, 9 a description of the data collection is included, but reporting directions on methodologic and data-analysis triangulation are missing. Therefore, the purpose of this study was to examine the process of methodologic and data-analysis triangulation in case studies. Accordingly, we conducted a scoping review to elicit descriptions of and directions for triangulation methods and analysis, drawing on case studies of nurse practitioners (NPs) in primary health care as an example. Case studies are recommended to evaluate the implementation of new roles in (primary) health care, such as that of NPs. 1 , 5 Case studies on new role implementation can generate a unique and in-depth understanding of specific roles (individual), teams (smaller groups), family practices or similar institutions (organization), and social and political processes in health care systems. 1 , 10 The integration of NPs into health care systems is at different stages of progress around the world. 11 Therefore, studies are needed to evaluate this process.

The methodological framework by Arksey and O’Malley 12 guided this scoping review. We examined the current scientific literature on the use of methodologic and data-analysis triangulation in case studies on NPs in primary health care. The review process included the following stages: (1) establishing the research question; (2) identifying relevant studies; (3) selecting the studies for inclusion; (4) charting the data; (5) collating, summarizing, and reporting the results; and (6) consulting experts in the field. 12 Stage 6 was not performed due to a lack of financial resources. The reporting of the review followed the PRISMA-ScR (Preferred Reporting Items for Systematic Reviews and Meta-Analyses extension for Scoping Review) guideline by Tricco et al 13 (guidelines for reporting systematic reviews and meta-analyses [ Supplementary Table A ]). Scoping reviews are not eligible for registration in PROSPERO.

Stage 1: Establishing the Research Question

The aim of this scoping review was to examine the process of triangulating methods and analysis in case studies on NPs in primary health care to improve the reporting. We sought to answer the following question: How have methodologic and data-analysis triangulation been conducted in case studies on NPs in primary health care? To answer the research question, we examined the following elements of the selected studies: the research question, the study design, the case definition, the selected data sources, and the methodologic and data-analysis triangulation.

Stage 2: Identifying Relevant Studies

A systematic database search was performed in the MEDLINE (via PubMed) and CINAHL (via EBSCO) databases between July and September 2020 to identify relevant articles. The following terms were used as keyword search strategies: (“Advanced Practice Nursing” OR “nurse practitioners”) AND (“primary health care” OR “Primary Care Nursing”) AND (“case study” OR “case studies”). Searches were limited to English- and German-language articles. Hand searches were conducted in the journals Nursing Inquiry , BMJ Open , and BioMed Central ( BMC ). We also screened the reference lists of the studies included. The database search was updated in August 2023. The complete search strategy for all the databases is presented in Supplementary Table B .

Stage 3: Selecting the Studies

Inclusion and exclusion criteria.

We used the inclusion and exclusion criteria reported in Table 1 . We included studies of NPs who had at least a master’s degree in nursing according to the definition of the International Council of Nurses. 14 This scoping review considered studies that were conducted in primary health care practices in rural, urban, and suburban regions. We excluded reviews and study protocols in which no data collection had occurred. Articles were included without limitations on the time period or country of origin.

Inclusion and Exclusion Criteria.

Screening process

After the search, we collated and uploaded all the identified records into EndNote v.X8 (Clarivate Analytics, Philadelphia, Pennsylvania) and removed any duplicates. Two independent reviewers (MCS and SA) screened the titles and abstracts for assessment in line with the inclusion criteria. They retrieved and assessed the full texts of the selected studies while applying the inclusion criteria. Any disagreements about the eligibility of studies were resolved by discussion or, if no consensus could be reached, by involving experienced researchers (MZ-S and RP).

Stages 4 and 5: Charting the Data and Collating, Summarizing, and Reporting the Results

The first reviewer (MCS) extracted data from the selected publications. For this purpose, an extraction tool developed by the authors was used. This tool comprised the following criteria: author(s), year of publication, country, research question, design, case definition, data sources, and methodologic and data-analysis triangulation. First, we extracted and summarized information about the case study design. Second, we narratively summarized the way in which the data and methodological triangulation were described. Finally, we summarized the information on within-case or cross-case analysis. This process was performed using Microsoft Excel. One reviewer (MCS) extracted data, whereas another reviewer (SA) cross-checked the data extraction, making suggestions for additions or edits. Any disagreements between the reviewers were resolved through discussion.

A total of 149 records were identified in 2 databases. We removed 20 duplicates and screened 129 reports by title and abstract. A total of 46 reports were assessed for eligibility. Through hand searches, we identified 117 additional records. Of these, we excluded 98 reports after title and abstract screening. A total of 17 reports were assessed for eligibility. From the 2 databases and the hand search, 63 reports were assessed for eligibility. Ultimately, we included 8 articles for data extraction. No further articles were included after the reference list screening of the included studies. A PRISMA flow diagram of the study selection and inclusion process is presented in Figure 1 . As shown in Tables 2 and ​ and3, 3 , the articles included in this scoping review were published between 2010 and 2022 in Canada (n = 3), the United States (n = 2), Australia (n = 2), and Scotland (n = 1).

PRISMA flow diagram.

Characteristics of Articles Included.

Overview of Within-Method, Between/Across-Method, and Data-Analysis Triangulation.

Research Question, Case Definition, and Case Study Design

The following sections describe the research question, case definition, and case study design. Case studies are most appropriate when asking “how” or “why” questions. 1 According to Yin, 1 how and why questions are explanatory and lead to the use of case studies, histories, and experiments as the preferred research methods. In 1 study from Canada, eg, the following research question was presented: “How and why did stakeholders participate in the system change process that led to the introduction of the first nurse practitioner-led Clinic in Ontario?” (p7) 19 Once the research question has been formulated, the case should be defined and, subsequently, the case study design chosen. 1 In typical case studies with mixed methods, the 2 types of data are gathered concurrently in a convergent design and the results merged to examine a case and/or compare multiple cases. 10

Research question

“How” or “why” questions were found in 4 studies. 16 , 17 , 19 , 22 Two studies additionally asked “what” questions. Three studies described an exploratory approach, and 1 study presented an explanatory approach. Of these 4 studies, 3 studies chose a qualitative approach 17 , 19 , 22 and 1 opted for mixed methods with a convergent design. 16

In the remaining studies, either the research questions were not clearly stated or no “how” or “why” questions were formulated. For example, “what” questions were found in 1 study. 21 No information was provided on exploratory, descriptive, and explanatory approaches. Schadewaldt et al 21 chose mixed methods with a convergent design.

Case definition and case study design

A total of 5 studies defined the case as an organizational unit. 17 , 18 - 20 , 22 Of the 8 articles, 4 reported multiple-case studies. 16 , 17 , 22 , 23 Another 2 publications involved single-case studies. 19 , 20 Moreover, 2 publications did not state the case study design explicitly.

Within-Method Triangulation

This section describes within-method triangulation, which involves employing at least 2 data-collection procedures within the same design approach. 6 , 7 This can also be called data source triangulation. 8 Next, we present the single data-collection procedures in detail. In 5 studies, information on within-method triangulation was found. 15 , 17 - 19 , 22 Studies describing a quantitative approach and the triangulation of 2 or more quantitative data-collection procedures could not be included in this scoping review.

Qualitative approach

Five studies used qualitative data-collection procedures. Two studies combined face-to-face interviews and documents. 15 , 19 One study mixed in-depth interviews with observations, 18 and 1 study combined face-to-face interviews and documentation. 22 One study contained face-to-face interviews, observations, and documentation. 17 The combination of different qualitative data-collection procedures was used to present the case context in an authentic and complex way, to elicit the perspectives of the participants, and to obtain a holistic description and explanation of the cases under study.

All 5 studies used qualitative interviews as the primary data-collection procedure. 15 , 17 - 19 , 22 Face-to-face, in-depth, and semi-structured interviews were conducted. The topics covered in the interviews included processes in the introduction of new care services and experiences of barriers and facilitators to collaborative work in general practices. Two studies did not specify the type of interviews conducted and did not report sample questions. 15 , 18

Observations

In 2 studies, qualitative observations were carried out. 17 , 18 During the observations, the physical design of the clinical patients’ rooms and office spaces was examined. 17 Hungerford et al 18 did not explain what information was collected during the observations. In both studies, the type of observation was not specified. Observations were generally recorded as field notes.

Public documents

In 3 studies, various qualitative public documents were studied. 15 , 19 , 22 These documents included role description, education curriculum, governance frameworks, websites, and newspapers with information about the implementation of the role and general practice. Only 1 study failed to specify the type of document and the collected data. 15

Electronic health records

In 1 study, qualitative documentation was investigated. 17 This included a review of dashboards (eg, provider productivity reports or provider quality dashboards in the electronic health record) and quality performance reports (eg, practice-wide or co-management team-wide performance reports).

Between/Across-Method Triangulation

This section describes the between/across methods, which involve employing both qualitative and quantitative data-collection procedures in the same study. 6 , 7 This procedure can also be denoted “methodologic triangulation.” 8 Subsequently, we present the individual data-collection procedures. In 3 studies, information on between/across triangulation was found. 16 , 20 , 21

Mixed methods

Three studies used qualitative and quantitative data-collection procedures. One study combined face-to-face interviews, documentation, and self-assessments. 16 One study employed semi-structured interviews, direct observation, documents, and service records, 20 and another study combined face-to-face interviews, non-participant observation, documents, and questionnaires. 23

All 3 studies used qualitative interviews as the primary data-collection procedure. 16 , 20 , 23 Face-to-face and semi-structured interviews were conducted. In the interviews, data were collected on the introduction of new care services and experiences of barriers to and facilitators of collaborative work in general practices.

Observation

In 2 studies, direct and non-participant qualitative observations were conducted. 20 , 23 During the observations, the interaction between health professionals or the organization and the clinical context was observed. Observations were generally recorded as field notes.

In 2 studies, various qualitative public documents were examined. 20 , 23 These documents included role description, newspapers, websites, and practice documents (eg, flyers). In the documents, information on the role implementation and role description of NPs was collected.

Individual journals

In 1 study, qualitative individual journals were studied. 16 These included reflective journals from NPs, who performed the role in primary health care.

Service records

Only 1 study involved quantitative service records. 20 These service records were obtained from the primary care practices and the respective health authorities. They were collected before and after the implementation of an NP role to identify changes in patients’ access to health care, the volume of patients served, and patients’ use of acute care services.

Questionnaires/Assessment

In 2 studies, quantitative questionnaires were used to gather information about the teams’ satisfaction with collaboration. 16 , 21 In 1 study, 3 validated scales were used. The scales measured experience, satisfaction, and belief in the benefits of collaboration. 21 Psychometric performance indicators of these scales were provided. However, the time points of data collection were not specified; similarly, whether the questionnaires were completed online or by hand was not mentioned. A competency self-assessment tool was used in another study. 16 The assessment comprised 70 items and included topics such as health promotion, protection, disease prevention and treatment, the NP-patient relationship, the teaching-coaching function, the professional role, managing and negotiating health care delivery systems, monitoring and ensuring the quality of health care practice, and cultural competence. Psychometric performance indicators were provided. The assessment was completed online with 2 measurement time points (pre self-assessment and post self-assessment).

Data-Analysis Triangulation

This section describes data-analysis triangulation, which involves the combination of 2 or more methods of analyzing data. 6 Subsequently, we present within-case analysis and cross-case analysis.

Mixed-methods analysis

Three studies combined qualitative and quantitative methods of analysis. 16 , 20 , 21 Two studies involved deductive and inductive qualitative analysis, and qualitative data were analyzed thematically. 20 , 21 One used deductive qualitative analysis. 16 The method of analysis was not specified in the studies. Quantitative data were analyzed using descriptive statistics in 3 studies. 16 , 20 , 23 The descriptive statistics comprised the calculation of the mean, median, and frequencies.

Qualitative methods of analysis

Two studies combined deductive and inductive qualitative analysis, 19 , 22 and 2 studies only used deductive qualitative analysis. 15 , 18 Qualitative data were analyzed thematically in 1 study, 22 and data were treated with content analysis in the other. 19 The method of analysis was not specified in the 2 studies.

Within-case analysis

In 7 studies, a within-case analysis was performed. 15 - 20 , 22 Six studies used qualitative data for the within-case analysis, and 1 study employed qualitative and quantitative data. Data were analyzed separately, consecutively, or in parallel. The themes generated from qualitative data were compared and then summarized. The individual cases were presented mostly as a narrative description. Quantitative data were integrated into the qualitative description with tables and graphs. Qualitative and quantitative data were also presented as a narrative description.

Cross-case analyses

Of the multiple-case studies, 5 carried out cross-case analyses. 15 - 17 , 20 , 22 Three studies described the cross-case analysis using qualitative data. Two studies reported a combination of qualitative and quantitative data for the cross-case analysis. In each multiple-case study, the individual cases were contrasted to identify the differences and similarities between the cases. One study did not specify whether a within-case or a cross-case analysis was conducted. 23

Confirmation or contradiction of data

This section describes confirmation or contradiction through qualitative and quantitative data. 1 , 4 Qualitative and quantitative data were reported separately, with little connection between them. As a result, the conclusions on neither the comparisons nor the contradictions could be clearly determined.

Confirmation or contradiction among qualitative data

In 3 studies, the consistency of the results of different types of qualitative data was highlighted. 16 , 19 , 21 In particular, documentation and interviews or interviews and observations were contrasted:

• Confirmation between interviews and documentation: The data from these sources corroborated the existence of a common vision for an NP-led clinic. 19
• Confirmation among interviews and observation: NPs experienced pressure to find and maintain their position within the existing system. Nurse practitioners and general practitioners performed complete episodes of care, each without collaborative interaction. 21
• Contradiction among interviews and documentation: For example, interviewees mentioned that differentiating the scope of practice between NPs and physicians is difficult as there are too many areas of overlap. However, a clear description of the scope of practice for the 2 roles was provided. 21

Confirmation through a combination of qualitative and quantitative data

Both types of data showed that NPs and general practitioners wanted to have more time in common to discuss patient cases and engage in personal exchanges. 21 In addition, the qualitative and quantitative data confirmed the individual progression of NPs from less competent to more competent. 16 One study pointed out that qualitative and quantitative data obtained similar results for the cases. 20 For example, integrating NPs improved patient access by increasing appointment availability.

Contradiction through a combination of qualitative and quantitative data

Although questionnaire results indicated that NPs and general practitioners experienced high levels of collaboration and satisfaction with the collaborative relationship, the qualitative results drew a more ambivalent picture of NPs’ and general practitioners’ experiences with collaboration. 21

Research Question and Design

The studies included in this scoping review evidenced various research questions. The recommended formats (ie, how or why questions) were not applied consistently. Therefore, no case study design should be applied because the research question is the major guide for determining the research design. 2 Furthermore, case definitions and designs were applied variably. The lack of standardization is reflected in differences in the reporting of these case studies. Generally, case study research is viewed as allowing much more freedom and flexibility. 5 , 24 However, this flexibility and the lack of uniform specifications lead to confusion.

Methodologic Triangulation

Methodologic triangulation, as described in the literature, can be somewhat confusing as it can refer to either data-collection methods or research designs. 6 , 8 For example, methodologic triangulation can allude to qualitative and quantitative methods, indicating a paradigmatic connection. Methodologic triangulation can also point to qualitative and quantitative data-collection methods, analysis, and interpretation without specific philosophical stances. 6 , 8 Regarding “data-collection methods with no philosophical stances,” we would recommend using the wording “data source triangulation” instead. Thus, the demarcation between the method and the data-collection procedures will be clearer.

Within-Method and Between/Across-Method Triangulation

Yin 1 advocated the use of multiple sources of evidence so that a case or cases can be investigated more comprehensively and accurately. Most studies included multiple data-collection procedures. Five studies employed a variety of qualitative data-collection procedures, and 3 studies used qualitative and quantitative data-collection procedures (mixed methods). In contrast, no study contained 2 or more quantitative data-collection procedures. In particular, quantitative data-collection procedures—such as validated, reliable questionnaires, scales, or assessments—were not used exhaustively. The prerequisites for using multiple data-collection procedures are availability, the knowledge and skill of the researcher, and sufficient financial funds. 1 To meet these prerequisites, research teams consisting of members with different levels of training and experience are necessary. Multidisciplinary research teams need to be aware of the strengths and weaknesses of different data sources and collection procedures. 1

Qualitative methods of analysis and results

When using multiple data sources and analysis methods, it is necessary to present the results in a coherent manner. Although the importance of multiple data sources and analysis has been emphasized, 1 , 5 the description of triangulation has tended to be brief. Thus, traceability of the research process is not always ensured. The sparse description of the data-analysis triangulation procedure may be due to the limited number of words in publications or the complexity involved in merging the different data sources.

Only a few concrete recommendations regarding the operationalization of the data-analysis triangulation with the qualitative data process were found. 25 A total of 3 approaches have been proposed 25 : (1) the intuitive approach, in which researchers intuitively connect information from different data sources; (2) the procedural approach, in which each comparative or contrasting step in triangulation is documented to ensure transparency and replicability; and (3) the intersubjective approach, which necessitates a group of researchers agreeing on the steps in the triangulation process. For each case study, one of these 3 approaches needs to be selected, carefully carried out, and documented. Thus, in-depth examination of the data can take place. Farmer et al 25 concluded that most researchers take the intuitive approach; therefore, triangulation is not clearly articulated. This trend is also evident in our scoping review.

Mixed-methods analysis and results

Few studies in this scoping review used a combination of qualitative and quantitative analysis. However, creating a comprehensive stand-alone picture of a case from both qualitative and quantitative methods is challenging. Findings derived from different data types may not automatically coalesce into a coherent whole. 4 O’Cathain et al 26 described 3 techniques for combining the results of qualitative and quantitative methods: (1) developing a triangulation protocol; (2) following a thread by selecting a theme from 1 component and following it across the other components; and (3) developing a mixed-methods matrix.

The most detailed description of the conducting of triangulation is the triangulation protocol. The triangulation protocol takes place at the interpretation stage of the research process. 26 This protocol was developed for multiple qualitative data but can also be applied to a combination of qualitative and quantitative data. 25 , 26 It is possible to determine agreement, partial agreement, “silence,” or dissonance between the results of qualitative and quantitative data. The protocol is intended to bring together the various themes from the qualitative and quantitative results and identify overarching meta-themes. 25 , 26

The “following a thread” technique is used in the analysis stage of the research process. To begin, each data source is analyzed to identify the most important themes that need further investigation. Subsequently, the research team selects 1 theme from 1 data source and follows it up in the other data source, thereby creating a thread. The individual steps of this technique are not specified. 26 , 27

A mixed-methods matrix is used at the end of the analysis. 26 All the data collected on a defined case are examined together in 1 large matrix, paying attention to cases rather than variables or themes. In a mixed-methods matrix (eg, a table), the rows represent the cases for which both qualitative and quantitative data exist. The columns show the findings for each case. This technique allows the research team to look for congruency, surprises, and paradoxes among the findings as well as patterns across multiple cases. In our review, we identified only one of these 3 approaches in the study by Roots and MacDonald. 20 These authors mentioned that a causal network analysis was performed using a matrix. However, no further details were given, and reference was made to a later publication. We could not find this publication.

Case Studies in Nursing Research and Recommendations

Because it focused on the implementation of NPs in primary health care, the setting of this scoping review was narrow. However, triangulation is essential for research in this area. This type of research was found to provide a good basis for understanding methodologic and data-analysis triangulation. Despite the lack of traceability in the description of the data and methodological triangulation, we believe that case studies are an appropriate design for exploring new nursing roles in existing health care systems. This is evidenced by the fact that case study research is widely used in many social science disciplines as well as in professional practice. 1 To strengthen this research method and increase the traceability in the research process, we recommend using the reporting guideline and reporting checklist by Rodgers et al. 9 This reporting checklist needs to be complemented with methodologic and data-analysis triangulation. A procedural approach needs to be followed in which each comparative step of the triangulation is documented. 25 A triangulation protocol or a mixed-methods matrix can be used for this purpose. 26 If there is a word limit in a publication, the triangulation protocol or mixed-methods matrix needs to be identified. A schematic representation of methodologic and data-analysis triangulation in case studies can be found in Figure 2 .

Schematic representation of methodologic and data-analysis triangulation in case studies (own work).

Limitations

This study suffered from several limitations that must be acknowledged. Given the nature of scoping reviews, we did not analyze the evidence reported in the studies. However, 2 reviewers independently reviewed all the full-text reports with respect to the inclusion criteria. The focus on the primary care setting with NPs (master’s degree) was very narrow, and only a few studies qualified. Thus, possible important methodological aspects that would have contributed to answering the questions were omitted. Studies describing the triangulation of 2 or more quantitative data-collection procedures could not be included in this scoping review due to the inclusion and exclusion criteria.

Conclusions

Given the various processes described for methodologic and data-analysis triangulation, we can conclude that triangulation in case studies is poorly standardized. Consequently, the traceability of the research process is not always given. Triangulation is complicated by the confusion of terminology. To advance case study research in nursing, we encourage authors to reflect critically on methodologic and data-analysis triangulation and use existing tools, such as the triangulation protocol or mixed-methods matrix and the reporting guideline checklist by Rodgers et al, 9 to ensure more transparent reporting.

Supplemental Material

Acknowledgments.

The authors thank Simona Aeschlimann for her support during the screening process.

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding: The author(s) received no financial support for the research, authorship, and/or publication of this article.

Supplemental Material: Supplemental material for this article is available online.

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

• View all journals
• Explore content
• About the journal
• Publish with us
• Sign up for alerts
• Open access
• Published: 02 September 2024

Green spaces provide substantial but unequal urban cooling globally

• Yuxiang Li 1 ,
• Jens-Christian Svenning   ORCID: orcid.org/0000-0002-3415-0862 2 ,
• Weiqi Zhou   ORCID: orcid.org/0000-0001-7323-4906 3 , 4 , 5 ,
• Kai Zhu   ORCID: orcid.org/0000-0003-1587-3317 6 ,
• Jesse F. Abrams   ORCID: orcid.org/0000-0003-0411-8519 7 ,
• Timothy M. Lenton   ORCID: orcid.org/0000-0002-6725-7498 7 ,
• William J. Ripple 8 ,
• Zhaowu Yu   ORCID: orcid.org/0000-0003-4576-4541 9 ,
• Shuqing N. Teng 1 ,
• Robert R. Dunn 10 &
• Chi Xu   ORCID: orcid.org/0000-0002-1841-9032 1  

Nature Communications volume  15 , Article number:  7108 ( 2024 ) Cite this article

154 Altmetric

Metrics details

• Climate-change mitigation
• Urban ecology

Climate warming disproportionately impacts countries in the Global South by increasing extreme heat exposure. However, geographic disparities in adaptation capacity are unclear. Here, we assess global inequality in green spaces, which urban residents critically rely on to mitigate outdoor heat stress. We use remote sensing data to quantify daytime cooling by urban greenery in the warm seasons across the ~500 largest cities globally. We show a striking contrast, with Global South cities having ~70% of the cooling capacity of cities in the Global North (2.5 ± 1.0 °C vs. 3.6 ± 1.7 °C). A similar gap occurs for the cooling adaptation benefits received by an average resident in these cities (2.2 ± 0.9 °C vs. 3.4 ± 1.7 °C). This cooling adaptation inequality is due to discrepancies in green space quantity and quality between cities in the Global North and South, shaped by socioeconomic and natural factors. Our analyses further suggest a vast potential for enhancing cooling adaptation while reducing global inequality.

Similar content being viewed by others

Global climate-driven trade-offs between the water retention and cooling benefits of urban greening

Water, energy and climate benefits of urban greening throughout Europe under different climatic scenarios

Greenery as a mitigation and adaptation strategy to urban heat

Introduction.

Heat extremes are projected to be substantially intensified by global warming 1 , 2 , imposing a major threat to human mortality and morbidity in the coming decades 3 , 4 , 5 , 6 . This threat is particularly concerning as a majority of people now live in cities 7 , including those cities suffering some of the hottest climate extremes. Cities face two forms of warming: warming due to climate change and warming due to the urban heat island effect 8 , 9 , 10 . These two forms of warming have the potential to be additive, or even multiplicative. Climate change in itself is projected to result in rising maximum temperatures above 50 °C for a considerable fraction of the world if 2 °C global warming is exceeded 2 ; the urban heat island effect will cause up to >10 °C additional (surface) warming 11 . Exposures to temperatures above 35 °C with high humidity or above 40 °C with low humidity can lead to lethal heat stress for humans 12 . Even before such lethal temperatures are reached, worker productivity 13 and general health and well-being 14 can suffer. Heat extremes are especially risky for people living in the Global South 15 , 16 due to warmer climates at low latitudes. Climate models project that the lethal temperature thresholds will be exceeded with increasing frequencies and durations, and such extreme conditions will be concentrated in low-latitude regions 17 , 18 , 19 . These low-latitude regions overlap with the major parts of the Global South where population densities are already high and where population growth rates are also high. Consequently, the number of people exposed to extreme heat will likely increase even further, all things being equal 16 , 20 . That population growth will be accompanied by expanded urbanization and intensified urban heat island effects 21 , 22 , potentially exacerbating future Global North-Global South heat stress exposure inequalities.

Fortunately, we know that heat stress can be buffered, in part, by urban vegetation 23 . Urban green spaces, and especially urban forests, have proven an effective means through which to ameliorate heat stress through shading 24 , 25 and transpirational cooling 26 , 27 . The buffering effect of urban green spaces is influenced by their area (relative to the area of the city) and their spatial configuration 28 . In this context, green spaces become a kind of infrastructure that can and should be actively managed. At broad spatial scales, the effect of this urban green infrastructure is also mediated by differences among regions, whether in their background climate 29 , composition of green spaces 30 , or other factors 31 , 32 , 33 , 34 . The geographic patterns of the buffering effects of green spaces, whether due to geographic patterns in their areal extent or region-specific effects, have so far been poorly characterized.

On their own, the effects of climate change and urban heat islands on human health are likely to become severe. However, these effects will become even worse if they fall disproportionately in cities or countries with less economic ability to invest in green space 35 or in other forms of cooling 36 , 37 . A number of studies have now documented the so-called ‘luxury effect,’ wherein lower-income parts of cities tend to have less green space and, as a result, reduced biodiversity 38 , 39 . Where the luxury effect exists, green space and its benefits become, in essence, a luxury good 40 . If the luxury effect holds among cities, and lower-income cities also have smaller green spaces, the Global South may have the least potential to mitigate the combined effects of climate warming and urban heat islands, leading to exacerbated and rising inequalities in heat exposure 41 .

Here, we assess the global inequalities in the cooling capability of existing urban green infrastructure across urban areas worldwide. To this end, we use remotely sensed data to quantify three key variables, i.e., (1) cooling efficiency, (2) cooling capacity, and (3) cooling benefit of existing urban green infrastructure for ~500 major cities across the world. Urban green infrastructure and temperature are generally negatively and relatively linearly correlated at landscape scales, i.e., higher quantities of urban green infrastructure yield lower temperatures 42 , 43 . Cooling efficiency is widely used as a measure of the extent to which a given proportional increase in the area of urban green infrastructure leads to a decrease in temperature, i.e., the slope of the urban green infrastructure-temperature relationship 42 , 44 , 45 (see Methods for details). This simple metric allows quantifying the quality of urban green infrastructure in terms of ameliorating the urban heat island effect. Meanwhile, the extent to which existing urban green infrastructure cools down an entire city’s surface temperatures (compared to the non-vegetated built-up areas) is referred to as cooling capacity. Hence, cooling capacity is a function of the total quantity of urban green infrastructure and its cooling efficiency (see Methods).

As a third step, we account for the spatial distributions of urban green infrastructure and populations to quantify the benefit of cooling mitigation received by an average urban inhabitant in each city given their location. This cooling benefit is a more direct measure of the cooling realized by people, after accounting for the within-city geography of urban green infrastructure and population density. We focus on cooling capacity and cooling benefit as the measures of the cooling capability of individual cities for assessing their global inequalities. We are particularly interested in linking cooling adaptation inequality with income inequality 40 , 46 . While this can be achieved using existing income metrics for country classifications 47 , here we use the traditional Global North/South classification due to its historical ties to geography which is influential in climate research.

Results and discussion

Our analyses indicate that existing green infrastructure of an average city has a capability of cooling down surface temperatures by ~3 °C during warm seasons. However, a concerning disparity is evident; on average Global South cities have only two-thirds the cooling capacity and cooling benefit compared to Global North cities. This inequality is attributable to the differences in both quantity and quality of existing urban green infrastructure among cities. Importantly, we find that there exists considerable potential for many cities to enhance the cooling capability of their green infrastructure; achieving this potential could dramatically reduce global inequalities in adaptation to outdoor heat stress.

Quantifying cooling inequality

Our analyses showed that both the quantity and quality of the existing urban green infrastructure vary greatly among the world’s ~500 most populated cities (see Methods for details, and Fig.  1 for examples). The quantity of urban green infrastructure measured based on remotely sensed indicators of spectral greenness (Normalized Difference Vegetation Index, NDVI, see Methods) had a coefficient of variation (CV) of 35%. Similarly, the quality of urban green infrastructure in terms of cooling efficiency (daytime land surface temperatures during peak summer) had a CV of 37% (Supplementary Figs.  1 , 2 ). The global mean value of cooling capacity is 2.9 °C; existing urban green infrastructure ameliorates warm-season heat stress by 2.9 °C of surface temperature in an average city. In truth, however, the variation in cooling capacity was great (global CV in cooling capacity as large as ~50%), such that few cities were average. This variation is strongly geographically structured. Cities closer to the equator - tropical and subtropical cities - tend to have relatively weak cooling capacities (Fig.  2a, b ). As Global South countries are predominantly located at low latitudes, this pattern leads to a situation in which Global South cities, which tend to be hotter and relatively lower-income, have, on average, approximately two-thirds the cooling capacity of the Global North cities (2.5 ± 1.0 vs. 3.6 ± 1.7°C, Wilcoxon test, p  = 2.7e-12; Fig.  2c ). The cities that most need to rely on green infrastructure are, at present, those that are least able to do so.

a , e , i , m , q Los Angeles, US. b , f , j , n , r Paris, France. c , g , k , o , s Shanghai, China. d , h , l , p , t Cairo, Egypt. Local cooling efficiency is calculated for different local climate zone types to account for within-city heterogeneity. In densely populated parts of cities, local cooling capacity tends to be lower due to reduced green space area, whereas local cooling benefit (local cooling capacity multiplied by a weight term of local population density relative to city mean) tends to be higher as more urban residents can receive cooling amelioration.

a Global distribution of cooling capacity for the 468 major urbanized areas. b Latitudinal pattern of cooling capacity. c Cooling capacity difference between the Global North and South cities. The cooling capacity offered by urban green infrastructure evinces a latitudinal pattern wherein lower-latitude cities have weaker cooling capacity ( b , cubic-spline fitting of cooling capacity with 95% confidence interval is shown), representing a significant inequality between Global North and South countries: city-level cooling capacity for Global North cities are about 1.5-fold higher than in Global South cities ( c ). Data are presented as box plots, where median values (center black lines), 25th percentiles (box lower bounds), 75th percentiles (box upper bounds), whiskers extending to 1.5-fold of the interquartile range (IQR), and outliers are shown. The tails of the cooling capacity distributions are truncated at zero as all cities have positive values of cooling capacity. Notice that no cities in the Global South have a cooling capacity greater than 5.5 °C ( c ). This is because no cities in the Global South have proportional green space areas as great as those seen in the Global North (see also Fig.  4b ). A similar pattern is found for cooling benefit (Supplementary Fig.  3 ). The two-sided non-parametric Wilcoxon test was used for statistical comparisons.

When we account for the locations of urban green infrastructure relative to humans within cities, the cooling benefit of urban green infrastructure realized by an average urban resident generally becomes slightly lower than suggested by cooling capacity (see Methods; Supplementary Fig.  3 ). Urban residents tend to be densest in the parts of cities with less green infrastructure. As a result, the average urban resident experiences less cooling amelioration than expected. However, this heterogeneity has only a minor effect on global-scale inequality. As a result, the geographic trends in cooling capacity and cooling benefit are similar: mean cooling benefit for an average urban resident also presents a 1.5-fold gap between Global South and North cities (2.2 ± 0.9 vs. 3.4 ± 1.7 °C, Wilcoxon test, p  = 3.2e-13; Supplementary Fig.  3c ). Urban green infrastructure is a public good that has the potential to help even the most marginalized populations stay cool; unfortunately, this public benefit is least available in the Global South. When walking outdoors, the average person in an average Global South city receives only two-thirds the cooling amelioration from urban green infrastructure experienced by a person in an average Global North city. The high cooling amelioration capacity and benefit of the Global North cities is heavily influenced by North America (specifically, Canada and the US), which have both the highest cooling efficiency and the largest area of green infrastructure, followed by Europe (Supplementary Fig.  4 ).

One way to illustrate the global inequality of cooling capacity or benefit is to separately look at the cities that are most and least effective in ameliorating outdoor heat stress. Our results showed that ~85% of the 50 most effective cities (with highest cooling capacity or cooling benefit) are located in the Global North, while ~80% of the 50 least effective are Global South cities (Fig.  3 , Supplementary Fig.  5 ). This is true without taking into account the differences in the background temperatures and climate warming of these cities, which will exacerbate the effects on human health; cities in the Global South are likely to be closer to the limits of human thermal comfort and even, increasingly, the limits of the temperatures and humidities (wet-bulb temperatures) at which humans can safely work or even walk, such that the ineffectiveness of green spaces in those cities in cooling will lead to greater negative effects on human health 48 , work 14 , and gross domestic product (GDP) 49 . In addition, Global South cities commonly have higher population densities (Fig.  3 , Supplementary Fig.  5 ) and are projected to have faster population growth 50 . This situation will plausibly intensify the urban heat island effect because of the need of those populations for housing (and hence tensions between the need for buildings and the need for green spaces). It will also increase the number of people exposed to extreme urban heat island effects. Therefore, it is critical to increase cooling benefit via expanding urban green spaces, so that more people can receive the cooling mitigation from a given new neighboring green space if they live closer to each other. Doing so will require policies that incentivize urban green spaces as well as architectural innovations that make innovations such as plant-covered buildings easier and cheaper to implement.

The axes on the right are an order of magnitude greater than those on the left, such that the cooling capacity of Charlotte in the United States is about 37-fold greater than that of Mogadishu (Somalia) and 29-fold greater than that of Sana’a (Yemen). The cities presenting lowest cooling capacities are most associated with Global South cities at higher population densities.

Of course, cities differ even within the Global North or within the Global South. For example, some Global South cities have high green space areas (or relatively high cooling efficiency in combination with moderate green space areas) and hence high cooling capacity. These cities, such as Pune (India), will be important to study in more detail, to shed light on the mechanistic details of their cooling abilities as well as the sociopolitical and other factors that facilitated their high green area coverage and cooling capabilities (Supplementary Figs.  6 , 7 ).

We conducted our primary analyses using a spatial grain of 100-m grid cells and Landsat NDVI data for quantifying spectral greenness. Our results, however, were robust at the coarser spatial grain of 1 km. We find a slightly larger global cooling inequality (~2-fold gap between Global South and North cities) at the 1-km grain using MODIS data (see Methods and Supplementary Fig.  17 ). MODIS data have been frequently used for quantifying urban heat island effects and cooling mitigation 44 , 45 , 51 . Our results reinforce its robustness for comparing urban thermal environments between cities across broad scales.

Influencing factors

The global inequality of cooling amelioration could have a number of proximate causes. To understand their relative influence, we first separately examined the effects of quality (cooling efficiency) and quantity (NDVI as a proxy indicator of urban green space area) of urban green infrastructure. The simplest null model is one in which cooling capacity (at the city scale) and cooling benefit (at the human scale) are driven primarily by the proportional area in a city dedicated to green spaces. Indeed, we found that both cooling capacity and cooling benefit were strongly correlated with urban green space area (Fig.  4 , Supplementary Fig.  8 ). This finding is useful with regards to practical interventions. In general, cities that invest in saving or restoring more green spaces will receive more cooling benefits from those green spaces. By contrast, differences among cities in cooling efficiency played a more minor role in determining the cooling capacity and benefit of cities (Fig.  4 , Supplementary Fig.  8 ).

a Relationship between cooling efficiency and cooling capacity. b Relationship between green space area (measured by mean Landsat NDVI in the hottest month of 2018) and cooling capacity. Note that the highest level of urban green space area in the Global South cities is much lower than that in the Global North (dashed line in b ). Gray bands indicate 95% confidence intervals. Two-sided t-tests were conducted. c A piecewise structural equation model based on assumed direct and indirect (through influencing cooling efficiency and urban green space area) effects of essential natural and socioeconomic factors on cooling capacity. Mean annual temperature and precipitation, and topographic variation (elevation range) are selected to represent basic background natural conditions; GDP per capita is selected to represent basic socioeconomic conditions. The spatial extent of built-up areas is included to correct for city size. A bi-directional relationship (correlation) is fitted between mean annual temperature and precipitation. Red and blue solid arrows indicate significantly negative and positive coefficients with p  ≤ 0.05, respectively. Gray dashed arrows indicate p  > 0.05. The arrow width illustrates the effect size. Similar relationships are found for cooling benefits realized by an average urban resident (see Supplementary Fig.  8 ).

A further question is what shapes the quality and quantity of urban green infrastructure (which in turn are driving cooling capacity)? Many inter-correlated factors are possibly operating at multiple scales, making it difficult to disentangle their effects, especially since experiment-based causal inference is usually not feasible for large-scale urban systems. From a macroscopic perspective, we test the simple hypothesis that the background natural and socioeconomic conditions of cities jointly affect their cooling capacity and benefit in both direct and indirect ways. To this end, we constructed a minimal structural equation model including only the most essential variables reflecting background climate (mean annual temperature and precipitation), topographic variation (elevation range), as well as gross domestic product (GDP) per capita and city area (see Methods; Fig.  4c ).

We found that the quantity of green spaces in a city (again, in proportion to its size) was positively correlated with GDP per capita and city area; wealthier cities have more green spaces. It is well known that wealth and green spaces are positively correlated within cities (the luxury effect) 40 , 46 ; our analysis shows that a similar luxury effect occurs among them at a global scale. In addition, larger cities often have proportionally more green spaces, an effect that may be due to the tendency for large cities (particularly in the US and Canada) to have lower population densities. Cities that were hotter and had more topographic variation tended to have fewer green spaces and those that were more humid tended to have more green spaces. Given that temperature and humidity are highly correlated with the geography of the Global South and Global North, it is difficult to know whether these effects are due to the direct effects of temperature and precipitation, for example, on the growth rate of vegetation and hence the transition of abandoned lots into green spaces, or are associated with historical, cultural and political differences that via various mechanisms correlate to climate. Our structural equation model explained only a small fraction of variation among cities in their cooling efficiency, which is to say the quality of their green space. Cooling efficiency was modestly influenced by background temperature and precipitation—the warmer a city, the greater the cooling efficiency in that city; conversely, the more humid a city the less the cooling efficiency of that city.

Our analyses suggested that the lower cooling adaptation capabilities of Global South cities can be explained by their lower quantity of green infrastructure and, to a much lesser extent, their weaker cooling efficiency (quality; Supplementary Fig.  2 ). These patterns appear to be in part structured by GDP, but are also associated with climatic conditions 39 , and other factors. A key question, unresolved by our work, is whether the climatic correlates of the size of green spaces in cities are due to the effects of climate per se or if they, instead, reflect correlates between contemporary climate and the social, cultural, and political histories of cities in the Global South 52 . Since urban planning has much inertia, especially in big cities, those choices might be correlated with climate because of the climatic correlates of political histories. It is also possible that these dynamics relate, in part, to the ways in which climate influences vegetation structure. However, this seems less likely given that under non-urban conditions vegetation cover (and hence cooling capacity) is normally positively correlated with mean annual temperature across the globe, opposite to our observed negative relationships for urban systems (Supplementary Fig.  9g ). Still, it is possible that increased temperatures in cities due to the urban heat island effects may lead to temperature-vegetation cover-cooling capacity relationships that differ from those in natural environments 53 , 54 . Indeed, a recent study found that climate warming will put urban forests at risk, and the risk is disproportionately higher in the Global South 55 .

Our model serves as a starting point for unraveling the mechanisms underlying global cooling inequality. We cannot rule out the possibility that other unconsidered factors correlated with the studied variables play important roles. We invite systematic studies incorporating detailed sociocultural and ecological variables to address this question across scales.

Potential of enhancing cooling and reducing inequality

Can we reduce the inequality in cooling capacity and benefits that we have discovered among the world’s largest cities? Nuanced assessments of the potential to improve cooling mitigation require comprehensive considerations of socioeconomic, cultural, and technological aspects of urban management and policy. It is likely that cities differ greatly in their capacity to implement cooling through green infrastructure, whether as a function of culture, governance, policy or some mix thereof. However, any practical attempts to achieve greater cooling will occur in the context of the realities of climate and existing land use. To understand these realities, we modeled the maximum additional cooling capacity that is possible in cities, given existing constraints. We assume that this capacity depends on the quality (cooling efficiency) and quantity of urban green infrastructure. Our approach provides a straightforward metric of the cooling that could be achieved if all parts of a city’s green infrastructure were to be enhanced systematically.

The positive outlook is that our analyses suggest a considerable potential of improving cooling capacity by optimizing urban green infrastructure. An obvious way is through increases in urban green infrastructure quantity. We employ an approach in which we consider each local climate zone 56 to have a maximum NDVI and cooling efficiency (see Methods). For a given local climate zone, the city with the largest NDVI values or cooling efficiency sets the regional upper bounds for urban green infrastructure quantities or quality that can be achieved. Notably, these maxima are below the maxima for forests or other non-urban spaces for the simple reason that, as currently imagined, cities must contain gray (non-green) spaces in the form of roads and buildings. In this context, we conduct a thought experiment. What if we could systematically increase NDVI of all grid cells in each city, per local climate zone type, to a level corresponding to the median NDVI of grid cells in that upper bound city while keeping cooling efficiency unchanged (see Methods). If we were able to achieve this goal, the cooling capacity of cities would increase by ~2.4 °C worldwide. The increase would be even greater, ~3.8°C, if the 90th percentile (within the reference maximum city) was reached (Fig.  5a ). The potential for cooling benefit to the average urban resident is similar to that of cooling capacity (Supplementary Fig.  10a ). There is also potential to reduce urban temperatures if we can enhance cooling efficiency. However, the benefits of increases in cooling efficiency are modest (~1.5 °C increases at the 90th percentile of regional upper bounds) when holding urban green infrastructure quantity constant. In theory, if we could maximize both quantity and cooling efficiency of urban green infrastructure (to 90th percentiles of their regional upper bounds respectively), we would yield increases in cooling capacity and benefit up to ~10 °C, much higher than enhancing green space area or cooling efficiency alone (Fig.  5a , Supplementary Fig.  10a ). Notably, such co-maximization of green space area and cooling efficiency would substantially reduce global inequality to Gini <0.1 (Fig.  5b , Supplementary Fig.  10b ). Our analyses thus provide an important suggestion that enhancing both green space quantity and quality can yield a synergistic effect leading to much larger gains than any single aspect alone.

a The potential of enhancing cooling capacity via either enhancing urban green infrastructure quality (i.e., cooling efficiency) while holding quantity (i.e., green space area) fixed (yellow), or enhancing quantity while holding quality fixed (blue) is much lower than that of enhancing both quantity and quality (green). The x-axis indicates the targets of enhancing urban green infrastructure quantity and/or quality relative to the 50–90th percentiles of NDVI or cooling efficiency, see Methods). The dashed horizontal lines indicate the median cooling capacity of current cities. Data are presented as median values with the colored bands corresponding to 25–75th percentiles. b The potential of reducing cooling capacity inequality is also higher when enhancing both urban green infrastructure quantity and quality. The Gini index weighted by population density is used to measure inequality. Similar results were found for cooling benefit (Supplementary Fig.  10 ).

Different estimates of cooling capacity potential may be reached based on varying estimates and assumptions regarding the maximum possible quantity and quality of urban green infrastructure. There is no single, simple way to make these estimates, especially considering the huge between-city differences in society, culture, and structure across the globe. Our example case (above) begins from the upper bound city’s median NDVI, taking into account different local climate zone types and background climate regions (regional upper bounds). This is based on the assumption that for cities within the same climate regions, their average green space quantity may serve as an attainable target. Still, urban planning is often made at the level of individual cities, often only implemented to a limited extent and made with limited consideration of cities in other regions and countries. A potentially more realistic reference may be taken from the existing green infrastructure (again, per local climate zone type) within each particular city itself (see Methods): if a city’s sparsely vegetated areas was systematically elevated to the levels of 50–90th percentiles of NDVI within their corresponding local climate zones within the city, cooling capacity would still increase, but only by 0.5–1.5 °C and with only slightly reduced inequalities among cities (Supplementary Fig.  11 ). This highlights that ambitious policies, inspired by the greener cities worldwide, are necessary to realize the large cooling potential in urban green infrastructure.

In summary, our results demonstrate clear inequality in the extent to which urban green infrastructure cools cities and their denizens between the Global North and South. Much attention has been paid to the global inequality of indoor heat adaptation arising from the inequality of resources (e.g., less affordable air conditioning and more frequent power shortages in the Global South) 36 , 57 , 58 , 59 . Our results suggest that the inequality in outdoor adaptation is particularly concerning, especially as urban populations in the Global South are growing rapidly and are likely to face the most severe future temperature extremes 60 .

Previous studies have been focusing on characterizing urban heat island effects, urban vegetation patterns, resident exposure, and cooling effects in particular cities 26 , 28 , 34 , 61 , regions 22 , 25 , 62 , or continents 32 , 44 , 63 . Recent studies start looking at global patterns with respect to cooling efficiency or green space exposure 35 , 45 , 64 , 65 . Our approach is one drawn from the fields of large-scale ecology and macroecology. This approach is complementary to and, indeed, can, in the future, be combined with (1) mechanism driven biophysical models 66 , 67 to predict the influence of the composition and climate of green spaces on their cooling efficiency, (2) social theory aimed at understanding the factors that govern the amount of green space in cities as well as the disparity among cities 68 , (3) economic models of the effects of policy changes on the amount of greenspace and even (4) artist-driven projects that seek to understand the ways in which we might reimagine future cities 69 . Our simple explanatory model is, ultimately, one lens on a complex, global phenomenon.

Our results convey some positive outlook in that there is considerable potential to strengthen the cooling capability of cities and to reduce inequalities in cooling capacities at the same time. Realizing this nature-based solution, however, will be challenging. First, enhancing urban green infrastructure requires massive investments, which are more difficult to achieve in Global South cities. Second, it also requires smart planning strategies and advanced urban design and greening technologies 37 , 70 , 71 , 72 . Spatial planning of urban green spaces needs to consider not only the cooling amelioration effect, but also their multifunctional aspects that involve multiple ecosystem services, mental health benefits, accessibility, and security 73 . In theory, a city can maximize its cooling while also maximizing density through the combination of high-density living, ground-level green spaces, and vertical and rooftop gardens (or even forests). In practice, the current cities with the most green spaces tend to be lower-density cities 74 (Supplementary Fig.  12 ). Still, innovation and implementation of new technologies that allow green spaces and high-density living to be combined have the potential to reduce or disconnect the negative relationship between green space area and population density 71 , 75 . However, this development has yet to be realized. Another dimension of green spaces that deserves more attention is the geography of green spaces relative to where people are concentrated within cities. A critical question is how best should we distribute green spaces within cities to maximize cooling efficiency 76 and minimize within-city cooling inequality towards social equity 77 ? Last but not least, it is crucial to design and manage urban green spaces to be as resilient as possible to future climate stress 78 . For many cities, green infrastructure is likely to remain the primary means people will have to rely on to mitigate the escalating urban outdoor heat stress in the coming decades 79 .

We used the world population data from the World’s Cities in 2018 Data Booklet 80 to select 502 major cities with population over 1 million people (see Supplementary Data  1 for the complete list of the studied cities). Cities are divided into the Global North and Global South based on the Human Development Index (HDI) from the Human Development Report 2019 81 . For each selected city, we used the 2018 Global Artificial Impervious Area (GAIA) data at 30 m resolution 82 to determine its geographic extent. The derived urban boundary polygons thus encompass a majority of the built-up areas and urban residents. In using this approach, rather than urban administrative boundaries, we can focus on the relatively densely populated areas where cooling mitigation is most needed, and exclude areas dominated by (semi) natural landscapes that may bias the subsequent quantifications of the cooling effect. Our analyses on the cooling effect were conducted at the 100 m spatial resolution using Landsat data and WorldPop Global Project Population Data of 2018 83 . In order to test for the robustness of the results to coarser spatial scales, we also repeated the analyses at 1 km resolution using MODIS data, which have been extensively used for quantifying urban heat island effects and cooling mitigation 44 , 45 , 51 . We discarded the five cities with sizes <30 km 2 as they were too small for us to estimate their cooling efficiency based on linear regression (see section below for details). We combined closely located cities that form contiguous urban areas or urban agglomerations, if their urban boundary polygons from GAIA merged (e.g., Phoenix and Mesa in the United States were combined). Our approach yielded 468 polygons, each representing a major urbanized area that were the basis for all subsequent analyses. Because large water bodies can exert substantial and confounding cooling effects, we excluded permanent water bodies including lakes, reservoirs, rivers, and oceans using the Copernicus Global Land Service (CGLS) Land Cover data for 2018 at 10 m resolution 84 .

Quantifying the cooling effect

As a first step, we calculated cooling efficiency for each studied city within the GAIA-derived urban boundary. Cooling efficiency quantifies the extent to which a given area of green spaces in a city can reduce temperatures. It is a measure of the effectiveness (quality) of urban green spaces in terms of heat amelioration. Cooling efficiency is typically measured by calculating the slope of the relationship between remotely-sensed land surface temperature (LST) and vegetation cover through ordinary least square regression 42 , 44 , 45 . It is known that cooling efficiency varies between cities. Influencing factors might include background climate 29 , species composition 30 , 85 , landscape configuration 28 , topography 86 , proximity to large water bodies 33 , 87 , urban morphology 88 , and city management practices 31 . However, the mechanism underlying the global pattern of cooling efficiency remains unclear.

We used Landsat satellite data provided by the United States Geological Survey (USGS) to calculate the cooling efficiency of each studied city. We used the cloud-free Landsat 8 Level 2 LST and NDVI data. For each city we calculated the mean LST in each month of 2018 to identify the hottest month, and then derived the hottest month LST; we used the cloud-free Landsat 8 data to calculate the mean NDVI for the hottest month correspondingly.

We quantified cooling efficiency for different local climate zones 56 separately for each city, to account for within-city variability of thermal environments. To this end, we used the Copernicus Global Land Service data (CGLS) 84 and Global Human Settlement Layers (GHSL) Built-up height data 89 of 2018 at the 100 m resolution to identify five types of local climate zones: non-tree vegetation (shrubs, herbaceous vegetation, and cultivated vegetation according to the CGLS classification system), low-rise buildings (built up and bare according to the CGLS classification system, with building heights ≤10 m according to the GHSL data), medium-high-rise buildings (built up and bare areas with building heights >10 m), open tree cover (open forest with tree cover 15–70% according to the CGLS system), and closed tree cover (closed forest with tree cover >70%).

For each local climate zone type in each city, we constructed a regression model with NDVI as the predictor variable and LST as the response variable (using the ordinary least square method). We took into account the potential confounding factors including topographic elevation (derived from MERIT DEM dataset 90 ), building height (derived from the GHSL dataset 89 ), and distance to water bodies (derived from the GSHHG dataset 91 ), the model thus became: LST ~ NDVI + topography + building height + distance to water. Cooling efficiency was calculated as the absolute value of the regression coefficient of NDVI, after correcting for those confounding factors. To account for the multi-collinearity issue, we conducted variable selection based on the variance inflation factor (VIF) to achieve VIF < 5. Before the analysis, we discarded low-quality Landsat pixels, and filtered out the pixels with NDVI < 0 (normally less than 1% in a single city). Cooling efficiency is known to be influenced by within-city heterogeneity 92 , 93 , and, as a result, might sometimes better fit non-linear relationships at local scales 65 , 76 . However, our central aim is to assess global cooling inequality based on generalized relationships that fit the majority of global cities. Previous studies have shown that linear relationships can do this job 42 , 44 , 45 , therefore, here we used linear models to assess cooling efficiency.

As a second step, we calculated the cooling capacity of each city. Cooling capacity is a positive function of the magnitude of cooling efficiency and the proportional area of green spaces in a city and is calculated based on NDVI and the derived cooling efficiency (Eq.  1 , Supplementary Fig.  13 ):

where CC lcz and CE lcz are the cooling capacity and cooling efficiency for a given local climate zone type in a city, respectively; NDVI i is the mean NDVI for 100-m grid cell i ; NDVI min is the minimum NDVI across the city; and n is the total number of grid cells within the local climate zone. Local cooling capacity for each grid cell i (Fig.  1 , Supplementary Fig.  7 ) can be derived in this way as well (Supplementary Fig.  13 ). For a particular city, cooling capacity may be dependent on the spatial configuration of its land use/cover 28 , 94 , but here we condensed cooling capacity to city average (Eq.  2 ), thus did not take into account these local-scale factors.

where CC is the average cooling capacity of a city; n lcz is the number of grid cells of the local climate zone; m is the total number of grid cells within the whole city.

As a third step, we calculated the cooling benefit realized by an average urban resident (cooling benefit in short) in each city. Cooling benefit depends not only on the cooling capacity of a city, but also on where people live within a city relative to greener or grayer areas of the city. For example, cooling benefits in a city might be low even if the cooling capacity is high if the green parts and the dense-population parts of a city are inversely correlated. Here, we are calculating these averages while aware that in any particular city the exposure of a particular person will depend on the distribution of green spaces in a city, and the occupation, movement trajectories of a person, etc. On the scale of a city, we calculated cooling benefit following a previous study 35 , that is, simply adding a weight term of population size per 100-m grid cell into cooling capacity in Eq. ( 1 ):

Where CB lcz is the cooling benefit of a given local climate zone type in a specific city, pop i is the number of people within grid cell i , \(\overline{{pop}}\) is the mean population of the city.

Where CB is the average cooling benefit of a city. The population data were obtained from the 100-m resolution WorldPop Global Project Population Data of 2018 83 . Local cooling benefit for a given grid cell i can be calculated in a similar way, i.e., local cooling capacity multiplied by a weight term of local population density relative to mean population density. Local cooling benefits were mapped for example cities for the purpose of illustrating the effect of population spatial distribution (Fig.  1 , Supplementary Fig.  7 ), but their patterns were not examined here.

Based on the aforementioned three key variables quantified at 100 m grid cells, we conducted multivariate analyses to examine if and to what extent cooling efficiency and cooling benefit are shaped by essential natural and socioeconomic factors, including background climate (mean annual temperature from ECMWF ERA5 dataset 95 and precipitation from TerraClimate dataset 96 ), topography (elevation range 90 ), and GDP per capita 97 , with city size (geographic extent) corrected for. We did not include humidity because it is strongly correlated with temperature and precipitation, causing serious multi-collinearity problems. We used piecewise structural equation modeling to test the direct effects of these factors and indirect effects via influencing cooling efficiency and vegetation cover (Fig.  4c , Supplementary Fig.  8c ). To account for the potential influence of spatial autocorrelation, we used spatially autoregressive models (SAR) to test for the robustness of the observed effects of natural and socioeconomic factors on cooling capacity and benefit (Supplementary Fig.  14 ).

Testing for robustness

We conducted the following additional analyses to test for robustness. We obtained consistent results from these robustness analyses.

(1) We looked at the mean hottest-month LST and NDVI within 3 years (2017-2019) to check the consistency between the results based on relatively short (1 year) vs. long (3-year average) time periods (Supplementary Fig.  15 ).

(2) We carried out the approach at a coarser spatial scale of 1 km, using MODIS-derived NDVI and LST, as well as the population data 83 in the hottest month of 2018. In line with our finer-scale analysis of Landsat data, we selected the hottest month and excluded low-quality grids affected by cloud cover and water bodies 98 (water cover > 20% in 1 × 1 km 2 grid cells) of MODIS LST, and calculated the mean NDVI for the hottest month. We ultimately obtained 441 cities (or urban agglomerations) for analysis. At the 1 km resolution, some local climate zone types would yield insufficient samples for constructing cooling efficiency models. Therefore, instead of identifying local climate zone explicitly, we took an indirect approach to account for local climate confounding factors, that is, we constructed a multiple regression model for a whole city incorporating the hottest-month local temperature 95 , precipitation 96 , and humidity (based on NASA FLDAS dataset 99 ), albedo (derived from the MODIS MCD43A3 product 100 ), aerosol loading (derived from the MODIS MCD19A2 product 101 ), wind speed (based on TerraClimate dataset 96 ), topography elevation 90 , distance to water 91 , urban morphology (building height 102 ), and human activity intensity (VIIRS nighttime light data as a proxy indicator 103 ). We used the absolute value of the linear regression coefficient of NDVI as the cooling efficiency of the whole city (model: LST ~ NDVI + temperature + precipitation + humidity + distance to water + topography + building height + albedo + aerosol + wind speed + nighttime light), and calculated cooling capacity and cooling benefit based on the same method. Variable selection was conducted using the criterion of VIF < 5.

Our results indicated that MODIS-based cooling capacity and cooling benefit are significantly correlated with the Landsat-based counterparts (Supplementary Fig.  16 ); importantly, the gap between the Global South and North cities is around two-fold, close to the result from the Landsat-based result (Supplementary Fig.  17 ).

(3) For the calculation of cooling benefit, we considered different spatial scales of human accessibility to green spaces: assuming the population in each 100 × 100 m 2 grid cell could access to green spaces within neighborhoods of certain extents, we calculated cooling benefit by replacing NDVI i in Eq. ( 3 ) with mean NDVI within the 300 × 300 m 2 and 500 × 500 m 2 extents centered at the focal grid cell (Supplementary Fig.  18 ).

(4) Considering cities may vary in minimum NDVI, we assessed if this variation could affect resulting cooling capacity patterns. To this end, we calculated the cooling capacity for each studied city using NDVI = 0 as the reference (i.e., using NDVI = 0 instead of minimum NDVI in Supplementary Fig.  13b ), and correlated it with that using minimum NDVI as the reference (Supplementary Fig.  19 ).

Quantifying between-city inequality

Inequalities in access to the benefits of green spaces in cities exist within cities, as is increasingly well-documented 104 . Here, we focus instead on the inequalities among cities. We used the Gini coefficient to measure the inequality in cooling capacity and cooling benefit between all studied cities across the globe as well as between Global North or South cities. We calculated Gini using the population-density weighted method (Fig.  5b ), as well as the unweighted and population-size weighted methods (Supplementary Fig.  20 ).

Estimating the potential for more effective and equal cooling amelioration

We estimated the potential of enhancing cooling amelioration based on the assumptions that urban green space quality (cooling efficiency) and quantity (NDVI) can be increased to different levels, and that relative spatial distributions of green spaces and population can be idealized (so that their spatial matches can maximize cooling benefit). We assumed that macro-climate conditions act as the constraints of vegetation cover and cooling efficiency. We calculated the 50th, 60th, 70th, 80th, and 90th percentiles of NDVI within each type of local climate zone of each city. For a given local climate zone type, we obtained the city with the highest NDVI per percentile value as the regional upper bounds of urban green infrastructure quantity. The regional upper bounds of cooling efficiency are derived in a similar way. For each local climate zone in a city, we generated a potential NDVI distribution where all grid cells reach the regional upper bound values for the 50th, 60th, 70th, 80th, or 90th percentile of urban green space quantity or quality, respectively. NDVI values below these percentiles were increased, whereas those above these percentiles remained unchanged. The potential estimates are essentially dependent on the references, i.e., the optimal cooling efficiency and NDVI that a given city can reach. However, such references are obviously difficult to determine, because complex natural and socioeconomic conditions could play important roles in determining those cooling optima, and the dominant factors are unknown at a global scale. We employed the simplifying assumption that background climate could act as an essential constraint according to our results. We therefore used the Köppen climate classification system 105 to determine the reference separately in each climate region (tropical, arid, temperate, and continental climate regions were involved for all studied cities).

We calculated potential cooling capacity and cooling benefit based on these potential NDVI maps (Fixed cooling efficiency in Fig.  5 ). We then calculated the potentials if cooling efficiency of each city can be enhanced to 50–90th percentile across all urban local climate zones within the corresponding biogeographic region (Fixed green space area in Fig.  5 ). We also calculated the potentials if both NDVI and cooling efficiency were enhanced (Enhancing both in Fig.  5) to a certain corresponding level (i.e., i th percentile NDVI +  i th percentile cooling efficiency). We examined if there are additional effects of idealizing relative spatial distributions of urban green spaces and humans on cooling benefits. To this end, the pixel values of NDVI or population amount remained unchanged, but their one-to-one correspondences were based on their ranking: the largest population corresponds to the highest NDVI, and so forth. Under each scenario, we calculated cooling capacity and cooling benefit for each city, and the between-city inequality was measured by the Gini coefficient.

We used the Google Earth Engine to process the spatial data. The statistical analyses were conducted using R v4.3.3 106 , with car v3.1-2 107 , piecewiseSEM v2.1.2 108 , and ineq v0.2-13 109 packages. The global maps of cooling were created using the ArcGIS v10.3 software.

Reporting summary

Further information on research design is available in the  Nature Portfolio Reporting Summary linked to this article.

Data availability

City population statistics data is collected from the Population Division of the Department of Economic and Social Affairs of the United Nations ( https://www.un.org/development/desa/pd/content/worlds-cities-2018-data-booklet ). Global North-South division is based on Human Development Report 2019 which from United Nations Development Programme ( https://hdr.undp.org/content/human-development-report-2019 ). Global urban boundaries from GAIA data are available from Star Cloud Data Service Platform ( https://data-starcloud.pcl.ac.cn/resource/14 ) . Global water data is derived from 2018 Copernicus Global Land Service (CGLS 100-m) data ( https://developers.google.com/earth-engine/datasets/catalog/COPERNICUS_Landcover_100m_Proba-V-C3_Global ), European Space Agency (ESA) WorldCover 10 m 2020 product ( https://developers.google.com/earth-engine/datasets/catalog/ESA_WorldCover_v100 ), and GSHHG (A Global Self-consistent, Hierarchical, High-resolution Geography Database) at https://www.soest.hawaii.edu/pwessel/gshhg/ . Landsat 8 LST and NDVI data with 30 m resolution are available at  https://developers.google.com/earth-engine/datasets/catalog/LANDSAT_LC08_C02_T1_L2 . Land surface temperature (LST) data with 1 km from MODIS Aqua product (MYD11A1) is available at https://developers.google.com/earth-engine/datasets/catalog/MODIS_061_MYD11A1 . NDVI (1 km) dataset from MYD13A2 is available at https://developers.google.com/earth-engine/datasets/catalog/MODIS_061_MYD13A2 . Population data (100 m) is derived from WorldPop ( https://developers.google.com/earth-engine/datasets/catalog/WorldPop_GP_100m_pop ). Local climate zones are also based on 2018 CGLS data ( https://developers.google.com/earth-engine/datasets/catalog/COPERNICUS_Landcover_100m_Proba-V-C3_Global ), and built-up height data is available from Global Human Settlement Layers (GHSL, 100 m) ( https://developers.google.com/earth-engine/datasets/catalog/JRC_GHSL_P2023A_GHS_BUILT_H ). Temperature data is calculated from ERA5-Land Monthly Aggregated dataset ( https://developers.google.com/earth-engine/datasets/catalog/ECMWF_ERA5_LAND_MONTHLY_AGGR ). Precipitation and wind data are calculated from TerraClimate (Monthly Climate and Climatic Water Balance for Global Terrestrial Surfaces, University of Idaho) ( https://developers.google.com/earth-engine/datasets/catalog/IDAHO_EPSCOR_TERRACLIMATE ). Humidity data is calculated from Famine Early Warning Systems Network (FEWS NET) Land Data Assimilation System ( https://developers.google.com/earth-engine/datasets/catalog/NASA_FLDAS_NOAH01_C_GL_M_V001 ). Topography data from MERIT DEM (Multi-Error-Removed Improved-Terrain DEM) product is available at https://developers.google.com/earth-engine/datasets/catalog/MERIT_DEM_v1_0_3 . GDP from Gross Domestic Product and Human Development Index dataset is available at https://doi.org/10.5061/dryad.dk1j0 . VIIRS nighttime light data is available at https://developers.google.com/earth-engine/datasets/catalog/NOAA_VIIRS_DNB_MONTHLY_V1_VCMSLCFG . City building volume data from Global 3D Building Structure (1 km) is available at https://doi.org/10.34894/4QAGYL . Albedo data is derived from the MODIS MCD43A3 product ( https://developers.google.com/earth-engine/datasets/catalog/MODIS_061_MCD43A3 ), and aerosol data is derived from the MODIS MCD19A2 product ( https://developers.google.com/earth-engine/datasets/catalog/MODIS_061_MCD19A2_GRANULES ). All data used for generating the results are publicly available at https://doi.org/10.6084/m9.figshare.26340592.v1 .

Code availability

The codes used for data collection and analyses are publicly available at https://doi.org/10.6084/m9.figshare.26340592.v1 .

Dosio, A., Mentaschi, L., Fischer, E. M. & Wyser, K. Extreme heat waves under 1.5 °C and 2 °C global warming. Environ. Res. Lett. 13 , 054006 (2018).

Article   ADS   Google Scholar  

Suarez-Gutierrez, L., Müller, W. A., Li, C. & Marotzke, J. Hotspots of extreme heat under global warming. Clim. Dyn. 55 , 429–447 (2020).

Article   Google Scholar  

Guo, Y. et al. Global variation in the effects of ambient temperature on mortality: a systematic evaluation. Epidemiology 25 , 781–789 (2014).

Article   PubMed   PubMed Central   Google Scholar  

Mora, C. et al. Global risk of deadly heat. Nat. Clim. Chang. 7 , 501–506 (2017).

Ebi, K. L. et al. Hot weather and heat extremes: health risks. Lancet 398 , 698–708 (2021).

Article   PubMed   Google Scholar  

Lüthi, S. et al. Rapid increase in the risk of heat-related mortality. Nat. Commun. 14 , 4894 (2023).

Article   ADS   PubMed   PubMed Central   Google Scholar  

United Nations Department of Economic Social Affairs, Population Division. in World Population Prospects 2022: Summary of Results (United Nations Fund for Population Activities, 2022).

Sachindra, D., Ng, A., Muthukumaran, S. & Perera, B. Impact of climate change on urban heat island effect and extreme temperatures: a case‐study. Q. J. R. Meteorol. Soc. 142 , 172–186 (2016).

Guo, L. et al. Evaluating contributions of urbanization and global climate change to urban land surface temperature change: a case study in Lagos, Nigeria. Sci. Rep. 12 , 14168 (2022).

Article   ADS   CAS   PubMed   PubMed Central   Google Scholar  

Liu, Z. et al. Surface warming in global cities is substantially more rapid than in rural background areas. Commun. Earth Environ. 3 , 219 (2022).

Mentaschi, L. et al. Global long-term mapping of surface temperature shows intensified intra-city urban heat island extremes. Glob. Environ. Change 72 , 102441 (2022).

Asseng, S., Spänkuch, D., Hernandez-Ochoa, I. M. & Laporta, J. The upper temperature thresholds of life. Lancet Planet. Health 5 , e378–e385 (2021).

Zander, K. K., Botzen, W. J., Oppermann, E., Kjellstrom, T. & Garnett, S. T. Heat stress causes substantial labour productivity loss in Australia. Nat. Clim. Chang. 5 , 647–651 (2015).

Flouris, A. D. et al. Workers’ health and productivity under occupational heat strain: a systematic review and meta-analysis. Lancet Planet. Health 2 , e521–e531 (2018).

Xu, C., Kohler, T. A., Lenton, T. M., Svenning, J.-C. & Scheffer, M. Future of the human climate niche. Proc. Natl Acad. Sci. USA 117 , 11350–11355 (2020).

Lenton, T. M. et al. Quantifying the human cost of global warming. Nat. Sustain. 6 , 1237–1247 (2023).

Harrington, L. J. et al. Poorest countries experience earlier anthropogenic emergence of daily temperature extremes. Environ. Res. Lett. 11 , 055007 (2016).

Bathiany, S., Dakos, V., Scheffer, M. & Lenton, T. M. Climate models predict increasing temperature variability in poor countries. Sci. Adv. 4 , eaar5809 (2018).

Alizadeh, M. R. et al. Increasing heat‐stress inequality in a warming climate. Earth Future 10 , e2021EF002488 (2022).

Tuholske, C. et al. Global urban population exposure to extreme heat. Proc. Natl Acad. Sci. USA 118 , e2024792118 (2021).

Article   CAS   PubMed   PubMed Central   Google Scholar  

Manoli, G. et al. Magnitude of urban heat islands largely explained by climate and population. Nature 573 , 55–60 (2019).

Article   ADS   CAS   PubMed   Google Scholar  

Wang, J. et al. Anthropogenic emissions and urbanization increase risk of compound hot extremes in cities. Nat. Clim. Chang. 11 , 1084–1089 (2021).

Article   ADS   CAS   Google Scholar  

Bowler, D. E., Buyung-Ali, L., Knight, T. M. & Pullin, A. S. Urban greening to cool towns and cities: a systematic review of the empirical evidence. Landsc. Urban Plan. 97 , 147–155 (2010).

Armson, D., Stringer, P. & Ennos, A. The effect of tree shade and grass on surface and globe temperatures in an urban area. Urban For. Urban Green. 11 , 245–255 (2012).

Wang, C., Wang, Z. H. & Yang, J. Cooling effect of urban trees on the built environment of contiguous United States. Earth Future 6 , 1066–1081 (2018).

Pataki, D. E., McCarthy, H. R., Litvak, E. & Pincetl, S. Transpiration of urban forests in the Los Angeles metropolitan area. Ecol. Appl. 21 , 661–677 (2011).

Konarska, J. et al. Transpiration of urban trees and its cooling effect in a high latitude city. Int. J. Biometeorol. 60 , 159–172 (2016).

Article   ADS   PubMed   Google Scholar  

Li, X., Zhou, W., Ouyang, Z., Xu, W. & Zheng, H. Spatial pattern of greenspace affects land surface temperature: evidence from the heavily urbanized Beijing metropolitan area, China. Landsc. Ecol. 27 , 887–898 (2012).

Yu, Z., Xu, S., Zhang, Y., Jørgensen, G. & Vejre, H. Strong contributions of local background climate to the cooling effect of urban green vegetation. Sci. Rep. 8 , 6798 (2018).

Richards, D. R., Fung, T. K., Belcher, R. & Edwards, P. J. Differential air temperature cooling performance of urban vegetation types in the tropics. Urban For. Urban Green. 50 , 126651 (2020).

Winbourne, J. B. et al. Tree transpiration and urban temperatures: current understanding, implications, and future research directions. BioScience 70 , 576–588 (2020).

Schwaab, J. et al. The role of urban trees in reducing land surface temperatures in European cities. Nat. Commun. 12 , 6763 (2021).

Vo, T. T. & Hu, L. Diurnal evolution of urban tree temperature at a city scale. Sci. Rep. 11 , 10491 (2021).

Wang, J. et al. Comparing relationships between urban heat exposure, ecological structure, and socio-economic patterns in Beijing and New York City. Landsc. Urban Plan. 235 , 104750 (2023).

Chen, B. et al. Contrasting inequality in human exposure to greenspace between cities of Global North and Global South. Nat. Commun. 13 , 4636 (2022).

Pavanello, F. et al. Air-conditioning and the adaptation cooling deficit in emerging economies. Nat. Commun. 12 , 6460 (2021).

Turner, V. K., Middel, A. & Vanos, J. K. Shade is an essential solution for hotter cities. Nature 619 , 694–697 (2023).

Hope, D. et al. Socioeconomics drive urban plant diversity. Proc. Natl Acad. Sci. USA 100 , 8788–8792 (2003).

Leong, M., Dunn, R. R. & Trautwein, M. D. Biodiversity and socioeconomics in the city: a review of the luxury effect. Biol. Lett. 14 , 20180082 (2018).

Schwarz, K. et al. Trees grow on money: urban tree canopy cover and environmental justice. PloS ONE 10 , e0122051 (2015).

Chakraborty, T., Hsu, A., Manya, D. & Sheriff, G. Disproportionately higher exposure to urban heat in lower-income neighborhoods: a multi-city perspective. Environ. Res. Lett. 14 , 105003 (2019).

Wang, J. et al. Significant effects of ecological context on urban trees’ cooling efficiency. ISPRS J. Photogramm. Remote Sens. 159 , 78–89 (2020).

Marando, F. et al. Urban heat island mitigation by green infrastructure in European Functional Urban Areas. Sust. Cities Soc. 77 , 103564 (2022).

Cheng, X., Peng, J., Dong, J., Liu, Y. & Wang, Y. Non-linear effects of meteorological variables on cooling efficiency of African urban trees. Environ. Int. 169 , 107489 (2022).

Yang, Q. et al. Global assessment of urban trees’ cooling efficiency based on satellite observations. Environ. Res. Lett. 17 , 034029 (2022).

Yin, Y., He, L., Wennberg, P. O. & Frankenberg, C. Unequal exposure to heatwaves in Los Angeles: Impact of uneven green spaces. Sci. Adv. 9 , eade8501 (2023).

Fantom N., Serajuddin U. The World Bank’s Classification of Countries by Income (The World Bank, 2016).

Iungman, T. et al. Cooling cities through urban green infrastructure: a health impact assessment of European cities. Lancet 401 , 577–589 (2023).

He, C. et al. The inequality labor loss risk from future urban warming and adaptation strategies. Nat. Commun. 13 , 3847 (2022).

Kii, M. Projecting future populations of urban agglomerations around the world and through the 21st century. npj Urban Sustain 1 , 10 (2021).

Paschalis, A., Chakraborty, T., Fatichi, S., Meili, N. & Manoli, G. Urban forests as main regulator of the evaporative cooling effect in cities. AGU Adv. 2 , e2020AV000303 (2021).

Hunte, N., Roopsind, A., Ansari, A. A. & Caughlin, T. T. Colonial history impacts urban tree species distribution in a tropical city. Urban For. Urban Green. 41 , 313–322 (2019).

Kabano, P., Harris, A. & Lindley, S. Sensitivity of canopy phenology to local urban environmental characteristics in a tropical city. Ecosystems 24 , 1110–1124 (2021).

Frank, S. D. & Backe, K. M. Effects of urban heat islands on temperate forest trees and arthropods. Curr. Rep. 9 , 48–57 (2023).

Esperon-Rodriguez, M. et al. Climate change increases global risk to urban forests. Nat. Clim. Chang. 12 , 950–955 (2022).

Stewart, I. D. & Oke, T. R. Local climate zones for urban temperature studies. Bull. Am. Meteorol. Soc. 93 , 1879–1900 (2012).

Biardeau, L. T., Davis, L. W., Gertler, P. & Wolfram, C. Heat exposure and global air conditioning. Nat. Sustain. 3 , 25–28 (2020).

Davis, L., Gertler, P., Jarvis, S. & Wolfram, C. Air conditioning and global inequality. Glob. Environ. Change 69 , 102299 (2021).

Colelli, F. P., Wing, I. S. & Cian, E. D. Air-conditioning adoption and electricity demand highlight climate change mitigation–adaptation tradeoffs. Sci. Rep. 13 , 4413 (2023).

Sun, L., Chen, J., Li, Q. & Huang, D. Dramatic uneven urbanization of large cities throughout the world in recent decades. Nat. Commun. 11 , 5366 (2020).

Liu, D., Kwan, M.-P. & Kan, Z. Analysis of urban green space accessibility and distribution inequity in the City of Chicago. Urban For. Urban Green. 59 , 127029 (2021).

Hsu, A., Sheriff, G., Chakraborty, T. & Manya, D. Disproportionate exposure to urban heat island intensity across major US cities. Nat. Commun. 12 , 2721 (2021).

Zhao, L., Lee, X., Smith, R. B. & Oleson, K. Strong contributions of local background climate to urban heat islands. Nature 511 , 216–219 (2014).

Wu, S., Chen, B., Webster, C., Xu, B. & Gong, P. Improved human greenspace exposure equality during 21st century urbanization. Nat. Commun. 14 , 6460 (2023).

Zhao, J., Zhao, X., Wu, D., Meili, N. & Fatichi, S. Satellite-based evidence highlights a considerable increase of urban tree cooling benefits from 2000 to 2015. Glob. Chang. Biol. 29 , 3085–3097 (2023).

Article   CAS   PubMed   Google Scholar  

Nice, K. A., Coutts, A. M. & Tapper, N. J. Development of the VTUF-3D v1. 0 urban micro-climate model to support assessment of urban vegetation influences on human thermal comfort. Urban Clim. 24 , 1052–1076 (2018).

Meili, N. et al. An urban ecohydrological model to quantify the effect of vegetation on urban climate and hydrology (UT&C v1. 0). Geosci. Model Dev. 13 , 335–362 (2020).

Nesbitt, L., Meitner, M. J., Sheppard, S. R. & Girling, C. The dimensions of urban green equity: a framework for analysis. Urban For. Urban Green. 34 , 240–248 (2018).

Hedblom, M., Prévot, A.-C. & Grégoire, A. Science fiction blockbuster movies—a problem or a path to urban greenery? Urban For. Urban Green. 74 , 127661 (2022).

Norton, B. A. et al. Planning for cooler cities: a framework to prioritise green infrastructure to mitigate high temperatures in urban landscapes. Landsc. Urban Plan 134 , 127–138 (2015).

Medl, A., Stangl, R. & Florineth, F. Vertical greening systems—a review on recent technologies and research advancement. Build. Environ. 125 , 227–239 (2017).

Chen, B., Lin, C., Gong, P. & An, J. Optimize urban shade using digital twins of cities. Nature 622 , 242–242 (2023).

Pamukcu-Albers, P. et al. Building green infrastructure to enhance urban resilience to climate change and pandemics. Landsc. Ecol. 36 , 665–673 (2021).

Haaland, C. & van Den Bosch, C. K. Challenges and strategies for urban green-space planning in cities undergoing densification: a review. Urban For. Urban Green. 14 , 760–771 (2015).

Shafique, M., Kim, R. & Rafiq, M. Green roof benefits, opportunities and challenges—a review. Renew. Sust. Energ. Rev. 90 , 757–773 (2018).

Wang, J., Zhou, W. & Jiao, M. Location matters: planting urban trees in the right places improves cooling. Front. Ecol. Environ. 20 , 147–151 (2022).

Lan, T., Liu, Y., Huang, G., Corcoran, J. & Peng, J. Urban green space and cooling services: opposing changes of integrated accessibility and social equity along with urbanization. Sust. Cities Soc. 84 , 104005 (2022).

Wood, S. & Dupras, J. Increasing functional diversity of the urban canopy for climate resilience: Potential tradeoffs with ecosystem services? Urban For. Urban Green. 58 , 126972 (2021).

Wong, N. H., Tan, C. L., Kolokotsa, D. D. & Takebayashi, H. Greenery as a mitigation and adaptation strategy to urban heat. Nat. Rev. Earth Environ. 2 , 166–181 (2021).

United Nations. Department of economic and social affairs, population division. in The World’s Cities in 2018—Data Booklet (UN, 2018).

United Nations Development Programme (UNDP). Human Development Report 2019: Beyond Income, Beyond Averages, Beyond Today: Inequalities in Human Development in the 21st Century (United Nations Development Programme (UNDP), 2019)

Li, X. et al. Mapping global urban boundaries from the global artificial impervious area (GAIA) data. Environ. Res. Lett. 15 , 094044 (2020).

Stevens, F. R., Gaughan, A. E., Linard, C. & Tatem, A. J. Disaggregating census data for population mapping using random forests with remotely-sensed and ancillary data. PloS ONE 10 , e0107042 (2015).

Buchhorn, M. et al. Copernicus global land cover layers—collection 2. Remote Sens 12 , 1044 (2020).

Gillerot, L. et al. Forest structure and composition alleviate human thermal stress. Glob. Change Biol. 28 , 7340–7352 (2022).

Article   CAS   Google Scholar  

Hamada, S., Tanaka, T. & Ohta, T. Impacts of land use and topography on the cooling effect of green areas on surrounding urban areas. Urban For. Urban Green. 12 , 426–434 (2013).

Sun, X. et al. Quantifying landscape-metrics impacts on urban green-spaces and water-bodies cooling effect: the study of Nanjing, China. Urban For . Urban Green. 55 , 126838 (2020).

Zhang, Q., Zhou, D., Xu, D. & Rogora, A. Correlation between cooling effect of green space and surrounding urban spatial form: Evidence from 36 urban green spaces. Build. Environ. 222 , 109375 (2022).

Pesaresi, M., Politis, P. GHS-BUILT-H R2023A - GHS building height, derived from AW3D30, SRTM30, and Sentinel2 composite (2018) . European Commission, Joint Research Centre (JRC) https://doi.org/10.2905/85005901-3A49-48DD-9D19-6261354F56FE (2023).

Yamazaki, D. et al. A high‐accuracy map of global terrain elevations. Geophys. Res. Lett. 44 , 5844–5853 (2017).

Wessel, P. & Smith, W. H. A global, self‐consistent, hierarchical, high‐resolution shoreline database. J. Geophys. Res. Solid Earth 101 , 8741–8743 (1996).

Ren et al. climatic map studies: a review. Int. J. Climatol. 31 , 2213–2233 (2011).

Zhou, X. et al. Evaluation of urban heat islands using local climate zones and the influence of sea-land breeze. Sust. Cities Soc. 55 , 102060 (2020).

Zhou, W., Huang, G. & Cadenasso, M. L. Does spatial configuration matter? Understanding the effects of land cover pattern on land surface temperature in urban landscapes. Landsc. Urban Plan 102 , 54–63 (2011).

Muñoz Sabater, J. ERA5-Land monthly averaged data from 1981 to present . Copernicus Climate Change Service (C3S) Climate Data Store (CDS) https://doi.org/10.24381/cds.68d2bb30 (2019).

Abatzoglou, J. T., Dobrowski, S. Z., Parks, S. A. & Hegewisch, K. C. TerraClimate, a high-resolution global dataset of monthly climate and climatic water balance from 1958–2015. Sci. Data 5 , 1–12 (2018).

Kummu, M., Taka, M. & Guillaume, J. H. Gridded global datasets for gross domestic product and Human Development Index over 1990–2015. Sci. Data 5 , 1–15 (2018).

Zanaga, D. et al. ESA WorldCover 10 m 2020 v100. https://doi.org/10.5281/zenodo.5571936 (2021).

McNally, A. et al. A land data assimilation system for sub-Saharan Africa food and water security applications. Sci. Data 4 , 1–19 (2017).

Schaaf C., & Wang Z. MODIS/Terra+Aqua BRDF/Albedo Daily L3 Global - 500m V061 . NASA EOSDIS Land Processes Distributed Active Archive Center. https://doi.org/10.5067/MODIS/MCD43A3.061 (2021).

Lyapustin A., & Wang Y. MODIS/Terra+Aqua Land Aerosol Optical Depth Daily L2G Global 1km SIN Grid V061 . NASA EOSDIS Land Processes Distributed Active Archive Center. https://doi.org/10.5067/MODIS/MCD19A2.061 (2022).

Li, M., Wang, Y., Rosier, J. F., Verburg, P. H. & Vliet, J. V. Global maps of 3D built-up patterns for urban morphological analysis. Int. J. Appl. Earth Obs. Geoinf. 114 , 103048 (2022).

Google Scholar  

Elvidge, C. D., Baugh, K., Zhizhin, M., Hsu, F. C. & Ghosh, T. VIIRS night-time lights. Int. J. Remote Sens. 38 , 5860–5879 (2017).

Zhou, W. et al. Urban tree canopy has greater cooling effects in socially vulnerable communities in the US. One Earth 4 , 1764–1775 (2021).

Beck, H. E. et al. Present and future Köppen-Geiger climate classification maps at 1-km resolution. Sci. Data 5 , 1–12 (2018).

R. Core Team. R: A Language and Environment for Statistical Computing . R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/ (2023).

Fox J., & Weisberg S. An R Companion to Applied Regression 3rd edn (Sage, 2019). https://socialsciences.mcmaster.ca/jfox/Books/Companion/ .

Lefcheck, J. S. piecewiseSEM: Piecewise structural equation modelling in r for ecology, evolution, and systematics. Methods Ecol. Evol. 7 , 573–579 (2016).

Zeileis, A. _ineq: Measuring Inequality, Concentration, and Poverty_ . R package version 0.2-13. https://CRAN.R-project.org/package=ineq (2014).

Download references

Acknowledgements

We thank all the data providers. We thank Marten Scheffer for valuable discussion. C.X. is supported by the National Natural Science Foundation of China (Grant No. 32061143014). J.-C.S. was supported by Center for Ecological Dynamics in a Novel Biosphere (ECONOVO), funded by Danish National Research Foundation (grant DNRF173), and his VILLUM Investigator project “Biodiversity Dynamics in a Changing World”, funded by VILLUM FONDEN (grant 16549). W.Z. was supported by the National Science Foundation of China through Grant No. 42225104. T.M.L. and J.F.A. are supported by the Open Society Foundations (OR2021-82956). W.J.R. is supported by the funding received from Roger Worthington.

Author information

Authors and affiliations.

School of Life Sciences, Nanjing University, Nanjing, China

Yuxiang Li, Shuqing N. Teng & Chi Xu

Center for Ecological Dynamics in a Novel Biosphere (ECONOVO), Department of Biology, Aarhus University, Aarhus, Denmark

Jens-Christian Svenning

State Key Laboratory of Urban and Regional Ecology, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing, China

University of Chinese Academy of Sciences, Beijing, China

Beijing Urban Ecosystem Research Station, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing, China

School for Environment and Sustainability, University of Michigan, Ann Arbor, MI, USA

Global Systems Institute, University of Exeter, Exeter, UK

Jesse F. Abrams & Timothy M. Lenton

Department of Forest Ecosystems and Society, Oregon State University, Corvallis, OR, USA

William J. Ripple

Department of Environmental Science and Engineering, Fudan University, Shanghai, China

Department of Applied Ecology, North Carolina State University, Raleigh, NC, USA

Robert R. Dunn

You can also search for this author in PubMed   Google Scholar

Contributions

Y.L., S.N.T., R.R.D., and C.X. designed the study. Y.L. collected the data, generated the code, performed the analyses, and produced the figures with inputs from J.-C.S., W.Z., K.Z., J.F.A., T.M.L., W.J.R., Z.Y., S.N.T., R.R.D. and C.X. Y.L., S.N.T., R.R.D. and C.X. wrote the first draft with inputs from J.-C.S., W.Z., K.Z., J.F.A., T.M.L., W.J.R., and Z.Y. All coauthors interpreted the results and revised the manuscript.

Corresponding authors

Correspondence to Shuqing N. Teng , Robert R. Dunn or Chi Xu .

Ethics declarations

Competing interests.

The authors declare no competing interests.

Peer review

Peer review information.

Nature Communications thanks Chris Webster and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. A peer review file is available.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Supplementary information, peer review file, description of additional supplementary files, supplementary data 1, reporting summary, rights and permissions.

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/ .

Reprints and permissions

About this article

Cite this article.

Li, Y., Svenning, JC., Zhou, W. et al. Green spaces provide substantial but unequal urban cooling globally. Nat Commun 15 , 7108 (2024). https://doi.org/10.1038/s41467-024-51355-0

Download citation

Received : 06 December 2023

Accepted : 05 August 2024

Published : 02 September 2024

DOI : https://doi.org/10.1038/s41467-024-51355-0

Share this article

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines . If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Quick links

• Explore articles by subject
• Guide to authors
• Editorial policies

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.