Atmospheric Column Research Articles

A One‐Dimensional Model of Atmospheric Sputtering at Io Driven by S++ and O+

AbstractIo, the closest of Jupiter's four Galilean moons, suffers from intense ion bombardment from Jupiter's magnetosphere. The constant atmospheric erosion by energetic ion precipitation, referred to atmospheric sputtering, serves as an important mechanism of Io's atmospheric escape. This study is devoted to a state‐of‐the‐art study of atmospheric sputtering at Io, with the aid of constantly accumulated understandings of Io's space environment and atmospheric photochemistry, as well as the updated laboratory measurements. A Monte Carlo model is constructed to track the energy degradation of incident S++ and O+ and atmospheric recoils from which the sputtering yields of different atmospheric species are determined. Our calculations suggest a total escape rate of 3 × 1029 atom s−1 on Io, and SO2 is the dominant sputtered species. Further investigations reveal that S++ is the most efficient species for atmospheric sputtering on Io, and sputtering yields increase substantially with increasing incident ion mass, energy, and incidence angle. The model sensitivity to different influence factors is also discussed, including scattering angle distribution, atmospheric column density, proton precipitation, inelastic process, and surface sputtering, of which the former two dominate.

High-resolution US methane emissions inferred from an inversion of 2019 TROPOMI satellite data: contributions from individual states, urban areas, and landfills

Abstract. We quantify 2019 annual mean methane emissions in the contiguous US (CONUS) at 0.25° × 0.3125° resolution by inverse analysis of atmospheric methane columns measured by the Tropospheric Monitoring Instrument (TROPOMI). A gridded version of the US Environmental Protection Agency (EPA) Greenhouse Gas Emissions Inventory (GHGI) serves as the basis for the prior estimate for the inversion. We optimize emissions and quantify observing system information content for an eight-member inversion ensemble through analytical minimization of a Bayesian cost function. We achieve high resolution with a reduced-rank characterization of the observing system that optimally preserves information content. Our optimal (posterior) estimate of anthropogenic emissions in CONUS is 30.9 (30.0–31.8) Tg a−1, where the values in parentheses give the spread of the ensemble. This is a 13 % increase from the 2023 GHGI estimate for CONUS in 2019. We find emissions for livestock of 10.4 (10.0–10.7) Tg a−1, for oil and gas of 10.4 (10.1–10.7) Tg a−1, for coal of 1.5 (1.2–1.9) Tg a−1, for landfills of 6.9 (6.4–7.5) Tg a−1, for wastewater of 0.6 (0.5–0.7), and for other anthropogenic sources of 1.1 (1.0–1.2) Tg a−1. The largest increase relative to the GHGI occurs for landfills (51 %), with smaller increases for oil and gas (12 %) and livestock (11 %). These three sectors are responsible for 89 % of posterior anthropogenic emissions in CONUS. The largest decrease (28 %) is for coal. We exploit the high resolution of our inversion to quantify emissions from 70 individual landfills, where we find emissions are on median 77 % larger than the values reported to the EPA's Greenhouse Gas Reporting Program (GHGRP), a key data source for the GHGI. We attribute this underestimate to overestimated recovery efficiencies at landfill gas facilities and to under-accounting of site-specific operational changes and leaks. We also quantify emissions for the 48 individual states in CONUS, which we compare to the GHGI's new state-level inventories and to independent state-produced inventories. Our posterior emissions are on average 27 % larger than the GHGI in the largest 10 methane-producing states, with the biggest upward adjustments in states with large oil and gas emissions, including Texas, New Mexico, Louisiana, and Oklahoma. We also calculate emissions for 95 geographically diverse urban areas in CONUS. Emissions for these urban areas total 6.0 (5.4–6.7) Tg a−1 and are on average 39 (27–52) % larger than a gridded version of the 2023 GHGI, which we attribute to underestimated landfill and gas distribution emissions.

INTEGRATION OF STUMPF'S RATIO MODEL AND RANDOM FOREST FOR SATELLITE-DERIVED BATHYMETRY ESTIMATION

Abstract. The development of remote sensing for coastal and marine environment mapping has significantly enhanced our understanding of these ecosystems, enabling improved mitigation strategies against the impacts of human activities. However, remote sensing must consider the complex interplay of the atmosphere and water column. Ongoing research focuses on refining water column correction techniques, including Depth Invariant Indices (DII), Radiative Transfer models, and bathymetry models. This study specifically aims to enhance the Stumpf's Ratio model (SRM) for bathymetry by employing the Random Forest (RF) machine learning regression algorithm. The resulting bathymetry model, which incorporates visible bands from a Sentinel-2 MSI image, and their Stumpf's ratios, outperforms other methods, yielding the highest accuracy with RMSE and R2 values of 1.25 m and 0.854, respectively. This was followed by the multivariate SRM with RMSE and R2 values of 2.196 m and 0.554 respectively These findings demonstrate the promising potential of using RF machine learning regression with SRM for bathymetry modelling.

XCO2 Super-Resolution Reconstruction Based on Spatial Extreme Random Trees

Carbon dioxide (CO2) is currently the most harmful greenhouse gas in the atmosphere. Obtaining long-term, high-resolution atmospheric column CO2 concentration (XCO2) datasets is of great practical significance for mitigating the greenhouse effect, identifying and controlling carbon emission sources, and achieving carbon cycle management. However, mainstream satellite observations provide XCO2 datasets with coarse spatial resolution, which is insufficient to support the needs of higher-precision research. To address this gap, in this study, we integrate spatial information with the extreme random trees model and develop a new machine learning model called spatial extreme random trees (SExtraTrees) to reconstruct a 1 km spatial resolution XCO2 dataset for China from 2016 to 2020. The results indicate that the predictive ability of spatial extreme random trees is more stable and has higher fitting accuracy compared to other methods. Overall, XCO2 in China shows an increasing trend year by year, with the spatial distribution revealing significantly higher XCO2 levels in eastern coastal regions compared to western inland areas. The contributions of this study are primarily in the following areas: (1) Considering the spatial heterogeneity of XCO2 and combining spatial features with the advantages of machine learning, we construct the spatial extreme random trees model, which is verified to have high predictive accuracy. (2) Using the spatial extreme random trees model, we reconstruct high-resolution XCO2 datasets for China from 2016 to 2020, providing data support for carbon emission reduction and related decision making. (3) Based on the generated dataset, we analyze the spatiotemporal distribution patterns of XCO2 in China, thereby improving emission reduction policies and sustainable development measures.

A Fast Forward Modelling Method for Simulating Satellite Observations Using Observing Path Tracking

The higher the atmosphere is, the larger the deviations in atmospheric temperature and humidity are between the vertical column atmosphere above the cross-section of a satellite instrument and a ray’s trajectory from the cross-section to the satellite. In general, satellite instruments that observe using cross-orbit scanning result in the difference between the observed radiance and the simulations using this method becoming incrementally larger and larger as the cross-section moves to the edge of the satellite’s orbit. The deviations depend on the distance from the column to the ray trajectory and on the horizontal gradient of variables in the distance. In fact, the horizontal gradient of water vapour is larger than the gradient of temperature in clear scenarios, which could introduce an impact of temperature and water vapour on the simulated radiance of a satellite. In this study, a new method to simulate upgoing and downgoing radiation synchronously was developed, using the observing path tracking method. The conventional vertical initial atmospheric profile (Exp.1) and the profiles along the upgoing and downgoing rays of the satellite’s observation (Exp.2) were established, in order to simulate the observed radiance of MWHS-II of FY-3D using global numerical forecasts with resolutions of 15 km and 25 km. The results showed that, for channels in the oxygen and water vapour absorption line on the microwave spectrum, deviations of the two atmospheric profiles were larger at the scan edge (0.01 K) than those at the nadir (0.001 K), and were larger in the upper atmosphere than in the lower atmosphere. The deviation was usually negative in low-latitude regions and was positive in southern high-latitude regions. Such results were obtained in experiments using both the numerical forecast method with 15 km grids and the forecast method with 25 km grids. Deviations were analysed for representative channels at 118 GHz and 183 GHz. Then, the results indicated that bigger deviations between the two experiments were observed in the water vapour absorption line than in the oxygen absorption line in the microwave spectrum. In conclusion, this indicates that, because of the greater horizontal gradient of water vapour, the stronger localisation of water vapour makes the atmospheric profile along the satellite’s observing ray have more increments in the simulated radiance at the scan edge, compared to the atmospheric column profile.

Ocean-driven interannual variability in atmospheric CO2 quantified using OCO-2 observations and atmospheric transport simulations

Interannual variability (IAV) in the atmospheric CO2 growth rate is caused by variation in the balance between uptake by land and ocean and accumulation of anthropogenic emissions in the atmosphere. While variations in terrestrial fluxes are thought to drive most of the observed atmospheric CO2 IAV, the ability to characterize ocean impacts has been limited by the fact that most sites in the surface CO2 monitoring network are located on coasts or islands or within the continental interior. NASA’s Orbiting Carbon-Observatory 2 (OCO-2) mission has observed the atmospheric total column carbon dioxide mole fraction (XCO2) from space since September 2014. With a near-global coverage, this dataset provides a first opportunity to directly observe IAV in atmospheric CO2 over remote ocean regions. We assess the impact of ocean flux IAV on the OCO-2 record using atmospheric transport simulations with underlying gridded air-sea CO2 fluxes from observation-based products. We use three observation-based products to bracket the likely range of ocean air-sea flux contributions to XCO2 variability (over both land and ocean) within the GEOS-Chem atmospheric transport model. We find that the magnitude of XCO2 IAV generated by the whole ocean is between 0.08-0.12 ppm throughout the world. Depending on location and flux product, between 20-80% of the IAV in the simulations is caused by IAV in air-sea CO2 fluxes, with the remainder due to IAV in atmospheric winds, which modulate the atmospheric gradients that arise from climatological ocean fluxes. The Southern Hemisphere mid-latitudes and low-latitudes are the dominant ocean regions in generating the XCO2 IAV globally. The simulation results based on all three flux products show that even within the Northern Hemisphere atmosphere, Southern Hemisphere ocean fluxes are the dominant source of variability in XCO2. Nevertheless, the small magnitude of the air-sea flux impacts on XCO2 presents a substantial challenge for detection of ocean-driven IAV from OCO-2. Although the IAV amplitude arising from ocean fluxes and transport is 20 to 50% of the total observed XCO2 IAV amplitude of 0.4 to 1.6 ppm in the Southern Hemisphere and the tropics, ocean-driven IAV represents only 10% of the observed amplitude in the Northern Hemisphere. We find that for all three products, the simulated ocean-driven XCO2 IAV is weakly anti-correlated with OCO-2 observations, although these correlations are not statistically significant (p>0.05), suggesting that even over ocean basins, terrestrial IAV obscures the ocean signal.

Impact of Saharan dust outbreaks on short‐range weather forecast errors in Europe

AbstractMineral dust, the most abundant atmospheric aerosol by mass, interacts with radiation directly and alters cloud properties indirectly. Many operational numerical weather prediction models account for aerosol direct effects by using climatological mean concentrations and neglect indirect effects. This simplification may lead to shortcomings in model forecasts during outbreaks of Saharan dust towards Europe, when climatological mean dust concentrations deviate strongly from actual concentrations. This study investigates errors in model analyses and short‐range forecasts during such events. We investigate a pronounced dust event in March 2021 using the pre‐operational ICOsahedral Nonhydrostatic weather and climate model with Aerosols and Reactive Trace gases (ICON‐ART) with prognostic calculation of dust and the operational European Centre for Medium‐Range Weather Forecasts (ECMWF) Integrated Forecasting System (IFS) model, which deploys a dust climatology. We compare model analysis and forecast with measurements from satellite and in situ instruments. We find that inclusion of prognostic aerosol and direct radiative effects from dust improves forecasts of surface radiation during clear‐sky conditions. However, dust‐induced cirrus clouds are strongly underestimated, highlighting the importance of representing indirect effects adequately. These findings are corroborated by systematic quantification of forecast errors against satellite measurements. For this we construct an event catalogue with 49 dust days over Central Europe between January 2018 and March 2022. We classify model cells by simulated and observed cloudiness and simulated dustiness in the total atmospheric column. We find significant overestimations of brightness temperature for cases with dust compared with cases without dust. For surface shortwave radiation, we find median overestimations of 6.2% during cloudy conditions with dust optical depth greater than 0.1, however these are not significant compared with cloudy conditions without dust. Our findings show that the pre‐operational ICON‐ART and the operational IFS model still do not reproduce cloudiness adequately during events with Saharan dust over Central Europe. Missing implementations of prognostic dust, particularly of indirect effects on cloud formation, lead to significant underestimations of cloudiness and potentially overestimations of surface radiation.

Combining Cn2 models to forecast the optical turbulence at Paranal

ABSTRACT We applied a combination of models to improve the forecasts of refractive index structure coefficient ($C^{2}_{n}$) profiles and seeing at Paranal using high-temporal and spatial resolution simulations with the Weather Research and Forecasting model. We assessed our method with Stereo−SCIntillation Detection And Ranging observations from several nights between 2016 and 2018. The combined approach consists of the turbulent kinetic energy-based model to estimate the $C^{2}_{n}$ profile within the boundary layer and another model for the free atmosphere. We tested the Dewan, Jackson−Dewan, and Gladstone models. The implementation of the combined method gives better results than those obtained using each model separately for the whole atmospheric column. However, a much better agreement with observations is obtained when we use a calibration method to improve the results. Calibrated seeing forecasts at Paranal showed a root mean squared error of 0.30 arcsec and a bias around −0.1 arcsec for all the nights of 2017 and 2018, which are similar to previous results obtained at Paranal during the same nights. Due to its performance and rapid execution, the proposed methodology could be implemented as an operational tool to forecast the $C^{2}_{n}$ profiles and the seeing at Paranal and potentially over other astronomical sites around the world.

Parameterization of the Elevated Convection With a Unified Convection Scheme (UNICON) and Its Impacts on the Diurnal Cycle of Precipitation

AbstractTo improve the simulation of nocturnal precipitation, we develop a parameterization for the elevated convection that is launched from the level of maximum grid‐mean moist static energy, when grid‐mean vertical flow is upward at the launching interface. The parameterized elevated convection is forced by both cold pool‐driven and partially resolved external mesoscale organized flows. Properties of the external mesoscale flow are estimated from grid‐mean vertical velocity and three dimensional advection tendencies of temperature and moisture. The new parameterization is implemented into a unified convection scheme (UNICON) and tested for both the single‐column case at the Southern Great Plain (SGP) site in US and in global simulations. At the SGP site, the parameterized elevated convection strengthens nocturnal convection, increases nocturnal precipitation, and better simulates the observed diurnal cycle of precipitation. It appears that the elevated and surface‐based convections interact with each other by stabilizing the atmospheric column, which affects subsequent convection. Global simulation shows that the elevated convection mostly occurs over the continents during the night, and also over the oceanic mid‐latitude warm air advection and storm track regions during summer. Without degrading global mean climate, the elevated convection improves the simulation of nocturnal precipitation in the mid‐US but simulates somewhat strong nocturnal precipitation in northeastern Asia. It seems that a key to simulate observed nocturnal precipitation is to appropriately parameterize the impacts of external organized flow on the elevated convection.

Lessons Learned from the Updated GEWEX Cloud Assessment Database

AbstractSince the first Global Energy and Water Exchanges cloud assessment a decade ago, existing cloud property retrievals have been revised and new retrievals have been developed. The new global long-term cloud datasets show, in general, similar results to those of the previous assessment. A notable exception is the reduced cloud amount provided by the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) Science Team, resulting from an improved aerosol–cloud distinction. Height, opacity and thermodynamic phase determine the radiative effect of clouds. Their distributions as well as relative occurrences of cloud types distinguished by height and optical depth are discussed. The similar results of the two assessments indicate that further improvement, in particular on vertical cloud layering, can only be achieved by combining complementary information. We suggest such combination methods to estimate the amount of all clouds within the atmospheric column, including those hidden by clouds aloft. The results compare well with those from CloudSat-CALIPSO radar–lidar geometrical profiles as well as with results from the International Satellite Cloud Climatology Project (ISCCP) corrected by the cloud vertical layer model, which is used for the computation of the ISCCP-derived radiative fluxes. Furthermore, we highlight studies on cloud monitoring using the information from the histograms of the database and give guidelines for: (1) the use of satellite-retrieved cloud properties in climate studies and climate model evaluation and (2) improved retrieval strategies.

Comparative Analysis and High−Precision Modeling of Tropospheric CH4 in the Yangtze River Delta of China Obtained from the TROPOMI and GOSAT

Remote sensing satellite monitoring involving the use of shortwave infrared (SWIR) solar backscatter radiation to measure atmospheric CH4 column concentrations provides wide−ranging and accurate data for quantitatively determining atmospheric CH4 emissions and is highly important for human studies of atmospheric composition and environmental protection. The ESA−launched Sentinel−2 satellite equipped with a tropospheric monitoring instrument (TROPOMI) can provide the concentration of CH4 columns in every piece of the global atmosphere every day. However, these data may be affected by surface albedo, SWIR, aerosols, cirrus cloud scattering, and other factors. The greenhouse gas observing satellite (GOSAT) launched by Japan has fairly accurate data that are minimally affected by the aforementioned factors; however, its data density is much less than that of the TROPOMI. In this study, we propose a CH4 model that combines the TROPOMI and GOSAT data. We construct the model by analyzing the data from the TROPOMI and GOSAT at the same location at the same time. Then, we apply the proposed model to a certain location at a certain time with TROPOMI data but without GOSAT data to obtain a large range of high−precision CH4 data. The most developed urban agglomeration in the Yangtze River Delta in China was selected for model construction and the correlations between the TROPOMI and GOSAT data and their spatial and temporal trends were analyzed. First, we analyzed the CH4 concentrations in the same area measured by both models. The results revealed a high degree of temporal and spatial correlation in the YRD region. The correlation coefficient reached 0.71 in the metropolitan area of the YRD. At the small−city scale, the correlation is much more significant, with the correlation reaching 0.80, 0.79, and 0.71 for Nanjing, Shanghai, and Ningbo, respectively. The most accurate model was screened through comparative construction to calibrate the TROPOMI data and high−precision and high−coverage CH4 concentration information was obtained for the study area. Five models (linear model, quadratic term model, cubic term model, lognormal model, and logistic model) were used to select the best−fitting model. The magnitudes of the differences in the CH4 concentrations calculated by each model were compared. The final results showed that the linear model, as the prediction model, had the highest accuracy, with a coefficient of determination (R22) of 0.542. To avoid the specificity of the constructed model, we used the same method in several simulations to validate. The coefficient of determination of the model constructed with different stochastic data was greater than 0.5. Subsequently, we used Nanjing as the study area and applied the same method to construct the model. The coefficient of determination of the model (R22) was approximately 0.601. The model constructed in this research can be used not only for data conversion between the same products from different sensors to obtain high−precision data products but also for calibrating newly developed satellite data products that utilize mature data products.

An Analytic Theory of Near-Surface Relative Humidity over Land

Abstract There is no simple explanation for the spatial structure of near-surface relative humidity over land. We present a diagnostic theory for zonally and temporally averaged near-surface relative humidity (RH) over land based on energy budgets of an atmospheric column in radiative–convective equilibrium. The theory analytically relates RH to the surface evaporative fraction (EF), has no calibrated parameters, and is quantitatively accurate when compared with RH from a reanalysis, and with cloud-permitting simulations over an idealized land surface. The theory is used to answer two basic questions. First, why is RH never especially low (e.g., 1%)? The theory shows that established lower bounds on EF over land and ocean are equivalent to lower bounds on RH that preclude particularly low values, at least for conditions typical of the modern Earth. Second, why is the latitudinal profile of RH over land shaped like the letter W, when both specific humidity and saturation specific humidity essentially decline monotonically from the equator to the poles? The theory predicts that the latitudinal profile of RH should look more like that of water stored in the soil (which also exhibits a W-shaped profile) than in the air (which does not).

Relative Contributions of Large‐Scale Atmospheric Circulation Dynamics and Anthropogenic Warming to the Unprecedented 2022 Yangtze River Basin Heatwave

AbstractThe unprecedented 2022 Yangtze River Basin (YRB) heatwave is a threat to human society and natural ecology, so the understanding of its underlying drivers is critical to regional climate adaptation and resilience. Here we conducted a multi‐method attribution analysis on the contribution of atmospheric circulation change and anthropogenic impacts to the occurrence probability and intensity of this extreme heatwave. Based on the nonstationary statistical analysis, the 2022 YRB heatwave is a 1‐in‐900‐year event and a 1‐in‐110‐year event with and without considering the 2022 YRB heatwave in the fitting, respectively. The large‐scale meteorological condition analysis shows that the 2022 YRB heatwave is featured with an anomalous high‐pressure system that favors a hot and dry atmospheric column, overlaid by anomalous subsidence and clear skies which leads to warming and greater solar heating. The ensemble constructed circulation analogue analyses show that the circulation anomaly fails to explain the observed 2022 YRB SAT anomalies fully. Specifically, 46% (0.132 ± 0.027°C decade−1) of the observed SAT trend during 1979–2022 (0.290 ± 0.048°C decade−1) is caused by anthropogenic warming and the associated thermodynamic feedback, while the remaining 54% (0.157 ± 0.038°C decade−1) of the trend is caused by changes in the large‐scale atmospheric circulation. Our findings on changes in atmospheric circulation patterns associated with YRB heatwave and anthropogenic contributions to YRB heatwave could provide valuable information for climate adaptation and mitigation strategies in the context of a warming climate.

A 7-year record of vertical profiles of radar measurements and precipitation estimates at Dumont d'Urville, Adélie Land, East Antarctica

Abstract. Studying precipitation falling over Antarctica is crucial as snowfall represents the main water input term for the polar cap. However, precipitation observations still remain scarce – and, more particularly, in the atmospheric column – due to numerous experimental issues related to the white continent. This paper aims at helping to close this observation gap by presenting 7 years of Micro Rain Radar (Metek MRR-2) data at the Dumont d'Urville station in coastal Adélie Land, East Antarctica. Statistics are calculated on three radar variables (equivalent reflectivity, mean Doppler velocity and signal-to-noise ratio (SNR)) to outline the main characteristics of the radar dataset. Seasonal and interannual variabilities are also investigated, but no significant temporal trends are detected, except for the seasonal mean Doppler velocity, which is higher in summer and lower in winter. We then use the snowfall rate (S) data from a collocated snow gauge to estimate the MRR precipitation profile from the radar equivalent reflectivity (Ze) through a locally derived Ze–S relation. We find the relation Ze=43.3S0.88. The processing method used to obtain this relation, data quality and uncertainty considerations are discussed in the paper. In order to give an example of application of the dataset, a brief statistical comparison of the MRR precipitation rate along the vertical with model data from the ERA5 reanalysis and the LMDZ climate model is performed, which notably shows that models underestimate heavy precipitation events. All datasets are available on the PANGAEA database with the associated DOI: https://doi.org/10.1594/PANGAEA.962727 (Wiener et al., 2023).

Comparison of Short‐Term Cloud Feedbacks at Top of the Atmosphere and the Surface in Observations and AMIP6 Models

AbstractWe compared short‐term cloud feedback, defined at the top of the atmosphere (TOA), the atmospheric column (ATM), and the surface (SFC), between observations and models participating in Atmospheric Model Intercomparison Project Phase 6 (AMIP6) for the period 2000–2014. The globally averaged net cloud feedbacks observed at TOA, ATM, and SFC are −0.06 ± 0.63, −0.17 ± 0.70, and 0.11 ± 0.81 W m−2 K−1, respectively. While most models produced TOA cloud feedbacks that agreed with the observations within uncertainty ranges, the intermodel spread at SFC and within ATM was relatively larger. This demonstrates that models are diverse in how their TOA feedback is distributed between ATM and SFC. Because short‐term cloud feedback is mainly driven by El Niño–Southern Oscillation (ENSO), the global‐mean cloud feedback was further decomposed into components from the ENSO and non‐ENSO regions. Results show that cloud feedback in these two regions tends to be inversely related. Compared to observations, almost all models overestimated the longwave cloud feedback in the ENSO region due to the overestimation of cloud amount changes for high‐topped clouds. For these models, it is the offset between deviations in ENSO and non‐ENSO regions that leads to the overall agreement of global mean with observations. Sensitivity tests show that the main conclusions still hold when alternative kernels are used in estimating cloud feedback.

Assessing the Monthly Trends in Precipitable Water Vapor over the Indian Subcontinent

This study estimates the trends in precipitable water vapor (PWV), atmospheric moisture budget (AMB), and the factors influencing them: air temperature, evapotranspiration (ET), convective available potential energy (CAPE), and vertical velocity (Omega) over the Indian subcontinent using ERA5 reanalysis data sets between 1980 and 2020. PWV is examined across three atmospheric layers (1000–850 hPa: lower layer; 850–500 hPa: middle layer; 500–300 hPa: upper layer), and the entire atmospheric column (EAC; 1000–300 hPa). The observed PWV trends exhibit variability within the EAC, ranging from −0.53 to 1.25 mm/decade across the study area, with the middle layer showing the most pronounced variation (–0.44 to 0.83 mm/decade), followed by the lower layer (0.10 to 0.45 mm/decade), and the upper layer (–0.02 to 0.23 mm/decade). These fluctuations in PWV are attributed to changes in air temperature, ET, CAPE, and Omega. This investigation, however, underscores the necessity of delving into the impacts of these influencing factors on PWV using finer resolution data, to enhance our comprehension of its spatial and temporal dynamics in the region. Furthermore, the annual AMB analysis reveals a declining trend in the study region. These findings collectively contribute to the understanding of regional water-energy cycles and the recent shifts in atmospheric dynamics.

Characterizing atmospheric aerosols using polarimetry and shadow hiding

Unpolarized sunlight is scattered by aerosols acquiring partial linear polarization. By aiming a ground-based detector vertically upward, it can record the polarimetric response of aerosols that are illuminated by the Sun. As the Sun sets, a portion of the sky is shadowed and the polarimetric response of the aerosols in the unshadowed region can be measured. This provides a means of scanning different portions of the atmospheric column with time. By comparing the measured polarimetric response with that of model agglomerated debris particles we can place constraints on the sizes and chemical composition of the aerosols in different portions of this column. We conducted a survey over 24 different epochs from April 2021 to December 2022, consisting of approximately 600 measurements of polarization of the atmosphere in twilight at the Ussuriysk Astrophysical Observatory. We found that most of the measurements correspond with water-ice particles or dust. However, on some occasions organic carbon dominated the measurements. These epochs correspond with increased fire seats in the region.

Ozone Trend Analysis in Natal (5.4°S, 35.4°W, Brazil) Using Multi-Linear Regression and Empirical Decomposition Methods over 22 Years of Observations

Ozone plays an important role in the Earth’s atmosphere. It is mainly formed in the tropical stratosphere and is transported by the Brewer–Dobson Circulation to higher latitudes. In the stratosphere, ozone can filter the incoming solar ultraviolet radiation, thus protecting life at the surface. Although tropospheric ozone accounts for only ~10%, it is a powerful GHG and pollutant, harmful to the health of the environment and living beings. Several studies have highlighted biomass burning as a major contributor to the tropospheric ozone budget. Our study focuses on the Natal site (5.40°S, 35.40°W, Brazil), one of the oldest ozone-observing stations in Brazil, which is expected to be influenced by fire plumes in Africa and Brazil. Many studies that examined ozone trends used the total atmospheric columns of ozone, but it is important to assess ozone separately in the troposphere and the stratosphere. In this study, we have used radiosonde ozone profiles and daily TCO measurements to evaluate the variability and changes of both tropospheric and stratospheric ozone separately. The dataset in this study comprises daily total columns of colocalized ozone and weekly ozone profiles collected between 1998 and 2019. The tropospheric columns were estimated by integrating ozone profiles measured by ozone sondes up to the tropopause height. The amount of ozone in the stratosphere was then deduced by subtracting the tropospheric ozone amount from the total amount of ozone measured by the Dobson spectrometer. It was assumed that the amount of ozone in the mesosphere is negligible. This produced three distinct time series of ozone: tropospheric and stratospheric columns as well as total columns. The present study aims to apply a new decomposition method named Empirical Adaptive Wavelet Decomposition (EAWD) that is used to identify the different modes of variability present in the analyzed signal. This is achieved by summing up the most significant Intrinsic Mode Functions (IMF). The Fourier spectrum of the original signal is broken down into spectral bands that frame each IMF obtained by the Empirical Modal Decomposition (EMD). Then, the Empirical Wavelet Transform (EWT) is applied to each interval. Unlike other methods like EMD and multi-linear regression (MLR), the EAWD technique has an advantage in providing better frequency resolution and thus overcoming the phenomenon of mode-mixing, as well as detecting possible breakpoints in the trend mode. The obtained ozone datasets were analyzed using three methods: MLR, EMD, and EAWD. The EAWD algorithm exhibited the advantage of retrieving ~90% to 95% of ozone variability and detecting possible breakpoints in its trend component. Overall, the MRL and EAWD methods showed almost similar trends, a decrease in the stratosphere ozone (−1.3 ± 0.8%) and an increase in the tropospheric ozone (+4.9 ± 1.3%). This study shows the relevance of combining data to separately analyze tropospheric and stratospheric ozone variability and trends. It highlights the advantage of the EAWD algorithm in detecting modes of variability in a geophysical signal without prior knowledge of the underlying forcings.

Aerosol optical depth disaggregation: toward global aerosol vertical profiles

Abstract Aerosol-cloud interactions constitute the largest source of uncertainty in assessments of anthropogenic climate change. This uncertainty arises in part from the difficulty in measuring the vertical distributions of aerosols, and only sporadic vertically resolved observations are available. We often have to settle for less informative vertically aggregated proxies such as aerosol optical depth (AOD). In this work, we develop a framework for the vertical disaggregation of AOD into extinction profiles, that is, the measure of light extinction throughout an atmospheric column, using readily available vertically resolved meteorological predictors such as temperature, pressure, or relative humidity. Using Bayesian nonparametric modeling, we devise a simple Gaussian process prior over aerosol vertical profiles and update it with AOD observations to infer a distribution over vertical extinction profiles. To validate our approach, we use ECHAM-HAM aerosol-climate model data which offers self-consistent simulations of meteorological covariates, AOD, and extinction profiles. Our results show that, while very simple, our model is able to reconstruct realistic extinction profiles with well-calibrated uncertainty, outperforming by an order of magnitude the idealized baseline which is typically used in satellite AOD retrieval algorithms. In particular, the model demonstrates a faithful reconstruction of extinction patterns arising from aerosol water uptake in the boundary layer. Observations however suggest that other extinction patterns, due to aerosol mass concentration, particle size, and radiative properties, might be more challenging to capture and require additional vertically resolved predictors.

Inversion and Validation of Atmospheric CO2 Column Concentration Inversion of Spaceborne IPDA Lidar Based on Atmospheric Environment Monitoring Satellite

Inversion and Validation of Atmospheric CO2 Column Concentration Inversion of Spaceborne IPDA Lidar Based on Atmospheric Environment Monitoring Satellite