ALBERT

Description: The purpose of the chemistry component of the model comparison is to assess to what extent differences in the formulation of chemical processes explain the variance between model results. Observed concentrations of chemical compounds are used to estimate to what degree the various models represent realistic situations. For readability, the materials for the chemistry experiment are reported in three separate sections. This section discussed the data used to evaluate the models in their simulation of the source gases and the Nitrogen compounds (NO(y)) and Chlorine compounds (Cl(y)) species.

Description: Space-based and airborne coherent Doppler lidars designed for measuring global tropospheric wind profiles in cloud-free air rely on backscatter, beta from aerosols acting as passive wind tracers. Aerosol beta distribution in the vertical can vary over as much as 5-6 orders of magnitude. Thus, the design of a wave length-specific, space-borne or airborne lidar must account for the magnitude of 8 in the region or features of interest. The SPAce Readiness Coherent Lidar Experiment under development by the National Aeronautics and Space Administration (NASA) and scheduled for launch on the Space Shuttle in 2001, will demonstrate wind measurements from space using a solid-state 2 micrometer coherent Doppler lidar. Consequently, there is a critical need to understand variability of aerosol beta at 2.1 micrometers, to evaluate signal detection under varying aerosol loading conditions. Although few direct measurements of beta at 2.1 micrometers exist, extensive datasets, including climatologies in widely-separated locations, do exist for other wavelengths based on CO2 and Nd:YAG lidars. Datasets also exist for the associated microphysical and chemical properties. An example of a multi-parametric dataset is that of the NASA GLObal Backscatter Experiment (GLOBE) in 1990 in which aerosol chemistry and size distributions were measured concurrently with multi-wavelength lidar backscatter observations. More recently, continuous-wave (CW) lidar backscatter measurements at mid-infrared wavelengths have been made during the Multicenter Airborne Coherent Atmospheric Wind Sensor (MACAWS) experiment in 1995. Using Lorenz-Mie theory, these datasets have been used to develop a method to convert lidar backscatter to the 2.1 micrometer wavelength. This paper presents comparison of modeled backscatter at wavelengths for which backscatter measurements exist including converted beta (sub 2.1).

Description: The typical fair weather electric field at the ground is between -100 and -300 V/m. At the NASA Kennedy Space Center and US Air Force Cape Canaveral Air Station (KSC) the electric field at the ground sometimes reaches -400 to -1200 V/m within an hour or two after sunrise on days that otherwise seem to be fair weather. We refer to the enhanced negative electric fields as the "sunrise enhancement." To investigate the sunrise enhancement at KSC we measured the electric field (E) in the first few hundred meters above the ground before and during several sunrise enhancements. From these E soundings we can infer the presence of charge layers and determine their thickness and charge density.

Description: This section contains a number of special diagnostics that are designed to examine certain mechanisms. Section 1 reports on the method used to test the photochemical partitioning in the models. Sections 2 and 3 represent efforts to examine the model calculated production and removal rates for ozone and how the values are combined with transport rates in the models to produce the simulated ozone distributions. Sections 4 and 5 concentrate on polar processes including the dynamics aspect of vortex confinement and the chemical aspects of chlorine activation.

Description: MM II defined a series of experiments to better understand and characterize model transport and to assess the realism of this transport by comparison to observations. Measurements from aircraft, balloon, and satellite, not yet available at the time of MM I [Prather and Remsberg, 1993], provide new and stringent constraints on model transport, and address the limits of our transport modeling abilities. Simulations of the idealized tracers the age spectrum, and propagating boundary conditions, and conserved HSCT-like emissions probe the relative roles of different model transport mechanisms, while simulations of SF6 and C02 make the connection to observations. Some of the tracers are related, and transport diagnostics such as the mean age can be derived from more than one of the experiments for comparison to observations. The goals of the transport experiments are: (1) To isolate the effects of transport in models from other processes; (2) To assess model transport for realistic tracers (such as SF6 and C02) for comparison to observations; (3) To use certain idealized tracers to isolate model mechanisms and relationships to atmospheric chemical perturbations; (4) To identify strengths and weaknesses of the treatment of transport processes in the models; (5) To relate evaluated shortcomings to aspects of model formulation. The following section are included:Executive Summary, Introduction, Age Spectrum, Observation, Tropical Transport in Models, Global Mean Age in Models, Source-Transport Covariance, HSCT "ANOY" Tracer Distributions, and Summary and Conclusions.

Description: A global lightning model that includes diurnal and annual lightning variation, and total flash density versus latitude for each major land and ocean, has been used as the basis for simulating the global electric circuit charging rate. A particular objective has been to reconcile the difference in amplitude ratios [AR=(max-min)/mean] between global lightning diurnal variation (AR approx. = 0.8) and the diurnal variation of typical atmospheric potential gradient curves (AR approx. = 0.35). A constraint on the simulation is that the annual mean charging current should be about 1000 A. The global lightning model shows that negative ground flashes can contribute, at most, about 10-15% of the required current. For the purpose of the charging rate simulation, it was assumed that each ground flash contributes 5 C to the charging process. It was necessary to assume that all electrified clouds contribute to charging by means other than lightning, that the total flash rate can serve as an indirect indicator of the rate of charge transfer, and that oceanic electrified clouds contribute to charging even though they are relatively inefficient in producing lightning. It was also found necessary to add a diurnally invariant charging current component. By trial and error it was found that charging rate diurnal variation curves in Universal time (UT) could be produced with amplitude ratios and general shapes similar to those of the potential gradient diurnal variation curves measured over ocean and arctic regions during voyages of the Carnegie Institute research vessels.

Description: Estimates of the effect of pulse stretching on satellite laser altimetry in particular the Geoscience Laser Altimeter System (GLAS), by cloud multiple scattering were made from an analytical method and from Monte Carlo calculations. The path delay of the return pulse was found to be largest for low-level clouds with particle radii (3-20 microns). The magnitude of the path delay was affected by several factors including cloud height, cloud optical depth, cloud particle size, particle shape, and receiver field of view. Polar aerosol and Rayleigh scattering usually added less than 1 cm to the overall path delay. Path delay estimates for all cloud conditions would be less if a simple Gaussian fit of the return pulse peak were used to measure the pulse's centroid. However, care must be taken in determining the centroid as factors such as pulse width, surface slope and the fitting method used will affect the estimate. A planned application for laser altimetry is high precision monitoring of the height change of polar ice sheets. In the absence of a correction algorithm, the required GLAS altimetry accuracies will not be achieved unless atmospheric channel information is used to remove profiles that are likely to be heavily contaminated by multiple scattering.

Description: The zonal mean eddy heat flux is directly proportional to the wave activity that propagates from the troposphere into the stratosphere. This quantity is a simple eddy diagnostic which is easily calculated from conventional meteorological analyses. Because this "wave driving" of the stratosphere has a strong impact on the stratospheric temperature, it is necessary to compare the impact of the flux with respect to stratospheric radiative changes caused by greenhouse gas changes. Hence, we must understand the precision and accuracy of the heat flux derived from our global meteorological analyses. Herein, we quantify the stratospheric heat flux using five different meteorological analyses, and show that there are 30% differences between these analyses during the disturbed conditions of the northern hemisphere winter. Such large differences result from the planetary differences in the stationary temperature and meridional wind fields. In contrast, planetary transient waves show excellent agreement amongst these five analyses, and this transient heat flux appears to have a long term downward trend.

Description: We examine concurrent upper tropospheric measurements of CN (diameter greater than 4 nm). NO, and NO(Y) during the SONEX Experiment over the North Atlantic (Oct.-Nov., 1997). Elevated CN and NO(Y) concentrations observed in the upper troposphere are attributed largely to enhancements in convective outflows. We estimate that less than 7% of observed high-CN plumes (greater than 10000 /cc) may be attributed to aircraft emissions. Dilution of high-CN convective and aircraft plumes appears to be much more rapid than losses of NO(X) and CN by oxidation and coagulation, respectively, and accounts for much of observed CN concentrations. When taking into account of different time scales against dilution for observable aircraft and convective high-CN plumes (estimated to be 1:4), the contribution by aircraft emissions to CN concentrations is significant, about 20% of the convective source. We find no evidence that particle formation in convective plumes is limited by OH oxidation of SO2.

Description: This study examines a unique data set returned by IMP8 and Geotail on January 29, 1995 during a substorm which resulted in the ejection of a plasmoid. The two spacecraft (s/c) were situated in the north lobe of the tail and both observed a traveling compression region (TCR). From single s/c observations only the length of the plasmoid in X and an estimate of its height in Z can be determined. However, we show that dual s/c measurements of TCRs can be used to model all three dimensions of the underlying plasmoid and to estimate of its rate of expansion or contraction. For this event plasmoid dimensions of Delta(X) approximates 18, Delta(Y) approximates 30, and Delta(Z) approximates 10 R(sub e) are inferred from the IMP8 and Geotail lobe magnetic field measurements. The earthward end of the plasmoid was inferred to be near the mean location of the near-earth neutral line, X approximates -26 R(sub e). Its center was underneath IMP 8 at X approximates -34 R(sub e) and its tailward end appeared to be near X approximates -44 R(sub e). Furthermore, a factor of approximately 2 increase in the amplitude of the TCR occurred in the 1.5 min it took to move from IMP 8 to Geotail. Modeled using conservation of the magnetic flux, this increase in lobe compression implies that the underlying plasmoid was expanding at a rate of approximately 140 km/s. Such an expansion is comparable to recently reported V(sub y) speeds in "young" plasmoids in this region of the tail. Finally, the Geotail measurements indicate that a reconfiguration of the lobe magnetic field closely followed the ejection of the plasmoid which moved magnetic flux tubes into the wake behind the plasmoid where they would convect into the near-earth neutral line and reconnect.

Description: Prior to 1991, major warmings (defined by increasing zonal mean temperatures and zonal mean easterly winds from 60degN to the pole at 10 hPa) typically occurred approximately once every two Arctic winters; a major warming in mid-Dec. 1998 was the first since Feb. 1991. The Dec. 1998 warming was also the second earliest on record. The earliest, and the only other major warming on record before the end of Dec. was in early Dec 1987; prior to that, the earliest was in late Dec./early Jan. 1984-85. The 1984-85 and 1987 warmings resulted in the warmest and weakest lower stratospheric polar vortices in the 20 years before 1998-99. Fig. 1 compares temperatures and vortex strength in 1998-99 with those in the previous 20 years, using the US National Center for Environmental Prediction (NCEP) record; 1987-88 and 1984-85 are also highlighted. The Dec. 1998 warming had a more pronounced effect on mid-stratospheric temperatures than the Dec. 1987 warming (Fig. 1a), although smaller than that of warmings later in winter (e.g., 1984-85). 10-hPa temperatures fell well below average again in late Jan. 1999 and remained unusually low until an early final warming began in late Feb. 840 K PV gradients (Fig. 1c) set a record minimum in Jan. 1999, but were near average in Feb before the final warming. The effect of the Dec. 1998 warming on lower stratospheric temperatures was comparable to that of other major warmings; there was a brief period of record-high minimum 46-hPa temperatures in early Jan 1999 (Fig. 1b), and temperatures then fell to near average for a short period in mid-Feb. Lower stratospheric PV gradients were the weakest on record during the 1998-99 winter (Fig. 1d). The evolution of the vortex and minimum temperatures during 1998-99 was remarkably similar to that during 1987-88, the only previous year when a major warming was observed before the end of Dec.

Description: Two different methods for retrieving Upper Tropospheric Humidities (UTH) from the TOVS (TIROS Operational Vertical Sounder) instruments aboard NOAA polar orbiting satellites are presented and compared. The first one, from the Environmental Technology Laboratory, computed by J. Bates and D. Jackson (hereafter BJ method), estimates UTH from a simplified radiative transfer analysis of the upper tropospheric infrared water vapor channel at wavelength measured by HIRS (6.3 micrometer). The second one results from a neural network analysis of the TOVS (HIRS and MSU) data developed at, the Laboratoire de Meteorologie Dynamique (hereafter the 3I (Improved Initialization Inversion) method). Although the two methods give very similar retrievals in temperate regions (30-60 N and S), an absolute bias up to 16% appears in the convective zone of the tropics. The two datasets have also been compared with UTH retrievals from infrared radiance measurements in the 6.3 micrometer channel from the geostationary satellite METEOSAT (hereafter MET method). The METEOSAT retrievals are systematically drier than the TOVS-based results by an absolute bias between 5 and 25%. Despite the biases, the spatial and temporal correlations are very good. The purpose of this study is to explain the deviations observed between the three datasets. The sensitivity of UTH to air temperature and humidity profiles is analysed as are the clouds effects. Overall, the comparison of the three retrievals gives an assessment of the current uncertainties in water vapor amounts in the upper troposphere as determined from NOAA and METEOSAT satellites.

Description: Carbonados are porous polycrystalline (with crystal sizes up to 100 micrometer) diamonds. Carbonado is found only in alluvial deposits in Bahia, Brazil and in the Central African Republic (CAR). Alluvial deposit host is 1.5 Ga while the carbonados are between 2.6 - 3.8 Ga. The process of fusing the carbonado microcrystals together is not fully understood, partly due to fact that the origin of these carbonado, is not known. Several modes of origins are proposed for carbonado. First, a crustal origin. Carbonados have a light carbon and helium isotopic signature. They contain an enrichment of the rare-earth elements (REE). Carbonados have tightly trapped atmospheric noble gases and contain an evidence of high He content despite the carbonado expected depletion of He at mantle temperatures. Carbonados have high porosity incompatible with high pressure mantle conditions. Second, a mantle origin is proposed based on similar REE pattern to kimberlites. The presence of nitrogen platelet (by IR spectra) indicates high temperature origin and syngenetic inclusions of rutile, ilmenite, and magnetite indicates high pressure and high temperature conditions consistent with mantle origin as well. Third, it is proposed that carbonado diamonds are a result of early impacts into crustal rocks. This is indicated by the rare and controversial occurrence of high pressure diamond polymorph, londsdaleite, frequently found in diamonds from meteorite impact sites, and by observation of planar deformation features, possibly indicating shock events. Finally, it is suggested that carbonados have an extraterrestrial origin, as indicated by a long term annealing based on observation of a zero-phonon line, attributed to paired nitrogen atoms in association with a vacancy.

Description: Vertical electron-density (N (sub e)) profiles, deduced from newly-available ISIS-II digital ionospheric topside-sounder data, are used to investigate the "polar-hole" region within the winter, nighttime polar cap ionosphere during solar minimum. The hole region is located around 0200 MLT near the poleward side of the auroral oval. Earlier investigations had revealed very low N (sub e) values in this region (down to 200/cu cm near 300 km). In the present study, such low N, values (approx. 100/cu cm) were only found near the ISIS (International Satellite for Ionospheric Study)-II altitude of 1400 km. The peak ionospheric concentration below the spacecraft remained fairly constant (approx. 10 (exp 5)/cu cm across the hole region but the altitude of the peak dropped dramatically. This peak dropped, surprisingly, to the vicinity of 100 km. These observations suggest that the earlier satellite in situ measurements, interpreted as deep holes in the ionospheric F-region concentration, could have been made during conditions of an extreme decrease in the altitude of the ionospheric N (sub e) peak. The observations, in combination with other data, indicate that the absence of an F-layer peak may be a frequent occurrence at high latitudes.

Description: We have examined the sea level height tide records at seven tide gauge sites in the region of southcentral Alaska that were affected by the 1964 Prince William Sound earthquake to determine the history of crustal uplift subsequent to the earthquake. There is considerable variation in the behavior depending on the location of the site relative to the 1964 rupture. At Seward, on the eastern side of the Kenai Peninsula we find a slow uplift that is consistent with elastic strain accumulation while at Seldovia and Nikiski on the western side of the Kenai we find a persistent rapid uplift of about 1 cm/yr that most likely represents a long term transient response to the earthquake, but which cannot be sustained over the expected recurrence interval for a great earthquake of several hundred years. Further to the southwest, at Kodiak, we find evidence that the rate of uplift, which is still several mm/yr, has slowed significantly over the past three and a half decades. To the east of the Kenai Peninsula we find subsidence at Cordova and an uncertain behavior at Valdez. At both of these sites there is a mathematically significant time-dependence to the uplift behavior, but the data confirming this time dependence are not as convincing as at Kodiak. At Anchorage, to the north there is little evidence of vertical motion since the earthquake. We compare these long term tide gauge records to recent GPS observations. In general there is reasonable consistency except at Anchorage and Cordova where the GPS measurement indicate somewhat more rapid uplift and subsidence, respectively.

Description: The geographic limits of the South Atlantic Anomaly, as defined by radiation damage, are compared to contours of geomagnetic total field intensity, as defined by the 1995 IGRF, for the present and recent past. The most likely secular variation of the geomagnetic field is estimated and used to extrapolate the geomagnetic field to the year 2100. This indicates that radiation damage to spacecraft and humans in space will likely increase and to cover a much larger geographic area.

Description: We have adopted the transport scenarios used in Part 1 to examine the sensitivity of stratospheric aircraft perturbations to transport changes in our 2-D model. Changes to the strength of the residual circulation in the upper troposphere and stratosphere and changes to the lower stratospheric K(sub zz) had similar effects in that increasing the transport rates decreased the overall stratospheric residence time and reduced the magnitude of the negative perturbation response in total ozone. Increasing the stratospheric K(sub yy) increased the residence time and enhanced the global scale negative total ozone response. However, increasing K(sub yy) along with self-consistent increases in the corresponding planetary wave drive, which leads to a stronger residual circulation, more than compensates for the K(sub yy)-effect, and results in a significantly weaker perturbation response, relative to the base case, throughout the stratosphere. We found a relatively minor model perturbation response sensitivity to the magnitude of K(sub yy) in the tropical stratosphere, and only a very small sensitivity to the magnitude of the horizontal mixing across the tropopause and to the strength of the mesospheric gravity wave drag and diffusion. These transport simulations also revealed a generally strong correlation between passive NO(sub y) accumulation and age of air throughout the stratosphere, such that faster transport rates resulted in a younger mean age and a smaller NO(y) mass accumulation. However, specific variations in K(sub yy) and mesospheric gravity wave strength exhibited very little NO(sub y)-age correlation in the lower stratosphere, similar to 3-D model simulations performed in the recent NASA "Models and Measurements" II analysis. The base model transport, which gives the most favorable overall comparison with inert tracer observations, simulated a global/annual mean total ozone response of -0.59%, with only a slightly larger response in the northern compared to the southern hemisphere. For transport scenarios which gave tracer simulations within some agreement with measurements, the annual/globally averaged total ozone response ranged from -0.45% to -0.70%. Our previous 1995 model exhibited overly fast transport rates, resulting in a global/annually averaged perturbation total ozone response of -0.25%, which is significantly weaker compared to the 1999 model. This illustrates how transport deficiencies can bias model simulations of stratospheric aircraft.

Description: According to the slip partitioning concept, the trench parallel component of relative plate motion in regions of oblique convergence is accommodated by strike-slip faulting in the overriding continental lithosphere. The pattern of postseismic surface deformation due to viscoelastic flow in the lower crust and asthenosphere following a major earthquake on such a fault is modified from that predicted from the conventual elastic layer over viscoelastic halfspace model by the presence of the subducting slab. The predicted effects, such as a partial suppression of the postseismic velocities by 1 cm/yr or more immediately following a moderate to great earthquake, are potentially detectable using contemporary geodetic techniques.

Description: The subject of neotectonics, covering the structures and structural activity of the last 5 million years (i.e., post-Miocene) is a well-recognized field, including "active tectonics," focussed on the last 500,000 years in a 1986 National Research Council report of that title. However, there is a cartographic gap between tectonic maps, generally showing all features regardless of age, and maps of current seismic or volcanic activity. We have compiled a map intended to bridge this gap, using modern data bases and computer-aided cartographic techniques. The maps presented here are conceptually descended from an earlier map showing tectonic and volcanic activity of the last one million years. Drawn by hand with the National Geographic Society's 1975 "The Physical World" map as a base, the 1981 map in various revisions has been widely reproduced in textbooks and various technical publications. However, two decades of progress call for a completely new map that can take advantage of new knowledge and cartographic techniques. The digital tectonic activity map (DTM), presented in shaded relief (Fig. 1) and schematic (Fig. 2) versions, is the result. The DTM is intended to show tectonism and volcanism of the last one million years, a period long enough to be representative of global activity, but short enough that features such as fault scarps and volcanos are still geomorphically recognizable. Data Sources and Cartographic Methods The DTM is based on a wide range of sources, summarized in Table 1. The most important is the digital elevation model, used to construct a shaded relief map. The bathymetry is largely from satellite altimetry, specifically the marine gravity compilations by Smith and Sandwell (1996). The shaded relief map was designed to match the new National Geographic Society world physical map (1992), although drawn independently, from the digital elevation model. The Robinson Projection is used instead of the earlier Van der Grinten one. Although neither conformal nor equal-area, the Robinson Projection provides a reasonable compromise and retains useful detail at high latitudes.

Description: There are procedures and methods for verification of coding algebra and for validations of models and calculations that are in use in the aerospace computational fluid dynamics (CFD) community. These methods would be efficacious if used by the glacier dynamics modelling community. This paper is a presentation of some of those methods, and how they might be applied to uncertainty management supporting code verification and model validation for glacier dynamics. The similarities and differences between their use in CFD analysis and the proposed application of these methods to glacier modelling are discussed. After establishing sources of uncertainty and methods for code verification, the paper looks at a representative sampling of verification and validation efforts that are underway in the glacier modelling community, and establishes a context for these within overall solution quality assessment. Finally, an information architecture and interactive interface is introduced and advocated. This Integrated Cryospheric Exploration (ICE) Environment is proposed for exploring and managing sources of uncertainty in glacier modelling codes and methods, and for supporting scientific numerical exploration and verification. The details and functionality of this Environment are described based on modifications of a system already developed for CFD modelling and analysis.

Description: The behavior of mesospheric odd nitrogen species during and following relativistic and diffuse auroral precipitation events is simulated, Below 75 km nitric oxide is enhanced in proportion to the ion pair production function associated with the electron precipitation and the length of the event. Nitrogen dioxide and nitric acid are also enhanced. At 65 km the percentage of odd nitrogen for N is 0.1%, HNO3 is 1.6%, NO2 is 15%, and NO is 83.3%. Between 75 and 85 km NO is depleted during particle events due to the faster destruction of NO by N relative to the production of NO by N reacting with O2. Recovery of NO depends on transport from the lower thermosphere, where NO is produced in abundant amounts during particle events.

Description: Long-lived tropospheric source gases, such as nitrous oxide, enter the stratosphere through the tropical tropopause, are transported throughout the stratosphere by the Brewer-Dobson circulation, and are photochemically destroyed in the upper stratosphere. These chemical constituents, or "tracers" can be used to track mixing and transport by the stratospheric winds. Much of our understanding about the stratospheric circulation is based on large scale gradients and other spatial features in tracer fields constructed from satellite measurements. The point of view presented in this paper is different, but complementary, in that transport is described in terms of tracer probability distribution functions (PDFs). The PDF is computed from the measurements, and is proportional to the area occupied by tracer values in a given range. The flavor of this paper is tutorial, and the ideas are illustrated with several examples of transport-related phenomena, annotated with remarks that summarize the main point or suggest new directions. One example shows how the multimodal shape of the PDF gives information about the different branches of the circulation. Another example shows how the statistics of fluctuations from the most probable tracer value give insight into mixing between different regions of the atmosphere. Also included is an analysis of the time-dependence of the PDF during the onset and decline of the winter circulation, and a study of how "bursts" in the circulation are reflected in transient periods of rapid evolution of the PDF. The dependence of the statistics on location and time are also shown to be important for practical problems related to statistical robustness and satellite sampling. The examples illustrate how physically-based statistical analysis can shed some light on aspects of stratospheric transport that may not be obvious or quantifiable with other types of analyses. An important motivation for the work presented here is the need for synthesis of the large and growing database of observations of the atmosphere and the vast quantities of output generated by atmospheric models.

Description: Magnetic remanence of crustal rocks can reside in three common rock-forming magnetic minerals: magnetite, pyrrhotite, and hematite. Thermoremanent magnetization (TRM) of magnetite and pyrrhotite is carried mostly by single domain (SD) grains. The TRM of hematite grains, however, is carried mostly by multidomain (NM) grains. This characteristic is illustrated by TRM acquisition curves for hematite of variable grainsizes. The transition between truly NM behavior and tendency towards SD behavior his been established between hematite grainsizes of 0. 1 and 0.05 mm. Coarse grainsize of lower crustal rocks and the large sensitivity of MD hematite grains to acquire TRM indicates that hematite could be a significant contributor to long-wavelength magnetic anomalies.

Description: We present an extension for the 2D (zonal mean) version of our Numerical Spectral Mode (NSM) that incorporates Hines' Doppler spread parameterization (DSP) for small scale gravity waves (GW). This model is applied to describe the seasonal variations and the semi-annual and quasi-biennial oscillations (SAO and QBO). Our earlier model reproduced the salient features of the mean zonal circulation in the middle atmosphere, including the QBO extension into the upper mesosphere inferred from UARS measurements. In the present model we incorporate also tropospheric heating to reproduce the upwelling at equatorial latitudes associated with the Brewer-Dobson circulation that affects significantly the dynamics of the stratosphere as Dunkerton had pointed out. Upward vertical winds increase the period of the QBO observed from the ground. To compensate for that, one needs to increase the eddy diffusivity and the GW momentum flux, bringing the latter closer to values recommended in the DSP. The QBO period in the model is 30 months (mo), which is conducive to synchronize this oscillation with the seasonal cycle of solar forcing. Multi-year interannual oscillations are generated through wave filtering by the solar driven annual oscillation in the zonal circulation. Quadratic non-linearities generate interseasonal variations to produce a complicated pattern of variability associated with the QBO. The computed temperature amplitudes for the SAO and QBO are in substantial agreement with observations at equatorial and extratropical latitudes. At high latitudes, however, the observed QBO amplitudes are significantly larger, which may be a signature of propagating planetary waves not included in the present model. The assumption of hydrostatic equilibrium not being imposed, we find that the effects from the vertical Coriolis force associated with the equatorial oscillations are large for the vertical winds and significant for the temperature variations even outside the tropics but are relatively small for the zonal winds.

Description: This paper reviews observations of stratospheric temperatures that have been made over a period of several decades. Those observed temperatures have been used to assess variations and trends in stratospheric temperatures. A wide range of observation datasets have been used, comprising measurements by radiosonde (1940s to the present), satellite (1979 - present), lidar (1979 - present) and rocketsonde (periods varying with location, but most terminating by about the mid-1990s). In addition, trends have also been assessed from meteorological analyses, based on radiosonde and/or satellite data, and products based on assimilating observations into a general circulation model. Radiosonde and satellite data indicate a cooling trend of the annual-mean lower stratosphere since about 1980. Over the period 1979-1994, the trend is 0.6K/decade. For the period prior to 1980, the radiosonde data exhibit a substantially weaker long-term cooling trend. In the northern hemisphere, the cooling trend is about 0.75K/decade in the lower stratosphere, with a reduction in the cooling in mid-stratosphere (near 35 km), and increased cooling in the upper stratosphere (approximately 2 K per decade at 50 km). Model simulations indicate that the depletion of lower stratospheric ozone is the dominant factor in the observed lower stratospheric cooling. In the middle and upper stratosphere both the well-mixed greenhouse gases (such as CO) and ozone changes contribute in an important manner to the cooling.

Description: 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.

Description: Tropospheric column ozone (TCO) and stratospheric column ozone (SCO) gridded data in the tropics for 1979-present are now available from NASA Goddard Space Flight Center via either direct ftp, world-NN,ide-NN,eb, or electronic mail. This note provides a brief overview of the method used to derive the data set including validation and adjustments.

Description: In this study, we examine the sensitivity of long lived tracers to changes in the base transport components in our 2-D model. Changes to the strength of the residual circulation in the upper troposphere and stratosphere and changes to the lower stratospheric K(sub zz) had similar effects in that increasing the transport rates decreased the overall stratospheric mean age, and increased the rate of removal of material from the stratosphere. Increasing the stratospheric K(sub yy) increased the mean age due to the greater recycling of air parcels through the middle atmosphere, via the residual circulation, before returning to the troposphere. However, increasing K(sub yy) along with self-consistent increases in the corresponding planetary wave drive, which leads to a stronger residual circulation, more than compensates for the K(sub yy)-effect, and produces significantly younger ages throughout the stratosphere. Simulations with very small tropical stratospheric K(sub yy) decreased the globally averaged age of air by as much as 25% in the middle and upper stratosphere, and resulted in substantially weaker vertical age gradients above 20 km in the extratropics. We found only very small stratospheric tracer sensitivity to the magnitude of the horizontal mixing across the tropopause, and to the strength of the mesospheric gravity wave drag and diffusion used in the model. We also investigated the transport influence on chemically active tracers and found a strong age-tracer correlation, both in concentration and calculated lifetimes. The base model transport gives the most favorable overall comparison with a variety of inert tracer observations, and provides a significant improvement over our previous 1995 model transport. Moderate changes to the base transport were found to provide modest agreement with some of the measurements. Transport scenarios with residence times ranging from moderately shorter to slightly longer relative to the base case simulated N2O lifetimes that were within the observational estimates of Volk et al. [1997]. However, only scenarios with rather fast transport rates were comparable with the Volk et al. estimates of CFCl3 lifetimes. This is inconsistent with model-measurement comparisons of mean age in which the base model or slightly slower transport rates compared the most favorably with balloon SF6 data. For all comparisons shown, large transport changes away from the base case resulted in simulations that were outside the range of measurements, and in many cases, far outside this range.

Description: An intensive field campaign involving measurement of various aerosol physical, chemical, and radiative properties was conducted at Sde Boker in the Negev Desert of Israel, from 18 February to 15 March 1997. Nephelometer measurements gave average background scattering coefficient values of about 25 M/m at 550 nm wavelength, but strong dust events caused the value of this parameter to rise up to about 800 M/m Backscattering fractions did not depend on aerosol loading, and generally fell in the range of 0.1 to 0.25, comparable to values reported for marine and Arctic environments. Chemical analysis of the aerosol revealed that, in the coarse size range (2 - 10 micrometer equivalent aerodynamic diameter (EAD)), calcium (Ca) was by far the most abundant element followed by silicon (Si), both of which are indicators for mineral dust. In the fine size fraction (〈 2 micrometers EAD), sulfur (S) generally was the dominant element, except during high dust episodes when Ca and Si were again the most abundant. Furthermore, fine black carbon (BC) correlates with S, suggesting that they may have originated from the same sources or source regions. An indication of the short-term effect of aerosol loading on radiative forcing was provided by measurements of global and diffuse solar radiation, which showed that during high turbidity periods (strong dust events) almost all of the solar radiation reaching the area is scattered or absorbed.

Description: This presentation provides a brief overview of the scientific activities of the International GPS Service (IGS). This was an approved activity of the International Association of Geodesy (IAG) with official start of service on 1 Jan 1994. The mission of the IGS is "To provide a service to support geodetic and geophysical research activities, through GPS data and data products." The presentation explains the concept of the IGS working group, and pilot projects, and reviews the current working groups and pilot projects.

Description: A model which simulates the sensitivity of satellite-measured solar radiation to changes in relative humidity was developed by AER under contract to the Office of Naval Research. AER uses the model to support meteorological analysis and evaluate the use of satellites in relation to marine aerosols/optical instrument reduction. During the project, AER utilized RADTMO, a COSMIC set of programs which computes scattered atmospheric radiation caused by aerosols.

Description: Ground-based solar absorption spectra were measured from Fairbanks, Alaska (65 degrees North, 148 degrees West) and from Ny Alesund, Spitsbergen, (79 degrees North, 12 degrees East) from March to September 1997 by Fourier Transform Infra-Red (FTIR) spectrometers.

Description: Remote sensing of chlorine monoxide (CIO) by the Microwave Limb Sounder (MLS) experiment aboard the Upper Atmosphere Research Satellite (UARS) has provided global measurements of variations in stratospheric free chlorine (for 1991 to 1997).

Description: There has been a question on whether it is possible to measure the surface wind speeds of tropical cyclones with a spaceborne Ku-band scatterometer.

Description: We examine energetic particle transport at or near Corotating Interaction Regions (CIRs).

Description: In global geodetic solutions vertical rates of site motion are usually estimated relative to the geocenter (center of figure) of the solid earth. The velocity of the geocenter is estimated assuming that the plates are rigid, that the velocities of the plates equal those in NUVEL-1A (DeMets et al. 1990, 1994) and that the uplift, subsidence, and intraplate deformation due to glacial isostatic adjustment is negligible.

Description: With aircraft-mounted in-situ and remote sensing instruments for dynamical, thermal, and chemical measurements, we studied two cases of tropopause folding.

Description: An accurate method of surveying the soil was developed by NASA and the Department of Agriculture. The method involves using ground penetrating radar to produce subsurface graphs. By examining printouts from the system's recorder, scientists can determine whether a site is appropriate for building, etc.

Description: A redesigned version of a soil/rock analyzer developed by Martin Marietta under a Langley Research Center contract is being marketed by Aurora Tech, Inc. Known as the Aurora ATX-100, it has self-contained power, an oscilloscope, a liquid crystal readout, and a multichannel spectrum analyzer. It measures energy emissions to determine what elements in what percentages a sample contains. It is lightweight and may be used for mineral exploration, pollution monitoring, etc.

Description: Target areas for sinking base holes, underground pipelines, etc., can be identified with the assistance of NASA Ames developed technology, by Airborne Pipeline Services, Inc. Subsurface features are computer processed; the system can cover 250 miles a day and was first developed by Applied Science, Inc.

Description: Current passive-microwave rain-retrieval methods are largely based on databases built off-line using cloud models.

Description: Geodetic measurements using the Global Positioning System and other techniques show north-south shortening near Los Angeles to be fastest across the northern part of the metropolitan area, where an ESE-striking, 5- to 40-km-wide belt lying to the south of San Gabriel Mountains and to the north of downtown and West Los Angeles is shortening at 5 mm/yr.

Description: Data from terrestrial GPS receivers are being used in growing numbers of applications requiring precise tropospheric sensing. One emerging application is the calibration of water vapor measurements from spaceborne microwave radiometers.

Description: Trace gases measured by the Atmospheric Trace Molecule Spectroscopy (ATMOS) instrument during the Mar/Apr 1992(AT-1), Apr 1993(AT-2), and Nov 1994(AT-3) space-shuttle missions have been mapped into equivalent latitude/potential temperature (EqL/0) coordinates.

Description: The late Holocene glacial moraine chronology in the southernmost South American Andes includes four (sup 14)C dated Neoglacial advances and retreats.

Description: Airborne measurements of NO, total reactive nitrogen(NO(sub y)), O(sub 3), and condensation nuclei (CN) were made, together with other tracers such as CO and H(sub 2)O, within air traffic corridors over the U.S. and North Atlantic regions in October and November 1997.

Description: We use WIND solar wind data and POLAR UV imaging data to study magnetospheric responses and substorm triggering mechanisms during and after interplanetary (IP) shock events.

Description: The M7.1 Hector Mine earthquake ruptured the Lavic Lake fault near Twentynine Palms, CA at 09:46 UTC October 16, 1999.

Description: During the 1998-1999 Antarctic field season, we installed three autonomous GPS stations in Marie Byrd Land, West Antarctica to measure glacio-isostatic rebound and rates of spreading across the West Antartic Rift System.

Description: Because of questions regarding the validity of proposed accreted terranes, the tectonostratigraphic evolution of the southern margin of North America plate in Mexico remains controversial.

Description: Sea surface height anomaly signals from satellite altimeter data are used to estimate heat storage. Since variability in sea surface height is mostly due to expansion and contraction of the water column it can be correlated with variations in the heat and salt content.

Description: We present box model calculations of ozone loss rates corresponding to the results of the Match experiment 1991/92.

Description: Many long-lived trace gas profiles observed by the Atmospheric Trace Molecule Spectroscopy (ATMOS) instrument around the developing polar vortex (the proto-vortex) during early November 1994 show distinctive minimum/maximum pairs (laminae).

Description: Dust is but one of many aerosols that are analyzed at the Jet Propulsion Laboratory in Pasadena. The purpose of this paper is to describe the process in analyzing and digitizing dust within a source region to better explain the work achieved by my internship. This paper will go over how to view collected data by Multi-angle Imaging SpectroRadiometer (MISR) [1] and the procedure of downloading data to be analyzed. With this data, one can digitize dust plumes using two methods called plume lines and plume polygons with the help of the software MISR INteractive eXplorer (MINX)[3]; thus, the theory of MINX's[3] algorithm and these methods are discussed in detail. Research was gathered from these techniques and emphasis is also focused on the obtained data and results.

Description: No abstract available

Description: Thin cirrus clouds (optical depth (OD) 〈 03) are often undetected by standard cloud masking in satellite aerosol retrieval algorithms. However, the Mu]tiangle Imaging Spectroradiometer (MISR) aerosol retrieval has the potential to discriminate between the scattering phase functions of cirrus and aerosols, thus separating these components. Theoretical tests show that MISR is sensitive to cirrus OD within Max{0.05 1 20%l, similar to MISR's sensitivity to aerosol OD, and MISR can distinguish between small and large crystals, even at low latitudes, where the range of scattering angles observed by MISR is smallest. Including just two cirrus components in the aerosol retrieval algorithm would capture typical MISR sensitivity to the natural range of cinus properties; in situations where cirrus is present but the retrieval comparison space lacks these components, the retrieval tends to underestimate OD. Generally, MISR can also distinguish between cirrus and common aerosol types when the proper cirrus and aerosol optical models are included in the retrieval comparison space and total column OD is 〉-0.2. However, in some cases, especially at low latitudes, cirrus can be mistaken for some combinations of dust and large nonabsorbing spherical aerosols, raising a caution about retrievals in dusty marine regions when cirrus is present. Comparisons of MISR with lidar and Aerosol Robotic Network show good agreement in a majority of the cases, but situations where cirrus clouds have optical depths 〉0.15 and are horizontally inhomogeneous on spatial scales shorter than 50 km pose difficulties for cirrus retrieval using the MISR standard aerosol algorithm..

Description: The 2010 Arctic Mars Analog Svalbard Expedition (AMASE) investigated two geologic settings using methodologies and techniques being developed or considered for future Mars missions, such as the Mars Science Laboratory (MSL), ExoMars, and Mars Sample Return. The Sample Analysis at Mars (SAM) [1] instrument suite, which will be on MSL, consists of a quadrupole mass spectrometer (QMS), a gas chromatograph (GC), and a tunable laser mass spectrometer (TLS); all will be applied to analyze gases created by pyrolysis of samples. During AMASE, a Hiden Evolved Gas Analysis-Mass Spectrometer (EGA-MS) system represented the EGA-MS capability of SAM. Another MSL instrument, CheMin, will use x-ray diffraction (XRD) and x-ray fluorescence (XRF) to perform quantitative mineralogical characterization of samples [e.g., 2]. Field-portable versions of CheMin were used during AMASE. AMASE 2010 focused on two sites that represented biotic and abiotic analogs. The abiotic site was the basaltic Sigurdfjell vent complex, which contains Mars-analog carbonate cements including carbonate globules which are excellent analogs for the globules in the ALH84001 martian meteorite [e.g., 3, 4]. The biotic site was the Knorringfjell fossil methane seep, which featured carbonates precipitated in a methane-supported chemosynthetic community [5]. This contribution focuses on EGA-MS analyses of samples from each site, with mineralogy comparisons to CheMin team results. The results give insight into organic content and organic-mineral associations, as well as some constraints on the minerals present.

Description: Due to gravitation the main reflector of the Effelsberg 100-m telescope of the Max Planck Institute for Radio Astronomy is deformed whenever it is tilted from zenith to arbitrary elevation angles. However, the resulting shape always is a paraboloid again, though with different parameters, a phenomenon which is called homologous deformation. In summer 2008, we have carried out measurements with a total station to determine the magnitude of these deformations in order to evaluate existing assumptions provided by the manufacturer from the telescope's design phase. The measurements are based on a newly developed approach with a Leica TCRP 1201 total station mounted head down near the subreflector. Mini-retro-reflectors are placed at various locations on the paraboloid itself and on the subreflector support structure. The results indicate that the measurement setup is suitable for the purpose and provides the information needed for a determination of elevation dependent delay corrections. The focal length changes only by about 8 mm when the telescope is tilted from 90. to 7.5. elevation angle.

Description: We describe the experience gained in a project to continuously monitor the local tie at the Geodetic Observatory Ny-Alesund. A PC-controlled robotic total station was used to monitor survey prisms that were attached to survey pillars of the local network and the monuments used for geodetic VLBI and GNSS measurements. The monitoring lasted for seven days and had a temporal resolution of six minutes. The raw angle and distance measurements show clear sinusoidal signatures with a daily period, most strongly for a four-day period with 24 hours of sunshine. The derived topocentric coordinates of the survey prisms attached to the GNSS monument and the VLBI radio telescope act as approximation for the local tie. We detect clear signatures at the mm-level. With the current approach we cannot distinguish between real motion of the prisms and potential thermal influences on the instrument used for the observations. However, the project shows that continuous local tie monitoring is feasible today and in the future can and should be used for all geodetic co-location stations.

Description: Temporal statistical position stability is required for VLBI sources to define a stable Celestial Reference Frame (CRF) and has been studied in many recent papers. This study analyzes the sources from the latest realization of the International Celestial Reference Frame (ICRF2) with the Allan variance, in addition to taking into account the apparent linear motions of the sources. Focusing on the 295 defining sources shows how they are a good compromise of different criteria, such as statistical stability and sky distribution, as well as having a sufficient number of sources, despite the fact that the most stable sources of the entire ICRF2 are mostly in the Northern Hemisphere. Nevertheless, the selection of a stable set is not unique: studying different solutions (GSF005a and AUG24 from GSFC and OPA from the Paris Observatory) over different time periods (1989.5 to 2009.5 and 1999.5 to 2009.5) leads to selections that can differ in up to 20% of the sources. Observing, recording, and network improvement are some of the causes, showing better stability for the CRF over the last decade than the last twenty years. But this may also be explained by the assumption of stationarity that is not necessarily right for some sources.

Description: GEOSAT is a multi-technique geodetic analysis software developed at Forsvarets Forsknings Institutt (Norwegian defense research establishment). The Norwegian Mapping Authority has now installed the software and has, together with Forsvarets Forsknings Institutt, adapted the software to deliver datum-free normal equation systems in SINEX format. The goal is to be accepted as an IVS Associate Analysis Center and to provide contributions to the IVS EOP combination on a routine basis. GEOSAT is based on an upper diagonal factorized Kalman filter which allows estimation of time variable parameters like the troposphere and clocks as stochastic parameters. The tropospheric delays in various directions are mapped to tropospheric zenith delay using ray-tracing. Meteorological data from ECMWF with a resolution of six hours is used to perform the ray-tracing which depends both on elevation and azimuth. Other models are following the IERS and IVS conventions. The Norwegian Mapping Authority has submitted test SINEX files produced with GEOSAT to IVS. The results have been compared with the existing IVS combined products. In this paper the outcome of these comparisons is presented.

Description: To achieve the goals of the VLBI2010 project and the Global Geodetic Observing System (GGOS), an automated monitoring of the reference points of the various geodetic space techniques, including Very Long Baseline Interferometry (VLBI), is desirable. The resulting permanent monitoring of the local-tie vectors at co-location stations is essential to obtain the sub-millimeter level in the combinations. For this reason a monitoring system was installed at the Geodetic Observatory Wettzell by the Geodetic Institute of the University of Karlsruhe (GIK) to observe the 20m VLBI radio telescope from May to August 2009. A specially developed software from GIK collected data from automated total station measurements, meteorological sensors, and sensors in the telescope monument (e.g., Invar cable data). A real-time visualization directly offered a live view of the measurements during the regular observation operations. Additional scintillometer measurements allowed refraction corrections during the post-processing. This project is one of the first feasibility studies aimed at determining significant deformations of the VLBI antenna due to, for instance, changes in temperature.

Description: Since 2008 the VLBI group at the Institute of Geodesy and Geophysics at TU Vienna has focused on the development of a new VLBI data analysis software called VieVS (Vienna VLBI Software). One part of the program, currently under development, is a unit for parameter estimation in so-called global solutions, where the connection of the single sessions is done by stacking at the normal equation level. We can determine time independent geodynamical parameters such as Love and Shida numbers of the solid Earth tides. Apart from the estimation of the constant nominal values of Love and Shida numbers for the second degree of the tidal potential, it is possible to determine frequency dependent values in the diurnal band together with the resonance frequency of Free Core Nutation. In this paper we show first results obtained from the 24-hour IVS R1 and R4 sessions.

Description: The ICRF2 catalog was constructed by the IERS/IVS Working Group with oversight by the IAU Working Group. Derived using data from August 1979 through March 2009, it is a great improvement over the original ICRF with 3414 extragalactic radio source positions, a noise floor of 40 microarcsec, and axis stability of 10 microarcsec. Significant refinements were made in the selection of defining sources, modeling, and the integration of CRF, TRF, and EOP. The adoption of the ICRF2 was approved by the IAU in Resolution B3 at the XXVII IAU General Assembly and became effective 1 January 2010.

Description: Three 12-meter radio telescopes are being built in Australia (the AuScope project) and one in New Zealand. These facilities will be fully-equipped for undertaking S and X-band geodetic VLBI observations and correlation will take place on a software correlator (part of the AuScope project). All sites are equipped with permanent GPS receivers to provide co-location of several space geodetic techniques. The following scientific tasks of geodesy and astrometry are considered. 1. Improvement and densification of the International Celestial Reference Frame in the southern hemisphere; 2. Improvement of the International Terrestrial Reference Frame in the region; 3. Measurement of intraplate deformation of the Australian tectonic plate.

Description: VLBI2010 holds out promise for greatly increased precision in measuring geodetic and Earth rotation parameters. As a by-product there will be a wealth of interesting new astronomical data. At the same time, astronomical knowledge may be needed to disentangle the astronomical and geodetic contributions to the measured delays and phases. This presentation explores this astro-geo link.

Description: Within almost all space-geodetic techniques, contributions of different analysis centers (ACs) are combined in order to improve the robustness of the final product. So far, the contributing series are assumed to be independent as each AC processes the observations in different ways. However, the series cannot be completely independent as each analyst uses the same set of original observations and many applied models are subject to conventions used by each AC. In this paper, it is shown that neglecting correlations between the contributing series yields too optimistic formal errors and small, but insignificant, errors in the estimated parameters derived from the adjustment of the combined solution.

Description: We present preliminary results of proof-of-concept studies for an automatic monitoring system of local site ties. The system is based on the usage of robotic total stations. A set of tests were performed with a Leica TCA2003 total station on the local network of Goddard s Geophysical and Astronomical Observatory (GGAO) and the 5-m VLBI antenna at this site. Both the TCA2003 and the VLBI antenna are controlled from a Matlab-coded control program. Running specific observational programs, data were collected that indicate that the reference point of the VLBI antenna can be automatically determined with an accuracy of 1 mm or better.

Description: In order to establish the position of the center of mass of the Earth in the International Celestial Reference Frame, observations of the Global Positioning Satellite (GPS) constellation using the IVS network are important. With a good frame-tie between the coordinates of the IVS telescopes and nearby GPS receivers, plus a common local oscillator reference signal, it should be possible to observe and record simultaneously signals from the astrometric calibration sources and the GPS satellites. The standard IVS solution would give the atmospheric delay and clock offsets to use in analysis of the GPS data. Correlation of the GPS signals would then give accurate orbital parameters of the satellites in the ICRF reference frame, i.e., relative to the positions of the astrometric sources. This is particularly needed to determine motion of the center of mass of the earth along the rotation axis.

Description: We have been developing a state-of-the-art tool to estimate the atmospheric path delays by raytracing through mesoscale analysis (MANAL) data, which is operationally used for numerical weather prediction by the Japan Meteorological Agency (JMA). The tools, which we have named KAshima RAytracing Tools (KARAT)', are capable of calculating total slant delays and ray-bending angles considering real atmospheric phenomena. The KARAT can estimate atmospheric slant delays by an analytical 2-D ray-propagation model by Thayer and a 3-D Eikonal solver. We compared PPP solutions using KARAT with that using the Global Mapping Function (GMF) and Vienna Mapping Function 1 (VMF1) for GPS sites of the GEONET (GPS Earth Observation Network System) operated by Geographical Survey Institute (GSI). In our comparison 57 stations of GEONET during the year of 2008 were processed. The KARAT solutions are slightly better than the solutions using VMF1 and GMF with linear gradient model for horizontal and height positions. Our results imply that KARAT is a useful tool for an efficient reduction of atmospheric path delays in radio-based space geodetic techniques such as GNSS and VLBI.

Description: The IVS scheduled a special astrometric VLBI session for the International Year of Astronomy 2009 (IYA09) commemorating 400 years of optical astronomy and 40 years of VLBI. The IYA09 session is the most ambitious geodetic session to date in terms of network size, number of sources, and number of observations. We describe the process of designing, coordinating, scheduling, pre-session station checkout, correlating, and analyzing this session.

Description: Recent investigations have shown significant shortcomings in the model which is proposed by the IERS to account for the variations in the Earth s rotation with periods around one day and less. To overcome this, an empirical model can be estimated more or less directly from the observations of space geodetic techniques. The aim of this paper is to evaluate the quality and reliability of such a model based on VLBI observations. Therefore, the impact of the estimation method and the analysis options as well as the temporal stability are investigated. It turned out that, in order to provide a realistic accuracy measure of the model coefficients, the formal errors should be inflated by a factor of three. This coincides with the noise floor and the repeatability of the model coefficients and it captures almost all of the differences that are caused by different estimation techniques. The impact of analysis options is small but significant when changing troposphere parameterization or including harmonic station position variations.

Description: The IVS runs two tropospheric products: The IVS tropospheric parameter rapid combination monitors the zenith wet delay (ZWD) and zenith total delay (ZTD) of the rapid turnaround sessions R1 and R4. Goal of the combination is the identification and the exclusion of outliers by comparison and the assessment of the precision of current VLBI solutions in terms of tropospheric parameters. The rapid combination is done on a weekly basis four weeks after the observation files are released on IVS Data Centers. Since tropospheric and geodetic parameters, such as vertical station components, can significantly correlate, the consistency of the ZTD can be a measure of the consistency of the corresponding TRF as well. The ZWD mainly rely on accurate atmospheric pressure data. Thus, besides estimation techniques, modeling and analyst s noise, ZWD reflects differences in the atmospheric pressure data applied to the VLBI analysis. The second product, called tropospheric parameter long-term combination, aims for an accurate determination of climatological signals, such as trends of the atmospheric water vapor observed by VLBI. Therefore, the long-term homogeneity of atmospheric pressure data plays a crucial role for this product. The paper reviews the methods applied and results achieved so far and describes the new maintenance through DGFI.

Description: The processing of the GPS observations of the site survey at Sheshan 25-m radio telescope in August 2008 is reported. Because each session in this survey is only about six hours, not allowing the subdaily high frequency variations in the station coordinates to be reasonably smoothed, and because there are serious cycle slips in the observations and a large volume of data would be rejected during the software automatic adjustment of slips, the ordinary solution settings of GAMIT needed to be adjusted by loosening the constraints in the a priori coordinates to 10 m, adopting the "quick" mode in the solution iteration, and combining Cview manual operation with GAMIT automatic fixing of cycle slips. The resulting coordinates of the local control polygon in ITRF2005 are then compared with conventional geodetic results. Due to large rotations and translations in the two sets of coordinates (geocentric versus quasi-topocentric), the seven transformation parameters cannot be solved for directly. With various trial solutions it is shown that with a partial pre-removal of the large parameters, high precision transformation parameters can be obtained with post-fit residuals at the millimeter level. This analysis is necessary to prepare the follow-on site and transformation survey of the VLBI and SLR telescopes at Sheshan

Description: While the overall architecture is similar for the station hardware in current S/X systems and in the VLBI2010 systems under development, various functions are implemented differently. Some of these differences, and the reasons behind them, are described here.

Description: A permanent geodetic VLBI station with a 22-m diameter antenna will be newly constructed in Korea by the National Geographic Information Institute (NGII) under the project Korea VLBI system for Geodesy (KVG) that aims at maintaining the Korean geodetic datum accurately on the International Terrestrial Reference Frame (ITRF). KVG can receive 2, 8, 22, and 43 GHz bands simultaneously in order to conduct geodetic and astronomical VLBI observations with Korea astronomical VLBI stations along with geodetic observations with IVS stations. This simultaneous four-band receiving capability is a unique feature of the KVG system. The KVG has started officially in October 2008. A new geodetic VLBI station will be constructed at Sejong city (about 120 km south of Seoul and about 20 km north-northwest of Daejeon) and construction of all systems will be completed in 2011.

Description: Automated monitoring of operational system parameters for the geodetic space techniques is becoming more important in order to improve the geodetic data and to ensure the safety and stability of automatic and remote-controlled observations. Therefore, the Wettzell group has developed the system monitoring software, SysMon, which is based on a reliable, remotely-controllable hardware/software realization. A multi-layered data logging system based on a fanless, robust industrial PC with an internal database system is used to collect data from several external, serial, bus, or PCI-based sensors. The internal communication is realized with Remote Procedure Calls (RPC) and uses generative programming with the interface software generator idl2rpc.pl developed at Wettzell. Each data monitoring stream can be configured individually via configuration files to define the logging rates or analog-digital-conversion parameters. First realizations are currently installed at the new laser ranging system at Wettzell to address safety issues and at the VLBI station O Higgins as a meteorological data logger. The system monitoring concept should be realized for the Wettzell radio telescope in the near future.

Description: Automating and remotely controlling observations are important for future operations in a Global Geodetic Observing System (GGOS). At the Geodetic Observatory Wettzell, in cooperation with the Max-Planck-Institute for Radio Astronomy in Bonn, a software extension to the existing NASA Field System has been developed for remote control. It uses the principle of a remotely accessible, autonomous process cell as a server extension for the Field System. The communication is realized for low transfer rates using Remote Procedure Calls (RPC). It uses generative programming with the interface software generator idl2rpc.pl developed at Wettzell. The user interacts with this system over a modern graphical user interface created with wxWidgets. For security reasons the communication is automatically tunneled through a Secure Shell (SSH) session to the telescope. There are already successful test observations with the telescopes at O Higgins, Concepcion, and Wettzell. At Wettzell the software is already used routinely for weekend observations. Therefore the first public release of the software is now available, which will also be useful for other telescopes.

Description: The next generation geodetic VLBI instrument is being developed with a goal of 1 mm position uncertainty in twenty-four hours. We have implemented a proof-of-concept system for a possible VLBI2010 signal chain, from feed through recorder, on the Westford (Massachusetts, USA) 18-m and MV-3 (Maryland, USA) 5-m antennas. Data have been obtained in four 512 MHz bands spanning the range 3.5 to 11 GHz to investigate the sensitivity and phase delay capability of the system. Using a new phase cal design, the phases have been aligned across four bands spanning 2 GHz with an RMS deviation of approximately eight degrees. Several components of the system will be improved for the prototype version of VLBI2010, including the feed, digital backend, and recorder, and these will be installed on a 12-m antenna that has been purchased and is ready for installation at the Goddard Space Flight Center outside of Washington, D.C., USA, site of the MV-3 antenna.

Description: Astrometric observations of radio source occultations by solar system bodies may be of large interest for testing gravity theories, dynamical astronomy, and planetary physics. In this paper, we present an updated list of the occultations of astrometric radio sources by planets expected in the coming years. Such events, like solar eclipses, generally speaking can only be observed in a limited region. A map of the shadow path is provided for the events that will occurr in regions with several VLBI stations and hence will be the most interesting for radio astronomy experiments.

Description: The purpose of the Russian VLBI "Quasar" Network is to carry out astrometrical and geodynamical investigations. Since 2006 purely domestic observational programs with data processing at the IAA correlator have been carried out. To maintain these geodynamical programs e-VLBI technology is being developed and tested. This paper describes the IAA activity of developing a real-time VLBI system using high-speed digital communication links.

Description: The construction project of Korea Geodetic VLBI officially started in October 2008. The construction of all systems will be completed by the end of 2011. The project was named Korea VLBI system for Geodesy (KVG), and its main purpose is to maintain the Korea Geodetic Datum. In case of the KVG system, an observation room with an H-maser frequency standard is located in a building separated from the antenna by several tens of meters. Therefore KVG system will adopt a so-called round-trip system to transmit reference signals to the antenna with reduction of the effect of path length variations. KVG s round-trip system is designed not only to use either metal or optical fiber cables, but also to measure path length variations directly. We present this unique round trip system for KVG.

Description: Zenith Tropospheric Delay (ZTD) due to water vapor derived from space geodetic techniques and numerical weather prediction simulated-reanalysis data exhibits non-linear and non-stationary properties akin to those in the crucial geophysical signals of interest to the research community. These time series, once decomposed into additive (and stochastic) components, have information about the long term global change (the trend) and other interpretable (quasi-) periodic components such as seasonal cycles and noise. Such stochastic component(s) could be a function that exhibits at most one extremum within a data span or a monotonic function within a certain temporal span. In this contribution, we examine the use of the combined Ensemble Empirical Mode Decomposition (EEMD) and Independent Component Analysis (ICA): the EEMD-ICA algorithm to extract the independent local oscillatory stochastic components in the tropospheric delay derived from the European Centre for Medium-Range Weather Forecasts (ECMWF) over six geodetic sites (HartRAO, Hobart26, Wettzell, Gilcreek, Westford, and Tsukub32). The proposed methodology allows independent geophysical processes to be extracted and assessed. Analysis of the quality index of the Independent Components (ICs) derived for each cluster of local oscillatory components (also called the Intrinsic Mode Functions (IMFs)) for all the geodetic stations considered in the study demonstrate that they are strongly site dependent. Such strong dependency seems to suggest that the localized geophysical signals embedded in the ZTD over the geodetic sites are not correlated. Further, from the viewpoint of non-linear dynamical systems, four geophysical signals the Quasi-Biennial Oscillation (QBO) index derived from the NCEP/NCAR reanalysis, the Southern Oscillation Index (SOI) anomaly from NCEP, the SIDC monthly Sun Spot Number (SSN), and the Length of Day (LoD) are linked to the extracted signal components from ZTD. Results from the synchronization analysis show that ZTD and the geophysical signals exhibit (albeit subtle) site dependent phase synchronization index.