Lagrangian Trajectory Model Research Articles

The quasi-stochastic nature of hail growth: Hail trajectory clusters in simulations of the Kingfisher, OK hailstorm

Abstract Current methods that incorporate synoptic, meso-alpha, and even some storm-scale processes to distinguish among storms that produce significantly severe versus merely severe hail have not been highly successful. These results suggest that smaller-scale processes, such as meso-beta, meso-gamma, and microscale, have a larger impact on hail production than previously theorized. Relatedly, previous research, often using hail trajectories driven by steady-state model fields, has also frequently painted hail trajectories as having simple, relatively uniform motions within a supercell. This study uses a Lagrangian hail trajectory model embedded in two idealized supercell simulations and updated every model timestep, in conjunction with a trajectory clustering technique, to examine the variability of hail trajectories both within and across the two simulations. The clustering method identified 11 distinct clusters in both supercell simulations. Clusters sourced embryos from almost all quadrants around the updraft core. Trajectories in both simulations exhibited complex and highly variable behavior, occasionally including multiple vertical loops. Cluster trajectories could follow similar paths but produce different distributions of mass at the surface, or conversely, could follow different paths but result in similar mass distributions at the surface. An aggregate review of all trajectories showed the larger hail in each simulation was somewhat more likely to be produced in regions with less strong peak updraft speeds, but other environmental profile information was more mixed, particularly for the simulation producing the largest hail. In sum, storm-scale information is highly impactful on hail trajectories, and needs to be considered in forecasting applications. Future hail trajectory modeling studies also need to incorporate more realistic storm-scale information.

Individual-based numerical experiment to describe the distribution of floating kelp within the Southern Benguela Upwelling System

Abstract Kelps are resilient organisms, capable of thriving in high-energy wave environments. However, when hydrodynamic drag forces exerted by the wave environment exceed the kelps’ structural limits, individuals become dislodged. Floating kelps generally follow ocean currents, traveling long distances until air-filled structures fail or the epibiont load becomes too great, causing them to sink to the seafloor. The ability of kelp to disperse over vast offshore and nearshore systems makes them important for organic subsidy and as a dispersal vector for marine organisms. Previous research on dislodged macroalgae focused on context-specific rafts, limiting insights into the broader ecological role of floating kelp. This study employed a site-specific Lagrangian trajectory model to describe the spatial distribution of floating Ecklonia maxima along the South African coastline. The model incorporated buoyancy and sinking using site-specific morphological data. Findings revealed that the distribution of floating E. maxima is influenced by oceanographic conditions, and seasonal patterns were also evident. Mesoscale features played a vital role in kelp accumulation on the surface and seafloor and acted as barriers to dispersal. This study offers essential insights into kelp’s role as an organic subsidy and provides numerical evidence for kelp’s potential as a carbon sink, contributing to a better understanding of kelp ecosystems and their ecological functions.

Stratospheric transport and tropospheric sink of solar geoengineering aerosol: a Lagrangian analysis

Stratospheric aerosol injection (SAI) aims to reflect solar radiation by increasing the stratospheric aerosol burden. To understand how the background circulation influences stratospheric transport of injected particles, we use a Lagrangian trajectory model (lacking numerical diffusion) to quantify particles’ number, flux, lifetime, and tropospheric sinks from a SAI injection strategy under present-day conditions. While particles are being injected, stratospheric particle number increases until reaching a steady-state. During the steady-state, the time series of particle number shows a dominant period of ~2 years (rather than a 1-year cycle), suggesting modulation by the quasi-biannual oscillation. More than half of particles, injected in the tropical lower stratosphere (15° S to 15° N, 65 hPa), undergo quasi-horizontal transport to the midlatitude. We find a zonal asymmetry of particles’ tropospheric sinks that are co-located with tropopause folding beneath the midlatitude jet stream, which can help predict tropospheric impacts of SAI (e.g., cirrus cloud thinning).

Insights into atmospheric trace gases, aerosols, and transport processes at a high-altitude station (2623 m a.s.l.) in Northeast Asia

As part of the efforts to understand long-range air pollution transport in Northeast Asia, a comprehensive field measurement campaign was conducted at the top of Changbai Mountain (2623 m a.s.l.) from 25th July to 21st August in 2021. The ambient concentrations of trace gases and chemical characteristics of fine particulate matter (PM2.5) were comparable to those at background or free tropospheric locations in China, but exceeded levels at corresponding sites in Europe and the United States. Well-defined diurnal variations were evident in most air pollutants, including carbon monoxide (CO), sulfur dioxide (SO2), and nitrogen oxide (NOx), reflecting the influence of boundary layer development and local anthropogenic activities. Ozone (O3) exhibited a unique diurnal distribution, with lower concentrations during the daytime and higher levels at nighttime, indicating a pronounced impact of downward transport of O3-rich air from the upper troposphere during the nighttime. During the campaign, organic carbon (OC) constituted the predominant component of PM2.5, comprising 43% of the total, while sulfate (SO42−) was the most abundant secondary inorganic component in PM2.5, accounting for 26%. Combining the analysis of air mass backward trajectories, we found that Changbai Mountain is primarily influenced by air masses originating from North Korea and the Sea of Japan, constituting 62% of the total, and local vegetation biological emissions and anthropogenic emissions in North Korea notably affect the observation site. Conversely, the secondary aerosol components were most abundant in marine air mass originating from the marginal seas of China and the strait between the Korean peninsula and southern Japan. Utilizing the Lagrangian photochemical trajectory model (LPTM), a significant increase in particulate SO42− concentrations was observed when air mass passed through the marginal seas of China, primarily due to SO2 oxidation reaction. This study contributes to an improved understanding of the atmospheric chemistry in the high-altitude regions of Northeast Asia and underscores the substantial role of secondary inorganic aerosol formation during the transport of marine air masses.

Moisture Sources and Pathways of Annual Maximum Precipitation in the Lancang‐Mekong River Basin

AbstractRecent extremely heavy precipitation has led to substantial economic losses and affected millions of residences in the Lancang‐Mekong River Basin (LMRB). This study analyzed the spatial‐temporal characteristics of the annual maximum precipitation (R1X) of the LMRB and identified the moisture sources and pathways conducive to R1Xs using a Lagrangian back trajectory model. Results show that India Ocean and Bay of Bengal (IO/BOB), local evapotranspiration, and West Pacific Ocean and East China (WP/EC) are the three main moisture transport pathways of the R1Xs in LMRB, contributing 68.3%, 20.4% and 11.3% of the trajectories, respectively. R1Xs in the downstream eastern area are affected by tropical cyclones bringing large amounts of moisture from the WP/EC. As tropical cyclones shifted northward under climate change impact, more extreme precipitation occurred over the LMRB due to moisture coming from WP/EC, but those from the IO/BOB had decreased because of the slowdown of flows across the Equator.

Heat and Moisture Anomalies During Crop Failure Events in the Southeastern Australian Wheat Belt

AbstractProlonged droughts and heatwaves are common causes of agricultural failure in Australia, yet the origins of these climate anomalies remain understudied. Here, we use a Lagrangian trajectory model driven by atmospheric reanalysis and constrained by satellite data to unravel the sources of precipitation and heat over the Southeastern Australia wheat belt. Furthermore, we assess the impact of local and upwind drought conditions on the moisture and heat imports to the region. Results indicate that the most extreme crop failure events over the wheat belt (i.e., 1994, 2002, and 2006) were associated with persistent high‐pressure systems. The ocean provided on average 72% of the moisture for precipitation and 39% of the heat arriving over the wheat belt, with the moisture sources substantially decreasing during crop failure events. Upwind drought further intensified rainfall deficits and heat stress during these events due to lower moisture and higher heat imports to the region. This was particularly clear during the initial phase of the Millennium Drought in 2002. Then, yield deficits exceeded 50%, and ∼4% of the precipitation originated from drought‐affected regions upwind, compared to the 9% that was expected climatologically from those regions. Simultaneously, the heat import from these regions upwind increased by ∼10 W m−2, from the climatological 23%–25%, during this event. While these results indicate a limited potential for upwind land management to mitigate downwind agricultural loss in the Southeastern Australia wheat belt, other agricultural regions with a higher climatic dependency on remote land may benefit from such strategies.

Impact of displacement ventilation and underfloor air distribution systems on concentrations of indoor particle in different seasons

The particle motion in a room is affected by the thermal environment and the air supply conditions caused by seasonal changes. In the study, a Lagrangian discrete trajectory model combined with the Eulerian fluid method was adopted to simulate the characteristics of indoor airflow and the diffusion of two typical indoor particles with considering different air supply conditions in different seasons. Specifically, the thermal comfort perception, draught rating (DR) of the indoor personnel, and the contaminant removal effectiveness (CRE, ε) of indoor particles under the influences of displacement ventilation (DV) and underfloor air distribution (UFAD) systems were investigated. The results revealed that when considering the indoor thermal comfort under summer and winter air supply conditions, UFAD was found be more effective than the DV system in protecting the temperature difference of indoor personnel. Additionally, the removal effect of PM1 in the breathing area was superior compared to that of the DV system. Regarding the removal effectiveness of PM2.5, the removal performance of both ventilation systems decreased, but the removal effect of the DV system on PM2.5 in the breathing area (with the largest εl of 1.55) was significantly better than that of the UFAD system under winter operating conditions. It is essential for the selection of the ventilation system to consider the thermal comfort of indoor occupants, as different ventilation systems have varying capabilities to remove particles of different sizes. Therefore, choosing an appropriate ventilation system is crucial for enhancing indoor air quality, and the findings of this study provide a basis for improving indoor air quality.

PM2.5 in Cape Town, South Africa: Chemical characterization and source apportionment using dispersion-normalised positive matrix factorization

Understanding fine particulate matter (PM2.5) composition and sources is beneficial to improving visibility, addressing climate change, and mitigating poor air quality and related public health effects. Source apportionment techniques have been instrumental in evaluating the impact of sources and secondary processes on the ambient PM2.5 concentrations in receptor areas. Positive Matrix Factorization (PMF) is now the most commonly used tool due to its ability to provide mixture resolution based on available PM2.5 compositional data. Sampling and analysis of PM2.5 was conducted in Cape Town, South Africa from April 2017 to April 2018. The resulting data were dispersion normalized to address the modifications of the source concentrations resulting from the varying dispersion conditions and thereby permit dispersion normalized PMF (DN-PMF) to be employed. DN-PMF quantified the 6 sources as 2-stroke vehicles/galvanizing industries (16.8%); soil/road dust (12.3%); sulphate/marine diesel (3.6%), traffic (15.7%), sea salt (21.8%), and heating/biomass burning/cooking (15.7%). In addition, air mass back trajectory analysis using the Hybrid Single Particle Lagrangian Trajectory (HYSPLIT) model identified long-range transport pathways to Cape Town. The HYSPLIT results showed air masses from the Atlantic SSW (6%), Atlantic SW (24%), Indian Ocean (31%), and Atlantic WSW (39%) influence air quality. The primary sources affected by the transport clusters were heating, 2-stroke vehicles/galvanizing, road and soil dust, and traffic emissions. These results show that reducing emissions from the local sources will improve air quality.

An ocean–wave–trajectory forecasting system for the eastern Baltic Sea: Validation against drifting buoys and implementation for oil spill modeling

We present the implementation and validation of OpenDrift, an open-source Lagrangian particle trajectory modeling framework for oil spill modeling in the coastal waters of Estonia in the Baltic Sea. The framework was coupled with ECMWF winds, NEMO-EST05 hydrodynamical model, and SWAN-EST wave model, and validated using six drift experiments from 2022. The sensitivity analysis revealed the importance of incorporating additional forcing factors, such as Stokes drift and currents, which generally improved the accuracy of the trajectory model compared to using wind alone. Nevertheless, the benefits of providing OpenDrift with, for example, the Stokes drift seemed to depend on whether currents are also included or not. The wind drift factors of the utilized drifters align closely with those commonly employed in oil spill modeling. Furthermore, the modeling results for hypothetical oil spills in severe weather conditions and high-risk regions emphasize the critical need for preparedness and rapid response strategies.

Characteristics of Water Vapor Transport for Extreme Summer Precipitation in the Eastern Southwest China and Its Impact Mechanism

To improve understanding of the characteristics of extreme summer rainfall and its water vapor transport in the eastern part of southwestern China (ESWC), this study analyzed data on daily precipitation from 118 meteorological stations in the ESWC from 1979 to 2020, as well as daily reanalysis data from ERA5 and daily reanalysis data from NCEP/NCAR. The study employed polynomial fitting, correlation, regression, clustering, and mixed single-particle Lagrangian trajectory (HYSPLITv5.0) modeling methods to simulate extreme summer precipitation and its water vapor transport characteristics in the ESWC and its possible formation mechanism. The results show that: (1) The contribution rate of extreme precipitation in the ESWC from 1979 to 2020 varied significantly on the interannual time scale. When the number of extreme precipitation days is high (low), the contribution rate of extreme precipitation is also high (low), while the contribution rate of general precipitation (the percentage of the sum of general precipitation to the total summer precipitation of that year) is often low (high). (2) When extreme precipitation occurs in the ESWC, compared with general precipitation, the high-level potential vortices are stronger, and the cold air from higher latitude is more likely to move southward. Meanwhile, the amount of water vapor input to the region is significantly larger than that of general precipitation. (3) There are four channels of water vapor sources in the ESWC during the period of extreme precipitation: the Bay of Bengal, the Arabian Sea, the western Pacific, and the northwest. The contribution of water vapor from the Bay of Bengal is the highest. The number of extreme summer precipitation days in the ESWC is significantly negatively correlated with the water vapor budget of the eastern boundary and positively correlated with Indian Ocean Basin-Wide (IOBW) index in the previous winter. (4) When the winter SST is high in the IOBW mode, it can cause the western Pacific subtropical high and the South Asian high to be stronger and shifted southward in summer, resulting in an increase in the number of extreme precipitation days in the ESWC.

An integrated monitoring system (IMS) for air quality: Observations of a unique ozone-exceedance event in Maryland

Ozone (O3) exceedance events in the mid-Atlantic states of the U.S. are typically observed during the summer months. They frequently result from local diurnal emissions of NO2 and VOCs where photochemical production of O3 is most prevalent. Occasionally, summertime O3-exceedances in the region are driven by long-range transport of pollutants or regional nocturnal dynamics. Discussed herein is the use of multiple observation platforms as an integrated monitoring system (IMS) to characterize the spatiotemporal confluence of long-range transport and nocturnal low-level jet induced O3-exceedances impacting air quality in the mid-Atlantic region. The IMS discussed is comprised of enhanced monitoring of vertically resolved observations of O3, atmospheric backscatter, winds, and state variables (e.g., temperature and relative humidity). This is in combination with surface-O3 observations, satellite-retrieved aerosol optical depth (by the moderate resolution imaging spectroradiometer (MODIS)), and simulations of air mass transport by NOAA's hybrid single-particle Lagrangian trajectory (HYSPLIT) model. Here we discuss an exceedance event that impacted the Mid-Atlantic region on May 20, 2021, where seven air monitoring sites in Maryland were driven by the combination of nocturnal low-level jet downmixing of pollutants and long-range transport of biomass-burning pollutants/O3-precursors. The most notable features of this case are this event's deviation from decadal O3-exceedance seasonality trends, the large spatial context of high surface-O3 impacting traditionally non-exceedance sites (sites that have not exceeded in years), the arrival of a well-mixed polluted air mass into the region, and nocturnal increases in surface-O3 induced by NLLJ.

Validation of OpenDrift-Based Drifter Trajectory Prediction Technique for Maritime Search and Rescue

Due to a recent increase in maritime activities in South Korea, the frequency of maritime distress is escalating and poses a significant threat to lives and property. The aim of this study was to validate a drift trajectory prediction technique to help mitigate the damages caused by maritime distress incidents. In this study, OpenDrift was verified using satellite drifter data from the Korea Hydrographic and Oceanographic Agency. OpenDrift is a Monte-Carlo-based Lagrangian trajectory modeling framework that allows for considering leeway, an important factor in predicting the movement of floating marine objects. The simulation results showed no significant differences in the performance of drift trajectory prediction when considering leeway using four evaluation methods (normalized cumulative Lagrangian separation, root mean squared error, mean absolute error, and Euclidean distance). However, leeway improved the performance in an analysis of location prediction conformance for maritime search and rescue operations. Therefore, the findings of this study suggest that it is important to consider leeway in drift trajectory prediction for effective maritime search and rescue operations. The results could help with future research on drift trajectory prediction of various floating objects, including marine debris, satellite drifters, and sea ice.

Optimizing Injection Locations Relaxes Altitude‐Lifetime Trade‐Off for Stratospheric Aerosol Injection

AbstractStratospheric Aerosol Injection (SAI) aims to offset some climate hazards by injecting aerosols into the stratosphere to reflect solar radiation. The lifetime of injected particles influences SAI's radiative efficacy—the ratio of radiative forcing to particle mass flux. We employ a Lagrangian trajectory model with particle sedimentation to simulate how background circulations influence the transport of passive particles (without microphysical growth) in the stratosphere and quantify sensitivities of particle lifetime to injection locations. At 20 km, optimizing injection locations can increase particle lifetime by >40%. Injection strategies can be constrained to maintain an interhemispheric balance of particle lifetime without significantly decreasing total lifetime. Generally, increasing injection altitude increases particle lifetime while also increasing costs and environmental impacts of deployment aircraft. Optimizing injection latitude and longitude can relax this altitude‐lifetime trade‐off by increasing lifetime without needing to increase altitude, which warrants further testing in global climate models with aerosol microphysics.

Tracking the pumice rafts from the Fukutoku-Okanoba submarine volcano with Satellites and a Lagrangian Particles trajectory model

On August 13th, 2021, the Fukutoku-Okanoba, a submarine volcano in the Northwest Pacific Ocean, erupted. Satellites detected various pumice rafts that had drifted westward to reach southern Japan over two months. To cope with the potential danger from pumice rafts, predicting their trajectories is crucial. Using a Lagrangian particle tracking model, the trajectories of the rafts were investigated. The model results showed strong sensitivity to the windage coefficient of pumice rafts, which is uncertain and could cause significant errors. An optimal windage coefficient was estimated by comparing the model results with satellite images using a skill score based on the distance between simulated particles and the nearest observed rafts divided by the travel distance of the particles. The optimal windage coefficients ranged between 2 and 3 % and produced pathways comparable to the observations from satellites. The simulation results showed that the pumice rafts moved from Fukutoku-Okanoba toward the Ryukyu Islands for approximately two months prior to being pushed by the north-easterly wind toward Taiwan against the Kuroshio. The methods presented here may become a valuable tool in managing coastal hazards due to diverse marine debris.

Effects of long-distance transport on O3 and secondary inorganic aerosols formation in Qingdao, China

This work presents a three-year continuous observation of ozone (O3) and fine particle (PM2.5), as well as their precursors in Qingdao, China, from September 2018 to August 2021. The annual concentrations of O3 and PM2.5 were measured as 35.9–44.5 ppbv and 31.6–34.2 μg m−3, respectively. Analysis of the interannual variations of O3 and PM2.5 concentration indicated the effective control measures for PM2.5 and O3 pollution in Qingdao in recent years. Nevertheless, we still observed 85 O3 episodes and 80 PM2.5 episodes during the whole observation. And it was found that the fraction of secondary inorganic aerosols (SIA) in PM2.5 mass significantly increased during PM2.5 episodes. By employing the Lagrangian photochemical trajectory model (LPTM), we investigated the roles of local production and regional transport (i.e., short-distance transport and long-distance transport) that play in the O3 and SIA formation in these episodes. The contribution of local production, short-distance transport, and long-distance transport to O3 concentration was calculated as 36%, 25%, and 39%, respectively. Long-distance transport seemed to play a more significant role in the SIA formation, accounting for 56% of SIA concentration. Moreover, the calculation results of the relative incremental reactivity (RIR) showed that both O3 and SIA formation were mainly affected by NOx and VOCs emissions during the process of long-distance transport. The observed unexpectedly high contribution of long-distance transport to O3 and SIA formation suggests that more extensive regional joint prevention and control policies on NOx and VOCs emissions are warranted to mitigate secondary air pollution.

Wind-driven advection across temperature gradients enhances iron-induced phytoplankton blooms in the Antarctic Polar Front

We demonstrate how the wind-driven Ekman transport enhances the advection and mixing of cells from the colder waters of the Surface Antarctic Waters from the south to the warmer waters of the northern Polar Front (PF) belt. This mechanism provides cells a mean ambient temperature near optimum levels for specific species and, ultimately, for community growth rates high enough to develop blooms under non-light limiting macronutrients and iron conditions. A Lagrangian trajectory model was constructed for tracking plankton cells as tracers forced by winds and surface currents. Depending on the region along the circumpolar front, increased winds can enhance this process across temperature gradients, and further accelerate such temperature-controlled growth. These results indicate that favorable temperature may enhance the growth rate even further when iron is sufficiently available, and thus have far-reaching implications for increased productivity in a future warming climate.

Improving trajectory calculations by FLEXPART 10.4+ using single-image super-resolution

Abstract. Lagrangian trajectory or particle dispersion models as well as semi-Lagrangian advection schemes require meteorological data such as wind, temperature and geopotential at the exact spatiotemporal locations of the particles that move independently from a regular grid. Traditionally, these high-resolution data have been obtained by interpolating the meteorological parameters from the gridded data of a meteorological model or reanalysis, e.g., using linear interpolation in space and time. However, interpolation is a large source of error for these models. Reducing them requires meteorological input fields with high space and time resolution, which may not always be available and can cause severe data storage and transfer problems. Here, we interpret this problem as a single-image super-resolution task. That is, we interpret meteorological fields available at their native resolution as low-resolution images and train deep neural networks to upscale them to a higher resolution, thereby providing more accurate data for Lagrangian models. We train various versions of the state-of-the-art enhanced deep residual networks for super-resolution (EDSR) on low-resolution ERA5 reanalysis data with the goal to upscale these data to an arbitrary spatial resolution. We show that the resulting upscaled wind fields have root-mean-squared errors half the size of the winds obtained with linear spatial interpolation at acceptable computational inference costs. In a test setup using the Lagrangian particle dispersion model FLEXPART and reduced-resolution wind fields, we find that absolute horizontal transport deviations of calculated trajectories from “true” trajectories calculated with un-degraded 0.5∘ × 0.5∘ winds are reduced by at least 49.5 % (21.8 %) after 48 h relative to trajectories using linear interpolation of the wind data when training on 2∘ × 2∘ to 1∘ × 1∘ (4∘ × 4∘ to 2∘ × 2∘) resolution data.

A Lagrangian Study of Tidal Advection of Mesospheric Water Vapor

AbstractPrevious studies have demonstrated that the diurnal variability of mesospheric water vapor involves a non‐linear interplay between dynamics, photochemistry, and catalytic processes. The inherent difficulty in measurement of winds and water vapor in this region makes the study of mesospheric water vapor complex. Here we study the sources of diurnal variability of mesospheric water vapor with the help of the specified dynamics configuration of the Whole Atmosphere Community Climate Model (WACCM) with thermosphere‐ionosphere eXtension (SD‐WACCM‐X) for typical equinox conditions during the solar minimum year 2009. Special emphasis is given to the advective transport by the migrating diurnal tide and its impact on the mean water vapor in the mesosphere. To isolate the contribution from advective and non‐advective processes, a Lagrangian parcel trajectory model is employed. Using this model, we incorporated the effect of both vertical and horizontal tidal advection on the diurnal variability of mesospheric water vapor. Our results indicate that the vertical tidal advection dominates the transport of water vapor in the upper mesosphere and that there is a significant contribution from horizontal tidal advection below 70 km. Among non‐advective processes, chemical loss of water vapor has the largest contribution. The Lagrangian parcel trajectory approach is further extended using the recently released version 2.07 of water vapor observations from the Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) instrument and dynamical tides from the empirical Climatological Tidal Model for the Thermosphere (CTMT). Differences in SD‐WACCM‐X and CTMT + SABER results are discussed in light of the veracity of SABER water vapor measurements.

Short Path Wind-Field Distance-Based Lagrangian Trajectory Model for Enhancing Atmospheric Dispersion Prediction Accuracy

Short Path Wind-Field Distance-Based Lagrangian Trajectory Model for Enhancing Atmospheric Dispersion Prediction Accuracy

Effects of spalled particles thermal degradation on a hypersonic flow field environment

The presence of spalled particles might affect the flow field and the material, thereby influencing the aerodynamic heating rates of the thermal protection system. In order to study the impact of particles on the flowfield, a two-way coupling is performed between a Lagrangian particle trajectory model and a hypersonic aerothermodynamics flow solver. Time-accurate solutions are computed for argon and air flowfields. Single-particle and multiple-particle simulations are performed and results are studied. The studies in the argon environment indicated that the particles start sublimating soon after they cross the shock and keep sublimating during the remainder of their time in the flow. In the air environment, the particles start releasing carbon vapor soon after their ejection, and the magnitudes of different carbon products vary based on the particle’s location relative to the shock. Similar behavior is observed for multiple-particle simulations but with an increase in the amount of released carbon vapor that diffuses and convects over a larger area. A single-particle simulation is also run by adding additional gas phase reactions, and it is found that the production rate of certain carbon products (C1 and CN) increases. A parametric study is conducted based on parameters that affect the motion of the particles. The results of the comprehensive study show that the carbon vapor released by spalled particles tends to change the composition of the flow field, particularly the upstream region of the shock, which affects the heat flux incident on the test sample.