FY-3气象卫星紫外臭氧总量探测仪辐亮度在轨定标与反演结果分析

FY-3气象卫星上搭载的紫外臭氧总量探测仪TOU是我国自主开发研制的首台用于全球臭氧总量定量测量的探测仪,自发射以来已成功在轨运行近两年.由于TOU发射前辐亮度定标存在偏差,为了得到高精度的产品,TOU必须进行在轨定标.本文介绍了基于辐射传输模式计算对TOU辐亮度进行在轨道定标的方法,定标过程中用于模拟辐亮度计算的臭氧总量由与TOU观测时刻相近的国外臭氧总量探测仪器MetOp/GOME-2提供.文章将在轨定标后TOU的反演结果与AURA/OMI以及地基的产品进行比较,结果表明,用辐射传输模式对TOU辐亮度进行在轨定标的方法是可行的,反演结果能够真实地反映臭氧的时空分布特性,在全球部分地基观测站所处的位置上对TOU, OMI以及地基的臭氧总量进行比较的结果表明,TOU与OMI的相对偏差均方根约为2.52%,TOU与地基以及OMI与地基观测结果之间的相对偏差均方根分别为4.45%和3.89%.     

Trends in total ozone and the effect of the equatorial zonal wind QBO

Trends in total column ozone have been analyzed in terms of the equatorial zonal wind. We used zonal monthly mean total ozone from Total Ozone Mapping Spectrometer (TOMS) and monthly mean zonal wind in the equatorial stratosphere at 30 hPa to define the phases of the quasi-biennial oscillation (QBO). Total column ozone trends have been assessed during the period 1979–2004, for both Hemispheres, and for each month, under three conditions considering, all the ozone dataset, ozone values during easterly phase and ozone values during westerly phase of the QBO. When the whole dataset is considered, negative trends are observed. From low to midlatitudes a zonal pattern is noticed with increasing negative values toward higher latitudes. When the data is filtered according to the QBO phase, statistically significant positive trends appear in the westerly case during January to May at low latitudes .The trend pattern in the case of the easterly phase presents more negative values.     

太湖沉积物再悬浮模拟方法

由于风浪对浅水湖泊底泥的强烈扰动作用,底泥再悬浮的生态效应研究成为目前国际研究热点之一.但方法的缺乏极大地限制了有关底泥再悬浮环境效应的深入研究.本研究应用一种新的沉积物再悬浮装置,模拟了太湖南部长兜港水域不同风浪影响下悬浮物的垂向分布,获得水柱总悬浮物量(T)与扰动频率(n)的定量关系:T=19.77×exp(n/2.61)-71.7.结合现场实测风情下的水柱总悬浮物量,建立了室内水动力条件、水柱总悬浮物量和现场风情的定量关系.研究2004年10月至2005年10月太湖全年代表性小风、中风及大风与再悬浮装置的扰动频率对应关系,估算了小风、中风和大风常规风情引起的总悬浮物量分别为159、230、425 g/m2.结果表明,模拟获得的太湖悬浮物垂向分布规律与现场状况较为接近,是目前较为适用的浅水水体沉积物再悬浮模拟方法.     

Remote sensing of the global distribution of total ozone and the inferred upper-tropospheric circulation from Nimbus IRIS experiments

The global distribution of total ozone is derived for the period April, May, June and July of 1969 from Nimbus-3 Infrared Interferometer Spectrometer (IRIS) experiment. Preliminary estimates of ozone amounts from Nimbus-4 IRIS for the same period of 1970 show similar results. The standard error of estimation of total ozone from both IRIS experiments is 6% with respect to Dobson Spectrophotometer measurements. A systematic variation in the ozone distribution from April to July in the tropical, middle and polar latitudes is observed indicating the changes in the lower stratospheric circulation.The total ozone measurements show a strong correlation with the upper tropospheric geopotential height in the extratropical latitudes. From this relationship total ozone is used as a quasi-stream function to deduce geostrophic winds at the 200 mb level over extratropical regions of the northern and southern hemispheres. These winds reveal the subtropical and polar jet streams over the globe.Allied research associates.     

A quarter century of TOMS total column ozone measurements over Brazil

The Earth's atmosphere has been probed since 1978 by the total ozone mapping spectrometer (TOMS) experiment. In polar orbits around the Earth, Nimbus 7, Meteor 3, and Earth Probe satellites have carried TOMS to measure the backscattered solar ultraviolet radiation, which is processed by Version 8 TOMS algorithm for the retrieval of column ozone values. These TOMS data and their behavior are reported and investigated in this work for some sites in the equatorial and tropical zones of Brazil. Comparison with measurements from Brewer and Dobson spectrophotometers in some cases is provided. An average difference as low as 1.5% between TOMS and Brewer measurements was obtained. In addition, more than two decades of TOMS monthly average column ozone measurements were applied to the calculation of ozone trends. Although small and in general negative, these trends were statistically significant for all sites in the winter season.     

Statistical interpolation of ozone measurements from satellite data (TOMS, SBUV and SAGE II) using the kriging method

This study demonstrates that ordinary kriging in spherical coordinates using experimental semi-variograms provides highly usable results, especially near the pole in winter and/or where there could be data missing over large areas. In addition, kriging allows display of the spatial variability of daily ozone measurements at different pressure levels. Three satellite data sets were used: Total Ozone Mapping Spectrometer (TOMS) data, Solar Backscattered UltraViolet (SBUV), and the Stratospheric Aerosol and Gas Experiment (SAGE II) ozone profiles. Since SBUV is a nadir-viewing instrument, measurements are only taken along the sun-synchronous polar orbits of the satellite. SAGE II is a limb-viewing solar occultation instrument, and measurements have high vertical resolution but poor daily coverage. TOMS has wider coverage with equidistant distribution of data (resolution 1° × 1.25°) but provides no vertical information. Comparisons of the resulting SBUV-interpolated (column-integrated) ozone field with TOMS data are strongly in agreement, with a global correlation of close to 98%. Comparisons of SBUV-interpolated ozone profiles with daily SAGE II profiles are relatively good, and comparable to those found in the literature. The interpolated ozone layers at different pressure levels are shown.     

Evidence for vertical ozone redistribution since 1967

Long-term measurements of the ozone concentration in the vicinity of the city of Berlin have been performed with ground based Dobson spectrophotometers and balloon borne systems. The respective experiments cover the past 24 years. All data have been reevaluated and corrected towards uniform calibration standards, leading to the longest European data set of total column density, altitude-dependent ozone partial pressures and the corresponding temperatures. Smoothing algorithms unravel significant long-term trends.The analysis shows an increase of ozone concentration within the middle stratosphere (below 31 km height) as well as in the troposphere over the past 24 years. On the contrary, ongoing ozone depletion in the lower stratosphere has been found.The large scale vertical redistribution of atmospheric ozone in the troposphere and the lower stratosphere seems to be in agreement with model calculations and trend predictions that have their roots in changes of the chemical composition and the ozone photochemistry due to anthropogenically induced trace gas concentrations.Deutscher Wetterdienst, Meteorologisches Observatorium Potsdam.Deutscher Wetterdienst, Meteorologisches Observatorium Lindenberg.     

Evidence of global-scale waves with zonal wave number zero in the stratosphere

Coherency spectra derived from time series of stratospheric quantities indicate oscillations in the frequency range below 0.5 d^–1 which are correlated on a global scale. Satellite observations of total ozone and stratospheric radiance (BUV and SIRS, Nimbus4, April–November 1970) have been used to derive phase relationships of such oscillations. As an example, an oscillation of total ozone with a period of 7.5 d and zonal wave number zero is analyzed in detail. The basic assumption is made and tested, that the oscillation reflects stratospheric planetary waves as obtained from Laplace's tidal equations. The observed latitudinal phase shifts for the total ozone oscillation are in good agreement with theoretical predictions. It is concluded from the observations of ozone and radiance that mainly divergence effects related to global-scale waves are responsible for the 7.5 d oscillations of total ozone at high and middle latitudes and at the equator whereas in the latitude range 10°S–20°S predominantly temperature effects are important. Meridional wind amplitudes of some 10 cm/s are sufficient to explain the high and mid-latitude ozone oscillations. At low latitudes vertical wind amplitudes of about 0.2 mm/s corresponding to height changes of the ozone layer of roughly ±20 m are obtained.     

Synoptic-scale fluctuations of total ozone in the atmosphere

A model, based on ozone-concentration tendency equation, is developed to study synoptic ozone-column variations. The application is referred to a middle-latitude site and to an atmospheric layer extending from the surface up to about 35-km altitude. Photochemical effects at the considered location for synoptic time scales are considered negligible. The data input consists of umkehr ozone profile, total ozone (obtained by Brewer No. 067, located at Rome) and horizontal wind at various levels. Analysis of several cases indicates that meridional advection is the main factor responsible for the observed synoptic-scale ozone fluctuations.     

Estimation of tropical cyclone parameters and wind fields from SAR images

The traditional method of Synthetic Aperture Radar(SAR)wind field retrieval is based on an empirical relation between the near surface winds and the normalized radar backscatter cross section to estimate wind speeds,where this relation is called the geophysical model function(GMF).However,the accuracy rapidly decreases due to the impact of rainfall on the measurement of SAR and the saturation of backscattered intensity under the condition of tropical cyclone.Because of no available instrument synchronously monitoring rain rate on the satellite platform of SAR,we have to derive the precipitation of the SAR observation time from non-simultaneous passive microwave observations of rain in combination with geostationary IR images,and then use the model of rain correction to remove the impact of rain on SAR wind field measurements.For the saturation of radar backscatter cross section in high wind speed conditions,we develop an approach to estimate tropical cyclone parameters and wind fields based on the improved Holland model and the SAR image features of tropical cyclone.To retrieve the low-to-moderate wind speed,the wind direction of tropical cyclone is estimated from the SAR image using wavelet analysis.And then the maximum wind speed and the central pressure of tropical cyclone are calculated by a least square minimization of the difference between the improved Holland model and the low-to-moderate wind speed retrieved from SAR.In addition,wind fields are estimated from the improved Holland model using the above-mentioned parameters of tropical cyclone as input.To evaluate the accuracy of our approach,the SAR images of typhoon Aere,typhoon Khanun,and hurricane Ophelia are used to estimate tropical cyclone parameters and wind fields,which are compared with the best track data and reanalyzed wind fields of the Joint Typhoon Warning Center(JTWC)and the Hurricane Research Division(HRD).The results indicate that the tropical cyclone center,maximum wind speed,and central pressure are generally consistent with the best track data,and wind fields agree well with reanalyzed data from HRD.     

The statistical and numerical study of the global distribution of coronal plasma and magnetic field near 2.5 Rs over a 10-year period

The basic characteristics of the global distribution for the corona plasma and magnetic field near 2.5 Rs are analyzed with the statistical and numerical methods for 136 Carrington Rotations (CRs) covering four different phases of solar activity. By using the observational data and the velocity distribution model in the corona, the statistical average distribution of the magnetic field, density and the coronal mass outputs are analyzed for the four different phases. Then, a numerical study of the global distribution near 2.5 Rs has been made by solving a self-consistent MHD system. Finally, the solar wind speed at 1 AU is given by mapping the speed at 2.5 Rs to that near 1 AU, and the comparison of the numerical results with the observational measurements and the simulation result of the Wang–Sheeley–Arge (WSA) model are made during more than 5 years. The numerical results indicate that the global distributions on the source surface of 2.5 Rs at different phases of solar activity could be used to predict the change of the solar wind in interplanetary space.     

The Nimbus-4 Backscatter Ultraviolet (BUV) atmospheric ozone experiment — tow years' operation

Summary In April 1970 the Backscatter Ultraviolet (BUV) experiment was placed into orbit aboard the Nimbus-4 satellite. This double monochromator experiment measures ultraviolet terrestrial radiance at twelve discrete wavelengths between 2550 Å and 3400 Å. Approximately 100 scans covering a 230 km square are made between terminator crossings on the daylight side of the earth. A colinear photometer channel with the same field of view is used to derive the Lambert reflectivity of the lower boundary of the scattering atmosphere. The extraterrestrial solar irradiance is measured at the northern terminator. The instrument has currently produced almost three years of nearly continuous data which are being used to infer the high-level ozone distribution and total ozone on a global basis. The high-level ozone data have been verified by independent coincident rocket ozone soundings, and the total ozone values show good agreement with Dobson spectrophotometer determinations as well as those made with the Infrared Interferometer Spectrometer also on Nimbus-4. An increase has been observed in equatorial radiance at 2550 Å relative to 2900 Å, which seems to indicate that the amount of ozone in the upper stratosphere is related to the eleven-year solar cycle.     

GPU-accelerated computing of three-dimensional solar wind background

High-performance computational models are required to make the real-time or faster than real-time numerical prediction of adverse space weather events and their influence on the geospace environment. The main objective in this article is to explore the application of programmable graphic processing units (GPUs) to the numerical space weather modeling for the study of solar wind background that is a crucial part in the numerical space weather modeling. GPU programming is realized for our Solar-Interplanetary-CESE MHD model (SIP-CESE MHD model) by numerically studying the solar corona/interplanetary solar wind. The global solar wind structures are obtained by the established GPU model with the magnetic field synoptic data as input. Meanwhile, the time-dependent solar surface boundary conditions derived from the method of characteristics and the mass flux limit are incorporated to couple the observation and the three-dimensional (3D) MHD model. The simulated evolution of the global structures for two Carrington rotations 2058 and 2062 is compared with solar observations and solar wind measurements from spacecraft near the Earth. The MHD model is also validated by comparison with the standard potential field source surface (PFSS) model. Comparisons show that the MHD results are in good overall agreement with coronal and interplanetary structures, including the size and distribution of coronal holes, the position and shape of the streamer belts, and the transition of the solar wind speeds and magnetic field polarities.     

Numerical modelling of pollutant dispersion in sea breeze conditions

The characteristics of air pollution in Tarragona (Spain) were investigated. Tarragona has an important petrochemical industry in a coastal region with a complex terrain. The numerical study was made in sea breeze conditions with a three-dimensional mesoscale model. Temporal and spatial variations of the wind fields have been used in the Eulerian equation for a non-reactive pollutant. The results of this study reveal the complexity of the dispersion patterns due to the combined effects of the sea breeze circulation and the orography. This work presents a comparison between the model output and the observed wind data by sodar and surface wind measurements. The evaluation shows that the model is capable of providing very realistic wind fields within this domain.     

High-altitude data assimilation system experiments for the northern summer mesosphere season of 2007

A global numerical weather prediction system is extended to the mesosphere and lower thermosphere (MLT) and used to assimilate high-altitude satellite measurements of temperature, water vapor and ozone from MLS and SABER during May–July 2007. Assimilated temperature and humidity from 100 to 0.001 hPa show minimal biases compared to satellite data and existing analysis fields. Saturation ratios derived diagnostically from these assimilated temperature and water vapor fields at PMC altitudes and latitudes compare well with seasonal variations in PMC frequency measured from the aeronomy of ice in the mesosphere (AIM) satellite. Synoptic maps of these diagnostic saturation ratios correlate geographically with three independent transient mesospheric cloud events observed at midlatitudes by SHIMMER on STPSat-1 and by ground observers during June 2007. Assimilated temperatures and winds reveal broadly realistic amplitudes of the quasi 5-day wave and migrating tides as a function of latitude and height. For example, analyzed winds capture the dominant semidiurnal MLT wind patterns at 55°N in June 2007 measured independently by a meteor radar. The 5-day wave and migrating diurnal tide also modulate water vapor mixing ratios in the polar summer MLT. Possible origins of this variability are discussed.     

Ozone content of the atmosphere. 2. Variability of the total ozone content in the atmosphere and the horizontal wind in the lower thermosphere (Central Europe)

The correlation between total ozone content lower thermosphere horizontal wind parameters, and standard indices of solar activity and geomagnetic activity has been studied. The satellite measurements of TOC for five observatories in Central Europe and the lower thermosphere wind measurements for Collm observatory (Germany) were used for 1996–2003. The quasi-periodic structure of these variations and the correlation between the corresponding periodograms were also studied. The quantitative evaluation of the statistically significant correlations and common periodicities were revealed.     

Temperature and ozone variations in the stratosphere

Examined are temperature and ozone variations in the Northern Hemisphere stratosphere during the period 1958–77, as estimated from radiosondes rocketsondes, ozonesondes, and Umkehr measurements. The temperature variation in the low tropical stratosphere is a combination of the variation associated with the quasi-biennial oscillation, and a variation nearly out of phase with the pronounced 3-yearly temperature oscillation (Southern Oscillation) present in the tropical troposphere since 1963. Based on radiosonde and rocketsonde data, the quasibiennial temperature oscillation can be traced as high as the stratopause, the phase varying with both height and latitude. However, the rocketsonde-derived temperature decrease of several degrees Celsius in the 25–55 km layer of the Western Hemisphere between 1969 (sunspot maximum) and 1976 (sunspot minimum) is not apparent in high-level radiosonde data, so that caution is advised with respect to a possible solar-terrestrial relation.There has been a strong quasi-biennial oscillation in ozone in the 8–16 km layer of the north polar region, with ozone minimum near the time of quasi-biennial west wind maximum at a height of 20 km in the tropics. A quasi-biennial oscillation in ozone (of similar phase) is also apparent from both ozonesonde data and Umkehr measurements in 8–16 and 16–24 km layers of north temperate latitudes, but not higher up. Both measurement techniques also suggest a slight overall ozone decrease in the same layers between 1969 and 1976, but no overall ozone change in the 24–32 km layer. Umkehr measurements indicate a significant 6–8% increase in ozone amount in all stratospheric layers between 1964 and 1970, and in 1977 the ozone amount in the 32–46 km layer was still 4% above average despite the predicted depletion due to fluorocarbon emissions. The decrease in ozone in the 32–46 km, layer of mid latitudes following the volcanic eruptions of Agung and Fuego is believed to be mostly fictitious and due to the bias introduced into the Umkehr technique by stratospheric aerosols of volcanic origin. Above-average water vapor amounts in the low stratosphere at Washington, DC, appear closely related to warm tropospheric temperatures in the tropics, presumably reflecting variations in strength of the Hadley circulation.     

Quantitative models of the fallout and dispersal of tephra from volcanic eruption columns

A theoretical model of clast fallout from convective eruption columns has been developed which quantifies how the maximum clast size dispersal is determined by column height and wind strength. An eruption column consists of a buoyant convecting region which rises to a heightH _B where the column density equals that of the atmosphere. AboveH _B the column rises further to a heightH _T due to excess momentum. BetweenH _T andH _B the column is forced laterally into the atmosphere to form an upper umbrella region. Within the eruption column, the vertical and horizontal velocity fields can be calculated from exprimental and theoretical studies and consideration of mass continuity. The centreline vertical velocity falls as a nearly linear function over most of the column's height and the velocity decreases as a gaussian function radially away from the centreline. Both column height and vertical velocity are strong functions of magma discharge rate. From calculations of the velocity field and the terminal fall velocity of clasts, a series of particle support envelopes has been constructed which represents positions where the column vertical velocity and terminal velocity are equal for a clast of specific size and density. The maximum range of a clast is determined in the absence of wind by the maximum width of the clast support envelope.The trajectories of clasts leaving their relevant support envelope at its maximum width have been modelled in columns from 6 to 43 km high with no wind and in a wind field. From these calculations the shapes and areas of maximum grain size contours of the air-fall deposit have been predicted. For the no wind case the theoretical isopleths show good agreement with the Fogo A plinian deposit in the Azores. A diagram has been constructed which plots, for a particular clast size, the maximum range normal to the dispersal axis against the downward range. From the diagram the column height (and hence magma discharge rate) and wind velocity can be determined. Historic plinian eruptions of Santa Maria (1902) and Mount St. Helens (1980) give maximum heights of 34 and 19 km respectively and maximum wind speeds at the tropopause of m/s and 30 m/s respectively. Both estimates are in good agreement with observations. The model has been applied to a number of other plinian deposits, including the ultraplinian phase of theA.D. 180 Taupo eruption in New Zealand which had an estimated column height of 51 km and wind velocity of 27 m/s.     

探讨利用SCIAMACHY数据反演温室气体二氧化碳

自工业革命以来CO2已经成为引起全球气候变暖的一种最主要的人为温室气体,由此引起的全球气候变化是当今社会经济可持续发展面临的最严峻挑战之一.为了更好的研究二氧化碳在时间与空间尺度上的格局及变化,综合分析CO2与太阳辐射、地面反照率、气溶胶等要素之间的关系,本文回顾了以往利用新型卫星传感器获取多源遥感数据的特点,分析了多种反演计算方法的优缺点,并结合地面实测数据进行对比分析的方法与理论,全面阐述国际上有关CO2柱浓度反演的思想.作者认为目前在地面观测站较少的情况下,针对先进的SCIAMACHY传感器数据利用DOAS相关算法及BESD算法进行CO2柱体浓度反演的数据连续,技术成熟,并由此提出了下一步工作中应解决的3个问题.