Search

We report the recovery and re-calibration of an extensive dataset of vertical profile measurements of the 1963/64 stratospheric aerosol layer measured from a two-site searchlight measurement facility at White Sands missile base and Sacramento Peak observatory, in New Mexico, US. The recovered dataset comprises 105 profiles of 550nm aerosol extinction (βp(z)) and is part of a longer program of measurements with the US Air Force Cambridge Research Laboratories (AFCRL) searchlight facility that began in February 1963. The recovered series of βp(z) profiles span the 13-month period December 1963 to December 1964 and provide a unique record of the altitude and vertical extent of the Northern Hemisphere dispersed portion of the aerosol cloud from the March 1963 Agung volcanic eruption. The data recovery methodology involved first re-digitizing the 105 original βp(z) profiles (βp(z)) from individual Figures within an AFCRL research report (Elterman, 1966a). The re-calibration involves inverting the original equation used to compute βp(z) in Elterman (1966a; 1966b) to retrieve a normalized detector response (Er(z) Er(35)) profile for each of the 105 re-digitized βp 0(z) profiles. An iterative procedure was then used to compute the re-calibration βp(z) profiles (βp Ra(z)), with the molecular extinction profile calculated with the corresponding daily molecular extinction profile calculated from local soundings, rather than the US Standard Atmosphere 1962 in the original dataset. Two-way molecular and aerosol transmittance corrections are applied using the MODTRAN (MODerate resolution atmospheric TRANsmission) code in transmission mode, applying a best-estimate aerosol phase function calculated from measurements, applied for the entire 2.76 to 35.2 km column. For the tropospheric aerosol transmittance, the AERONET aerosol phase function from White Sands High Energy Laser Systems Test Facility (HELSTF) was applied (2.76 to 10.7 km), a separate stratospheric phase function applied between 11.2 and 35.2 km, calculated from a set of particle size distributions measured by the U-2 high-altitude aircraft over a US region in the vicinity of White Sands in early 1964. Errors were estimated taking as a reference the errors determined in the computation of βp 0(z). Using available tabulated data from the original procedure the errors in the re-digitalization of βp 0(z) and in the retrieval of the Er(z) Er(35) procedures were calculated and later added to the original estimates. Both the βp R(z) and the stratospheric aerosol optical depth magnitudes showed higher magnitudes than βp 0(z) and the original stratospheric aerosol optical depth, however their magnitudes show a reasonable agreement with other contemporary observations. Both the original and recalibration datasets are being submitted to PANGAEA data repository for its storage and public access. [ABSTRACT FROM AUTHOR]

We report the recovery and processing methodology of the first ever multi-year lidar dataset of the stratospheric aerosol layer. A Q-switched Ruby lidar measured 66 vertical profiles of 694nm attenuated backscatter at Lexington, Massachusetts between January 1964 and August 1965, with an additional 9 profile measurements conducted from College, Alaska during July and August 1964. We describe the processing of the recovered lidar backscattering ratio profiles to produce mid-visible (532nm) stratospheric aerosol extinction profiles (sAEP532) and stratospheric aerosol optical depth (sAOD532) measurements, utilizing a number of contemporary measurements of several different atmospheric variables. Stratospheric soundings of temperature, and pressure generate an accurate local molecular backscattering profile, with nearby ozone soundings determining the ozone absorption, those profiles then used to correct for two-way ozone transmittance. Two-way aerosol transmittance corrections were also applied based on nearby observations of total aerosol optical depth (across the troposphere and stratosphere) from sun photometer measurements . We show the two-way transmittance correction has substantial effects on the retrieved sAEP532 and sAOD532, calculated without the corrections resulting in substantially lower values of both variables, as it was not applied in the original processing producing the lid ar scattering ratio profiles we rescued. The combined transmittance corrections causes the aerosol extinction to increase by 67 % for Lexington and 27 % for Fairbanks, for sAOD532 the increases 66 % and 26 % respectively. Comparing the magnitudes of the aerosol extinction and sAOD with the few contemporary available measurements reported show a better agreement in the case of the two way transmittance corrected values. The sAEP and sAOD timeseries at Lexington show a surprisingly large degree of variability, three periods where the stratospheric aerosol layer had suddenly elevated optical thickness, the highest sAOD532 of 0.07 measured at the end of March 1965. The two other periods of enhanced sAOD532 are both two-month periods where the lidars show more than 1 night where ret rieved sAOD532 exceeded 0.05: in January and February 1964 and November and December 1964. Interactive stratospheric aerosol model simulations of the 1963 Agung cloud illustrate that although substantial variation in mid-latitude sAOD532 is expected from the seasonal cycle in the Brewer-Dobson circulation, the Agung cloud dispersion will have caused much slower increase than the more episodic variations observed, with also different timing, elevated optical thickness from Agung occurring in winter and spring. The abruptness and timing of the steadily increasing sAOD from January to July 1965 suggests this variation was from a different source than Agung, possibly from one or both of the two VEI3 eruptions that occurred in 1964/65: Trident, Alaska and Vestmannaeyjar, Heimey, south of Iceland. A detailed error analysis of the uncertainties in each of the variables involved in the processing chain was conducted, relative errors of 54 % for Fairbanks and 44 % Lexington for the uncorrected sAEP532, corrected sAEP532 of 61 % and 64 % respectively. The analysis of the uncertainties identified variables that, with additional data recovery and reprocessing could reduce these relative error levels. Data described in this work are available at https://doi.pangaea.de/10.1594/PANGAEA.922105 (Dataset in Review) (Antuña-Marrero et al., 2020a). [ABSTRACT FROM AUTHOR]

A key limitation of volcanic forcing datasets for the Pinatubo period, is the large uncertainty that remains with respect to the extent of the optical depth of the Pinatubo aerosol cloud in the first year after the eruption, the saturation of the SAGE-II instrument restricting it to only be able to measure the upper part of the aerosol cloud in the tropics. Here we report the recovery of stratospheric aerosol measurements from two ship-borne lidars, both of which measured the tropical reservoir of volcanic aerosol produced by the June 1991 Mount Pinatubo eruption. The lidars were on-board two Soviet vessels, each ship crossing the Atlantic, their measurement datasets providing unique observational transects of the Pinatubo cloud across the tropics from Europe to the Caribbean (~ 40° N to 8° N) from July to September 1991 (the Prof Zubov ship) and from Europe to south of the Equator (8° S to ~ 40° N) between January and February 1992 (the Prof Vizeship). Our philosophy with the data recovery is to follow the same algorithms and parameters appearing in the two peer-reviewed articles that presented these datasets in the same issue of GRL in 1993, and here we provide all 48 lidar soundings made from the Prof. Zubov, and 11 of the 20 conducted from the Prof. Vize, ensuring we have reproduced the aerosols backscatter and extinction values in the Figures of those two papers. These original approaches used thermodynamic properties from the CIRA-86 standard atmosphere to derive the molecular backscattering, vertically and temporally constant values applied for the aerosol backscatter to extinction ratio and the correction factor of the aerosols backscattering wavelength dependence. We demonstrate this initial validation of the recovered stratospheric aerosol extinction profiles, providing full details of each dataset in this paper's Supplement S1, the original text files of the backscatter ratio, the calculated aerosols backscatter and extinction profiles. We anticipate the data providing potential new observational case studies for modelling analyses, including a 1-week series of consecutive soundings (in September 1991) at the same location showing the progression of the entrainment of part of the Pinatubo plume into the upper troposphere and the formation of an associated cirrus cloud. The Zubov lidar dataset illustrates how the tropically confined Pinatubo aerosol cloud transformed from a highly heterogeneous vertical structure in August 1991, maximum aerosol extinction values around 19 km for the lower layer and 23-24 for the upper layer, to a more homogeneous and deeper reservoir of volcanic aerosol in September 1991. We encourage modelling groups to consider new analyses of the Pinatubo cloud, comparing to the recovered datasets, with the potential to increase our understanding of the evolution of the Pinatubo aerosol cloud and its effects. Data described in this work are available at https://doi.pangaea.de/10.1594/PANGAEA.912770 (Antuña-Marrero et al., 2020). [ABSTRACT FROM AUTHOR]

Accurate quantification of the effects of volcanic eruptions on climate is a key requirement for better attribution of anthropogenic climate change. Here we use the UM-UKCA composition-climate model to simulate the atmospheric evolution of the volcanic aerosol clouds from the three largest eruptions of the 20th century: 1963 Agung, 1982 El Chichón and 1991 Pinatubo. The model has interactive stratospheric chemistry and aerosol microphysics, with coupled aerosol-radiation interactions for realistic composition-dynamics feedbacks. Our simulations align with the design of the Interactive Stratospheric Aerosol Model Intercomparison (ISA-MIP) "Historical Eruption SO2 Emissions Assessment". For each eruption, we perform 3-member ensemble model experiments with upper, mid-point and lower estimates for SO2 emission, each initialised to a meteorological state to match the observed phase of the quasi-biennial oscillation (QBO) at the times of the eruptions. We assess how each eruption's emitted SO2 evolves into a tropical reservoir of volcanic aerosol and analyse the subsequent dispersion to mid-latitudes. We compare the simulations to the three volcanic forcing datasets used in historical integrations for the two most recent Coupled Model Intercomparison Project (CMIP) assessments: the Global Space-based Stratospheric Aerosol Climatology (GloSSAC) for CMIP6, and the Sato et al. (1993) and Ammann et al. (2003) datasets used in CMIP5. We also assess the vertical extent of the volcanic aerosol clouds by comparing simulated extinction to Stratospheric Aerosol and Gas Experiment II (SAGE-II) v7.0 satellite aerosol data (1985-1995) for Pinatubo and El Chichón, and to 1964-65 northern hemisphere ground-based lidar measurements for Agung. As an independent test for the simulated volcanic forcing after Pinatubo, we also compare to the shortwave (SW) and longwave (LW) Top-of-the-Atmosphere flux anomalies measured by the Earth Radiation Budget Experiment (ERBE) satellite instrument. For the Pinatubo simulations, an injection of 10 to 14 Tg SO2 gives the best match to the High Resolution Infrared Sounder (HIRS) satellite-derived global stratospheric sulphur burden, with good agreement also to SAGE-II mid-visible and near-infrared extinction measurements. This 10-14 Tg range of emission also generates a heating of the tropical stratosphere that is comparable with the temperature anomaly seen in the ERA-Interim reanalyses. For El Chichón the simulations with 5 Tg and 7 Tg SO2 emission give best agreement with the observations. However, these runs predict a much deeper volcanic cloud than present in the CMIP6 data, with much higher aerosol extinction than the GloSSAC data up to October 1984, but better agreement during the later SAGE-II period. For 1963 Agung, the 9 Tg simulation compares best to the forcing datasets with the model capturing the lidar-observed signature of peak extinction descending from 20 km in 1964 to 16 km in 1965. Overall, our results indicate that the downward adjustment to previous SO2 emission estimates for Pinatubo as suggested by several interactive modelling studies is also needed for the Agung and El Chichón aerosol clouds. This strengthens the hypothesis that interactive stratospheric aerosol models may be missing an important removal or redistribution process (e.g. effects of co-emitted ash) which changes how the tropical reservoir of volcanic aerosol evolves in the initial months after an eruption. Our analysis identifies potentially important inhomogeneities in the CMIP6 dataset for all three periods that are hard to reconcile with variations predicted by the interactive stratospheric aerosol model. We also highlight large differences between the CMIP5 and CMIP6 volcanic aerosol datasets for the Agung and El Chichón periods. Future research should aim to reduce this uncertainty by reconciling the datasets with additional stratospheric aerosol observations. [ABSTRACT FROM AUTHOR]

In the present study, we report the first comparison of the aerosol properties measured with sun photometer at Camagüey, Cuba, with the MODerate resolution Imaging Spectroradiometer (MODIS) instruments on Terra and Aqua satellites. We compared the aerosol optical depth at 550 nm (AOD) and the Ångström Exponent (AE) from the sun photometer for the period 2008 to 2014 with the same variables measured by both MODIS instruments, that are spatially and temporally coincident. The comparison includes AOD derived with both Deep Blue (DB) and Dark Target (DT) algorithms from MODIS Collection 6. The AOD derived with DT algorithm for Terra and Aqua agrees better with AOD from the sun photometer than the AOD derived with DB. Additionally there is little difference between AOD from both satellite instruments, when they are compared with sun photometer AOD, allowing to combine AOD from Terra and Aqua for more comprehensive climatological statistics. The comparison of the AE showed similar results with reports in the literature about the little skills of the current DT and DB algorithms for its retrieval. In addition, we report the comparison of the broadband AOD (BAOD) from pyrheliometer measurements located at Camagüey site and other three meteorological stations along Cuba, with AOD measurements from the sun photometer and from MODIS onboard Terra and Aqua. The comparison of the BAOD from the four sites as a whole with coincident AOD from MODIS onboard Terra and Aqua showed similar results than the ones of the comparison between the sun photometer AOD and the AOD from the two satellite instruments. In the comparison between the BAOD and the AOD at each one of the eight individual sun photometer wavelengths, the results improve in the spectral range 400 to 675 nm, with the best result at 500 nm. The BAOD typical uncertainty ranges from 0.04 to 0.06 at this band. The results from the BAOD comparisons demonstrate its reliability for characterizing AOD at sites with no sun photometer and for extending backward in time AOD estimates. [ABSTRACT FROM AUTHOR]

The molecular backscatter signal is necessary in aerosol lidars to determine the aerosol backscatter and extinction. It is derived using atmospheric temperature and pressure profiles from simultaneous soundings or climatological values from radiosonde measurements. In cases when no radiosonde measurements are available at or near the lidar stations the U.S. Standard Atmosphere is commonly used. Although the atmospheric reanalysis can also be used to obtain the atmospheric temperature and pressure profiles at any station, we found no reports of its use for this purpose. Therefore, the goal of this research is to compare molecular backscatter coefficients derived from radiosonde measurements with ones derived using the ERA-Interim reanalysis and the U.S. Standard Atmosphere 1976 (USSA-1976), for the locations of the lidar stations from the Latin America Lidar Network (LALINET). The differences between radiosonde measurements and ERA-Interim temperature profiles are smaller than the differences between radiosonde measurements and the USSA-1976 temperature profiles. In many stations, the differences between the temperature in radiosonde measurements and ERA-Interim are smaller than 2 K and the greatest differences are about 4 K. However, the differences between temperature in radiosonde measurements and in the USSA-1976 are smaller than 10 K only in three stations, and the greatest differences are about 20 K. The profiles of relative differences of the molecular backscatter coefficients between radiosonde measurements and ERA-Interim show negative values at most of the pressure levels. Positive values greater than 0.5 % are only observed above 8 km. However, the profiles of relative differences of the molecular backscatter coefficients between radiosonde measurements and the USSA-1976 shows positive values from surface to 12 km in all stations. Only in four stations are observed negative differences above 14 km with values smaller than -5 %. Comodoro Rivadavia presents the best behavior because is located near to 45 degree of latitude, but even in this station the differences between radiosonde measurements and the USSA-1976 are bigger than those observed between radiosonde measurements and ERA-Interim. [ABSTRACT FROM AUTHOR]

ABSTRACT The ocean occupies 95% of the Caribbean's area and plays a leading role in the region's climate, thus making the sea surface temperature ( SST) a very important regional climate index. This, in conjunction with the lack of a regionally consistent, quality-controlled surface temperature dataset increases the scientific value of using SST to characterize the regional climate and its trends. This study determines the magnitudes of the long-term SST trends in the Wider Caribbean ( WC) and the Antilles. We overcome the presence of discontinuity points in the SST time series using the change point statistical technique. Annual mean SST trends combining the subperiods 1906-1969 and 1972-2005 are 1.32 ± 0.41 °C per century for the Antilles and 1.08 ± 0.32 °C per century for the WC. For the same regions during the subperiod 1972-2005, the corresponding trends are 1.41 ± 0.68 °C per century and 1.18 ± 0.49 °C per century, illustrating the warming intensification during the last four decades. A significant correlation is found between the SSTs in the Caribbean and Atlantic Multidecadal Oscillation ( AMO) index, suggesting a potential link between Caribbean SSTs and the mechanisms governing the Atlantic basin-wide SSTs. Finally, the capabilities of two state-of-the-art coupled climate models, the Norwegian Earth System Model ( NorESM1-M) and the Bergen Climate Model ( BCM), to simulate SST in the Caribbean were tested. Both models produce an adequate simulation of the annual mean SST anomalies and SST seasonal cycle for the WC and the Antilles. The simulated annual and monthly mean SSTs are colder in the two models compared to the observations, a common feature among the majority of general circulation models participating in the Coupled Model Intercomparison Project Phase 5. However, despite these minor deficiencies both BCM and NorESM1-M are considered adequate for conducting SST simulations relevant for future climate change research in the Caribbean. [ABSTRACT FROM AUTHOR]

On 22 April 2015, the Calbuco volcano in Chile (Lat: 41.33 ∘ S, Long: 72.62 ∘ W) erupted after 43 years of inactivity followed by a great amount of aerosol injection into the atmosphere. The pyroclastic material dispersed into the atmosphere posed a potential threat to aviation traffic and air quality over affected a large area. The plumes and debris spread from its location to Patagonian and Pampean regions, reaching the Atlantic and Pacific Oceans and neighboring countries, such as Argentina, Brazil and Uruguay, driven by the westerly winds at these latitudes. The presence of volcanic aerosol layers could be identified promptly at the proximities of Calbuco and afterwards by remote sensing using satellites and lidars in the path of the dispersed aerosols. The Cloud-Aerosol Lidar and Pathfinder Satellite Observations (CALIPSO), Moderate Resolution Imaging Spectroradiometer (MODIS) on board of AQUA/TERRA satellites and Ozone Mapping and Profiler Suite (OMPS) on board of Suomi National Polar-orbiting Partnership (Suomi NPP) satellite were the space platforms used to track the injected layers and a multi-channel lidar system from Latin America Lidar Network (LALINET) SPU Lidar station in South America allowed us to get the spatial and temporal distribution of Calbuco ashes after its occurrence. The SPU lidar stations co-located Aerosol Robotic Network (AERONET) sunphotometers to help in the optical characterization. Here, we present the volcanic layer transported over São Paulo area and the detection of aerosol plume between 18 and 20 km. The path traveled by the volcanic aerosol to reach the Metropolitan Area of São Paulo (MASP) was tracked by CALIPSO and the aerosol optical and geometrical properties were retrieved at some points to monitor the plume evolution. Total attenuated backscatter profile at 532 nm obtained by CALIPSO revealed the height range extension of the aerosol plume between 18 and 20 km and are in agreement with SPU lidar range corrected signal at 532 nm. The daily evolution of Aerosol Optical Depth (AOD) at 532 and 355 nm, retrieved from AERONET sunphotometer, showed a substantial increasing on 27 April, the day of the volcanic plume detection at Metropolitan Area of São Paulo (MASP), achieving values of 0. 33 ± 0. 16 and 0. 22 ± 0. 09 at 355 and 532 nm, respectively. AERONET aerosol size distribution was dominated by fine mode aerosol over coarse mode, especially on 27 and 28 April. The space and time coincident aerosol extinction profiles from SPU lidar station and OMPS LP from the Calbuco eruption conducted on 27 April agreed on the double layer structure. The main objective of this study was the application of the transmittance method, using the Platt formalism, to calculate the optical and physical properties of volcanic plume, i.e., aerosol bottom and top altitude, the aerosol optical depth and lidar ratio. The aerosol plume was detected between 18 and 19.3 km, with AOD value of 0.159 at 532 nm and Ånsgtröm exponent of 0. 61 ± 0. 58 . The lidar ratio retrieved was 76 ± 27 sr and 63 ± 21 sr at 532 and 355 nm, respectively. Considering the values of these parameters, the Calbuco volcanic aerosol layers could be classified as sulfates with some ash type. [ABSTRACT FROM AUTHOR]