- Pre-registration nursing students
- No definition of master’s degree in nursing described in the publication
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).
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.
Author | Contandriopoulos et al | Flinter | Hogan et al | Hungerford et al | O’Rourke | Roots and MacDonald | Schadewaldt et al | Strachan et al |
---|---|---|---|---|---|---|---|---|
Country | Canada | The United States | The United States | Australia | Canada | Canada | Australia | Scotland |
How or why research question | No information on the research question | Several how or why research questions | What and how research question | No information on the research question | Several how or why research questions | No information on the research question | What research question | What and why research questions |
Design and referenced author of methodological guidance | Six qualitative case studies Robert K. Yin | Multiple-case studies design Robert K. Yin | Multiple-case studies design Robert E. Stake | Case study design Robert K. Yin | Qualitative single-case study Robert K. Yin Robert E. Stake Sharan Merriam | Single-case study design Robert K. Yin Sharan Merriam | Multiple-case studies design Robert K. Yin Robert E. Stake | Multiple-case studies design |
Case definition | Team of health professionals (Small group) | Nurse practitioners (Individuals) | Primary care practices (Organization) | Community-based NP model of practice (Organization) | NP-led practice (Organization) | Primary care practices (Organization) | No information on case definition | Health board (Organization) |
Overview of Within-Method, Between/Across-Method, and Data-Analysis Triangulation.
Author | Contandriopoulos et al | Flinter | Hogan et al | Hungerford et al | O’Rourke | Roots and MacDonald | Schadewaldt et al | Strachan et al |
---|---|---|---|---|---|---|---|---|
Within-method triangulation (using within-method triangulation use at least 2 data-collection procedures from the same design approach) | ||||||||
: | ||||||||
Interviews | X | x | x | x | x | |||
Observations | x | x | ||||||
Public documents | x | x | x | |||||
Electronic health records | x | |||||||
Between/across-method (using both qualitative and quantitative data-collection procedures in the same study) | ||||||||
: | ||||||||
: | ||||||||
Interviews | x | x | x | |||||
Observations | x | x | ||||||
Public documents | x | x | ||||||
Electronic health records | x | |||||||
: | ||||||||
Self-assessment | x | |||||||
Service records | x | |||||||
Questionnaires | x | |||||||
Data-analysis triangulation (combination of 2 or more methods of analyzing data) | ||||||||
: | ||||||||
: | ||||||||
Deductive | x | x | x | |||||
Inductive | x | x | ||||||
Thematic | x | x | ||||||
Content | ||||||||
: | ||||||||
Descriptive analysis | x | x | x | |||||
: | ||||||||
: | ||||||||
Deductive | x | x | x | x | ||||
Inductive | x | x | ||||||
Thematic | x | |||||||
Content | x |
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
“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.
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.
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.
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
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.
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
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).
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
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.
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.
In 1 study, qualitative individual journals were studied. 16 These included reflective journals from NPs, who performed the role in primary health care.
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.
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).
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.
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.
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.
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.
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
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.
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:
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.
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
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, 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.
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
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.
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.
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).
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.
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.
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.
Simply put, the case method is a discussion of real-life situations that business executives have faced.
On average, you'll attend three to four different classes a day, for a total of about six hours of class time (schedules vary). To prepare, you'll work through problems with your peers.
Often, executives are surprised to discover that the objective of the case study is not to reach consensus, but to understand how different people use the same information to arrive at diverse conclusions. When you begin to understand the context, you can appreciate the reasons why those decisions were made. You can prepare for case discussions in several ways.
In self-reflection.
The time you spend here is deeply introspective. You're not only working with case materials and assignments, but also taking on the role of the case protagonist—the person who's supposed to make those tough decisions. How would you react in those situations? We put people in a variety of contexts, and they start by addressing that specific problem.
The discussion group is a critical component of the HBS experience. You're working in close quarters with a group of seven or eight very accomplished peers in diverse functions, industries, and geographies. Because they bring unique experience to play you begin to see that there are many different ways to wrestle with a problem—and that’s very enriching.
The faculty guides you in examining and resolving the issues—but the beauty here is that they don't provide you with the answers. You're interacting in the classroom with other executives—debating the issue, presenting new viewpoints, countering positions, and building on one another's ideas. And that leads to the next stage of learning.
Once you leave the classroom, the learning continues and amplifies as you get to know people in different settings—over meals, at social gatherings, in the fitness center, or as you are walking to class. You begin to distill the takeaways that you want to bring back and apply in your organization to ensure that the decisions you make will create more value for your firm.
Pioneered by HBS faculty, the case method puts you in the role of the chief decision maker as you explore the challenges facing leading companies across the globe. Learning to think fast on your feet with limited information sharpens your analytical skills and empowers you to make critical decisions in real time.
To get the most out of each case, it's important to read and reflect, and then meet with your discussion group to share your insights. You and your peers will explore the underlying issues, compare alternatives, and suggest various ways of resolving the problem.
There's more than one way to prepare for a case discussion, but these general guidelines can help you develop a method that works for you.
Read the professor's assignment or discussion questions.
The assignment and discussion questions help you focus on the key aspects of the case. Ask yourself: What are the most important issues being raised?
Each case begins with a text description followed by exhibits. Ask yourself: What is the case generally about, and what information do I need to analyze?
Put yourself in the shoes of the case protagonist, and own that person's problems. Ask yourself: What basic problem is this executive trying to resolve?
Sort out relevant considerations and do the quantitative or qualitative analysis. Ask yourself: What recommendations should I make based on my case data analysis?
The key to being an active listener and participant in case discussions—and to getting the most out of the learning experience—is thorough individual preparation.
We've set aside formal time for you to discuss the case with your group. These sessions will help you to become more confident about sharing your views in the classroom discussion.
Actively express your views and challenge others. Don't be afraid to share related "war stories" that will heighten the relevance and enrich the discussion.
If the content doesn't seem to relate to your business, don't tune out. You can learn a lot about marketing insurance from a case on marketing razor blades!
Actively apply what you're learning to your own specific management situations, both past and future. This will magnify the relevance to your business.
People with diverse backgrounds, experiences, skills, and styles will take away different things. Be sure to note what resonates with you, not your peers.
Being exposed to so many different approaches to a given situation will put you in a better position to enhance your management style.
What can i expect on the first day, what happens in class if nobody talks, does everyone take part in "role-playing".
You are accessing a machine-readable page. In order to be human-readable, please install an RSS reader.
All articles published by MDPI are made immediately available worldwide under an open access license. No special permission is required to reuse all or part of the article published by MDPI, including figures and tables. For articles published under an open access Creative Common CC BY license, any part of the article may be reused without permission provided that the original article is clearly cited. For more information, please refer to https://www.mdpi.com/openaccess .
Feature papers represent the most advanced research with significant potential for high impact in the field. A Feature Paper should be a substantial original Article that involves several techniques or approaches, provides an outlook for future research directions and describes possible research applications.
Feature papers are submitted upon individual invitation or recommendation by the scientific editors and must receive positive feedback from the reviewers.
Editor’s Choice articles are based on recommendations by the scientific editors of MDPI journals from around the world. Editors select a small number of articles recently published in the journal that they believe will be particularly interesting to readers, or important in the respective research area. The aim is to provide a snapshot of some of the most exciting work published in the various research areas of the journal.
Original Submission Date Received: .
Find support for a specific problem in the support section of our website.
Please let us know what you think of our products and services.
Visit our dedicated information section to learn more about MDPI.
Teach programming using task-driven case studies: pedagogical approach, guidelines, and implementation.
1.1. task-driven teaching, 1.2. case studies, 1.3. task-driven case studies, 2. background, 2.1. problem-based learning, 2.2. case studies in teaching, 2.3. task-driven teaching methodology, 2.4. games in teaching and automatic feedback, 3. task-driven case studies’ pedagogical framework.
3.4. task-driven case study lifecycle, 3.5. guidelines.
4.1. know the learning goals.
5.1. implement code first.
Click here to enlarge figure
6. writing study guide.
6.10. review the guide, 6.11. start with a minimal version and improve it, 7. course execution, 7.1. monitor progress continuously, 7.2. favor individual achievements.
8.1. teachers’ view, 8.2. students’ view, 8.2.1. objective, 8.2.2. method, 8.2.3. results, 9. potential drawbacks, 10. conclusions, author contributions, institutional review board statement, informed consent statement, data availability statement, conflicts of interest.
# | Question | Answer | No. | % |
---|---|---|---|---|
1 | Do you like learning by implementing a game? | 103 | 92 | |
No, I would rather implement something else | 9 | 8 | ||
2 | How do you like the implemented game? | Poor (is it still a game?) | 0 | 0 |
Below average (I would never play it) | 18 | 16 | ||
56 | 50 | |||
Above average (I would definitely try it) | 36 | 32 | ||
Excellent (from now on, I will play only this game) | 2 | 2 | ||
3 | Would you show your game to your friend or a family member? | 87 | 78 | |
No | 25 | 22 | ||
4 | In practical lessons, you prefer to implement: | 84 | 75 | |
Multiple simple independent tasks | 28 | 25 | ||
5 | Do you think you understood programming principles better by implementing one large project? | 91 | 81 | |
No | 21 | 19 | ||
6 | From the point of view of assignment organization, you prefer: | 97 | 87 | |
To get the assignment in the beginning and to solve it on my own | 15 | 13 | ||
7 | Did you have problems with dependencies between lessons, i.e., that you had to solve a previous lesson to be able to continue? | Yes | 48 | 43 |
64 | 57 | |||
8 | Was the difficulty of the tasks in the case study balanced? | Yes | 48 | 43 |
64 | 57 | |||
9 | Were the tasks too easy? | Yes, most of the time, I just needed to repeat what was written | 3 | 3 |
109 | 97 | |||
10 | Were the tasks described clearly enough? | 57 | 51 | |
No, often I had to ask the teacher or colleagues for help | 55 | 49 | ||
11 | Would you like learning with study guides for a case study in future courses? * | 103 | 92 | |
No | 8 | 7 | ||
12 | Do you think that working with a study guide: | Limits me because I cannot do what I want | 22 | 20 |
90 | 80 | |||
13 | When did you implement the tasks for a given lesson? | I programmed mostly before the lesson | 12 | 11 |
I programmed during the lesson | 24 | 21 | ||
76 | 68 | |||
14 | Which properties of studying with study guides for a case study do you consider most important (choose max. 3)? | I implemented a large project | 69 | 62 |
I implemented a game | 34 | 30 | ||
81 | 72 | |||
I had a study guide that lad me to good practices | 55 | 49 | ||
I worked incrementally, but the game was always playable | 42 | 38 | ||
15 | What did you like about practical lessons in the OOP course? | |||
16 | What did you dislike about practical lessons in the OOP course? | |||
17 | What would you change or improve about practical lessons in the OOP course? |
The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
Porubän, J.; Nosál’, M.; Sulír, M.; Chodarev, S. Teach Programming Using Task-Driven Case Studies: Pedagogical Approach, Guidelines, and Implementation. Computers 2024 , 13 , 221. https://doi.org/10.3390/computers13090221
Porubän J, Nosál’ M, Sulír M, Chodarev S. Teach Programming Using Task-Driven Case Studies: Pedagogical Approach, Guidelines, and Implementation. Computers . 2024; 13(9):221. https://doi.org/10.3390/computers13090221
Porubän, Jaroslav, Milan Nosál’, Matúš Sulír, and Sergej Chodarev. 2024. "Teach Programming Using Task-Driven Case Studies: Pedagogical Approach, Guidelines, and Implementation" Computers 13, no. 9: 221. https://doi.org/10.3390/computers13090221
Further information, mdpi initiatives, follow mdpi.
Subscribe to receive issue release notifications and newsletters from MDPI journals
The objective of this paper is to estimate the equivalent permeability of the rock surrounding the tailrace tunnel of the Azad Dam pumped storage power plant, using geostatistical methods. The permeability of the rock mass is a critical factor that influences the estimation of water flow rates. Since the tunnel passes through various geological units with different permeabilities, it is crucial to estimate the equivalent permeability for each unit in order to predict the water seepage from that unit into the tunnel. In order to estimate the permeability along the tunnel and underground structures, twelve exploratory boreholes were drilled, and water pressure tests were conducted. Due to the distribution of the exploratory boreholes, a study and statistical analysis are necessary to determine the permeability of the rock mass for each geological unit. Using geostatistical log kriging, multiple indicator kriging with four thresholds, and multiple indicator kriging with five thresholds, the permeability of the rock mass at the tunnel route was estimated. The results indicate that at least 40% of the rock mass has low permeability, while the remaining mass of the tunnel passes through rocks with moderate to high permeability. The accuracy of the estimated permeability values was evaluated by predicting the water inflow into the tunnel using the estimated permeability values and comparing it to the observed flow. Numerical models were generated for each geological unit to estimate the water inflow into the tunnel, based on the results of the geostatistical methods. Log kriging, multiple indicator kriging with four thresholds, and multiple indicator kriging with five thresholds were used to calculate the water inflow, resulting in 94.15, 94.15, and 127.5 L per second, respectively. The results of the modeling were compared to the observed water flow into the tunnel. Comparing the modeling results to both the statistical methods and observed values showed errors of 31.2%, 31.2%, and 6.9%, respectively. Of the three methods, the multiple indicator kriging computational method with five thresholds was found to be the most accurate, with the least amount of error and the closest approximation to the actual value. As a result, it was selected as the best method.
This is a preview of subscription content, log in via an institution to check access.
Subscribe and save.
Price includes VAT (Russian Federation)
Instant access to the full article PDF.
Rent this article via DeepDyve
Institutional subscriptions
No datasets were generated or analysed during the current study.
Anvari AA, Katibeh H, Sharifzade M (2010) A new approach for computing permeability of fault zones case study: the upper reservoir of Azad pumped-storage power station in Iran. Arch Min Sci 55(3):605–621
Google Scholar
Atilla Öztürk C, Nasuf E (2002) Geostatistical assessment of rock zones for tunneling. Tunn Undergr Space Technol 17(3):275–285. https://doi.org/10.1016/S0886-7798(02)00023-8
Article Google Scholar
Bastante FG et al (2005) Evaluation of the resources of a slate deposit using indicator kriging. Eng Geol 81(4):407–418
Doostmohammadi M, Jafari A, Asghari O (2014) Geostatistical modeling of uniaxial compressive strength along the axis of the Behesht-Abad tunnel in Central Iran. Bull Eng Geol Environ 74:789–802
Ewert F-K (2015) Evaluation and interpretation of water pressure tests. In: Bell AL (ed) Grouting in the ground: Proceedings of the conference organized by the Institution of Civil Engineers and held in London on 25–26 November 1992. Thomas Telford Publishing, pp 141–162. https://doi.org/10.1680/gitg.19287.0010
Chapter Google Scholar
Farhadian H, Aalianvari A, Katibeh H (2012) Optimization of analytical equations of groundwater seepage into tunnels: A case study of Amirkabir tunnel. J Geol Soc India 80(1):96–100
Ferrari F, Apuani T, Giani GP (2014) Rock Mass Rating spatial estimation by geostatistical analysis. Int J Rock Mech Min Sci 70:162–176
Goovaerts P (1997) Geostatistics for natural resources evaluation. Oxford University Press, New York
Book Google Scholar
HÉBERT, Étienne, Jean CÔTÉ, and Marc SMITH (2016) A geostatistical analysis of variations of permeability within a compacted dam core. Proceedings of the 17th Nordic Geotechnical Meeting, Reykjavik, Islande
Hezarkhani A (2006) Mineralogy and fluid inclusion investigations in the Reagan Porphyry System, Iran, the path to an uneconomic porphyry copper deposit. J Asian Earth Sci 27(5):598–612
Isaaks EH, Srivastava RM (1989) Applied Geostatistics. Oxford University Press, United Kingdom
Kasmaee S, Gholamnejad J, Yarahmadi A, Mojtahedzadeh H (2010) Reserve estimation of the high phosphorous stockpile at the Choghart iron mine of Iran using geostatistical modeling. MinSci Tech 20:0855–0860
Katibeh H, Aalianvari A (2012) Common Approximations to the water inflow into Tunnels. In: Javaid MS (ed) Drainage systems. InTech. https://doi.org/10.5772/36003
Kim Y, Moon J-S (2020) Change of groundwater inflow by cutoff grouting thickness and permeability coefficient. Geomech Eng 21(2):165
Lee H, Tae-Min Oh, Park C (2020) Analysis of permeability in rock fracture with effective stress at deep depth. Geomech Eng 22(5):2020
Liu C-W, Jang ChSh, Liao ChM (2004) Evaluation of arsenic contamination potential using indicator kriging in the Yun-Lin aquifer (Taiwan). Sci Total Environ 321:173–188
Pinheiro M et al (2016) Geostatistical simulation to map the spatial heterogeneity of geomechanical parameters: a case study with rock mass rating. Eng Geol 205(2016):93–103
Rezaei M, Habibi H (2023) Designing of the Beheshtabad water transmission tunnel based on the hybrid empirical method. Struct Eng Mech 86(5):621–633
Rezaei M, Latifi S (2018) Development of a new empirical relation between RMR and Q rock mass classification systems. Iranian J Eng Geol, 12(2)
Rezaei M, Davoodi PK, Najmoddini I (2019) Studying the correlation of rock properties with P-wave velocity index in dry and saturated conditions. J Appl Geophys 169:49–57
Sinha RS (1989) Underground structure desing and insprumentation, Elsevier, Technology & Engineering pp 480
Soltani-Mohammadi S, Erhan Tercan A (2012) Constrained multiple indicator kriging using sequential quadratic programming. Comput Geosci 48:211–219
Wang Y, Rezaei M, Abdullah RA, Hasanipanah M (2023) Developing two hybrid algorithms for predicting the elastic modulus of intact rocks. Sustainability 15(5):4230
Wei L, Zou J-F, An W (2020) “Solution for a circular tunnel in strain-softening rock with seepage forces. Geomech Eng 22(6):2020
Zou J-F, Wei An, Liang Li (2020) Analytical solution for steady seepage and groundwater inflow into an underwater tunnel. Geomech Eng 20(3):267
Download references
Author information, authors and affiliations.
Department of Mining Engineering, Faculty of Engineering, University of Kashan, Kashan, Islamic Republic of Iran
Sanaz Khoubani, Ali Aalianvari & Saeed Soltani-Mohammadi
You can also search for this author in PubMed Google Scholar
All of the authors have equal contributions.
Correspondence to Ali Aalianvari .
Conflict of interest.
The authors declare no competing interests.
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
Reprints and permissions
Khoubani, S., Aalianvari, A. & Soltani-Mohammadi, S. Estimation of Rock Mass Equivalent Permeability Around Tunnel Route Using the Geostatistical Methods: A Case Study. Iran J Sci Technol Trans Civ Eng (2024). https://doi.org/10.1007/s40996-024-01608-1
Download citation
Received : 27 April 2024
Accepted : 23 August 2024
Published : 05 September 2024
DOI : https://doi.org/10.1007/s40996-024-01608-1
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
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.
Scientific Reports volume 14 , Article number: 20483 ( 2024 ) Cite this article
Metrics details
In the lower atmosphere, CO 2 emissions impact human health and ecosystems, making data at this level essential for addressing carbon-cycle and public-health questions. The atmospheric concentration of CO 2 is crucial in urban areas due to its connection with air quality, pollution, and climate change, becoming a pivotal parameter for environmental management and public safety. In volcanic zones, geogenic CO 2 profoundly affects the environment, although hydrocarbon combustion is the primary driver of increased atmospheric CO 2 and global warming. Distinguishing geogenic from anthropogenic emissions is challenging, especially through air CO 2 concentration measurements alone. This study presents survey results on the stable isotope composition of carbon and oxygen in CO 2 and airborne CO 2 concentration in Naples’ urban area, including the Campi Flegrei caldera, a widespread hydrothermal/volcanic zone in the metropolitan area. Over the past 50 years, two major volcanic unrests (1969–72 and 1982–84) were monitored using seismic, deformation, and geochemical data. Since 2005, this area has experienced ongoing unrest, involving the pressurization of the underlying hydrothermal system as a causal factor of the current uplift in the Pozzuoli area and the increased CO 2 emissions in the atmosphere. To better understand CO 2 emission dynamics and to quantify its volcanic origin a mobile laboratory was used. Results show that CO 2 levels in Naples’ urban area exceed background atmospheric levels, indicating an anthropogenic origin from fossil fuel combustion. Conversely, in Pozzuoli's urban area, the stable isotope composition reveals a volcanic origin of the airborne CO 2 . This study emphasizes the importance of monitoring stable isotopes of atmospheric CO 2 , especially in volcanic areas, contributing valuable insights for environmental and public health management.
Introduction.
The equilibrium among natural CO 2 emissions, biotic uptake on land, and ocean absorption regulates long-term fluctuations in airborne CO 2 , establishing the greenhouse effect essential for the biosphere's existence on Earth. Human activities, particularly fossil fuel combustion, vehicle mobility, house heating, and waste management, disrupt the carbon cycle, leading to an increase in airborne CO 2 levels 1 , 2 , 3 , 4 , 5 . Disruption of this equilibrium worsens the effects of global warming and climate changes.
Global temperature data from Copernicus ( https://climate.copernicus.eu/ accessed on 2024, January 10), shows that the mean near-surface temperature in 2023 was ~ 1.4 ± 0.12 °C above the 1850–1900 average. This marked the warmest year in the 174-year observational record, surpassing the joint warmest years of 2016 and 2020. Notably, the last decade (2014–2023) encompasses the nine warmest years on record. Real-time data from specific locations reveals a continued increase in CO 2 levels in 2023, while consolidated concentration datasets of CO 2 , methane, and nitrous oxide reached their highest records in 2022.
Several causes contribute to global warming and climate change 6 . Since the eighteenth century the industrialization has led to the gradual abandonment of rural areas and the concentration of people in urbanized zones. Industries, mainly relying on electrical power generated by hydrocarbon combustion, settled in suburban areas contribute significantly to CO 2 emissions 5 , 7 . Urban growth, characterized by skyscrapers and increased vehicle mobility, results in continuous large-scale carbon dioxide release, predominantly concentrated in urban areas, significantly impacting the global atmospheric composition.
Earth degassing, driven by natural sources like soil respiration, volcanic degassing, and photosynthesis, contributes to atmospheric CO 2 concentrations 8 . Regions of active volcanism, responsible for a significant portion of natural gas emissions, release CO 2 of magmatic origin, particularly during eruptions, accounting for ~ 1% of global CO 2 emissions annually 9 , 10 , 11 . Although this percentage is modest on a global scale, locally, natural emissions may have a more substantial environmental impact, raising hazards for local populations 12 , 13 , 14 , 15 . For example, during the recent outgassing crisis at Vulcano, Italy 16 , 17 , gas hazards increased due to either diffuse degassing or crater plume emissions, though human health risk threshold value was not exceeded 18 , 19 , 20 .
Naples, with around 1 million residents, ranks third in population among Italian cities and is the most densely populated city in Europe. Its strategic location in Mediterranean shipping routes and heavy ship traffic in the harbour make it a potential major source of anthropogenic CO 2 . The city is located in a volcanic area with active volcanic and hydrothermal zones, making it an ideal study area to investigate the coexistence of human-related and natural CO 2 emissions.
This study presents the results of a spatial survey on airborne CO 2 in the metropolitan area of Naples. The survey aimed to collect measurements of airborne CO 2 concentration and stable isotopes of CO 2 to differentiate between volcanic and anthropogenic sources, identifying sources that elevate airborne CO 2 concentrations above the background. The study area includes Naples’ downtown and a broad urbanized zone extending from the western edge of Vesuvius volcano to Bacoli and Cuma in the east, and Agnano crater in the north, encompassing the active volcanic/hydrothermal zone of Campi Flegrei (Fig. 1 a). The Campi Flegrei area has experienced significant volcanic activity, including supereruptions, the oldest one dating back 40,000 years 21 , 22 . This area exhibits continuous degassing and seismic activity (i.e., Solfatara and Pisciarelli in the municipality of Pozzuoli). Anomalies in CO 2 emissions occur from soils via diffuse degassing and from fumaroles 23 , 24 , 25 , 26 , 27 , particularly in the Solfatara area (Zone A in Fig. 1 a). The most recent eruption originated from Monte Nuovo in 1538 A.D. Since then, this system has been in a state of persistent degassing and fluctuating seismic activity, leading to ground motion known as bradiseism. The study area has also increased the degassing since 2005 and is currently in unrest 28 , 29 , 30 , 31 , 32 , 33 , 34 , 35 , 36 . Human-related and geologic CO 2 emissions have distinct stable isotopic signatures, allowing differentiation in the air at the district scale through a combination of concentrations and isotopic measurements 12 , 13 , 18 , 19 , 20 , 37 , 38 , 39 . The results of the spatial survey enable a comparison between volcanic CO 2 emissions and those of anthropogenic origin in the urbanized area of Naples.
Study area, the route used during survey and dataset distribution. The survey was conducted in May 2023. ( a ) The blue line represents the route used during the survey. The selected subsets for Solfatara area (orange zone A), downtown Naples (green zone B), and airport (ice blue zone C) are shown. ( b ) Probability plot for concentration dataset. Global average value of airborne CO 2 concentration is reported as reference. (blue line indicates 423 ppm vol) for a comparison with the average CO 2 concentration over the target area (50% cumulative probability). ( c ) Histograms for both the oxygen isotope (δ 18 O–CO 2 ) and carbon isotope (δ 13 C–CO 2 ) compositions. ( d ) Three-dimensional view of the study area showing the atmospheric CO 2 concentration measurements at their respective locations. The height of the vertical bars is proportional to the concentration levels. The colour scale and bar height indicate that the highest CO 2 concentration was detected near the port. ( e ) Placement of the measurements within the study area. The colour scale is identical to that in subplot ( d ) and indicates the CO 2 concentration measured in the air. The maps ( a ), ( d ), and ( e ) were generated in Qgis 3.34 environment ( https://qgis.org/download/ ).
We developed a measurement program to detect and quantify the spatial variability of CO 2 concentration and its stable isotopes in the near-surface air of the Naples metropolitan area (Fig. 1 a). The dataset enables a better determination of the influence of meteorological factors and multiple greenhouse gas sources on the nature of the urban CO 2 dome 40 , 41 , 42 , 43 , which is considerably more challenging to identify than its mere presence. For this study, the wind direction was selected as the meteorological factor influencing CO 2 dispersal, while other meteorological factors (e.g., temperature, atmospheric pressure, relative humidity) can be averaged over the survey's completion time (11.4 h of acquisition during daytime over 24 h), as variations at the meteorological station from National Research Counsil (i.e., C.N.R. Long: 432,409; Lat: 4,520,399 UTM) are likely suitable for the entire Naples metropolitan area. Throughout the survey period, the weather remained consistently sunny. Table 1 presents the statistics of both environmental variables and atmospheric measurements.
Figure 1 b–c shows statistical distributions of measurements collected during the survey. The dataset collected in the Naples metropolitan area shows airborne CO 2 concentrations higher than 423 ppm vol (Fig. 1 b), which is the global reference for airborne CO 2 concentration for May 20, 2023 ( https://www.climate.gov/climatedashboard accessed on July 2, 2024). The probability plot 44 reveals three independent subsets of CO 2 concentrations. The 50% cumulative distribution indicates that the average value for the background CO 2 concentration in the urban area of Naples is 448.1 ± 1.0 ppm vol. The background population comprises more than 98.9% of the cumulative dataset, while the anomalous subset constitutes less than 0.1% of the cumulative dataset, with CO 2 concentrations exceeding 1300 ppm vol (Fig. 1 b).
Regarding stable isotopes, the carbon isotope composition of airborne CO 2 (reported in delta notation δ 13 C–CO 2 against the Vienna Pee Dee isotopic ratio-VPDB) shows values more 13 C-depleted than the theoretical background air (δ 13 C–CO 2 = − 8‰ vs VPDB). This result indicates that a source of CO 2 forces airborne CO 2 concentration above background values. This gas source has a 13 C-depleted isotopic signature and establishes an urban CO 2 dome in the Naples metropolitan area. Furthermore, the statistical parameters of the data distribution (skewness = − 2.27, kurtosis = 15.74) indicate that the dataset has a peak at δ 13 C–CO 2 = − 10.40‰, which is more 13 C-depleted than the theoretical atmospheric CO 2 value 45 .
The range of values for δ 18 O-CO 2 is wide compared to the spatial and temporal scales of the collected measurements. The oxygen isotope composition of airborne CO 2 depends on both the hydrology of the region and oxygen isotope fractionation in plant leaves during photosynthesis 46 , 47 , 48 . These factors change over spatial and temporal scales different from those of the measurement acquisition (i.e., ~ 10 4 –10 5 m and approximately 24 h, respectively). The oxygen isotope values are almost normally distributed (skewness = − 1.50, kurtosis = 7.3) throughout the study area (Fig. 2 b). Gaussian fitting of the oxygen values has a peak at δ 18 O–CO 2 = − 3.16‰ versus VPDB, which is more 18 O-depleted than the expected value for a coastal area of the Mediterranean region 49 , 50 , 51 , 52 .
Spatial variations of the CO 2 measurements collected during survey throughout the target area (cell size 10 m). The maps were generated in Qgis 3.34 environment ( https://qgis.org/download/ ) using Measurement interpolation generated in SAGA GIS environment ( https://saga-gis.sourceforge.io/en/index.html ). ( a ) Spatial variation of the airborne CO 2 concentration. ( b ) Spatial variation of the oxygen isotope composition of the airborne CO 2 . ( c ) Spatial variation of the carbon isotope composition of the airborne CO 2 . Traces of the concentration profiles are reported (black lines). See text for description.
The collected dataset was utilized to investigate the spatial variation of airborne CO 2 (Fig. 2 ). These data allow investigation of whether urban CO 2 sources affect atmospheric chemistry at a district scale or over the urban area (i.e., at the local scale, ~ 10 4 –10 5 m). The results illustrate a heterogeneous distribution of airborne CO 2 concentration over the Naples metropolitan area, with a concentration gradient from the coast to the inland, likely influenced by local atmospheric circulation. Granieri et al. 14 , who conducted detailed micrometeorological studies on atmospheric circulation in the Naples area for gas dispersal simulations, noted a diurnal sea breeze blowing from SW to NE, pushing clean air inland from the seaside during morning hours. A supply of clean air from the sea would dilute CO 2 concentration at relatively low levels. However, measured CO 2 concentrations in the urban area of Naples suggest that atmospheric circulation is insufficient to reduce atmospheric CO 2 concentrations to background levels, at least on days with similar weather conditions to the ones of the day of measurements. Further research should address the issue concerning the critical atmospheric circulation conditions that help to reduce the concentration level of CO 2 . The implementation of atmospheric CO 2 monitoring programs in urban areas, particularly when integrated with stable isotope composition analyses, is posited as an effective method for detecting anthropogenic or natural forcings influencing atmospheric CO 2 levels. Elevated atmospheric CO 2 concentrations are frequently correlated with increased levels of other pollutants, suggesting that these monitoring programs can significantly enhance public health management strategies. Additionally, in urbanized regions located within volcanic zones, atmospheric CO 2 monitoring is crucial for mitigating volcanic risks associated with gas emissions (i.e., the gas hazard) 19 . Examination of the dataset reveals areas with high airborne CO 2 concentrations, notably near Naples' harbour, where the highest CO 2 concentration was measured (Fig. 1 d,e), and the district of Museo square, among others (Zone B in Fig. 1 a).
Figure 2 illustrates additional zones with high concentrations of airborne CO 2 . The airborne CO 2 concentrations achieve 572 ppm vol in a zone situated in the northeastern sector of the investigated area. While this level does not surpass any established risk threshold for human health 53 , it exceeds the reference value recorded at NOAA Global Monitoring Laboratory for the investigated time frame (424 ppm by volume) by > 33%. A land use survey in the metropolitan area of Naples reveals the presence of the airport, particularly the runways and aircraft parking areas adjacent to the route used for data collection. Another zone exhibiting elevated concentrations of airborne CO 2 is identified on the western side of downtown Naples, within the municipality of Pozzuoli (i.e., transect B–B′ in Fig. 3 b). This area is renowned for its evidence of the underlying volcanic hydrothermal system of Campi Flegrei 26 , 27 , 28 , 29 , with airborne CO 2 concentrations reaching 567 ppm vol. The spatial distribution of airborne CO 2 concentrations in this zone appears more heterogeneous compared to other areas, attributable to the presence of several high concentration nuclei near Bagnoli and Baia (Figs. 2 a and 3 a), eastward and westward of Solfatara, respectively.
Transects through selected zones of the study area to inspect lateral variations of airborne CO 2 concentration (blue line), δ 18 O–CO 2 (red line), and δ 13 C-CO 2 (blue line). ( a ) A–A′ transect (Bacoli). ( b ) B–B′ transect (Solfatara). ( c ) C–C′ transect (Downtown). ( d ) D–D′ transect (Portici).
The δ 18 O–CO 2 has been recognized as a tracer of photosynthesis and the hydrologic cycle's effects on airborne CO 2 . These processes play a pivotal role in the fractionation of oxygen in airborne CO 2 at vastly different spatiotemporal scales. While the hydrologic cycle exhibits seasonal effects at the regional scale, notable changes in vegetation (e.g., transition from C3 to C4 or CAM plant dominant types) account for variations in the oxygen isotope composition due to differences in photosynthesis. Since the survey was completed in a few hours, the spatial variations in the oxygen isotope composition resulting from these processes are expected to have negligible effects on the spatial variations of δ 18 O–CO 2 , which constitutes an ancillary factor for identifying variations in the source of CO 2 at the district scale 12 , 13 , 18 , 20 , 51 , 54 .
The kriging interpolation of the δ 18 O-CO 2 dataset reveals a zone with slightly 18 O-depleted airborne CO 2 westward of downtown Naples, where the δ 18 O–CO 2 = ~ − 2‰. Near Baia, where high concentrations of CO 2 were measured (Fig. 2 a), the airborne CO 2 exhibits more 18 O-depleted values, reaching δ 18 O–CO 2 = − 5.38‰ through a steep isotopic gradient (e.g., transect A–A′ in Fig. 3 a). The δ 18 O–CO 2 abruptly increases to approximately − 2‰ northwestwardly along the transect A–A′ (Fig. 3 a). Airborne CO 2 shows less 18 O-depleted values near Solfatara. A concentration profile across the Pozzuoli area (Fig. 3 b) depicts the least 18 O-depleted CO 2 in the air, having δ 18 O–CO 2 = − 0.06‰ in the vicinity of Solfatara and toward the northeast (Fig. 3 b). The δ 18 O-CO 2 values decrease to an average of − 2.5‰ northeast of Astroni. Downtown Naples has been identified as an area where airborne CO 2 exhibits more 18 O-depleted values, although zones with δ 18 O–CO 2 < − 6.5‰ are heterogeneously distributed between Pianura and Capodimonte, where CO 2 exhibits more 18 O-depleted CO 2 (i.e., transect C–C′ in Fig. 3 c). In this zone, heavily 18 O-depleted CO 2 (δ 18 O–CO 2 < − 16.0‰) was measured in the harbour district.
Additionally, a wide zone elongated NW–SE exhibits δ 18 O–CO 2 < − 5.3‰, extending from the eastern edge of downtown Naples to the west of Torre del Greco, coinciding with a densely urbanized area and a widespread industrialized area (i.e., Area Est-Centro direzionale). Figure 3 d illustrates a step gradient of δ 18 O–CO 2 that separates the coastal zone where δ 18 O–CO 2 = ~ − 5.06‰ from the inland area where δ 18 O–CO 2 = ~ − 2.85‰. The ∆ 18 O–CO 2 = 2.21 represents an order of magnitude greater than the accuracy of the oxygen isotope determination (± 0.25‰). In summary, the spatial variations of the measurements show strong fluctuations of δ 18 O–CO 2 in different zones. Kriging interpolation of the δ 18 O–CO 2 dataset reveals areas with slightly 18 O-depleted airborne CO 2 westward of Naples' downtown, and more 18 O-depleted values eastwards of downtown Naples. Similarly, wide variations in δ 13 C–CO 2 values correspond to spatial variations in the carbon isotopic signature of airborne CO 2 (Fig. six). Dataset statistics indicate that airborne CO 2 is 13 C-depleted compared to standard air. Cross-sections show trends indicating potential CO 2 sources with 13 C-depleted or enriched signatures in different areas, with notable variations near Baia and downtown Naples. These results suggest considerable variability in emission sources at the scale of the urbanized zone, and a dominant source of CO 2 with a 13 C-depleted signature. This expectation arises because the carbon isotope signature of airborne CO 2 can track the source of the gas 55 , 56 .
The cross-section through the urbanized areas of Bacoli (Figs. 2 a and 3 a) shows an average value of δ 13 C–CO 2 = − 10.5‰, indicating airborne CO 2 to be more 13 C-depleted than theoretical air and global reference values recorded by NOAA ( https://www.climate.gov/climatedashboard accessed on July 2, 2024). A significant change in the carbon isotope composition of airborne CO 2 is evident at Baia, where a decrease to a value of δ 13 C–CO 2 < − 14‰ was measured, coinciding with concentration values higher than those measured at Bacoli (Fig. 2 c). Low values of δ 13 C-CO 2 indicates that a heavily 13 C-depleted source of CO 2 is responsible for forcing airborne CO 2 above background levels and is the main contributor to increased CO 2 concentration. The carbon isotope composition increases to less 13 C-depleted values north of Cuma and achieves δ 13 C–CO 2 = ~ − 9‰ in the northern zone of the target area. The B–B′ cross-section shows a different trend compared to the A–A′ profile (Fig. 3 a,b respectively). Specifically, δ 13 C–CO 2 decreases from approximately − 9 to − 11‰. Continuing along the Solfatara profile (Fig. 3 b), a sudden increase in δ 13 C–CO 2 value is observed, reaching values of approximately − 8‰ at the highest concentration values observed along the same profile.
This trend appears to be clearly opposite to that observed in the A–A′ profile, suggesting the presence of potential CO 2 sources with a less 13 C-depleted signature compared to those forcing airborne CO 2 concentration in adjacent areas. The alternative hypothesis, suggesting that clean air with δ 13 C–CO 2 = − 8‰ produces the observed values, can be dismissed based on the evidence that 13 C-enrichment correlates with an increase in CO₂ concentration. This trend contradicts the expectation of a decrease in CO₂ concentration, which would be consistent with the clean air hypothesis. Furthermore, the A–A′, C–C′, and D–D′ profiles demonstrate that CO₂ concentrations exhibit opposite trends in comparison with δ 13 C–CO 2 . Specifically, these transects reveal that increases in CO₂ concentration coincide spatially with decreases in δ 13 C–CO 2 , indicating that the effective source of CO₂ in these zones is more 13 C-depleted. In the surrounding area, δ 13 C–CO 2 values average around − 10‰ regardless of CO 2 concentration in the air. These 13 C-depleted values reduce evidences of spatial 13 C-enrichment in airborne CO 2 . Therefore, the gas source which causes rise in CO 2 concentration above background levels in the area of Solfatara has a carbon isotope composition only slightly 13 C-depleted compared to the VPDB standard. Accordingly, to the northeast of the Astroni crater, δ 13 C–CO 2 decrease sharply to values ranging between − 10 and − 11‰. Moreover, zones with high airborne CO 2 concentrations near both Bagnoli and Posillipo also show heavily 13 C-depleted isotopic composition (i.e., δ 13 C–CO 2 = − 14.69‰ and δ 13 C–CO 2 = − 13.85‰, respectively).
The C–C′ profile (Fig. 3 c) crosses Naples’ downtown (Fig. 2 ), which is busiest by vehicle during morning hours. The airborne CO 2 has δ 13 C–CO 2 values from − 17.65 to − 8.54‰ with an average δ 13 C–CO 2 = ~ − 11‰. High CO 2 concentrations along this profile occur at the harbour district (Fig. 3 c and Fig. 2 a), which coincides with the zone having the most 13 C-depleted values of airborne CO 2 (Fig. 2 c). A comparison with other profiles reveals that a 13 C-depleted source of CO 2 forces the airborne CO 2 concentration in downtown Naples more efficiently than in peripheral zones to the west (i.e., Bacoli, Baia, and Posillipo). This source is less effective in forcing CO 2 concentration in the zone near Pozzuoli (B–B′ profile), where the source of CO 2 has a less 13 C-depleted carbon isotope composition. This northwest-oriented profile shows a zone with less 13 C-depleted values of airborne CO 2 to northwest (Fig. 2 c), consistent with a decrease in airborne CO 2 concentrations (Fig. 2 a).
Remarkable variations in the stable isotope composition of airborne CO 2 can be identified east of the urban area of Naples (Fig. 2 c). In concordance with δ 18 O–CO 2 , δ 13 C–CO 2 shows remarkable variations along the seaside compared to the inland along the D–D′ profile (Fig. 3 d). Airborne CO 2 concentrations fluctuate, superimposed on a decrease from the seaside to the inland. According to this trend, the carbon isotope composition shows an opposite trend from the most 13 C-depleted values in the coastal zone to the less 13 C-depleted CO 2 inland, revealing that potential sources of CO 2 with heavily 13 C-depleted signatures force airborne CO 2 concentration in the coastal zone near Portici and Torre del Greco. These sources are less effective in forcing CO 2 concentration inland, near San Giorgio a Cremano.
Measurements of CO 2 concentration, combined with stable isotope compositions of airborne CO 2 , provide relevant data for distinguishing between natural and anthropogenic CO 2 emissions in the atmosphere, and potentially tracking the gas dispersal from various sources of greenhouse gases at the urban spatial scale (i.e., 10 4 –10 5 m). This method overcomes the inherent difficulty of studying CO 2 dispersion caused by its high background level and subtle spatial variations of airborne CO 2 concentration. Indeed, various sources of CO 2 have different isotopic signatures for both carbon and oxygen.
There are several methods for tracking the dispersion of gases emitted from a source into the atmosphere. The methods commonly used to track gas dispersion are based on models that require a priori knowledge of the source, the amount of gas emitted, and the geometry of the dispersion area. Isotopic studies combined with atmospheric chemistry follow a different paradigm. The data collected from field measurements underwent analysis utilizing the Keeling plot method mass balance models for oxygen and carbon isotopes 49 , 50 , 57 . The Keeling plot method facilitates the determination of the primary CO 2 source at the local level using observational data.
At the same time, the mass balance model for oxygen and carbon isotopes allows an assessment of the influences of the individual CO 2 sources on the local air composition. The mathematical expressions governing this model were developed within the framework of previous studies 20 and are expounded concisely upon in the method section dedicated to assessing additional CO 2 in the atmosphere. This method allows for detecting the forcing effects introduced by the gas sources on the composition of the atmosphere. The measurements utilized in the theoretical model results (see Eq. ( 11 ) in this study) furnish point-by-point estimates of additional CO 2 concentration (i.e., the C fs ) along the trajectory.
Subsequently, the interpolation of C fs values employing the Kriging algorithm model facilitates the simulation of CO 2 dispersion. This algorithm generates a predictive layer for δ 13 C–CO 2 , δ 18 O–CO 2 , CO 2 concentration, and C fs. This method has been successfully applied to detect chemical and isotopic effects on the air in the La Fossa caldera on the island of Vulcano, both during periods of quiescent outgassing and during the recent period of increased volcanic outgassing in 2021 20 .
The Keeling plot illustrates a correlation between the carbon isotope composition of CO 2 and the inverse of airborne CO 2 concentration. Figure 4 shows the concentration dataset normalized by the global reference for airborne CO 2 concentration (i.e., 423 ppm vol). Each straight line on this plot represents binary mixing between the atmospheric background and an additional CO 2 source. The intercept on the isotopic axis provides the carbon isotopic signatures, facilitating the identification of the CO 2 emission source.
The correlation between δ 13 C–CO 2 and the inverse of airborne CO 2 concentration (i.e., Keeling plot). Data were normalized against the Global reference values recorded by NOAA (a https://www.climate.gov/climatedashboard accessed on July 2, 2024. ( a ) Dataset collected over the target area. ( b ) Urbanized areas of Naples. Green circles distinguish the subset of measurement collected near the airport (zone C in Fig. 1 a) from those collected in downtown Naples (zone B in Fig. 1 a). ( c ) Pozzuoli–Solfatara–Agnano area (zone A in Fig. 1 a) Yellow circles distinguish the subset of magmatic origin from that of anthropogenic origin in the area (blue circles).
Figure 4 a displays several mixing lines between background air and various potential sources of CO 2 , including natural (e.g., soil and plant respiration or volcanic degassing) and anthropogenic origins (e.g., combustion of fossil fuels or natural gas and landfill CO 2 emissions), whose isotopic signatures were retrieved from previous studies 37 . A geometric mean regression is recommended for the analysis of a scattered dataset (i.e., R 2 < 0.980) in the Keeling plot due to the inherent bias associated with determining the carbon isotopic signature through the utilization of a linear regression model 58 . The line representing the isotopic signature of the forcing source can be derived by applying a standard regression and subsequently dividing by the r-coefficient. This corrective approach aims to approximate the geometric mean regression through the utilization of a standard estimate obtained from a linear regression model.
The dataset collected over the target area reveals a variety of mixing lines, highlighting the inherent complexity of identifying a single CO 2 source. The alignments of δ 13 C–CO 2 in the Keeling plot suggests that fossil fuel combustion is a significant source of greenhouse gases, resulting in airborne CO 2 concentrations ranging from > 600 to ~ 1410 ppm vol (i.e., normalized values are from 0.7 to 0.3, respectively). However, multiple CO 2 sources can influence airborne CO 2 concentrations in the target area, especially at low to intermediate values (i.e., from 423 to 600 ppm vol, corresponding to normalized values ranging 0.7–1). These results support findings that human-related activities, such as urban mobility by vehicles and household heating, predominantly based on the combustion of fossil fuels, contribute significantly to rise the airborne CO 2 concentration. Nonetheless, natural CO 2 emissions, such as those from volcanic outgassing, which is estimated on the synoptic scale to account for approximately 1% of total annual emissions, can locally play a pivotal role in the amount of CO 2 injected into the atmosphere.
A sector in Naples’ downtown (i.e., Zone B in Fig. 1 a), distinct from Zone A, which includes the Campi Flegrei volcanic/hydrothermal zone and the western suburbs of Naples (i.e., Bagnoli and Posillipo), can serve as a test site to quantify the specific contribution to increasing airborne CO 2 concentrations caused by human-related emissions. Figure 4 b illustrates δ 13 C–CO 2 against CO 2 concentrations, showing good agreement with the mixing line between background air and CO 2 produced by fossil fuel combustion, characterized by a heavily 13 C-depleted signature (i.e., δ 13 C–CO 2 = − 29.94‰). Furthermore, data collected in the airport zone (i.e., Zone C in Fig. 1 a), where high levels of airborne CO 2 concentrations have been measured, indicate that the CO 2 source affecting both concentration and isotope composition of airborne CO 2 is of anthropogenic origin (i.e., δ 13 C–CO 2 = − 29.31 ‰).
Figure 4 c illustrates the complex distribution of concentration and carbon isotope composition values detected in the study area, predominantly located in the urban area of Pozzuoli, in the western suburbs of Naples. Results of cluster analysis applied to a subset of measurements collected in the Zone A (Fig. 1 a) reveal that multiple CO 2 sources play an almost equivalent role in elevating the concentration of airborne CO 2 above background levels. One subset of measurements, with CO 2 concentrations in the range 423–700 ppm vol, exhibits an isotopic signature in good agreement with the mixing line between background air and CO 2 produced by the combustion of fossil fuels (i.e., δ 13 C–CO 2 = − 32.93‰). Another subset of measurements indicates that δ 13 C–CO 2 of the air increases as CO 2 concentrations rise due to the influence of a less 13 C-depleted CO 2 source, with δ 13 C–CO 2 ≈ − 1.97‰. Although slightly lower, this value aligns with the carbon isotopic signature of CO 2 emitted from Pisciarelli and Bocca Grande fumaroles.
Those data retrieved from application of laboratory techniques to condensed fumarolic fluids have accuracy ± 0.1‰ 26 . Differences in the range Δ 13 C < 0.4‰ can be neglected because of the accuracy of the measurements with Deltaray (i.e., ± 0.25‰ according to 12 , 13 , 18 , 37 , 58 ).
Equation ( 11 ), included in the method section, facilitates the calculation of additional CO 2 in the air owing to either natural (i.e., volcanic/hydrothermal CO 2 ) or anthropogenic (i.e., produced by the combustion of fossil fuels) emissions. This calculation is based on input parameters in a theoretical model and measurements of airborne CO 2 concentration, δ 13 C–CO 2 , and δ 18 O–CO 2 in the field. C fs provides the concentration of the forcing source of CO 2 , exceeding local background levels in the atmosphere. A combination of the positioning of the endogenous sources of CO 2 and results of the Keeling plot helps distinguish the application of the mass balance model to the dispersal of volcanic CO 2 in the zone Solfatara (i.e., Zone A in Fig. 1 a) and the dispersal of CO 2 produced by the combustion of fossil fuels downtown Naples (i.e., Zone B in Fig. 1 a).
Figure 5 shows dispersions of CO 2 from anthropogenic origin in Naples’ downtown. In particular, the excess CO 2 concentrations in air produced by hydrocarbon combustion, which has a 13 C-depleted isotope composition compared to standard air (Fig. 4 b). For the calculation of the additional amount of CO 2 in the air, an anthropogenic source of CO 2 with the isotopic signature δ 13 C = − 31.00‰ and δ 18 O = − 16.00‰ has been adopted as the model parameters. Figure 5 a shows the fossil fuel-derived CO 2 has a heterogeneous distribution across the target area. A CO 2 dome 40 , 41 , 42 , 43 appears irregular and has numerous lobulations. The dome encloses islands where hydrocarbon combustion forces the CO 2 above the atmospheric background and generates concentration peaks even greater than + 300 ppm above the airborne CO 2 background (Fig. 5 b). One such island of high CO 2 concentration is well delineated in the harbour area, which is renowned for being among the Mediterranean's major harbours. In fact, the burning of hydrocarbons sustains the majority of the ship traffic in these areas. Another area with high CO 2 concentrations is located in the western downtown, near one of the most densely populated areas of Naples.
Dispersal of anthropogenic CO 2 in downtown Naples (cell size 10 m). The maps were generated in Qgis 3.34 environment ( https://qgis.org/download/ ) using Measurement interpolation generated in SAGA GIS environment ( https://saga-gis.sourceforge.io/en/index.html ). ( a ) CO 2 concentration map that shows the CO 2 concentration excess above the reference background. The concentration excess value of 83 ppm vol has been set as the threshold for transparency. ( b ) Vertical profile (black line in subplot a) of the excess CO 2 concentration across downtown Naples. ( c ) Wind vectors and speed recorded at C.N.R. station. ( d ) Wind direction frequency during survey.
The results of isotopic investigations prove the anthropogenic origin of atmospheric CO 2 . It is reasonable to assume that most of the anthropogenic CO 2 found in downtown Naples is the result of hydrocarbon combustion produced by urban mobility, given that the average air temperature during measurement collection was 22 °C (with an air temperature range of 19–23 °C). Within the one-day measurement acquisition timescale, variations in wind intensity and direction affecting the dispersion of CO 2 cannot be ruled out. This is particularly expected in Zone B, where the wind can influence the dispersion of emitted plumes near Naples' harbour. However, the data on wind direction (Fig. 5 c) and speed indicate that during the acquisition time window, the atmospheric circulation brought in SSW air, which is generally less enriched in anthropogenic CO 2 . Given the morphology of the study area and the local effects of densely built environments (Fig. 5 d), it is reasonable to assume a dilution effect of anthropogenic CO 2 due to the influence of less CO 2 —rich air from the sea. Accordingly, the anthropogenic CO 2 concentration along the C–C′ profile (Fig. 5 c) shows a notable increase in airborne CO 2 near the harbour and above a pedestrian area, suggesting that proximal sources of greenhouse gas emissions in the nearby areas are responsible for the increase in CO 2 above background levels.
Measurements collected at Pozzuoli (Zone A in Fig. 1 a) reveal multiple origins for CO 2 present in the air, namely volcanic and anthropogenic. Although human-related activities cause high concentrations of airborne CO 2 , a comparison with downtown can be made concerning the dispersal of geogenic CO 2 in the Pozzuoli area because the Campi Flegrei volcanic/hydrothermal system was in a state of unrest at the time of measurement collection 28 , 29 , 30 , 31 , 32 , 33 , 34 , 35 , 36 . Results of the cluster analysis provide a subset for calculating the amount of geogenic CO 2 that the main degassing zones at Campi Flegrei discharge into the atmosphere. The isotope composition and the airborne CO 2 concentration values of this subset were used in Eq. ( 11 ), with the values for the isotopic signatures of both carbon and oxygen for CO 2 emitted at Solfatara and Pisciarelli in May 2023 serving as model parameters (i.e., δ 13 C–CO 2 = − 1.67‰ vs VPDB and δ 18 O–CO 2 = − 7.85‰ vs VPDB) to obtain point-to-point calculations of volcanic CO 2 dispersal. These values provide insight into the dispersal of volcanic CO 2 from the main degassing vents at Campi Flegrei based on direct measurements and model parameters (Fig. 6 a). A comparison with the downtown map (Fig. 5 a) shows a more homogeneous dispersal of volcanic CO 2 along an N–S oriented dispersal zone. Furthermore, the volcanic CO 2 concentration is higher than 124 ppm vol above background levels in the area lying between Pozzuoli and Pisciarelli alone. Considering a background air CO 2 concentration of 423 ppm vol and the volcanic input calculated in the present study, this result is in good agreement with dispersal simulations averaged over a whole diurnal cycle obtained using the DISGAS software 14 . Measurements of concentration, corroborated by isotopic determinations, reveal volcanic CO 2 dispersion in the area of Bagnoli and eastward towards the urbanized area of Naples. In this area, measurements of airborne CO 2 concentrations alone are not able to track the dispersal of volcanic CO 2 because comparable absolute concentration values are found throughout the urban areas, where the additional CO 2 has anthropogenic origins. According to Granieri et al. 14 , inland air circulation prevails during nighttime in the Gulf of Naples, when volcanic CO 2 dispersal occurs towards the sea.
Dispersal of volcanic CO 2 in Pozzuoli-Solfatara zone (cell size 10 m). The maps were generated in Qgis 3.34 environment ( https://qgis.org/download/ ) using Measurement interpolation generated in SAGA GIS environment ( https://saga-gis.sourceforge.io/en/index.html ). ( a ) CO 2 concentration map that shows the CO 2 concentration excess above the reference background. The concentration excess value of 75 ppm vol has been set as the threshold for transparency. ( b ) Wind vectors and speed recorded at C.N.R. station. ( c ) Wind direction frequency during survey. ( d ) Vertical profile (black line in subplot a) of the excess CO 2 concentration across Pozzuoli-Solfatara area.
At the time of the survey, weather datasets reveal that NW winds blew from Solfatara towards the sea, even during the early morning (i.e., by ~ 8:00 UTC), after which sea breeze dominated air circulation from the SE throughout the daytime hours (Figura 6b,c). Arguably, the dispersal of volcanic CO 2 results from a combination of volcanic CO 2 dispersal and the residual layer that develops during nighttime and has not yet been disrupted by diurnal atmospheric turbulence. These results show that spatial surveys for studying airborne CO 2 helps in identifying multiple sources of greenhouse gases at the district scale of urban areas. Furthermore, stable isotope measurements allow an assessment of the impact of either volcanic degassing or anthropogenic emissions on airborne CO 2 concentrations.
The results of this study illustrate that integrating measurements of carbon and oxygen isotopic composition with those of CO 2 concentration aids in elucidating the genesis and development of CO 2 dome in urbanized areas. This represents a step forward in evaluating the impact of specific carbon dioxide sources, whether anthropogenic or natural, on the progression of climate change, as it facilitates the discernment of the underlying causes of urban domes through direct investigations.
The findings of this study also suggest that surveys conducted in urban areas such as Naples can be utilized to identify the primary regions for continuous monitoring of both natural and anthropogenic CO 2 emissions against global warming. Climate change has reached a global scale and threatens the stability of various vital sectors, including infrastructure, the economy, electricity production, international relations, biodiversity, and freshwater and food resources. Climate change affects all regions of the world, and its macroscopic effects manifest through extreme weather events, producing vast damage in cities and rural areas.
The international community is implementing a series of measures to combat ongoing climate change, which significantly impacts economic and social systems globally. For instance, several ambitious plans aim to reduce greenhouse gas emissions by 2050, mainly CO 2 . To achieve such ambitious goals, it is crucial to estimate and monitor CO 2 emissions, especially in urban areas where most CO 2 is produced through hydrocarbon combustion. Currently, no monitoring tools are available to detect near-real-time CO 2 emissions for individual countries. Therefore, efforts to monitor CO 2 in the air on a regional scale (synoptic ~ 10 6 m) with low latency (through the publication of hourly, daily, weekly, and annual data) via networks of stations installed in densely urbanized areas are becoming increasingly relevant. However, monitoring CO₂ in the atmosphere is not straightforward due to the high background concentration (approximately 400 ppm vol), which limits the potential for spatial variability. Consequently, monitoring the concentration alone may not always provide sufficient data for real-time estimation. Various studies demonstrate that integrating isotopic and concentration data provides information on the origin of CO 2 emissions 12 , 13 , 16 , 18 , 20 , 37 , 38 , 39 , 58 , 59 , 60 .
The δ 18 O–CO 2 largely depends on the CO 2 partitioning among the atmosphere, hydrosphere, lithosphere, and biosphere and can be deciphered through isotopic fractionation processes. Recent studies 12 , 18 , 19 , 20 show that it is possible to quantify atmospheric CO 2 emissions from natural and anthropogenic sources, isotopically characterized by δ 13 C–CO 2 and δ 18 O–CO 2 values, through integrated monitoring of atmospheric CO 2 concentration, isotopic composition, and meteorological data (direct investigations).
Therefore, the implementation of an active monitoring system is urgent and represents a paradigm shift in quantifying atmospheric CO 2 emissions at the scale of individual urbanized areas, compared to the currently applied methods based on statistical data at the national level for countries that are signatories to the United Nations Framework Convention on Climate Change 61 .
The instrument employed for data acquisition in this study is a Delta Ray–Thermo Fisher Scientific. It measures the concentration of the isotopologues 13 COO, 12 COO, and CO 18 O based on the adsorption strength of light in the mid-infrared region (~ 4.3 μm) following the Lambert–Beer law. The 13 C/ 12 C and 18 O/ 16 O ratios are calculated using different concentration ratios of the isotopologues, while the total CO 2 concentration is determined by summing the concentrations of the three CO 2 isotopologues. Stable isotope ratios are expressed in agreement with the VPDB scale using the δ-notation (i.e., δ 13 C–CO 2 and δ 18 O–CO 2 , respectively) within the CO 2 concentration range of 200–3500 ppm vol.
The Delta Ray instrument is equipped with the QTegra software. A specially designed template includes protocols for recording δ 13 C–CO 2 , δ 18 O–CO 2 , and CO 2 concentration values, along with information on the sample list, acquisition parameters, referencing, evaluation settings, and sample definition. Instrument calibration and referencing against two working standards ensure an accuracy of ± 0.25‰ for isotope determinations and ± 1 ppm vol for CO 2 concentration measurements.
The instrument records each measurement of δ 13 C–CO 2 and δ 18 O–CO 2 at a frequency of 1 Hz. Before data acquisition, the instrument conducts isotope ratio referencing on the working standards at a fixed CO 2 concentration (i.e., CO 2 = 400 ppm vol) approximating background airborne CO 2 . After purging the unknown air sample for 60 s, the instrument skips the purge and measures the concentration of CO 2 isotopologues in the air. Once the air has purged the gas inlet, the instrument calculates δ 13 C–CO 2 and δ 18 O–CO 2 , as well as CO 2 concentration.
An off-road vehicle housed the instrument, and the equipment for measuring δ 13 C–CO 2 , δ 18 O–CO 2 , and airborne CO 2 concentrations during the studies across the urbanized zone of Naples. The positioning of the vehicle was recorded by a global positioning system device (GARMIN GPSMAP® 64 s), time-synchronized with the Delta Ray's internal clock. In specific urban environments 12 , 37 , 38 , 39 , 55 , 56 , 62 and, more recently, in volcanic regions 18 , 20 , investigations have been conducted utilizing mobile laboratories to analyze the spatial variability of CO 2 .
An inverter (12 V input–output, pure sine wave) was connected to the car's electrical system, supplying power to the instrument (~ 300 W). A stainless-steel capillary (1/16 in.; Swagelok-typeTM, 3 m long) was connected to the instrument's inlet, with the other end attached to the front of the car roof (~ 2.3 m above the ground) to avoid potential contamination from the gasoline engine exhaust. The air passed through a filter (2 μm, 1/16 in, capillary aperture) to prevent contamination from dust on the roads. Considering the volume of the sampling capillary, the instrument's flow rate (approximately 100–110 ml min −1 ), and the average speed of the mobile laboratory (approximately 3.5 m s −1 ), the delay between measurements and their corresponding positions is approximately 25 m. This delay is comparable to the GPS positioning.
A route of approximately 170 km (Table 1 ) was designed in the laboratory to obtain a continuous, non-overlapping path, covering various environments in the wide urbanized area of Naples (Fig. 1 a). The route includes Miseno, Bacoli, Agnano, Campi Flegrei caldera, Pozzuoli, Capodimonte, Bagnoli, Posillipo to the east of Naples' downtown, and Portici, Ercolano, Torre del Greco, and San Giorgio a Cremano to the west, respectively. The route was planned to ensure that segments did not overlap, preventing an increase in the statistical weight of some route segments over others. The route was meticulously followed using a routing application (e.g., Google Maps). The survey was completed in thirteen hours at an average speed of 13 km h −1 , with the spatial density of measurements corresponding to the metric order (~ 4 m average distance between measurements). The dataset encompasses ~ 41,000 georeferenced measurements for δ 13 C–CO 2 , δ 18 O–CO 2 , and CO 2 concentration, respectively 63 . This method was already employed for a simultaneous airborne CO 2 spatial survey at Vulcano and revealed the dispersion of volcanic CO 2 through direct measurements 18 , 20 .
The data acquired from onsite measurements underwent processing utilizing the Keeling plot approach and mass balance models for oxygen and carbon isotopes. The Keeling plot enables the identification of the predominant CO 2 source at the local scale through observational data. The mass balance model for oxygen and carbon isotopes aims to quantify the impact of the CO 2 source on the local air. The algebraic equations for the model were developed as part of a previous study 18 and are detailed in the following paragraph of this paper, addressing the assessment of either volcanic or anthropogenic CO 2 in the air at Naples’ urban area. This methodology integrates measurements of stable CO 2 isotopes in the air with isotopic signatures of both the local CO 2 source, determined through the Keeling plot method 49 , 50 , and CO 2 in the background air. The theoretical outcomes of the model facilitate the partitioning of CO 2 in the air between the local background air and the CO 2 source.
The Keeling plot 49 , 50 , is the method broadly used to identify the isotopic signature of the gas source that increases CO 2 concentrations at the atmospheric background. The Keeling plot method facilitates the examination of the primary origin of atmospheric CO 2 by analyzing the δ 13 C–CO 2 against the reciprocal of CO 2 concentration. This method relies on mass balance principles, wherein a local CO 2 source alters the concentration from the atmospheric baseline. Mathematically, this is expressed by equations:
where C and δ 13 C denote CO 2 concentration and δ 13 C–CO 2 , respectively. Subscripts denote measured values (m), atmospheric background (a), and local source (fs). The linear combination of these equations generates a straight line in the δ 13 C versus 1/C plot 50 , 51 , 52 , 53 , 54 , 55 , 56 , 57 , 58 as delineated by equation
Equation ( 3 ) provides insight into the carbon isotope composition of the local CO 2 source under constant background and CO 2 source conditions.
To identify the main CO 2 source downtown Naples, a subset of measurements was selected. This subset encompasses measurements collected in an area of 24.60 km 2 centred in the Plebiscito square district of Naples (436,676.0 E and 4,520,817.0 N). Another subset, with its centre in the Pozzuoli area (Lat: 428,058.0 E; Long: 4,520,131.0 N, UTM), was selected for comparative purposes with data collected in Naples’ downtown. Specifically focusing on Pozzuoli (Zone A), the assessment focused on volcanic CO 2 as the primary source of CO 2 in the air. The measurements used in the theoretical model results (Eq. ( 11 ) reported below in this study) provide the concentration of the isotopically marked CO 2 source (e.g., volcanic or anthropogenic), causing the airborne CO 2 concentration to exceed the background concentration, point by point within the area. In the case of the Solfatara-Pisciarelli degassing area (Zone A in Fig. 1 ), the circular area is 47.28 km 2 , and the theoretical model provides the concentration of volcanic CO 2 (C V ). Following this calculation, the interpolation of C V values using the Kriging algorithm generates simulations of CO 2 from the forcing source (volcanic or anthropogenic for Solfatara-Pisciarelli and Naples' downtown, respectively). This algorithm produces the prediction layer for δ 13 C–CO 2 , δ 18 O–CO 2 , CO 2 concentration, and concentration (C V or C F , respectively) based on the assumption that each interpolating variable changes linearly with the distance between adjacent measurements. This assumption aligns with the expected homogeneity of spatial variations in atmospheric variables at the local scale 7 . Kriging interpolation is a geostatistical method used to estimate unknown values of each spatial variable based on known measurements of δ 13 C–CO 2 , δ 18 O–CO 2 , and CO 2 concentration at specific measurement points. The spatial correlation of the data is modeled using a Gaussian variogram, a standard variogram model defined by the equation:
where γ(h) is the semivariance at lag distance h, C 0 is the nugget, C is the partial sill, and a is the range. The kriging system of equations is set up using the Gaussian variogram model to determine the weights assigned to each known data point. These weights are calculated to minimize the estimation variance for the unknown points. The Gaussian model ensures smooth interpolation with continuous and differentiable transitions between estimated values, reflecting the assumed autocorrelation structure of the data.
Based on variogram analysis, CO 2 concentration measurements (Supplementary Fig. S1 online) are spatially dependent up to 700 m (i.e., the range), beyond which they become substantially independent. The range for δ 13 C–CO 2 , indicating the distance at which spatial correlation between carbon isotope measurements becomes negligible, has also been set to 700 m for kriging interpolation. For δ 18 O–CO 2 measurements, the range was determined to be 800 m. The partial sill was calculated as 1780 for CO₂ concentration, 1.62 for δ 13 C–CO 2 , and 1.65 for δ 18 O–CO 2 , indicating the variance attributable to the spatial structure for each variable. Simulations of stable isotope variables, airborne CO 2 concentration, and volcanic CO 2 dispersion were executed using the SAGA GIS software package ( https://saga-gis.sourceforge.io/en/index.html ).
An appropriate mass balance model for airborne CO 2 incorporates both isotopic parameters and concentration. Utilizing literature values for δ 13 C–CO 2 and δ 18 O–CO 2 of standard air (e.g., δ 13 C–CO 2 = − 8‰ and δ 18 O–CO 2 = − 0.1‰ 50 ) alongside values specific to CO 2 of external sources (e.g., either volcanic/hydrothermal or fossil fuel derived CO 2 ), an isotopic mass balance model incorporates four unknowns: background air CO 2 concentration, CO 2 concentration in the forcing source of gas, air CO 2 mixing fraction, and volcanic CO 2 mixing fraction.
The model is expressed by Eq. ( 5 ), which represents the CO 2 concentration in the air:
where C represents the CO 2 concentration and X denotes the mixing fraction between forcing source and atmospheric CO 2 , with subscripts m, a, and fs referring to measured, background, and local forcing source of CO 2 , respectively. This model operates under the assumptions that external source (i.e., volcanic or fossil fuel derived CO 2 ) significantly elevates CO 2 concentration relative to background levels.
The binary mixing equation to determine the relative weights of CO 2 from volcanic and atmospheric sources is given by Eq. ( 6 ):
Similarly, Eqs. ( 7 )
describe the isotopic mass balance models for carbon and oxygen isotopes of CO 2 , respectively. The combination of Eq. ( 6 ) and ( 7 ) provides Eq. ( 9 ), which allows for the calculation of X a
Using Eq. ( 9 ) in Eq. ( 8 ) yields Eq. ( 10 ), enabling the determination of X FS
By employing both Eqs. ( 9 ) and ( 10 ) and rearranging Eq. ( 5 ), we derive Eq. ( 11 )
which provides the concentration of CO 2 produced by the local effective gas source in the air C fs .
Cluster analysis was conducted to explore the relationships between airborne CO 2 concentrations and carbon isotope composition. Cluster analysis facilitates the classification of observational datasets into distinct classes based on specified similarity criteria. The objective of this analysis is to discern several groups of data that exhibit internal homogeneity (i.e., similarity criteria) while displaying heterogeneity among themselves concerning both CO 2 concentration and stable isotope compositions (i.e., δ 13 C–CO 2 and δ 18 O–CO 2 values). Various clustering methods are available for partitioning datasets (e.g., k-means, hierarchical, and two-way clustering), each differing in the requirement of preselecting the number of clusters, statistical properties of the dataset, or computational complexity.
Hierarchical clustering enables the grouping of objects such that those within a group are similar to each other and distinct from objects in other groups. Hierarchical clustering holds an advantage over alternative methods as it obviates the necessity of specifying the number of clusters a priori. The hierarchical structure of clusters can be formed using partitioning algorithms, initially considering all objects as individual clusters. Subsequently, through an iterative process, objects are assigned to different clusters based on principles maximizing the inter-cluster distances. One variant of hierarchical clustering is agglomerative clustering, where each object begins as its own cluster, and pairs of smaller clusters are successively merged until all data is encompassed within a single cluster. Essentially, hierarchical clustering assesses object similarity (i.e., distance) to form new clusters. Cluster merging is predicated on the Euclidean distance metric, reflecting the sum of squares of object coordinates in Euclidean space. Calculation of Euclidean distances leads to the updating of the distance matrix, with the iterative process culminating in the merging of the last two clusters into a final cluster encompassing the entire dataset.
Multiple approaches exist for computing inter-cluster distances and updating the proximity matrix, with some (e.g., single linkage or complete linkage) assessing minimum or maximum distances between objects from different clusters. In the cluster analysis of our dataset, we employed the Ward approach, which evaluates cluster variance rather than directly measuring distances, aiming to minimize variance among clusters. In Ward's method, the distance between two clusters is contingent upon the increase in the sum of squares when the clusters are combined. Ward's method implementation seeks to minimize the sum of squares distances of points from cluster centroids. In contrast to other distance-based methods, Ward's method exhibits less susceptibility to noise and outliers. Hence, in this paper, the Ward method is preferred over alternative methods for clustering.
This study presents findings from a spatial survey conducted in the metropolitan area of Naples, Italy, aimed at examining potential variations in atmospheric CO 2 sources. The urban zone of Naples was chosen due to its diverse CO 2 sources, including those from both geological (e.g., volcanic/hydrothermal emissions) and anthropogenic (e.g., combustion-related) origins. Situated within the extensive volcanic zone hosting Vesuvius, Campi Flegrei, and active volcano on Ischia, Naples provides a compelling location for such investigation owing to its dense urban population compared to other urban areas in the European continent.
Identification of CO 2 sources was facilitated through a combination of stable isotopic analysis and concentration measurements. Stable isotopic composition (i.e., carbon and oxygen isotopic ratios) and airborne CO 2 concentration were measured using a high-precision laser-based analyzer installed in an SUV vehicle. Measurements, recorded at 1 Hz, were synchronized with GPS data to ascertain spatial positioning, achieving a spatial resolution on a metric scale.
Spatial variations in both isotopic composition and concentration were derived from the dataset using the kriging algorithm with Gaussian autocorrelation. Resulting maps delineated three zones characterized by elevated CO 2 concentrations exhibiting distinct stable isotopic signatures. The zone with the highest CO 2 concentration encompassed Naples’ downtown and harbour district, while intermediate concentrations were observed inland across the urban area. Spatial simulations indicated lower CO 2 concentrations along the seaside to the west of downtown, consistent with local morning atmospheric circulation patterns oriented from SW to NE. Additional zones of heightened CO 2 concentrations were identified near the airport, situated northeast of downtown, and in proximity to inhabited areas such as Pozzuoli and Pisciarelli, near Solfatara to the west. These last areas (Pozzuoli and Pisciarelli) exhibit manifestations of a broad hydrothermal/magmatic system beneath the Campi Flegrei, constituting a geological source of airborne CO 2 . Anthropogenic CO 2 emissions, primarily from vehicular engine combustion, were found to elevate CO 2 concentrations above background levels in downtown Naples, near the airport, and in the vicinity of Solfatara.
A mixing model incorporating stable isotope composition and airborne CO 2 concentration allowed quantification of CO 2 contributions from different sources. Geochemical modeling based on this approach revealed spatial dispersal patterns of additional CO 2 near Solfatara and downtown Naples, with volcanic CO 2 dispersing northeastward under prevailing morning winds northeast oriented. This volcanic CO 2 extends beyond the hydrothermal zone, supplementing anthropogenic CO 2 emissions from vehicular traffic.
This study underscores the utility of combining isotopic and CO 2 concentration data for discerning the dispersion of both endogenous greenhouse CO 2 and emissions from anthropogenic activities. Particularly relevant in densely populated volcanic/hydrothermal regions, this methodology effectively distinguishes between natural and anthropogenic gas emissions in the atmosphere, overcoming challenges associated with high background levels and subtle spatial variations of the airborne CO 2 . Measurements in Naples were collected within a single day, during the diurnal phase of the planetary boundary layer (PBL) evolution, under turbulent conditions and mixing of the atmospheric layer closest to the ground. Consequently, the CO 2 dispersal maps represent average conditions for the urban area of Naples. Establishing monitoring programs for the concentration and isotopic composition of airborne CO 2 in Naples and other cities is crucial for studying the impact of the daily evolution of the PBL on potential variations in airborne CO 2 . This is particularly important in areas where geogenic sources (i.e., volcanic or hydrothermal) coexist with anthropogenic CO 2 emissions (e.g., from fossil fuel combustion) resulting from high population density.
The datasets generated during and/or analyzed during the current study are available in the ZENODO repository, https://zenodo.org/records/11300873 .
Calvert, J. G., Heywood, J. B., Sawyer, R. F. & Seinfeld, J. H. Achieving acceptable air quality: Some reflections on controlling vehicle emissions. Science 261 (5117), 37–45 (1993).
Article ADS CAS PubMed Google Scholar
Hasanbeigi, A., Price, L. & Lin, E. Emerging energy–efficiency and CO 2 emission–reduction technologies for cement and concrete production: A technical review. Renew. Sustain. Energy Rev. 16 (8), 6220–6238. https://doi.org/10.1016/j.rser.2012.07.019 (2012).
Article CAS Google Scholar
Andrew, R. M. Global CO 2 emissions from cement production, 1928–2017. Earth System Data Sci. 10 (4), 2213–2239. https://doi.org/10.5194/essd-10-2213-2018 (2018).
Article ADS Google Scholar
Paraschiv, S. & Paraschiv, L. S. Trends of carbon dioxide (CO 2 ) emissions from fossil fuels combustion (coal, gas and oil) in the EU member states from 1960 to 2018. Energy Rep. 6 (Supplement 8), 237–242. https://doi.org/10.1016/j.egyr.2020.11.116 (2020).
Article Google Scholar
Wimbadi, R. W., Djalante, R. & Mori, A. Urban experiments with public transport for low carbon mobility transitions in cities: A systematic literature review (1990–2020). Sustain. Cities Soc. 72 , 103023. https://doi.org/10.1016/j.scs.2021.103023 (2021).
Ahmadi, P., Dincer, I. & Rosen, M. A. Environmental Impact Assessments of Integrated Multigeneration Energy Systems. In Causes, Impacts and Solutions to Global Warming (eds Dincer, I. et al. ) 751–777 (Springer New York, New York, NY, 2013).
Chapter Google Scholar
Oke, T. R., Mills, G., Christen, A. & Voogt, J. A. Urban climates (Cambridge University Press, Cambridge, 2017).
Book Google Scholar
Umair, M., Kim, D. & Choi, M. Impact of climate, rising atmospheric carbon dioxide, and other environmental factors on water-use efficiency at multiple land cover types. Sci. Rep. 10 , 11644. https://doi.org/10.1038/s41598-020-68472-7 (2020).
Article CAS PubMed PubMed Central Google Scholar
Dasgupta, R. & Hirschmann, M. M. The deep carbon cycle and melting in Earth’s interior. Earth Planet. Sci. Lett. 298 , 1–13 (2010).
Article ADS CAS Google Scholar
Burton, M. R., Sawyer, G. M. & Granieri, D. Deep carbon emissions from volcanoes. Rev. Mineral Geochem. 75 , 323–354 (2013).
Aiuppa, A., Fischer, T. P., Plank, T. & Bani, P. CO 2 flux emissions from the Earth’s most actively degassing volcanoes, 2005–2015. Sci. Rep. 9 , 5442. https://doi.org/10.1038/s41598-019-41901-y (2019).
Article ADS CAS PubMed PubMed Central Google Scholar
Di Martino, R. M. R. & Capasso, G. On the complexity of anthropogenic and geological sources of carbon dioxide: Onsite differentiation using isotope surveying. Atm. Env. 2556 , 118446. https://doi.org/10.1016/j.atmosenv.2021.118446 (2021).
Di Martino, R. M. R. & Capasso, G. Fast and accurate both carbon and oxygen isotope determination in volcanic and urban gases using laser–based analyzer. Geophys. Res. Abstr. 21 , EGU2019-5095 (2019).
Google Scholar
Granieri, D., Costa, A., Macedonio, G., Bisson, M. & Chiodini, G. Carbon dioxide in the urban area of Naples: Contribution and effects of the volcanic source. J. Volcanol. Geoth. Res. 260 , 52–61 (2013).
Granieri, D. et al. Atmospheric dispersion of natural carbon dioxide emissions on Vulcano Island. Italy. J. Geophys Res. Solid Earth 119 , 5398–5413. https://doi.org/10.1002/2013JB010688 (2014).
Di Martino, R. M. R., Gurrieri, S., Camarda, M., Capasso, G. & Prano, V. Hazardous changes in soil CO 2 emissions at Vulcano, Italy, in 2021. J.G.R. Solid Earth 127 , e2022JB024516. https://doi.org/10.1029/2022JB024516 (2022).
Federico, C. et al. Inferences on the 2021 ongoing volcanic unrest at Vulcano Island (Italy) through a comprehensive multidisciplinary surveillance network. Remote Sens. 15 , 1405. https://doi.org/10.3390/rs15051405 (2023).
Di Martino, R. M. R. and Gurrieri, S. Stable isotope surveys reveal variations in the air CO 2 during the unrest event at Vulcano, Italy, in 2021. EGU General Assembly 2023, Vienna, Austria, 24–28 Apr 2023, EGU23-11463 (2023), https://doi.org/10.5194/egusphere-egu23-11463
Gurrieri, S., Di Martino, R. M. R., Camarda, M. & Francofonte, V. Monitoring CO 2 hazard of volcanic origin: a case study at the island of Vulcano (Italy) during 2021–2022. Geosciences 13 , 266. https://doi.org/10.3390/geosciences13090266 (2023).
Di Martino, R. M. R. & Gurrieri, S. Quantification of the volcanic carbon dioxide in the air of Vulcano Porto by stable isotope surveys. J. Geophys. Res. Atmos. 128 , e2022JD037706 (2023).
Orsi, G., D’Antonio, M., De Vita, S. & Gallo, G. The Neapolitan Yellow Tuff, a large-magnitude trachytic phreatoplinian eruption: Eruptive dynamics, magma withdrawal and caldera collapse. J. Volcanol. Geotherm. Res. 53 , 275–287. https://doi.org/10.1016/0377-0273(92)90086-S (1992).
Orsi, G., De Vita, S. & Di Vito, M. The restless, resurgent Campi Flegrei nested caldera (Italy): Constraints on its evolution and configuration. J. Volcanol. Geotherm. Res. 74 (3–4), 179–214. https://doi.org/10.1016/S0377-0273(96)00063-7 (1996).
Chiodini, G., Cioni, R., Guidi, M., Raco, B. & Marini, L. Soil CO 2 flux measurements in volcanic and geothermal areas. Appl. Geochem. 1 , 543–552 (1998).
Chiodini, G. & Marini, L. Hydrothermal gas equilibria the H 2 O-H 2 -CO 2 -CO-CH 4 system. Geochim. Cosmochim. Acta 15 , 2673–2687 (1998).
Chiodini, G. et al. CO 2 degassing and energy release at Solfatara Volcano, Campi Flegrei. Italy. J. Geophys. Res. 106 (B8), 16213–16221 (2001).
Caliro, S. et al. The origin of the fumaroles of La Solfatara (Campi Flegrei South Italy). Geochim. Cosmochim. Acta 71 , 3040–3055 (2007).
Chiodini, G. et al. Long term variations of the Campi Flegrei (Italy) volcanic system as revealed by the monitoring of hydrothermal activity. J. Geophys. Res. 115 (B03205), 531–542 (2010).
Caliro, S., Chiodini, G. & Paonita, A. Geochemical evidences of magma dynamics at Campi Flegrei (Italy). Geochim. Cosmochim. Acta 132 , 1–15. https://doi.org/10.1016/j.gca.2014.01.021 (2014).
Chiodini, G. et al. Hydrothermal pressure-temperature control on CO 2 emissions and seismicity at Campi Flegrei (Italy). J. Volcanol. Geoth. Res. 414 , 107245. https://doi.org/10.1016/j.jvolgeores.2021.107245 (2021).
Chiodini, G., Caliro, S., De Martino, P., Avino, R. & Gherardi, F. Early signals of new volcanic unrest at Campi Flegrei caldera? Insights from geochemical data and physical simulations. Geology 40 , 943–946. https://doi.org/10.1130/G33251.1 (2012).
Chiodini, G. et al. Magmas near the critical degassing pressure drive volcanic unrest towards a critical state. Nat. Com. 7 , 1–9. https://doi.org/10.1038/ncomms13712 (2016).
Chiodini, G., Pappalardo, L., Aiuppa, A. & Caliro, S. The geological CO 2 degassing history of a long-lived caldera. Geology 43 , 767–770. https://doi.org/10.1130/G36905.1 (2015).
Chiodini, G. et al. Magma degassing as a trigger of bradyseismic events: The case of Phlegrean fields (Italy). Geophys. Res. Lett. 30 (8), 1434. https://doi.org/10.1029/2002GL016790 (2003).
Chiodini, G. et al. Evidence of thermal-driven processes triggering the 2005–2014 unrest at Campi Flegrei caldera. Earth Planet. Sci. Lett. 414 , 58–67. https://doi.org/10.1016/j.epsl.2015.01.012 (2015).
De Siena, L. et al. Source and dynamics of a volcanic caldera unrest: Campi Flegrei, 1983–84. Sci. Rep. 7 , 8099. https://doi.org/10.1038/s41598-017-08192-7 (2017).
Tramelli, A. et al. Statistics of seismicity to investigate the Campi Flegrei caldera unrest. Sci. Rep. 11 , 7211. https://doi.org/10.1038/s41598-021-86506-6 (2021).
Di Martino, R. M. R. & Gurrieri, S. Theoretical principles and application to measure the flux of carbon dioxide in the air of urban zones. Atm. Env. 288 (4), 119302. https://doi.org/10.1016/j.atmosenv.2022.119302 (2022).
Carapezza, M. L., Di Martino, R. M. R., Capasso, G., Di Gangi, F., Ranaldi, M., Tarchini, L. Stable isotope composition and airborne concentration of CO 2 in Rome capital city (Italy). EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-11294, 10.5194/egusphere-egu24-11294, (2024).
Gurrieri, S. and Di Martino, R. M. R.: Spatial Isotopic Analysis of Airborne CO 2 : Insights from Etna Volcano and Madonie Mountains, Italy, Surveys, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-5556, 10.5194/egusphere-egu24-5556, (2024).
Idso, C. D., Idso, S. B. & Balling, R. C. Jr. The Urban CO 2 dome of Phoenix, Arizona. Phys. Geogr. 19 (2), 95–108. https://doi.org/10.1080/02723646.1998.10642642 (1998).
Idso, C. D., Idso, S. B. & Balling, R. C. An intensive two-week study of an urban CO 2 dome in Phoenix, Arizona. USA. Atm. Env. 35 , 995–1000 (2001).
Balling, R. C., Cerveny, R. S. & Idso, C. D. Does the urban CO 2 dome of Phoenix, Arizona contribute to its heat island?. Geophys. Res. Lett. 28 , 4599–4601. https://doi.org/10.1029/2000GL012632 (2001).
Jacobson, M. Z. Enhancement of local air pollution by urban CO 2 domes. Environ. Sci. Technol 44 , 2497–2502 (2010).
Sinclair, A. J. Selection of threshold values in geochemical data using probability graphs. J. Geochem. Explor. 3 , 129–149. https://doi.org/10.1016/0375-6742(74)90030-2 (1974).
Yakir, D. The stable isotopic composition of atmospheric CO2. In The Atmosphere, Treatise on Geochemistry (eds Holland, H. D. & Turekian, K. K.) 175–212 (Elsevier-Pergamon, The Netherlands, 2003).
Park, R. & Epstein, S. Carbon isotope fractionation during photosynthesis. Geochim. Cosmochim. Acta 44 , 5–15 (1960).
Flanagan, L. B., Brooks, J. R., Vanery, G. T. & Ehleringr, J. R. Discrimination against C 18 O 16 O during photosynthesis and the oxygen isotope ratio of respired CO 2 in boreal forest ecosystems. Global Biogeochem. Cy. 11 (1), 83–98 (1997).
Farquhar, G. D. et al. Vegetation effects on the isotope composition of oxygen in atmospheric CO 2 . Nature 363 , 439–443 (1993).
Keeling, C. D. The concentration and isotopic abundance of carbon dioxide in rural and marin areas. Geochim. Cosmochim. Acta 13 , 322–334 (1958).
Keeling, C. D. The concentration and isotopic abundance of carbon dioxide in rural areas. Geochim. Cosmochim. Acta 24 , 277–298 (1961).
Buenning, N., Noone, D., Randerson, J., Riley, W. J. & Still, C. The response of the 18 O/ 16 O composition of atmospheric CO 2 to changes in environmental conditions. J. Geoph. Res. Biogeosci. 119 , 55–79. https://doi.org/10.1002/2013JG002312 (2014).
Francey, R. J. & Tans, P. P. Latitudinal variation in oxygen-18 of atmospheric CO 2 . Nature 327 (6122), 495–497 (1987).
National Institute for Occupational Safety and Health (NIOSH). Documentation for Immediately Dangerous to Life or Health Concentrations (IDLHs) for Carbon Dioxide (1997). Available online: hCp:// www.cdc.gov/niosh/idlh/124389.html (Accessed 15 August 2023).
Shumacher, M. et al. Oxygen isotopic signature of CO 2 from combustion processes. Atmos. Chem. Phys. 11 , 1473–1490. https://doi.org/10.5194/acp-11-1473-2011 (2011).
Newman, S., Xu, X., Affek, H. P., Stolper, E. & Epstein, S. Changes in mixing ratio and isotopic composition of CO2 in urban air from the Los Angeles basin, California, between 1972 and 2003. J. Geophys. Res. 113 , D23304. https://doi.org/10.1029/2008JD009999 (2008).
Newman, S. et al. Toward consistency between trends in bottom-up CO2 emissions and top-down atmospheric measurements in the Los Angeles megacity. Atmos. Chem. Phys. 16 , 3843–3863. https://doi.org/10.5194/acp-16-3843-2016 (2016).
Pataki, D. E. et al. The application and interpretation of keeling plots in terrestrial carbon cycle research. Global Biogeochem. Cycles 17 (1), 1022 (2003).
Venturi, S. et al. A multi–instrumental geochemical approach to assess the environmental impact of CO2–rich gas emissions in a densely populated area: The case of Cava dei Selci (Latium, Italy). App. Geochem. 101 , 109–126 (2019).
Di Martino, R. M. R., Capasso, G. & Camarda, M. Spatial domain analysis of carbon dioxide from soils on Vulcano Island: Implications for CO 2 output evaluation. Chem. Geol. 444 , 59–70. https://doi.org/10.1016/j.chemgeo.2016.09.037 (2016).
Di Martino, R. M. R., Capasso, G., Camarda, M., De Gregorio, S. & Prano, V. Deep CO 2 release revealed by stable isotope and diffuse degassing surveys at Vulcano (Aeolian Islands) in 2015–2018. J. Volcanol. Geoth. Res. 401 , 106972. https://doi.org/10.1016/j.jvolgeores.2020.106972 (2020).
United Nations Framework Convention on Climate Change, May 9, 1992, S. Treaty Doc No. 102-38, 1771 U.N.T.S. 107.
Imasu, R. & Tanabe, Y. Diurnal and seasonal variations of carbon dioxide (CO 2 ) concentration in urban, suburban and rural areas around Tokyo. Atmosphere 9 , 367. https://doi.org/10.3390/atmos9100367 (2018).
Di Martino, R. M. R., Gurrieri, S., Paonita, A., Caliro, S., Santi, A. Dataset on Atmospheric CO 2 in the Metropolitan Area of Naples, Italy. Zenodo (2024).
Download references
Authors and affiliations.
Istituto Nazionale di Geofisica e Vulcanologia, Sezione di Palermo, Via Ugo La Malfa, 153, 90146, Palermo, Italy
Roberto M. R. Di Martino, Sergio Gurrieri & Antonio Paonita
Istituto Nazionale di Geofisica e Vulcanologia, Sezione di Napoli, Osservatorio Vesuviano, Via Diocleziano 328, 80124, Naples, Italy
Stefano Caliro & Alessandro Santi
You can also search for this author in PubMed Google Scholar
R.M.R.D.M.: Conceptualization, Methodology, Software, Validation, Formal analysis, Investigation, Writing Original Draft, Writing—Review & Editing, Visualization, Supervision; S.G.: Conceptualization, Investigation, Writing—Review & Editing; A.P.: Writing—Review & Editing; S.C.: Writing—Review & Editing; A.S.: Investigation, Writing—Review & Editing.
Correspondence to Roberto M. R. Di Martino .
Competing interests.
The authors declare no competing interests.
Publisher's note.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary figures., supplementary table s1., 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
Cite this article.
Di Martino, R.M.R., Gurrieri, S., Paonita, A. et al. Unveiling spatial variations in atmospheric CO 2 sources: a case study of metropolitan area of Naples, Italy. Sci Rep 14 , 20483 (2024). https://doi.org/10.1038/s41598-024-71348-9
Download citation
Received : 30 May 2024
Accepted : 27 August 2024
Published : 03 September 2024
DOI : https://doi.org/10.1038/s41598-024-71348-9
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.
Sign up for the Nature Briefing: Anthropocene newsletter — what matters in anthropocene research, free to your inbox weekly.
IMAGES
VIDEO
COMMENTS
A case study is one of the most commonly used methodologies of social research. This article attempts to look into the various dimensions of a case study research strategy, the different epistemological strands which determine the particular case study type and approach adopted in the field, discusses the factors which can enhance the effectiveness of a case study research, and the debate ...
A case study is a detailed study of a specific subject, such as a person, group, place, event, organization, or phenomenon. Case studies are commonly used in social, educational, clinical, and business research. A case study research design usually involves qualitative methods, but quantitative methods are sometimes also used.
The purpose of case study research is twofold: (1) to provide descriptive information and (2) to suggest theoretical relevance. Rich description enables an in-depth or sharpened understanding of the case. It is unique given one characteristic: case studies draw from more than one data source. Case studies are inherently multimodal or mixed ...
A case study is a research method that involves an in-depth examination and analysis of a particular phenomenon or case, such as an individual, organization, community, event, or situation. It is a qualitative research approach that aims to provide a detailed and comprehensive understanding of the case being studied. Case studies typically ...
Qualitative case study methodology provides tools for researchers to study complex phenomena within their contexts. When the approach is applied correctly, it becomes a valuable method for health ...
Case study method is the most widely used method in academia for researchers interested in qualitative research (Baskarada, 2014). Research students select the case study as a method without understanding array of factors that can affect the outcome of their research.
Purpose of case study methodology. Case study methodology is often used to develop an in-depth, holistic understanding of a specific phenomenon within a specified context. 11 It focuses on studying one or multiple cases over time and uses an in-depth analysis of multiple information sources. 16,17 It is ideal for situations including, but not limited to, exploring under-researched and real ...
A case study is a research approach that is used to generate an in-depth, multi-faceted understanding of a complex issue in its real-life context. It is an established research design that is used extensively in a wide variety of disciplines, particularly in the social sciences. A case study can be defined in a variety of ways (Table 5), the ...
A case study protocol outlines the procedures and general rules to be followed during the case study. This includes the data collection methods to be used, the sources of data, and the procedures for analysis. Having a detailed case study protocol ensures consistency and reliability in the study.
Case study methodology can entail the study of one or more "cases," that could be described as instances, examples, or settings where the problem or phenomenon can be examined. The researcher is tasked with defining the parameters of the case, that is, what is included and excluded. This process is called bounding the case, or setting boundaries.
Definitions of qualitative case study research. Case study research is an investigation and analysis of a single or collective case, intended to capture the complexity of the object of study (Stake, 1995).Qualitative case study research, as described by Stake (), draws together "naturalistic, holistic, ethnographic, phenomenological, and biographic research methods" in a bricoleur design ...
A Case study is: An in-depth research design that primarily uses a qualitative methodology but sometimes includes quantitative methodology. Used to examine an identifiable problem confirmed through research. Used to investigate an individual, group of people, organization, or event. Used to mostly answer "how" and "why" questions.
The definition above is an example of an all-inclusive descriptive definition of case study research represented by Yin (2003).According to the definition of case study research, there is no doubt that this research strategy is one of the most powerful methods used by researchers to realize both practical and theoretical aims.
This article defends case study methodology as an appropriate methodology, giving a description, the process and its strengths and weaknesses. The Case Study Approach. This article by Crowe et al gives a nice overview of case studies and includes several examples from health science research.
The case study is not a research method, but researchers select methods of data collection and analysis that will generate material suitable for case studies. Freud (1909a, 1909b) conducted very detailed investigations into the private lives of his patients in an attempt to both understand and help them overcome their illnesses.
Qualitative case study methodology provides tools for researchers to study complex phenomena within their contexts. When the approach is applied correctly, it becomes a valuable method for health science research to develop theory, evaluate programs, and develop interventions. The purpose of this paper is to
A case study is a detailed study of a specific subject, such as a person, group, place, event, organisation, or phenomenon. Case studies are commonly used in social, educational, clinical, and business research. A case study research design usually involves qualitative methods, but quantitative methods are sometimes also used.
A case study methodology that combines a real-time longitudinal three-year study with nine retrospective case studies about the same phenomenon and enhances three kinds of validity: construct, internal and external is described. Expand. 1,483. Highly Influential.
It's been 100 years since Harvard Business School began using the case study method. Beyond teaching specific subject matter, the case study method excels in instilling meta-skills in students.
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 ...
Qualitative case study methodology enables researchers to conduct an in-depth exploration of intricate phenomena within some specific context. By keeping in mind research students, this article presents a systematic step-by-step guide to conduct a case study in the business discipline. Research students belonging to said discipline face issues ...
Simply put, the case method is a discussion of real-life situations that business executives have faced. Harvard Business School. The Learning Experience. The Case Study Method. On average, you'll attend three to four different classes a day, for a total of about six hours of class time (schedules vary). To prepare, you'll work through problems ...
Despite the effort invested to improve the teaching of programming, students often face problems with understanding its principles when using traditional learning approaches. This paper presents a novel teaching method for programming, combining the task-driven methodology and the case study approach. This method is called a task-driven case study. The case study aspect should provide a real ...
The objective of this paper is to estimate the equivalent permeability of the rock surrounding the tailrace tunnel of the Azad Dam pumped storage power plant, using geostatistical methods. The permeability of the rock mass is a critical factor that influences the estimation of water flow rates. Since the tunnel passes through various geological units with different permeabilities, it is ...
Case study is a common methodology in the social sciences (management, psychology, science of education, political science, sociology). A lot of methodological papers have been dedicated to case study but, paradoxically, the question "what is a case?" has been less studied. Hence the fact that researchers conducting a case study are ...
Zhang et al. [18] utilized theoretical and semi-empirical methods to analyze the impact of operating backpressure on PEMFC reactions and performance. The study demonstrated that variations in operating backpressure can affect reversible thermodynamic potentials, open-circuit voltage (OCV), membrane conductivity, and mass transfer characteristics.
ELISpot response rates and magnitudes for case-control vaccine recipients were both substantially lower than observed previously with this vaccine regimen 3, 4, 6 and in the HVTN 705 pilot study (59.6% vs. 82.4% positive response frequency in vaccine non-cases, respectively; median readout 74 vs. 308, Figure S15), differences that were not ...
Abstract. The purpose of this article is to provide a comprehensive view of the case study process from the researcher's perspective, emphasizing methodological considerations. As opposed to other qualitative or quantitative research strategies, such as grounded theory or surveys, there are virtually no specific requirements guiding case research.
This methodology integrates measurements of stable CO 2 isotopes in the air with isotopic signatures of both the local CO 2 source, determined through the Keeling plot method 49,50, and CO 2 in ...
When conducting QLR, time is the lens used to inform the overall study design and processes of data collection and analysis. While QLR is an evolving methodology, spanning diverse disciplines (Holland et al., 2006), a key feature is the collection of data on more than one occasion, often described as waves (Neale, 2021).Thus, researchers embarking on designing a new study need to consider ...