Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cloud top phase”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

Diurnal variation of marine stratocumulus over San Nicolas Island during July 1987

A study is presented of surface measurements made at San Nicolas Island during the intensive field observation marine stratocumulus phase of the First International Satellite Cloud Climatology Program Regional Experiment, July 1987, to retrieve the average diurnal variation of marine stratocumulus and associated surface variables. The diurnal behavior of the cloud base appears to be consistent with model-predicted uncoupling of the cloud layer and the subcloud layer as the turbulent flux of moisture is limited by solar heating close to the cloud base. Also, variation in surface air temperature is consistent with the limitation of the turbulent flux of heat between the two layers, screening the surface from the effect of longwave cooling from the cloud top.

Blaskovic, Miriam↗

A FIRE-ACE/SHEBA Case Study of Mixed-Phase Arctic Boundary Layer Clouds: Entrainment Rate Limitations on Rapid Primary Ice Nucleation Processes

Observations of long-lived mixed-phase Arctic boundary layer clouds on 7 May 1998 during the First International Satellite Cloud Climatology Project (ISCCP) Regional Experiment (FIRE)Arctic Cloud Experiment (ACE)Surface Heat Budget of the Arctic Ocean (SHEBA) campaign provide a unique opportunity to test understanding of cloud ice formation. Under the microphysically simple conditions observed (apparently negligible ice aggregation, sublimation, and multiplication), the only expected source of new ice crystals is activation of heterogeneous ice nuclei (IN) and the only sink is sedimentation. Large-eddy simulations with size-resolved microphysics are initialized with IN number concentration N(sub IN) measured above cloud top, but details of IN activation behavior are unknown. If activated rapidly (in deposition, condensation, or immersion modes), as commonly assumed, IN are depleted from the well-mixed boundary layer within minutes. Quasi-equilibrium ice number concentration N(sub i) is then limited to a small fraction of overlying N(sub IN) that is determined by the cloud-top entrainment rate w(sub e) divided by the number-weighted ice fall speed at the surface v(sub f). Because w(sub c)< 1 cm/s and v(sub f)> 10 cm/s, N(sub i)/N(sub IN)<< 1. Such conditions may be common for this cloud type, which has implications for modeling IN diagnostically, interpreting measurements, and quantifying sensitivity to increasing N(sub IN) (when w(sub e)/v(sub f)< 1, entrainment rate limitations serve to buffer cloud system response). To reproduce observed ice crystal size distributions and cloud radar reflectivities with rapidly consumed IN in this case, the measured above-cloud N(sub IN) must be multiplied by approximately 30. However, results are sensitive to assumed ice crystal properties not constrained by measurements. In addition, simulations do not reproduce the pronounced mesoscale heterogeneity in radar reflectivity that is observed.

Fridlin, Ann↗

Wave Forcing of Saturn's Equatorial Oscillation

Ground-based measurements and Cassini data from CIRS thermal-infrared spectra and radio-occultation soundings have characterized the spatial structure and temporal behavior of a 15-year equatorial oscillation in Saturn's stratosphere. The equatorial region displays a vertical pattern of alternating warm and cold anomalies and, concomitantly, easterly and westerly winds relative to the cloud-top winds, with a peak-to-peak amplitude of 200 m/s. Comparison of the Cassini data over a four-year period has established that the pattern of mean zonal winds and temperatures descends at a rate of roughly I scale height over 4 years. This behavior is reminiscent of the equatorial oscillations in Earth's middle atmosphere. Here the zonal-mean spatial structure and descending pattern are driven by the absorption of vertically propagating waves. The maximum excursions in the pattern of easterly and westerly winds is determined by the limits of the zonal phase velocities of the waves. Here we report on the characterization of the waves seen in the temperature profiles retrieved from the Cassini radio-occultation soundings. The equatorial profiles exhibit a complex pattern of wavelike structure with dimensions one pressure scale height and smaller. We combine a spectral decomposition with a WKBJ analysis, where the vertical wavelength is assumed to vary slowly with the ambient static stability and doppler-shifted phase velocity of the wave. Use of the temperature and zonal wind maps from CIRS makes this approach viable. On Earth, the wave forcing associated with the equatorial oscillations generates secondary meridional circulations that affect the mean flow and planetary wave ducting well away from the equator. This may relate to the triggering of the recently reported mid-latitude storms on Saturn.

Flasar, F. M.↗

Remote Sensing Data from CLARET: A Prototype Cart Data Set

A data set containing radiation, meteorological, and cloud sensor observations is documented. It was prepared for use by the Department of Energy's Atmospheric Radiation Measurement (ARM) program and other interested scientists. These data are a precursor of the types of data that ARM Cloud And Radiation Testbed (CART) sites will provide. The data are from the Cloud Lidar And Radar Exploratory Test (CLARET) conducted by the Wave Propagation Laboratory during autumn 1989 in the Denver-Boulder area of Colorado primarily for the purpose of developing new cloud-sensing techniques on cirrus. After becoming aware of this experiment, ARM scientists requested archival of subsets or the data to assist in the developing ARM program. Five CLARET cases were selected: two with cirrus, one with stratus, one with mixed-phase clouds, and one with clear skies. The cases range from 2 to 9.5 h in length. A pyranometer, pyrgeometer, pyrheliometer, and an infrared radiometer constituted the ensemble of instruments that provided surface radiation data. A lidar, radar, and ceilometer observed the cloud geometrical structure, and visual reports and all-sky camera observations were assimilated to provide cloud cover data. Radiosondes, wind profiler, RASS (profiling virtual temperature), microwave radiometers (observing column integrated liquid water and water vapor), and standard surface measurements provided meteorological data. Satellite data from the stratus case and one cirrus case were analyzed for statistics on cloud cover and top height. The main body of the selected data are available on diskette from the Wave Propagation Laboratory or Los Alamos National Laboratory. In addition to documenting the data set, this report describes CLARET and gives a bibliography of publications associated with the project. Some preliminary results of CLARET' research are also summarized. Simultaneous CO 2 lidar and radar backscatter measurements were shown to provide estimates of the effective radius of ice particles. Simultaneous radar and infrared radiometer data appear useful for estimating column-integrated numbers and average sizes of ice cloud particles. Ice water content obtained with this method compared favorably with values from another empirical technique using radar data alone. Depolarization of the CO 2 lidar signal from ice clouds was surprisingly small, suggesting that calculation of backscatter from nonspherical particles for this lidar is a tractable problem. Examples are also cited of CO 2 lidar measurements of the effective radius of water cloud drop size distributions and of inference of the size of pristine ice crystals that assume a particular orientation in the air. These parameters are all important to radiative transfer through clouds.

Clouds (Meteorology)↗

Mariner 9 television limb observations of dust and ice hazes on Mars

Results are presented for a statistical analysis of over 3000 Mariner 9 television reflectance profiles of the planetary limb of Mars for the period from the decaying phase of a global dust storm (southern summer) to southern winter solstice. It is shown that during the decaying phase of the storm, the dust is maintained aloft to high elevations by vertical mixing primarily in the equatorial zone, apparently by the direct or indirect action of the diurnal tide; the result is that the dust decay is much faster in high-latitude regions than near the equator. Significant dust remains aloft in Martian equatorial regions long after the major dust storms have subsided. At the top of the dusty atmosphere, a region exists where temperatures relax toward values expected for a dust-free CO2 atmosphere; this is therefore a region of convective overturning topped by thin condensation hazes resembling terrestrial noctilucent clouds, thereby showing occasionally gravity wave structures.

Anderson, E.↗

Inhomogeneous models of the Venus clouds containing sulfur

Based on the suggestion that elemental sulfur is responsible for the yellow color of Venus, calculations are compared at 3.4 microns of the reflectivity phase function of two sulfur containing inhomogeneous cloud models with that of a homogeneous model. Assuming reflectivity observations with 25% or less total error, comparison of the model calculations leads to a minimum detectable mass of sulfur equal to 7% of the mass of sulfuric acid for the inhomogeneous drop model. For the inhomogeneous cloud model the comparison leads to a minimum detectable mass of sulfur between 17% and 38% of the mass of the acid drops, depending upon the actual size of the large particles. It is concluded that moderately accurate 3.4 microns reflectivity observations are capable of detecting quite small amounts of elemental sulfur at the top of the Venus clouds.

Smith, S. M.↗

International Satellite Cloud Climatology Project

The International Satellite Cloud Climatology Project (ISCCP) was established as the first project of the World Climate Research Programme (WCP-2) to collect and analyze satellite radiance measurements to infer the global distribution of cloud radiative properties and their diurnal and seasonal variations. The operational phase of ISCCP began in July 1983 and is currently planned to continue through June 2010. [Mission Objectives] To produce a global, reduced resolution, infrared and visible, calibrated and normalized radiance data set containing basic information on the radiative properties of the atmosphere from which cloud parameters can be derived. To stimulate and coordinate basic research on techniques for inferring the physical properties of clouds from the condensed radiance data set and to apply the resulting algorithms to derive and validate a global cloud climatology for improving the parameterization of clouds in climate models. To promote research using ISCCP data and contributing to improved understanding of the Earth's radiation budget (top of the atmosphere and surface) and hydrological cycle. [Temporal_Coverage: Start_Date=1983-07-01; Stop_Date=] [Spatial_Coverage: Southernmost_Latitude=-90; Northernmost_Latitude=90; Westernmost_Longitude=-180; Easternmost_Longitude=180].

Rossow, William B.↗

Dispersion and Aging of Volcanic Aerosols after the La Soufriere Eruption in April 2021

Volcanic aerosols change the atmospheric composition and thereby affect weather and climate. Aerosol dynamic processes such as nucleation, condensation, and coagulation modify the shape, size, and mass of aerosol particles, which influence their atmospheric lifetime and radiative properties. Nevertheless, most models omit these processes for ash particles. In this work, we explore the ash aerosol aging and sulfate production during the first 4 days following the 2021 La Soufrière (St. Vincent) eruption with the ICON-ART model (ICOsahedral Nonhydrostatic model with Aerosol and Reactive Trace gases). Online coupling of ICON-ART with a one-dimensional volcanic plume model calculates volcanic emission, which makes it possible to resolve the different eruption phases of the noncontinuous La Soufrière eruption. We compared our simulated aerosol distribution and composition with observations from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) instrument, the Multiangle Imaging SpectroRadiometer (MISR) Research Aerosol (RA) Algorithm, and the Barbados Cloud Observatory (BCO). We show that online coupling is essential to adequately model the emissions and plume development close to the volcano. The modeled aerosol aging is in very good agreement with observations from MISR near the emission source and with CALIOP at larger distances. Furthermore, particle aging occurs faster in the troposphere than in the stratosphere due to the availability of water vapor and OH, but a layer of coated ash appears at the plume top due to faster oxidation of SO2 and lofting by aerosol-radiation interaction. This paper gives the first direct comparison of aerosol aging in volcanic eruption plumes between simulations and observations

volcanic aerosols↗

Effect of Scattering Angle on DSCOVR/EPIC Observations

The Earth Polychromatic Imaging Camera (EPIC) on the Deep Space Climate Observatory (DSCOVR) routinely captures reflected radiation from the whole sunlit side of the Earth in the near backward direction to monitor the changing planet. The instrument had routinely operated until 27 June 2019, when the spacecraft was placed in an extended safe hold due to degradation of an inertial navigation unit. DSCOVR returned to full operations on 2 March 2020. Since then, the range of scattering angles between the incident sunlight and sensor direction has been larger than before and the largest scattering angle reaches ∼178°, only 2° from perfect backscattering, proving a unique opportunity to study the top-of-atmosphere (TOA) reflectance under such extreme conditions. In the paper, we compare EPIC global spectral reflectances in 2021–2016. We found that there are four occasions when the scattering angle reaches about 178° and associated with them enhanced global daily average spectral reflectances in 2021. The scattering angle related reflectance enhancements are not found in 2016 data when the maximum scattering angle is about 174.5°. CERES data do not show such occasions in global daily reflected shortwave flux. As a result, those enhanced reflectance occasions are primarily due to the change in scattering angle. The enhancement due to changes in scattering angle depends strongly on wavelength, primarily because of wavelength dependence of cloud scattering phase function. Radiative transfer calculations show that the change in scattering angles has the largest impact on reflectance in the red and NIR channels at 680 nm and 780 nm and the smallest influence on reflectance in the UV channel at 388 nm, consistent with EPIC observations. The change of global average cloud amount also plays an important role in the reflectance enhancement. The influence of the cloud effect depends on whether the change is in phase or not with the change of scattering angle.

DSCOVR EPIC↗

Forward and Inverse Models for Satellite Remote Sensors using Principal Component Analysis

Satellite remote sensors such as AIRS on Aqua, CrIS on S-NPP, NOAA20 and JPSS-2, IASI on Metop A, B, and C make millions of observations each day with thousands of spectral channels for each observation; this poses challenges for efficiently inversion of the inherently large dataset as needed to retrieve atmospheric and surface properties. This presentation will illustrate the use of Principal Component Analysis (PCA) to speed up radiative transfer forward model calculations and to stabilize the inversion algorithms. A Principal Component-based radiative transfer model (PCRTM) developed at NASA Langley Research Center can simulate top of atmosphere (TOA) radiance or reflectance spectra from 50 cm-1 to 50000 cm-1 (200 m to 0.20 m quickly and accurately. PCRTM demonstrated very high accuracy relative to reference line-by-line radiative transfer models and it saves orders of magnitude computational time. Examples of the PCRTM model developed for hyperspectral sensors such as AIRS, CrIS, IASI, NAST-I, SHIS, CPF, TEMPO, SBG, OMI, and SCIAMACHY will be presented. In addition to using the PCRTM as forward model, the NASA Langley developed inversion algorithm also uses PCA to compress the state vector into a compressed dimension to speed up and stabilize the inversion process. Examples of retrieved atmospheric temperature, water vapor, CO2, CO, CH4, N2O, and O3 profiles, cloud properties (optical depth, size, phase, and height), and surface properties (surface emissivity spectra and skin temperatures) will be presented. This algorithm is being transitioned to the NASA Sounder SIPS and NASA's Goddard Earth Sciences Data and Information Services Center (GES DISC).

forward model↗

The End of Comet Shoemaker-Levy 9

The pieces of the former Comet Shoemaker-Levy 9 struck Jupiter during the week of July 16-22 1994. We divide the general impact phenomena into four phases: (1) Atmospheric entry should have been accompanied by an entry flash, probably only detectable by the Galileo spacecraft. (2) Aerodynamic forces cause the impactor to explode. We show that the apparent altitude of the explosions is consistent with the analytical model proposed by Zahnle and Chyba et al, with numerical models calculated by Mac Low and Zahnle and Yabe et al., and with dynamically-based estimates of fragment masses Suggested by Scotti and Melosh and Asphaug and Benz. (3) The explosion produces a fireball that propagates most quickly up the wake, as predicted by Takata et al. and Crawford et al. The fireball's luminosity fades over a timescale of 30-60 seconds. (4) Each explosion produced a plume that reached thousands of kilometers above the cloud-tops. The plumes appear to have been dusty and SO illuminated when they reached sunlight. The plumes rose and fell over a period of ten or twenty minutes, governed by the simple laws of ballistics. When plume material (dust) fell back on the atmosphere it was reheated and glowed brightly in the thermal infrared for a period of order ten to twenty minutes; a large fraction of the impact energy was spent thus. Geometry introduces additional factors that complicate interpretation of the observed light curves.

Zahnle, K. J.↗

Cloud-Radiative Driving of the Madden-Julian Oscillation as Seen by the A-Train

Cloud and water vapor radiative heating anomalies associated with convection may be an effective source of moist static energy driving the Madden-Julian Oscillation (MJO). In this paper five years of radiative heating profiles derived from CloudSat radar and Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation data are analyzed to document radiative heating anomalies during the MJO. Atmospheric shortwave absorption and surface longwave radiation anomalies are of opposite sign and 10-20% as large as top-of-atmosphere outgoing longwave radiation (OLR) anomalies, confirming that OLR provides a useful estimate of the total column radiative heating anomaly. Positive anomalies generally peak about one week before the MJO peak and are smallest over the Indian Ocean. Anomalies over the Maritime Continent are strongest, and coincident with the MJO peak. Shortwave heating profile anomalies are about half as large as longwave anomalies in the active region of the MJO but generally of opposite sign; thus shortwave heating damps the longwave destabilization of the lower troposphere. The exception is the onset phase of the MJO, where shortwave and longwave heating anomalies due to thin cirrus are both positive in the upper troposphere and exert a stabilizing influence. Specific humidity anomalies in the middle troposphere reach 0.5 g kg(exp. -1), but the associated clear sky heating anomaly is very small. Radiative enhancement of column moist static energy becomes significant about 10 days before the MJO peak, when precipitation anomalies are still increasing, and then remains high after the MJO peak after precipitation has begun to decline.

symbols↗

The CHROMA Cloud-Top Pressure Retrieval Algorithm for the Plankton, Aerosol, Cloud, Ocean Ecosystem (PACE) Satellite Mission

This paper provides the theoretical basis and simulated retrievals for the Cloud Height Retrieval from O2 Molecular Absorption (CHROMA) algorithm. Simulations are performed for the Ocean Color Instrument (OCI), which is the primary payload on the forthcoming NASA Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission, and the Ocean Land Colour Instrument (OLCI) currently flying on the Sentinel 3 satellites. CHROMA is a Bayesian approach which simultaneously retrieves cloud optical thickness (COT), cloud-top pressure and height (CTP and CTH respectively), and (with a significant prior constraint) surface albedo. Simulated retrievals suggest that the sensor and algorithm should be able to meet the PACE mission goal for CTP error, which is ±60 mb for 65 % of opaque (COT ≥3) single-layer clouds on global average. CHROMA will provide pixel-level uncertainty estimates, which are demonstrated to have skill at telling low-error situations from high-error ones. CTP uncertainty estimates are well-calibrated in magnitude, although COT uncertainty is overestimated relative to observed errors. OLCI performance is found to be slightly better than OCI overall, demonstrating that it is a suitable proxy for the latter in advance of PACE's launch. CTP error is only weakly sensitive to correct cloud phase identification or assumed ice crystal habit/roughness. As with other similar algorithms, for simulated retrievals of multi-layer systems consisting of optically thin cirrus clouds above liquid clouds, retrieved height tends to be underestimated because the satellite signal is dominated by the optically thicker lower layer. Total (liquid plus ice) COT also becomes underestimated in these situations. However, retrieved CTP becomes closer to that of the upper ice layer for ice COT ≈3 or higher.

PACE↗

Aerosol Direct, Indirect, Semidirect, and Surface Albedo Effects from Sector Contributions Based on the IPCC AR5 Emissions for Preindustrial and Present-day Conditions

The anthropogenic increase in aerosol concentrations since preindustrial times and its net cooling effect on the atmosphere is thought to mask some of the greenhouse gas-induced warming. Although the overall effect of aerosols on solar radiation and clouds is most certainly negative, some individual forcing agents and feedbacks have positive forcing effects. Recent studies have tried to identify some of those positive forcing agents and their individual emission sectors, with the hope that mitigation policies could be developed to target those emitters. Understanding the net effect of multisource emitting sectors and the involved cloud feedbacks is very challenging, and this paper will clarify forcing and feedback effects by separating direct, indirect, semidirect and surface albedo effects due to aerosols. To this end, we apply the Goddard Institute for Space Studies climate model including detailed aerosol microphysics to examine aerosol impacts on climate by isolating single emission sector contributions as given by the Coupled Model Intercomparison Project Phase 5 (CMIP5) emission data sets developed for Intergovernmental Panel on Climate Change (IPCC) AR5. For the modeled past 150 years, using the climate model and emissions from preindustrial times to present-day, the total global annual mean aerosol radiative forcing is -0.6 W/m(exp 2), with the largest contribution from the direct effect (-0.5 W/m(exp 2)). Aerosol-induced changes on cloud cover often depends on cloud type and geographical region. The indirect (includes only the cloud albedo effect with -0.17 W/m(exp 2)) and semidirect effects (-0.10 W/m(exp 2)) can be isolated on a regional scale, and they often have opposing forcing effects, leading to overall small forcing effects on a global scale. Although the surface albedo effects from aerosols are small (0.016 W/m(exp 2)), triggered feedbacks on top of the atmosphere (TOA) radiative forcing can be 10 times larger. Our results point out that each emission sector has varying impacts by geographical region. For example, the single sector most responsible for a net positive radiative forcing is the transportation sector in the United States, agricultural burning and transportation in Europe, and the domestic emission sector in Asia. These sectors are attractive mitigation targets.

Bauer, Susanne E.↗

Seasonal and Interannual Variabilities in Tropical Tropospheric Ozone

This paper presents a detailed characterization of seasonal and interannual variability in tropical tropospheric column ozone (TCO). TCO time series are derived from 20 years (1979-1998) of total ozone mapping spectrometer (TOMS) data using the convective cloud differential (CCD) method. Our study identifies three regions in the tropics with distinctly different zonal characteristics related to seasonal and interannual variability. These three regions are the eastern Pacific, Atlantic, and western Pacific. Results show that in both the eastern and western Pacific seasonal-cycle variability of northern hemisphere (NH) TCO exhibits maximum amount during NH spring whereas largest amount in southern hemisphere (SH) TCO occurs during SH spring. In the Atlantic, maximum TCO in both hemispheres occurs in SH spring. These seasonal cycles are shown to be comparable to seasonal cycles present in ground-based ozonesonde measurements. Interannual variability in the Atlantic region indicates a quasi-biennial oscillation (QBO) signal that is out of phase with the QBO present in stratospheric column ozone (SCO). This is consistent with high pollution and high concentrations of mid-to-upper tropospheric O3-producing precursors in this region. The out of phase relation suggests a UV modulation of tropospheric photochemistry caused by the QBO in stratospheric O3. During El Nino events there is anomalously low TCO in the eastern Pacific and high values in the western Pacific, indicating the effects of convectively-driven transport of low-value boundary layer O3 (reducing TCO) and O3 precursors including H2O and OH. A simplified technique is proposed to derive high-resolution maps of TCO in the tropics even in the absence of tropopause-level clouds. This promising approach requires only total ozone gridded measurements and utilizes the small variability observed in TCO near the dateline. This technique has an advantage compared to the CCD method because the latter requires high-resolution footprint measurements of both reflectivity and total ozone in the presence of tropopause-level cloud tops.

Ziemke, J. R.↗

Ice-Over-Water Cloud Properties in an Artificial Neural Network Approach

Clouds are a crucial component of the atmospheric energy system, particularly the radiative balance within, above, and below the troposphere. The vertical distribution of cloud mass and phase determines layer heating rates, the loss of radiation to space, and the amount of radiative heating at the surface. Thus, it is important to know how clouds are distributed both vertically and horizontally at all times of day. Satellite remote sensing is the only approach available to monitor clouds day and night around the globe. In this paper several artificial neural network (ANN) algorithms, employing several Aqua MODIS infrared channels, profiles of relative humidity and temperature from GMAO numerical weather analyses, and the retrieved total cloud visible optical depth, are trained to detect multilayer ice-over-water cloud systems and to retrieve some of their properties as identified by a year of 2008 Aqua MODIS data matched with CloudSat and CALIPSO (CC) cloud profiles. The CC lidar and radar profiles provide the vertical structure that serves as output truth for the multilayer algorithm. The neural networks were trained using one year (2008) of cloud top height data from the CC dataset, with correlation around 0.94 (0.95) and MAE as low as 0.82 (0.81) km for nonpolar regions during the day (night). Applying the trained ANN to independent year 2009 MODIS data resulted in a combined ML and single layer hit rate of 86.4% (85.1%) for nonpolar regions during the day (night). Since the ANN is trained using near-nadir MODIS pixels, infrared radiance corrections were developed as a function of view zenith angle from MODIS and applied to off-nadir pixels when processing MODIS swath data. The multilayer amount derived with the ANN is relatively invariant with increasing view zenith angle compared to the multilayer amount without the corrections.

MODIS↗

Modeling of 2008 Kasatochi Volcanic Sulfate Direct Radiative Forcing: Assimilation of OMI SO2 Plume Height Data and Comparison with MODIS and CALIOP Observations

Volcanic SO2 column amount and injection height retrieved from the Ozone Monitoring Instrument (OMI) with the Extended Iterative Spectral Fitting (EISF) technique are used to initialize a global chemistry transport model (GEOS-Chem) to simulate the atmospheric transport and lifecycle of volcanic SO2 and sulfate aerosol from the 2008 Kasatochi eruption, and to subsequently estimate the direct shortwave, top-of-the-atmosphere radiative forcing of the volcanic sulfate aerosol. Analysis shows that the integrated use of OMI SO2 plume height in GEOS-Chem yields: (a) good agreement of the temporal evolution of 3-D volcanic sulfate distributions between model simulations and satellite observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) and Cloud-Aerosol Lidar with Orthogonal Polarisation (CALIOP), and (b) an e-folding time for volcanic SO2 that is consistent with OMI measurements, reflecting SO2 oxidation in the upper troposphere and stratosphere is reliably represented in the model. However, a consistent (approx. 25 %) low bias is found in the GEOS-Chem simulated SO2 burden, and is likely due to a high (approx.20 %) bias of cloud liquid water amount (as compared to the MODIS cloud product) and the resultant stronger SO2 oxidation in the GEOS meteorological data during the first week after eruption when part of SO2 underwent aqueous-phase oxidation in clouds. Radiative transfer calculations show that the forcing by Kasatochi volcanic sulfate aerosol becomes negligible 6 months after the eruption, but its global average over the first month is -1.3W/sq m, with the majority of the forcing-influenced region located north of 20degN, and with daily peak values up to -2W/sq m on days 16-17. Sensitivity experiments show that every 2 km decrease of SO2 injection height in the GEOS-Chem simulations will result in a approx.25% decrease in volcanic sulfate forcing; similar sensitivity but opposite sign also holds for a 0.03 m increase of geometric radius of the volcanic aerosol particles. Both sensitivities highlight the need to characterize the SO2 plume height and aerosol particle size from space. While more research efforts are warranted, this study is among the first to assimilate both satellite-based SO2 plume height and amount into a chemical transport model for an improved simulation of volcanic SO2 and sulfate transport.

Wang, J.↗

Interpretation of satellite-measured bidirectional reflectance from Cirrus cloudy atmospheres

The interpretation of the observed bidirectional reflectance from cirrus cloudy atmospheres is presented. A theoretical model was developed for the computation of the transfer of solar radiation in an anisotropic medium with particular applications to oriented ice crystals in cirrus clouds. In this model, the adding principle for radiative transfer was used with modifications to account for the anisotropy of scattering particles and the associated scattering phase matrix. The single-scattering properties, including the phase function, single-scattering albedo, and extinction cross sections, were used for randomly and horizontally oriented hexagonal ice crystals in radiative transfer computations. The radiative transfer model developed for the cirrus clouds was modified to account for the scattering contributions from the atmosphere and the surface. In order to test the relevance and significance of the ice crystal model for the interpretation of observed bidirectional reflectance from satellites, visible radiances collected on the half hour by the GOES series were selected. The data are calibrated and corrected with the proper filter functions, then navigated to match selected landmark data. A number of clear and cloudy cases during the cirrus IFO of the Fire experiment were chosen for theoretical analyses. The cloud particle shape and size distributions that were taken during satellite overpasses are used in radiative transfer calculations. The sensitivities of the shape, orientation, and size distribution of ice crystals on the reflected intensities at the top of the atmosphere are investigated. Finally, the relative importance of these cloud microphysical properties in the interpretation of satellite bidirectional reflectance are assessed and presented.

Takano, Y.↗