A numerical study of the contributions of dust source regions to the global dust budget

Contributions of the nine potential dust source regions (North and South Africa, the Arabian Peninsula, Central Asia, eastern and western China, North and South America, and Australia) to the global dust budget are investigated with a global dust transport model. A six-year simulation (1990 to 1995) indicates that the greatest contributor to the global dust budget is found to be North Africa (the Sahara Desert), which accounts for 58% of the total global dust emission and 62% of the total global dust load in the atmosphere. Australian dust dominates the southern hemisphere. The dust emission and atmospheric dust load originating from East Asia (eastern and western China) are estimated to be 214 Tg yr^− 1 and 1.1 Tg, respectively, which are 11% and 6% of the total global dust emission and dust load. Dust from East Asia dominates the atmospheric load over China and Mongolia (about 70%), Korea (60%), Japan (50%), and the North Pacific Ocean (40%). The contribution of dust originating from regions other than East Asia to the dust load over these East Asian countries and the North Pacific Ocean cannot be ignored. The simulated total dust deposition flux on Greenland suggests a possible overestimation of the Saharan dust and an underestimation of the East Asian dust in the Arctic region, which may be a common problem with global dust transport models. Possible reasons for the underestimation of the East Asian dust are discussed.     

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.     

Modeling mineral dust emissions from Chinese and Mongolian deserts

The present study investigates the frequency and intensity of mineral dust emissions over the deserts of eastern Asia from 1996 to 2001. Mineral dust emissions are simulated using a physical dust emission scheme over a region extending from 35.5°N to 47°N and from 73°E to 125°E. The input parameters required by the dust emission model are (1) surface features data including aerodynamic roughness length, soil dry size distribution and texture; and (2) meteorological surface data, mainly wind speed, soil moisture and snow cover. The way by which these surface features and meteorological data can be assessed is described and discussed. The influence of soil moisture and snow cover is taken into account and their effects on simulated dust emission are quantified.The simulations reproduce on a daily basis the location and intensity of the severe events of April 1998 and spring 2001 as recorded by the meteorological stations and/or described in various studies. Based on 6 yr of simulations, the main dust source regions are identified and their relative contributions to the total dust emissions are quantified.The seasonal cycle of the dust storms frequency is well reproduced with a maximum in spring. The simulations suggest that it is mainly controlled by the emissions occurring in the Taklimakan desert in latter spring and in summer, and by those occurring in the northern deserts of China in winter. The Taklimakan desert appears to be the most frequent and steady source of dust emissions during the studied period. On the other hand, in the Gobi desert, only a few dust emission events are simulated, but the dust amount emitted during each event is generally very large. In the northern deserts of China, dust emissions are frequent and their intensity is variable.These results show an important annual and inter-annual variability of the emitted dust (between 100 Mt yr^− 1 and 460 Mt yr^− 1), mainly controlled by the occurrence of severe events in the Gobi desert and in the northern deserts of China.     

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.     

High-resolution bispectrum speckle interferometry and two-dimensional radiative transfer modeling of the Red Rectangle

We present the first diffraction-limited K-band image of the Red Rectangle with 76 mas resolution, an H-band image with 75 mas resolution, and an RG 715 filter image ( 800 nm wavelength) with 78 mas resolution (corresponding to 25 AU for a distance of 330 pc). The H and K images were reconstructed from 6 m telescope speckle data and the RG 715 image from 2.2 m telescope data using the speckle masking bispectrum method. At all wavelengths the images show a compact, highly symmetric bipolar nebula, suggesting a toroidal density distribution of the circumstellar material. No direct light from the central binary can be seen as it is obscured by a dust disk or circumbinary torus. Our first high-resolution H−K color image of the nebula shows a broad red plateau of H−K≈ 2^m in the bright inner regions.The optical and near-infrared images and the available photometric continuum observations in a wide range of ultraviolet to centimeter wavelengths enabled us to model the Red Rectangle in detail using a two-dimensional radiative transfer code. Our model matches both the high-resolution images and the spectral energy distribution of this object very well, making the following picture much more certain. The central close binary system with a total luminosity of 3000 L is embedded in a very dense, compact circumbinary torus which has an average number density n_H ≈5×10¹² cm⁻³, an outer radius of the dense inner region of R≈30 AU (91 mas), and a ρ∝r⁻² density distribution. The full opening angle of the bipolar outflow cavities in our model is 70°. By comparing the observed and theoretical images, we derived an inclination angle of the torus to the line of sight of 7°±1°.The radiative transfer calculations show that the dust properties in the Red Rectangle are spatially inhomogeneous. The modeling confirms that the idea of large grains in the long-lived disk around the Red Rectangle (Jura et al., 1997 [ApJ, 474, 741]) is quantitatively consistent with the observations. In our models, unusually large, approximately millimeter-sized grains dominate the emission of the compact, massive torus. Models with smaller average grain sizes can possibly be found in future studies, for instance, if it turns out that the radio spectrum is not mainly caused by continuum dust emission. Therefore, the large grains suggested by our models require further confirmation by both new observations and radiative transfer calculations. Assuming a dust-to-gas ratio ρ_d/ρ_g of 0.005, the dense torus mass is 0.25 M. The model gives a lower limit of 0.0018 M, for the mass of the large particles, which produce a gray extinction of A≈ 28^m, towards the center. A much smaller mass of submicron-sized dust grains is presumably located in the polar outflow cavities, their conical surface layers, and in the outer low-density parts of the torus (where ρ∝r⁻⁴, in the region of 30 AUr 2000 AU corresponding to 0.^′′09–6^′′).     

Mid-infrared (5–100 μm) reflectance spectra and optical constants of ten phyllosilicate minerals

We have derived the real and imaginary indices of refraction for 10 phyllosilicate minerals—montmorillonite, beidellite, nontronite, hectorite, saponite, illite, illite–smectite (60/40 interlayered) kaolinite, halloysite, and serpentine—from 100–2000 cm⁻¹ (5–100 μm) at 2 cm⁻¹ spectral sampling using classical Lorentz–Lorenz dispersion theory. We present the real and imaginary indices and the oscillator parameters with which they were modeled. Use of these optical constants will aid in the modeling of thermal infrared spectra of planets, asteroids, interplanetary and interstellar dust, and protoplanetary disks around nearby stars.     

Identification of “small” dust impacts in the Ulysses dust detector data

Since October 1990, 3 weeks after the launch of the Ulysses spacecraft, the dust detector onboard recorded impacts of cosmic dust particles. Besides dust impacts, the detector recorded noise from a variety of sources. So far, a very rigid scheme had been applied to eliminate noise from impact data. The data labeled “big” dust impacts previously led to the identification of interstellar dust and of dust streams from Jupiter. The analysis presented here is concerned with data of signals of small amplitudes which are strongly contaminated by noise. Impacts identified in this data set are called “small” impacts. It is shown that dust impacts can be clearly distinguished from noise for most of the events due to the multi-coincidence characteristics of the instrument. 516 “small” impacts have been identified. For an additional 119 events, strong arguments can be given that they are probably small dust impacts. Thereby, the total number of dust impacts increases from 333 to 968 in the time period from 28 October 1990 to 31 December 1992. This increase permits a better statistical analysis, especially of the Jupiter dust streams which consist mostly of small and fast particles. Additional dust streams have been identified between the already known streams before and after Jupiter flyby. The dependence of the deflection from the Jupiter direction, the stream intensity and width on Jupiter distance support the assertion that they have been emitted from the Jovian system. The masses of the 635 “small” dust particles range from 6 × 10⁻¹⁷ to 3 × 10⁻¹⁰ g with a mean value of 1 × 10⁻¹² g, which compares to a range from 1 × 10⁻¹⁶ to 4 × 10⁻⁹ g with a mean value of 2 × 10⁻¹¹ g for the previously identified 333 “big” dust particles.     

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.     

Evidence for enhanced hydration on the northern flank of Olympus Mons, Mars

In this paper we report about a small region on the northern scarp of Olympus Mons showing an increase of the 3 μm hydration band in the OMEGA spectra, together with low superficial temperatures. Although water ice clouds can occurs on the flank of big martian volcanoes, radiative transfer modeling indicates that atmospheric water ice alone cannot justify the shape of the observed band. A fit of the 1.9–3 μm absorption features is obtained by hypothesizing that the study region consists of a mixture of dust and water ice covered by an optically thin (τ=0.08 at 3 μm) layer of dust. Thermal modeling also suggests that water ice in this region may be stable during most of the martian year due to the saturation of the atmosphere. If water ice is responsible for the observed spectral behavior, it might consist of a number of ice or snow patches possibly deposited in small depressions.     

Lidar atmospheric measurements on Mars and Earth

The LIDAR instrument operating from the surface of Mars on the Phoenix Mission measured vertical profiles of atmospheric dust and water ice clouds at temperatures around −65 °C. An equivalent lidar system was utilized for measurements in the atmosphere of Earth where dust and cloud conditions are similar to Mars. Coordinated aircraft in situ sampling provided a verification of lidar measurement and analysis methods and also insight for interpretation of lidar derived optical parameters in terms of the dust and cloud microphysical properties. It was found that the vertical distribution of airborne dust above the Australian desert is quite similar to what is observed in the planetary boundary layer above Mars. Comparison with the in situ sampling is used to demonstrate how the lidar derived optical extinction coefficient is related to the dust particle size distribution. The lidar measurement placed a constraint on the model size distribution that has been used for Mars. Airborne lidar measurements were also conducted to study cirrus clouds that form in the Earth’s atmosphere at a similar temperature and humidity as the clouds observed with the lidar on Mars. Comparison with the in situ sampling provides a method to derive the cloud ice water content (IWC) from the Mars lidar measurements.     

Brief description of bare essentials of surface transfer and results from simulations with the HAPEX-MOBILHY data

The Bare Essentials of Surface Transfer (BEST) land surface scheme is briefly described and the key physical parameterisations discussed. Results are then presented to illustrate how the model performs, with forcing data for HAPEX-MOBILHY, compared to a series of other schemes in the simulation of evaporation and sensible heat. The implications of the models performance, and some indications of the future development of the scheme are provided. The basic version of BEST was found to overestimate evaporation for the HAPEX-MOBILHY data, simulating 816 mm yr⁻¹ compared to a range of 550 to 816 mm yr⁻¹ for all models. A modification to the β parameterisation reduced the evaporation to 759 mm yr⁻¹ which, although an improvement, is still probably too high.     

The dust environment at Titan

To explain the observed abundances of CO₂ in Titan's atmosphere, a relatively high water deposition into the atmosphere needs to be invoked due to the importance of H₂O photolysis in CO₂ production. A likely source of H₂O is icy dust particles from space. This paper considers the direct dust input to Titan's atmosphere from the interplanetary environment, and also ejecta particles from micrometeoroid impacts with the icy satellites Hyperion, Iapetus and Phoebe. It is found that the likely mass influx to Titan is 10^–16 to 10^–15 kg m^–2 s^–1. This mass influx is an order of magnitude too low to explain the observed levels of CO₂ in Titan's atmosphere in the context of a recent photochemical model. This leads one to speculate as to the likelihood of one large impact to Titan in the recent past;i.e., that the atmosphere is not in equilibrium but is cnrrently losing CO₂.     

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.     

Chemical sources of haze formation in Titan's atmosphere

A prominent feature of Titan's atmosphere is a thick haze region that acts as the end product of hydrocarbon and nitrile chemistry. Using a one-dimensional photochemical model, an investigation into the chemical mechanisms responsible for the formation of this haze region is conducted. The model derives profiles for Titan's atmospheric constituents that are consistent with observations. Included is an updated benzene profile that matches more closely with—recent ISO observations (Icarus 161 (2003) 383), replacing the profile given in the benzene study of Wilson et al. (J. Geophys. Res. 108 (2003) 5014). Using these profiles, pathways from polyynes, aromatics, and nitriles are considered, as well as possible copolymerization among the pathways. The model demonstrates that the growth of polycyclic aromatic hydrocarbons throughout the lower stratosphere plays an important role in furnishing the main haze layer, with nitriles playing a secondary role. The peak chemical production of haze layer ranges from 140 to 300 km peaking at an altitude of 220 km, with a production rate of 3.2×10⁻¹⁴ gcm⁻² s⁻¹. Possible mechanisms for polymerization and copolymerization and suggestions for further kinetic study are discussed, along with the implications for the distribution of haze in Titan's atmosphere.     

Interannual Variability of Martian Global Dust Storms: Simulations with a Low-Order Model of the General Circulation

We present simulations of the interannual variability of martian global dust storms (GDSs) with a simplified low-order model (LOM) of the general circulation. The simplified model allows one to conduct computationally fast long-term simulations of the martian climate system. The LOM is constructed by Galerkin projection of a 2D (zonally averaged) general circulation model (GCM) onto a truncated set of basis functions. The resulting LOM consists of 12 coupled nonlinear ordinary differential equations describing atmospheric dynamics and dust transport within the Hadley cell. The forcing of the model is described by simplified physics based on Newtonian cooling and Rayleigh friction. The atmosphere and surface are coupled: atmospheric heating depends on the dustiness of the atmosphere, and the surface dust source depends on the strength of the atmospheric winds. Parameters of the model are tuned to fit the output of the NASA AMES GCM and the fit is generally very good.Interannual variability of GDSs is possible in the LOM, but only when stochastic forcing is added to the model. The stochastic forcing could be provided by transient weather systems or some surface process such as redistribution of the sand particles in storm generating zones on the surface. The results are sensitive to the value of the saltation threshold, which hints at a possible feedback between saltation threshold and dust storm activity. According to this hypothesis, erodable material builds up as a result of a local process, whose effect is to lower the saltation threshold until a GDS occurs. The saltation threshold adjusts its value so that dust storms are barely able to occur.     

DuneXpress

The DuneXpress observatory will characterize interstellar and interplanetary dust in-situ, in order to provide crucial information not achievable with remote sensing astronomical methods. Galactic interstellar dust constitutes the solid phase of matter from which stars and planetary systems form. Interplanetary dust, from comets and asteroids, represents remnant material from bodies at different stages of early solar system evolution. Thus, studies of interstellar and interplanetary dust with DuneXpress in Earth orbit will provide a comparison between the composition of the interstellar medium and primitive planetary objects. Hence DuneXpress will provide insights into the physical conditions during planetary system formation. This comparison of interstellar and interplanetary dust addresses directly themes of highest priority in astrophysics and solar system science, which are described in ESA’s Cosmic Vision. The discoveries of interstellar dust in the outer and inner solar system during the last decade suggest an innovative approach to the characterization of cosmic dust. DuneXpress establishes the next logical step beyond NASA’s Stardust mission, with four major advancements in cosmic dust research: (1) analysis of the elemental and isotopic composition of individual interstellar grains passing through the solar system, (2) determination of the size distribution of interstellar dust at 1 AU from 10^− 14 to 10^− 9 g, (3) characterization of the interstellar dust flow through the planetary system, (4) establish the interrelation of interplanetary dust with comets and asteroids. Additionally, in supporting the dust science objectives, DuneXpress will characterize dust charging in the solar wind and in the Earth’s magnetotail. The science payload consists of two dust telescopes of a total of 0.1 m² sensitive area, three dust cameras totaling 0.4 m² sensitive area, and a nano-dust detector. The dust telescopes measure high-resolution mass spectra of both positive and negative ions released upon impact of dust particles. The dust cameras employ different detection methods and are optimized for (1) large area impact detection and trajectory analysis of submicron sized and larger dust grains, (2) the determination of physical properties, such as flux, mass, speed, and electrical charge. A nano-dust detector searches for nanometer-sized dust particles in interplanetary space. A plasma monitor supports the dust charge measurements, thereby, providing additional information on the dust particles. About 1,000 grains are expected to be recorded by this payload every year, with 20% of these grains providing elemental composition. During the mission submicron to micron-sized interstellar grains are expected to be recorded in statistically significant numbers. DuneXpress will open a new window to dusty universe that will provide unprecedented information on cosmic dust and on the objects from which it is derived.     

Characterization of Asian Dust during ACE-Asia

ACE-Asia was a multi-national collaboration organized to investigate and understand the chemistry, radiative properties, and climatic effects of mineral dust and other aerosol particles in the East Asia/Northwest Pacific region. Studies conducted at the Gosan and Zhenbeitai surface supersites show striking variations in aerosol concentrations and properties that were affected by the occurrence and origins of the Asian dust storms, air mass pathways, and mixing during the transport. Investigations conducted on the R/V Ronald H. Brown (RHB) showed that dust had a pervasive influence on the chemical composition, size distribution, and optical properties of the aerosol. Analyses using an aerosol time-of-flight mass spectrometer on the RHB showed that most of the coarse-particle nitrate and sulfate in post-frontal air was associated with dust, and more remarkably, that competitive or exclusionary processes evidently are involved in the uptake or production of these substances. Studies conducted onboard the NCAR C-130 aircraft showed that coarse-mode dust was less absorbing and less hygroscopic than pollution aerosol and that there was little correlation in light scattering and absorption by the sub- vs. super-micrometer aerosol. Below 2 km, dust was commonly mixed with pollutants, and this had a stronger influence on the optical properties of the submicrometer particles than the coarse-mode dust; at higher altitudes, the dust was less affected by pollution. Single particle analyses of C-130 samples showed that the mixing of black carbon (BC) with dust was common, but only certain types of BC particles were aggregated. Models were used in the planning, execution and interpretative phases of ACE-Asia; and summaries of modeling results are presented to illustrate the progress being made in identifying new dust sources; in depicting the time-varying, three-dimensional structure of dust plumes; and in quantifying the production, transport, and deposition of Asian dust.     

Radiative cooling within illuminated layers of dust on (pre)-planetary surfaces and its effect on dust ejection

Temperature gradients in dust beds embedded in a low pressure gaseous environment induce a lift of particles under certain conditions. This effect can erode planetesimals and enables entrainment of dust into the martian atmosphere. Here, we consider a numerical model to calculate the temperature profile in a dust bed which is subject to illumination. We consider the situation when the illumination is switched on and heats the dust bed’s surface and when it is switched off again after a certain time. The calculations focus on the heat transfer by infrared radiation within the dust layer. We find that radiative transfer within the dust bed modifies the absolute temperatures and temperature gradients significantly. This is important for effects which are sensitive to absolute temperatures, i.e. ice sublimation or melting of solids. For low thermal conductivity dust beds of 0.001 W m⁻¹ K⁻¹ it determines the temperature structure of the dust. For higher thermal conductivities the modifications are moderate with respect to dust eruptions as the order of magnitude of temperature gradients stays the same.     

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.     

Comparison of Mediterranean sea level fields for the period 1961–2000 as given by a data reconstruction and a 3D model

Two Mediterranean sea level distributions spanning the last decades are examined. The first one is a reconstruction of sea level obtained by a reduced-space optimal interpolation applied to tide gauge and altimetry data. The second distribution is obtained from a 3D (baroclinic) regional circulation model. None of the two representations includes the mechanical atmospheric forcing. Results are presented for two different periods: 1993–2000 (for which altimetry data are available) and 1961–2000 (the longest period common to both distributions).The first period is examined as a test period for the model, since the reconstruction is very similar to altimetry observations. The modelled sea level is in fair agreement with the reconstruction in the Western Mediterranean and in the Aegean Sea (except in the early nineties), but in the Ionian Sea the model departs from observations. For the whole period 1961–2000 the main feature is a marked positive trend in the Ionian Sea (up to 1.8 mm yr^− 1), observed both in the reconstruction and in the model. Also the distribution of positive trends in the Western Mediterranean (mean value of 1.1 mm yr^− 1) and the smaller trends in the Aegean Sea (0.5 mm yr^− 1) are similar in the reconstruction and in the model, despite the first implicitly accounts for sea level variations due to remote sources such as ice melting and the second does not. The interannual sea level variability associated with key regional events such as the Eastern Mediterranean Transient is apparently captured by the reconstruction but not by the model (at least in its present configuration). Hence, the reconstruction can be envisaged as a useful tool to validate further long-term numerical simulations in the region.