Diagnosis and simulation of a rapidly developing cyclone related to a severe dust storm in East Asia

Mongolian cyclones are important to the outbreaks of severe dust storms in Northeast Asia. In this paper, we conduct a diagnostic study and a numerical simulation of a rapidly developing Mongolian cyclone that produced the April 5 to 7, 2000 severe dust storm in East Asia. The surface pressure of the cyclone decreased rapidly on April 5 to 6 causing strong surface winds. The diagnosis shows that the cyclone was a “dry cyclone”, as the convergence of moisture flux was weak and the baroclinic forcing was strong. The analysis of the Q vectors also reveals that ageostrophic wind was significant in the middle and upper part of the cyclone, resulting in strong vertical motion in the lower part. The three-dimensional structure of the cyclone is characterized by an ascending southern warm current and a descending northern cold current. This structure is favorable to the release of available potential energy and the intensification of the cyclone. In the second part of this study, we develop a coupled system of a dust emission scheme and a mesoscale numerical model (MM5V3) and applied the system to the simulation of the Mongolian cyclone. The main features of the cyclone identified from the diagnostic analysis are successfully reproduced. The predicted dust storm-affected area is found to be consistent with the meteorological observations and satellite remote sensing. It is shown that the coupled system is capable of predicting the spatial and temporal variations of the dust storm.     

The characteristics of Asian dust events in Northeast Asia during the springtime from 1993 to 2004

The characteristics of Asian dust events that occurred in Northeast Asia during the springtime from 1993 to 2004 are investigated using 3-hourly SYNOP reports (World Meteorological Organization). Occurrences of blowing sand and dust storm are low in 1997 and 1999, but have increased rapidly since 2000. The maximum occurrence was recorded in 2002. Wind velocity of 6.5 m s^− 1 as a threshold wind velocity is not so exactly consistent with the occurrence of blowing sand. However, wind velocity of 14 m s^− 1 as a strong wind causing dust storm had similar tendency to those of dust storm and Dust Storm Index.Source regions of Asian dust are divided into three regions (A: dry arid, B: semi-arid, and C: cultivated), based upon the occurrence of blowing sand and dust storm. Eight meteorological stations are selected in three regions, which have frequent occurrences of blowing sand. Source regions of Asian dust that affect the Korean peninsula are gradually extending eastward. Positive anomalies of NDVI occurred in 1994, 1995, and 1998 when temperature was high and precipitation was heavy. However, the frequent occurrence of the dust phenomena is not always consistent with lots of vegetation, high temperature, and much precipitation in this study.     

An investigation of sand–dust storm events and land surface characteristics in China using NOAA NDVI data

Observations from 560 weather stations in China show that sand–dust storms occur most frequently in April in north China. The region consists of Sub-dry Mid Temperate, Dry Mid Temperate, Sub-dry South Temperate and Dry South Temperate Zones and much of the land surface is desert or semi-desert: it is relatively dry with minimal rainfall and a high annual mean temperature. In most regions of China, the annual mean frequency of sand–dust events decreased sharply between 1980 and 1997 and then increased from 1997 to 2000. Statistical analyses demonstrate that the frequency of sand–dust storms correlates highly with wind speed, which in turn is strongly related to land surface features; on the other hand, a significant correlation between storm events and other atmospheric quantities such as precipitation and temperature was not observed. Accordingly, land surface cover characteristics (vegetation, snowfall and soil texture) may play a significant role in determining the occurrence of sand–dust storms in China. Analysis of Normalized Difference Vegetation Index derived from National Oceanic and Atmospheric Administration and Empirical Orthogonal Function show that since 1995 surface vegetation cover in large areas of Northern China has significantly deteriorated. Moreover, a high correlation is shown to exist among the annual occurrence of sand–dust storms, surface vegetation cover and snowfall. This suggests that the deterioration of surface vegetation cover may strongly influence the occurrence of sand–dust storms in China. Soils with coarse and medium textures are found to be more associated with sand–dust storms than other soils.     

Numerical simulation and synoptic analysis of dust emission and transport in East Asia

This study is concerned with the mechanisms of dust storm development in East Asia and the characteristics of the responsible synoptic systems. Two severe East Asian dust storms which occurred in spring 2002 are analyzed using synoptic and remote sensing data. The relationships between the formation and the movement of the dust storms and the evolution of the synoptic systems are examined. It is shown that a dust storm may develop when a synoptic system moves to the desert area of Northeast Asia with a surface wind speed exceeding 6 m s^− 1. Numerical simulations of the two dust storms are carried out using a dust storm forecasting model. The performance of the model is verified with observations. The dust sources are found to be consistent with the desert regions in Northeast Asia, but cover a somewhat larger area than the observations suggest. Finally, we present a conceptual model of dust storm generation and movement in East Asia on the basis of numerical modeling and synoptic analysis.     

An integrated dust storm prediction system suitable for east Asia and its simulation results

An integrated dust storm modeling system is developed for the prediction of dust storms. The system couples a wind erosion scheme, a dust transportation model and the Penn State/UCAR modeling system (MM5) with a geographic information database. The system can be used for the prediction of dust emission rate and dust concentration associated with individual dust storm events. Two severe dust storm events occurred in spring 2002, one on the 19th–22nd of March and the other on the 6th–9th of April. The integrated modeling system is used to simulate the two events. The numerical results are compared with surface weather records and satellite images and good agreement is found between the model results and observation in dust concentration distribution and evolutions. The Gobi Desert in southern Mongolia and the Badain Jaran Desert, Tengger Desert and Hunshandake sandy land in Inner Mongolia (China) are identified to be the dust sources for the two events. The dominant modes of dust particles over western Inner Mongolia and Mongolia are from 2 to 11 μm in size, and 2 to 22 μm over Beijing and its surrounding area. The emission of particles in the 2–11 μm size range is found to be most important for Northeast Asian dust storms.     

Albedo models for the residual south polar cap on Mars: Implications for the stability of the cap under near-perihelion global dust storm conditions

It is uncertain whether the residual (perennial) south polar cap on Mars is a transitory or a permanent feature in the current Martian climate. While there is no firm evidence for complete disappearance of the cap in the past, clearly observable changes have been documented. Observations suggest that the perennial cap lost more CO₂ material in the spring/summer season prior to the Mariner 9 mission than in those same seasons monitored by Viking and Mars Global Surveyor. In this paper we examine one process that may contribute to these changes—the radiative effects of a planet encircling dust storm that starts during late Martian southern spring on the stability of the perennial south polar cap. To approach this, we model the radiative transfer through a dusty planetary atmosphere bounded by a sublimating CO₂ surface.A critical parameter for this modeling is the surface albedo spectrum from the near-UV to the thermal-IR, which was determined from both space-craft and Earth-based observations covering multiple wavelength regimes. Such a multi-wavelength approach is highly desirable since one spectral band by itself cannot tightly constrain the three-parameter space for polar surface albedo models, namely photon “scattering length” in the CO₂ ice and the amounts of intermixed water and dust.Our results suggest that a planet-encircling dust storm with onset near solstice can affect the perennial cap's stability, leading to advanced sublimation in a “dusty” year. Since the total amount of solid CO₂ removed by a single storm may be less than the total CO₂ thickness, a series of dust storms would be required to remove the entire residual CO₂ ice layer from the south perennial cap.     

Identification and physical retrieval of dust storm using three MODIS thermal IR channels

In this paper, the dust event on 7 April 2001 in northern China is investigated with three MODIS thermal infrared (IR) bands. It is found that for the dust cloud, the observed 11 μm minus 12 μm brightness temperature difference (BTD) is always negative, while the BTD of 8.5 μm minus 11 μm varies from positive to negative depending on the dust concentration. Based on these distinguishing properties, we develop a dust mask algorithm to identify the dust storm occurrence and spatial extent. The algorithm can be used successfully in both the daytime and nighttime. Using the Mie spherical scattering theory, the thermal radiation transfer through the single dust layer is performed with the widely used forward model DISTORT. Our calculations show that the dust-like aerosols can well explain the observed BTD although both of the complex refractive index and particle size of aerosols will significantly influence the BTD. When the complex refractive index is fixed (dust-like aerosols in this paper), then the dust optical thickness and effective radii of dust particles can be retrieved from the brightness temperature (BT) of the 11 μm channel and the BTD of 11 μm minus 12 μm channels, respectively. The integral dust column density can also be derived from the retrieved dust optical thickness and effective radius.     

Aeolian dust experiment on climate impact: An overview of Japan–China joint project ADEC

The Aeolian Dust Experiment on Climate Impact (ADEC) was initiated in April 2000 as a joint five-year Japan–China project. The goal was to understand the impact of aeolian dust on climate via radiative forcing (RF). Field experiments and numerical simulations were conducted from the source regions in northwestern China to the downwind region in Japan in order to understand wind erosion processes temporal and spatial distribution of dust during their long-range transportation chemical, physical, and optical properties of dust and the direct effect of radiative forcing due to dust. For this, three intensive observation periods (IOP) were conducted from April 2002 to April 2004.The in situ and network observation results are summarized as follows: (1) In situ observations of the wind erosion process revealed that the vertical profile of moving sand has a clear size dependency with height and saltation flux and that threshold wind velocity is dependent on soil moisture. Results also demonstrated that saltation flux is strongly dependent on the parent soil size distribution of the desert surface. (2) Both lidar observations and model simulations revealed a multiple dust layer in East Asia. A numerical simulation of a chemical transport model, CFORS, illustrated the elevated dust layer from the Taklimakan Desert and the lower dust layer from the Gobi Desert. The global-scale dust model, MASINGAR, also simulated the dust layer in the middle to upper free troposphere in East Asia, which originated from North Africa and the Middle East during a dust storm in March 2003. Raman lidar observations at Tsukuba, Japan, found the ice cloud associated with the dust layer at an altitude of 6 to 9 km. Analysis from lidar and the radio-sonde observation suggested that the Asian dust acted as ice nuclei at the ice-saturated region. These results suggest the importance of dust's climate impact via the indirect effect of radiative forcing due to the activation of dust into ice nuclei. (3) Studies on the aerosol concentration indicated that size distributions of aerosols in downwind regions have bimodal peaks. One peak was in the submicron range and the other in the supermicron range. The main soluble components of the supermicron peak were Na⁺, Ca²⁺, NO₃⁻, and Cl⁻. In the downwind region in Japan, the dust, sea salt, and a mixture of the two were found to be dominant in coarse particles in the mixed boundary layer. (4) Observation of the optical properties of dust by sky-radiometer, particle shoot absorption photometer (PSAP), and Nephelometer indicated that unpolluted dust at source region has a weaker absorption than originally believed.A sensitivity experiment of direct RF by dust indicated that single scattering albedo is the most important of the optical properties of dust and that the sensitivity of instantaneous RF in the shortwave region at the top of the atmosphere to the refractive index strongly depends on surface albedo. A global scale dust model, MASINGAR, was used for evaluation of direct RF due to dust. The results indicated the global mean RF at the top and the bottom of the atmosphere were − 0.46 and − 2.13 W m^− 2 with cloud and were almost half of the RF with cloud-free condition.     

The impact of Eurasian dust storms and anthropogenic emissions on atmospheric nutrient deposition rates in forested Japanese catchments and adjacent regional seas

Bulk precipitation and stream water chemistry data from 1993 to 2005 are used to analyze the relationship between Eurasian dust storms and nutrient deposition rates in the Kutsuki experimental forest (near Lake Biwa). From 2000 to 2005, atmospheric deposition, total nitrogen (TN), total phosphorus (TP) and dissolved silica (DSi) deposition rates increased by 26%, 132%, and 38%, respectively in the Kutsuki experimental forest. These TN and TP increases are associated with three seasonal factors: the increasing frequency and intensity of Eurasian spring dust events (March/April); the annual typhoon period (late August/September); and autumn/early winter (October to December) monsoons. The annual typhoon and monsoon winter periods are drivers for atmospheric TP and DSi deposition due to the correlation between the deposition and precipitation. In addition, increased spring dust deposition is a primarily driver for TN deposition changes. Increased emissions from urbanized areas in China (and likely Korea) affect the chemical properties of aerosols reaching downwind Japanese regions. Aerosol processes are responsible for increasing TN in aerosols, which are affected primarily by anthropogenic emissions. From 2000 to 2005, coal burning emissions from East Asia have contributed to an increase in TP (and possibly DSi) deposition rates. The observed increase in nutrient deposition did not noticeably impact short-term (5 year) stream water fluxes in the Kutsuki experimental forest. Due to plant uptake, the forest ecosystem retained atmospherically deposited N and P. Finally, the observed increases in nutrient deposition rates over the East China Sea and the Sea of Japan may significantly influence intra-annual net primary production. It is recommended that earth system modelling incorporate changes in atmospheric nutrient deposition rates and their impacts on the regional carbon cycle as well as aquatic and terrestrial ecosystems.     

Monte-Carlo and multifluid modelling of the circumnuclear dust coma II. Aspherical-homogeneous, and spherical-inhomogeneous nuclei

We use our newly developed Dust Monte-Carlo (DMC) simulation technique [Crifo, J.F., Lukianov, G.A., Rodionov, A.V., Zakharov, V.V., 2005. Icarus 176, 192-219] to study the dynamics of dust grains in the vicinity of some of the benchmark aspherical, homogeneous cometary nuclei and of the benchmark spherical, inhomogeneous nuclei studied by us precedingly. We use the interim unrealistic simplifying assumptions of grain sphericity, negligible nucleus rotation rate, and negligible tidal force, but take accurately into account the nucleus gravitational force, gas coma aerodynamic force, and solar radiation pressure force, and consider the full mass range of ejectable spherical grains. The resulting complicated grain motions are described in detail, as well as the resulting complicated and often counter-intuitive dust coma structure. The results are used to answer several important questions: (1) When computing coma dust distributions, (a) is it acceptable to take into consideration only one or two of the above mentioned forces (as currently done)? (b) to which accuracy must these forces be known, in particular is it acceptable to represent the gravity of an aspherical nucleus by a spherically symmetric gravity? (c) how do the more efficient but less general Dust Multi-Fluid (DMF) computations compare with the DMC results? (2) Are there simple structural relationships between the dust coma of a nucleus at small heliocentric distance rh, and that of the same nucleus at large rh? (3) Are there similarities between the gas coma structures and the associated dust coma structures? (4) Are there dust coma signatures revealing non-ambiguously a spherical nucleus inhomogeneity or an homogeneous nucleus asphericity? (5) What are the implications of the apparently quite general process of grain fall-backs for the evolution of the nucleus surface, and for the survival of a landed probe?     

A miniature laser anemometer for measurement of wind speed and dust suspension on Mars

Suspended dust is a dominant component of the Martian environment. It has a major influence on atmospheric circulation and it is deposited widely over the planetary surface causing a serious hazard to instrumentation. In order to study dust transport, quantification of the wind flow and dust concentration are vital. A simple laser-based anemometer system is presented that is able to measure suspended dust grain velocity and turbulence from a landed spacecraft. This system has advantages over other techniques of wind speed determination in being insensitive to contamination or atmospheric conditions such as temperature, pressure or composition. For the first time it would allow direct measurement of the suspended dust concentration on Mars. A prototype instrument has been constructed and successfully tested in a wind tunnel facility under simulated Martian conditions. The optics are simple in design, light weight and the instrument has low power consumption. It is also robust and the output signal is easily interpreted, producing only a small data volume. Future improvements will be discussed, specifically modification to measure wind direction, the possibility of obtaining information about dust grain size and the construction of a flight model.     

A review of East Asian reports of aurorae and comets circa AD 775

Given that a strong ¹⁴C variation in AD 775 has recently been suggested to be due to the largest solar flare ever recorded in history, it is relevant to investigate whether celestial events observed around that time may have been aurorae, possibly even very strong aurorae, or otherwise related to the ¹⁴C variation (e.g. a suggested comet impact with Earth's atmosphere). We critically review several celestial observations from AD 757 to the end of the 770s, most of which were previously considered to be true, and in some cases, strong aurorae; we discuss in detail the East Asian records and their wording. We conclude that probably none among the events after AD 770 was actually an aurora, including the event in AD 776 Jan, which was misdated for AD 774 or 775; the observed white qi phenomenon that happened above the moon in the south‐east was most probably a halo effect near the full Moon – too late in any case to be related to the ¹⁴C variation in AD 774/5. There is another report of a similar (or identical) white qi phenomenon above the moon, reported just before a comet observation and dated to AD 776 Jan; the reported comet observed by the Chinese was misdated to AD 776, but actually sighted in AD 767. Our critical review of East Asian reports of aurorae circa AD 775 shows some very likely true Chinese auroral displays observed and reported for AD 762; there were also several events prior to AD 771 that may have been aurorae but are questionable. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Coral Sr/Ca and Mg/Ca records in Con Dao Island off the Mekong Delta: Assessment of their potential for monitoring ENSO and East Asian monsoon

The climate of the South China Sea (SCS) is dominated by the East Asian monsoon (EAM) and can be related to the El Niño-Southern Oscillation (ENSO) owing to the interaction between ENSO and the EAM. An annually-banded coral (Porites sp.) collected from Con Dao Island in the southern SCS was measured for Sr/Ca and Mg/Ca ratios at near-monthly resolution through the annual bands of AD 1948–1999. This island is only ~ 90 km from the Mekong Delta coast and thus significantly influenced by riverine discharge, suggesting relatively severe environmental stress on corals. The Sr/Ca time series shows a clear annual cyclicity chiefly modulated by sea-surface temperature (SST), whereas the Mg/Ca time series exhibits an indistinct annual cyclicity, indicating that the previously-proposed coral Mg/Ca thermometry is greatly disturbed. An instrumental SST record in Con Dao Island (since 1980) has been compared with the Sr/Ca time series to calibrate a Sr/Ca thermometer. The Sr/Ca vs. SST comparison shows that the Sr/Ca thermometer is sometimes disturbed by some factor and that almost all of the disturbances occur around the annual-maximum SST in the warm/wet season. The Sr/Ca data around the annual-minimum SST in the cool/dry season is almost free from the disturbance and thus useful as a SST proxy. The disturbances of the Sr/Ca and Mg/Ca thermometers may be ascribed to the Mekong River discharge and its accompanying phenomena (i.e., large freshwater input, suspended-sediment loads, intense phytoplankton blooms, etc.), which are likely to disturb coral physiological processes. Applying the Sr/Ca thermometer to the whole Sr/Ca time series provides a SST reconstruction from 1948 through 1999. Reconstructed annual-minimum SSTs show a clear quasi-biennial oscillation significantly correlated with ENSO, indicating that the annual-minimum SST in the southern SCS tends to be higher (lower) in El Niño (La Niña) phases. This is compatible with previous observations that the East Asian winter monsoon is weakened (strengthened) in El Niño (La Niña) phases. The reconstructed SST record suggests a warming of 1.0 °C for the latter half of the 20th century. The Sr/Ca and Mg/Ca time series exhibit similar decadal-to-bidecadal variations, which do not seem to be primarily due to SST variability but rather due to some other factor possibly related to disturbance or fluctuation of coral physiological processes. Although both of our Sr/Ca and Mg/Ca data are affected, to a greater or lesser extent, by some non-temperature factor, a part of the Sr/Ca data provides a useful SST proxy and suggests that coral-based SST reconstruction in the southern SCS may be an effective means for monitoring the EAM and ENSO.     

A reanalysis of the Apollo light scattering observations, and implications for lunar exospheric dust

Conspicuous excess brightness, exceeding that expected from coronal and zodiacal light (CZL), was observed above the lunar horizon in the Apollo 15 coronal photographic sequence acquired immediately after orbital sunset (surface sunrise). This excess brightness systematically faded as the Command Module moved farther into shadow, eventually becoming indistinguishable from the CZL background. These observations have previously been attributed to scattering by ultrafine dust grains (radius ∼0.1 microns) in the lunar exosphere, and used to obtain coarse estimates of dust concentration at several altitudes and an order-of-magnitude estimate of ∼10⁻⁹ g cm⁻² for the column mass of dust near the terminator, collectively referred to as model “0”.We have reanalyzed the Apollo 15 orbital sunset sequence by incorporating the known sightline geometries in a Mie-scattering simulation code, and then inverting the measured intensities to retrieve exospheric dust concentration as a function of altitude and distance from the terminator. Results are presented in terms of monodisperse (single grain size) dust distributions. For a grain radius of 0.10 microns, our retrieved dust concentration near the terminator (∼0.010 cm⁻³) is in agreement with model “0” at z=10 km, as is the dust column mass (∼3–6×10⁻¹⁰ g cm⁻²), but the present results indicate generally larger dust scale heights, and much lower concentrations near 1 km (<0.08 cm⁻³ vs. a few times 0.1 cm⁻³ for model “0"). The concentration of dust at high altitudes (z>50 km) is virtually unconstrained by the measurements. The dust exosphere extends into shadow a distance somewhere between 100 and 200 km from the terminator, depending on the uncertain contribution of CZL to the total brightness. These refined estimates of the distribution and concentration of exospheric dust above the lunar sunrise terminator should place new and more rigorous constraints on exospheric dust transport models, as well as provide valuable support for upcoming missions such as the Lunar Atmosphere and Dust Environment Explorer (LADEE).     

Effects of two-temperature electrons, external oblique magnetic field, concentration of charged dust grains and higher-order nonlinearity on dust ion-acoustic solitary waves in Saturn's E-ring

The purpose of the present work is to investigate some nonlinear properties of the dust ion-acoustic (DIA) solitary waves in a four-component hot-magnetized dusty plasma consisting of charged dust grains, positively charged ions and two-temperature isothermal electrons. Applying a reductive perturbation theory, a nonlinear Korteweg-de Vries (KdV) equation for the first-order perturbed potential and a linear inhomogeneous KdV-type equation for the second-order perturbed potentials are derived. Stationary solutions of these coupled equations are obtained using a renormalization method. A method based on energy consideration is used to obtain a condition for stable solitons. The effects of two different types of isothermal electrons, external oblique magnetic field, concentration of negatively (positively) charged dust grains and higher-order nonlinearity on the nature of the DIA solitary waves are discussed. The numerical results are applied to Saturn's E-ring.     

Morphological and ice-dynamical changes on the Tasman Glacier, New Zealand, 1990–2007

This paper presents data concerning recent (1990–2007) surface morphological and ice-dynamical changes on the Tasman Glacier, New Zealand. We use remote-sensing data to derive rates of lake growth, glacier velocities and rates of glacier surface lowering. Between 1990 and 2007, the glacier terminus receded ~ 3.5 km and a large ice-contact proglacial lake developed behind the outwash head. By 2007 the lake area was ~ 6 km² and had replaced the majority of the lowermost 4 km of the glacier tongue. There is evidence that lake growth is proceeding at increasing rates — the lake area doubled between 2000 and 2007 alone. Measured horizontal glacier velocities decline from 150 m a^− 1 in the upper glacier catchment to almost zero at the glacier terminus and there is a consequent down-glacier increase in surface debris cover. Surface debris mapping shows that a large catastrophic rockfall onto the glacier surface in 1991 is still evident as a series of arcuate debris ridges below the Hochstetter icefall. Calculated glacier surface lowering is most clearly pronounced around the terminal area of the glacier tongue, with down-wasting rates of 4.2 ± 1.4 m a^− 1 in areas adjacent to the lateral moraine ridges outside of the current lake extent. Surface lowering rates of approximately 1.9 ± 1.4 m a^− 1 are common in the upper areas of the glacier. Calculations of future lake expansion are dependent on accurate bathymetric and bed topography surveys, but published data indicate that a further 8–10 km of the glacier is susceptible to calving and further lake development in the future.     

Comment on “Comparison of the composition of the Tempel 1 ejecta to the dust in Comet C/Hale-Bopp 1995 O1 and YSO HD 100546” by C.M. Lisse, K.E. Kraemer, J.A. Nuth III, A. Li, D. Joswiak [2007. Icarus 187, 69-86]

Lisse et al. [Lisse, C.M., Kraemer, K.E., Nuth III, J.A., Li, A., Joswiak, D., 2007. Icarus 187, 69-86] recently presented a new analysis of an ISO spectrum of Comet C/1995 O1 (Hale-Bopp), from which they claimed the identification of many new dust species. Among them are PAHs, which were not found in our first analysis of the ISO spectra. We present here a re-examination of the ISO observations of Comet Hale-Bopp. From the absence of PAHs features in the 5.25-8.5 μm region, we infer that PAHs are at least twice less abundant than derived by Lisse et al. The carbonate feature at 7.00 μm is marginally present, but lower by a factor 2 to 3 than predicted by the model of Lisse et al.     

Rebuttal to “Comment on the paper “Comparison of the composition of the Tempel 1 ejecta to the dust in Comet C/Hale-Bopp 1995 O1 and YSO HD 100546” by C.M. Lisse, K.E. Kraemer, J.A. Nuth III, A. Li, and D. Joswiak"

This response is to address the comments made by Drs. J. Crovisier and D. Bockelee-Morvan concerning the spectral analysis of Lisse et al. [Lisse, C.M., Kraemer, K.E., Nuth, J.A., Li, A., Joswiak, D., 2007. Icarus 187, 69-86] of the mid-IR ISO SWS spectrum of Comet Hale-Bopp 1995 O1 taken on October 6, 1996, and to support the conclusions made in Lisse et al. concerning the positive detection of PAHs in this comet. We also present some additional information determined from the Deep Impact and STARDUST missions, demonstrating the presence of PAHs in other comets, to support the plausibility of the Hale-Bopp PAH detection.