Wallops Flight Facility

International audience

Abstract. We present Bedmap2, a new suite of gridded products describing surface elevation, ice-thickness and the seafloor and subglacial bed elevation of the Antarctic south of 60° S. We derived these products using data from a variety of sources, including many substantial surveys completed since the original Bedmap compilation (Bedmap1) in 2001. In particular, the Bedmap2 ice thickness grid is made from 25 million measurements, over two orders of magnitude more than were used in Bedmap1. In most parts of Antarctica the subglacial landscape is visible in much greater detail than was previously available and the improved data-coverage has in many areas revealed the full scale of mountain ranges, valleys, basins and troughs, only fragments of which were previously indicated in local surveys. The derived statistics for Bedmap2 show that the volume of ice contained in the Antarctic ice sheet (27 million km3) and its potential contribution to sea-level rise (58 m) are similar to those of Bedmap1, but the mean thickness of the ice sheet is 4.6% greater, the mean depth of the bed beneath the grounded ice sheet is 72 m lower and the area of ice sheet grounded on bed below sea level is increased by 10%. The Bedmap2 compilation highlights several areas beneath the ice sheet where the bed elevation is substantially lower than the deepest bed indicated by Bedmap1. These products, along with grids of data coverage and uncertainty, provide new opportunities for detailed modelling of the past and future evolution of the Antarctic ice sheets.

▪ Abstract We survey the newly developed Hilbert spectral analysis method and its applications to Stokes waves, nonlinear wave evolution processes, the spectral form of the random wave field, and turbulence. Our emphasis is on the inadequacy of presently available methods in nonlinear and nonstationary data analysis. Hilbert spectral analysis is here proposed as an alternative. This new method provides not only a more precise definition of particular events in time-frequency space than wavelet analysis, but also more physically meaningful interpretations of the underlying dynamic processes.

The confidence limit is a standard measure of the accuracy of the result in any statistical analysis. Most of the confidence limits are derived as follows. The data are first divided into subsections and then, under the ergodic assumption, the temporal mean is substituted for the ensemble mean. Next, the confidence limit is defined as a range of standard deviations from this mean. However, such a confidence limit is valid only for linear and stationary processes. Furthermore, in order for the ergodic assumption to be valid, the subsections have to be statistically independent. For non‐stationary and nonlinear processes, such an analysis is no longer valid. The confidence limit of the method here termed EMD/HSA (for empirical mode decomposition/Hilbert spectral analysis) is introduced by using various adjustable stopping criteria in the sifting processes of the EMD step to generate a sample set of intrinsic mode functions (IMFs). The EMD technique acts as a pre‐processor for HSA on the original data, producing a set of components (IMFs) from the original data that equal the original data when added back together. Each IMF represents a scale in the data, from smallest to largest. The ensemble mean and standard deviation of the IMF sample sets obtained with different stopping criteria are calculated, and these form a simple random sample set. The confidence limit for EMD/HSA is then defined as a range of standard deviations from the ensemble mean. Without evoking the ergodic assumption, subdivision of the data stream into short sections is unnecessary; hence, the results and the confidence limit retain the full‐frequency resolution of the full dataset. This new confidence limit can be applied to the analysis of nonlinear and non‐stationary processes by these new techniques. Data from length‐of‐day measurements and a particularly violent recent earthquake are used to demonstrate how the confidence limit is obtained and applied. By providing a confidence limit for this new approach, a stable range of stopping criteria for the decomposition or sifting phase (EMD) has been established, making the results of the final processing with HSA, and the entire EMD/HSA method, more definitive.

Determining trend and implementing detrending operations are important steps in data analysis. Yet there is no precise definition of "trend" nor any logical algorithm for extracting it. As a result, various ad hoc extrinsic methods have been used to determine trend and to facilitate a detrending operation. In this article, a simple and logical definition of trend is given for any nonlinear and nonstationary time series as an intrinsically determined monotonic function within a certain temporal span (most often that of the data span), or a function in which there can be at most one extremum within that temporal span. Being intrinsic, the method to derive the trend has to be adaptive. This definition of trend also presumes the existence of a natural time scale. All these requirements suggest the Empirical Mode Decomposition (EMD) method as the logical choice of algorithm for extracting various trends from a data set. Once the trend is determined, the corresponding detrending operation can be implemented. With this definition of trend, the variability of the data on various time scales also can be derived naturally. Climate data are used to illustrate the determination of the intrinsic trend and natural variability.

Interferometric synthetic‐aperture radar data collected by ERS‐1/2 and Radarsat‐1 satellites show that Antarctic Peninsula glaciers sped up significantly following the collapse of Larsen B ice shelf in 2002. Hektoria, Green and Evans glaciers accelerated eightfold between 2000 and 2003 and decelerated moderately in 2003. Jorum and Crane glaciers accelerated twofold in early 2003 and threefold by the end of 2003. In contrast, Flask and Leppard glaciers, further south, did not accelerate as they are still buttressed by an ice shelf. The mass loss associated with the flow acceleration exceeds 27 km 3 per year, and ice is thinning at rates of tens of meters per year. We attribute this abrupt evolution of the glaciers to the removal of the buttressing ice shelf. The magnitude of the glacier changes illustrates the importance of ice shelves on ice sheet mass balance and contribution to sea level change.

Instantaneous frequency (IF) is necessary for understanding the detailed mechanisms for nonlinear and nonstationary processes. Historically, IF was computed from analytic signal (AS) through the Hilbert transform. This paper offers an overview of the difficulties involved in using AS, and two new methods to overcome the difficulties for computing IF. The first approach is to compute the quadrature (defined here as a simple 90° shift of phase angle) directly. The second approach is designated as the normalized Hilbert transform (NHT), which consists of applying the Hilbert transform to the empirically determined FM signals. Additionally, we have also introduced alternative methods to compute local frequency, the generalized zero-crossing (GZC), and the teager energy operator (TEO) methods. Through careful comparisons, we found that the NHT and direct quadrature gave the best overall performance. While the TEO method is the most localized, it is limited to data from linear processes, the GZC method is the most robust and accurate although limited to the mean frequency over a quarter wavelength of temporal resolution. With these results, we believe most of the problems associated with the IF determination are resolved, and a true time–frequency analysis is thus taking another step toward maturity.

The biogeochemistry of continental shelf systems plays an important role in the global elemental cycling of nitrogen and carbon, but remains poorly quantified. We have developed a high‐resolution physical‐biological model for the U.S. east coast continental shelf and adjacent deep ocean that is nested within a basin‐wide North Atlantic circulation model in order to estimate nitrogen fluxes in the shelf area of the Middle Atlantic Bight (MAB). Our biological model is a relatively simple representation of nitrogen cycling processes in the water column and organic matter remineralization at the water‐sediment interface that explicitly accounts for sediment denitrification. Climatological and regionally integrated means of nitrate, ammonium, and surface chlorophyll are compared with its model equivalents and were found to agree within 1 standard deviation. We also present regional means of primary production and denitrification, and statistical measures of chlorophyll pattern variability. A nitrogen budget for the MAB shows that the sediment denitrification flux is quantitatively important in determining the availability of fixed nitrogen and shelf primary production (it was found to remove 90% of all the nitrogen entering the MAB). Extrapolation of nitrogen fluxes estimated for the MAB to the North Atlantic basin suggests that shelf denitrification removes 2.3 × 10 12 mol N annually; this estimate exceeds estimates of N 2 fixation by up to an order of magnitude. Our results emphasize the importance of representing shelf processes in biogeochemical models.

Biomass‐burning plumes and haze layers were observed during the ABLE 2A flights in July/August 1985 over the central Amazon Basin. The haze layers occurred at altitudes between 1000 and 4000 m and were usually only some 100 to 300‐m thick but extended horizontally over several 100 km. They could be traced by satellite imaging and trajectory studies to biomass burning at the southern perimeter of the Amazon Basin, with transport times estimated to be 1–2 days. These layers strongly influenced the chemical and optical characteristics of the atmosphere over the eastern Amazon Basin. The concentrations of CO, CO 2 , O 3 , and NO were significantly elevated in the plumes and haze layers relative to the regional background. The NO/CO ratio in fresh plumes was much higher than in the aged haze layers, suggesting that more than 80% of the NO x in the haze layers had been converted to nitrate and organic nitrogen species subsequent to emission. The haze aerosol was composed predominantly of organic material, NH 4 + , K + , NO 3 − , SO 4 = , and anionic organic species (formate, acetate, and oxalate). While the concentrations of most aerosol ions were substantially higher in the haze layers than in the regional background aerosol, the ratios between the aerosol ions in the haze layer aerosols were very similar to those in the boundary layer aerosol over the central Amazon region. Simultaneous measurements of trace gas and aerosol species in the haze layers made it possible to derive emission ratios for CO, NO x , NH 3 , sulfur oxides, and aerosol constituents relative to CO 2 . Regional and global emission estimates based on these ratios indicate that biomass burning is an important contributor in the global and regional cycles of carbon, sulfur, and nitrogen species. Similar considerations suggest that photochemical ozone production in the biomass‐burning plumes contributes significantly to the regional ozone budget.

The Moderate Resolution Imaging Spectroradiometer (MODIS) will add a significant new capability for investigating the 70% of the Earth's surface that is covered by oceans, in addition to contributing to the continuation of a decadal scale time series necessary for climate change assessment in the oceans. Sensor capabilities of particular importance for improving the accuracy of ocean products include high SNR and high stability for narrow or spectral bands, improved onboard radiometric calibration and stability monitoring, and improved science data product algorithms. Spectral bands for resolving solar-stimulated chlorophyll fluorescence and a split window in the 4-/spl mu/m region for SST will result in important new global ocean science products for biology and physics. MODIS will return full global data at 1-km resolution. The complete suite of Levels 2 and 3 ocean products is reviewed, and many areas where MODIS data are expected to make significant, new contributions to the enhanced understanding of the oceans' role in understanding climate change are discussed. In providing a highly complementary and consistent set of observations of terrestrial, atmospheric, and ocean observations, MODIS data will provide important new information on the interactions between Earth's major components.

Abstract A new method, the Hilbert–Huang Transform (HHT), developed initially for natural and engineering sciences has now been applied to financial data. The HHT method is specially developed for analysing non‐linear and non‐stationary data. The method consists of two parts: (1) the empirical mode decomposition (EMD), and (2) the Hilbert spectral analysis. The key part of the method is the first step, the EMD, with which any complicated data set can be decomposed into a finite and often small number of intrinsic mode functions (IMF). An IMF is defined here as any function having the same number of zero‐crossing and extrema, and also having symmetric envelopes defined by the local maxima, and minima respectively. The IMF also thus admits well‐behaved Hilbert transforms. This decomposition method is adaptive, and, therefore, highly efficient. Since the decomposition is based on the local characteristic time scale of the data, it is applicable to non‐linear and non‐stationary processes. With the Hilbert transform, the IMF yield instantaneous frequencies as functions of time that give sharp identifications of imbedded structures. The final presentation of the results is an energy–frequency–time distribution, which we designate as the Hilbert Spectrum. Comparisons with Wavelet and Fourier analyses show the new method offers much better temporal and frequency resolutions. The EMD is also useful as a filter to extract variability of different scales. In the present application, HHT has been used to examine the changeability of the market, as a measure of volatility of the market. Published in 2003 by John Wiley & Sons, Ltd.

Four distinct meteorological regimes in the Amazon basin have been examined to distinguish the contributions from boundary layer aerosol and convective available potential energy (CAPE) to continental cloud structure and electrification. The lack of distinction in the electrical parameters (peak flash rate, lightning yield per unit rainfall) between aerosol‐rich October and aerosol‐poor November in the premonsoon regime casts doubt on a primary role for the aerosol in enhancing cloud electrification. Evidence for a substantial role for the aerosol in suppressing warm rain coalescence is identified in the most highly polluted period in early October. The electrical activity in this stage is qualitatively peculiar. During the easterly and westerly wind regimes of the wet season, the lightning yield per unit of rainfall is positively correlated with the aerosol concentration, but the electrical parameters are also correlated with CAPE, with a similar degree of scatter. Here cause and effect are difficult to establish with available observations. This ambiguity extends to the “green ocean” westerly regime, a distinctly maritime regime over a major continent with minimum aerosol concentration, minimum CAPE, and little if any lightning.

Recent advances in the determination of the mass balance of polar ice sheets show that the Greenland Ice Sheet is losing mass by near-coastal thinning, and that the West Antarctic Ice Sheet, with thickening in the west and thinning in the north, is probably thinning overall. The mass imbalance of the East Antarctic Ice Sheet is likely to be small, but even its sign cannot yet be determined. Large sectors of ice in southeast Greenland, the Amundsen Sea Embayment of West Antarctica, and the Antarctic Peninsula are changing quite rapidly as a result of processes not yet understood.

Since 1996, quality assurance experiments of electrochemical concentration cell (ECC) ozonesondes of two different model types (SPC‐6A and ENSCI‐Z) have been conducted in the environmental simulation facility at the Research Centre Juelich within the framework of the Juelich Ozone Sonde Intercomparison Experiment (JOSIE). The experiments have shown that the performance characteristics of the two ECC‐sonde types can be significantly different, even when operated under the same conditions. Particularly above 20 km the ENSCI‐Z sonde tends to measure 5–10% more ozone than the SPC‐6A sonde. Below 20 km the differences are 5% or less, but appear to show some differences with year of manufacture. There is a significant difference in the ozone readings when sondes of the same type are operated with different cathode sensing solutions. Testing the most commonly used sensing solutions showed that for each ECC‐manufacturer type the use of 1.0% KI and full buffer gives 5% larger ozone values compared with the use of 0.5% KI and half buffer, and as much as 10% larger values compared with 2.0% KI and no buffer. For ozone sounding stations performing long term measurements this means that changing the sensing solution type or ECC‐sonde type can easily introduce a change of ±5% or more in their records, affecting determination of ozone trends. Standardization of operating procedures for ECC‐sondes yields a precision better than ±(3–5)% and an accuracy of about ±(5–10)% up to 30 km altitude.

A network of 10 southern hemisphere tropical and subtropical stations, designated the Southern Hemisphere Additional Ozonesondes (SHADOZ) project and established from operational sites, provided over 1000 ozone profiles during the period 1998–2000. Balloon‐borne electrochemical concentration cell (ECC) ozonesondes, combined with standard radiosondes for pressure, temperature, and relative humidity measurements, collected profiles in the troposphere and lower to midstratosphere at: Ascension Island; Nairobi, Kenya; Irene, South Africa; Réunion Island; Watukosek, Java; Fiji; Tahiti; American Samoa; San Cristóbal, Galapagos; and Natal, Brazil. The archived data are available at: 〈 http://croc.gsfc.nasa.gov/shadoz 〉. 1 In this paper, uncertainties and accuracies within the SHADOZ ozone data set are evaluated by analyzing: (1) imprecisions in profiles and in methods of extrapolating ozone above balloon burst; (2) comparisons of column‐integrated total ozone from sondes with total ozone from the Earth‐Probe/Total Ozone Mapping Spectrometer (TOMS) satellite and ground‐based instruments; and (3) possible biases from station to station due to variations in ozonesonde characteristics. The key results are the following: (1) Ozonesonde precision is 5%. (2) Integrated total ozone column amounts from the sondes are usually to within 5% of independent measurements from ground‐based instruments at five SHADOZ sites and overpass measurements from the TOMS satellite (version 7 data). (3) Systematic variations in TOMS‐sonde offsets and in ground‐based‐sonde offsets from station to station reflect biases in sonde technique as well as in satellite retrieval. Discrepancies are present in both stratospheric and tropospheric ozone. (4) There is evidence for a zonal wave‐one pattern in total and tropospheric ozone, but not in stratospheric ozone.

Repeated laser‐altimeter surveys and modelled snowfall/summer melt show average ice loss from Greenland between 1997 and 2003 was 80 ± 12 km 3 yr −1 , compared to about 60 km 3 yr −1 for 1993/4–1998/9. Half of the increase was from higher summer melting, with the rest caused by velocities of some glaciers exceeding those needed to balance upstream snow accumulation. Velocities of one large glacier almost doubled between 1997 and 2003, resulting in net loss from its drainage basin by about 20 km 3 of ice between 2002 and 2003.

Aircraft laser-altimeter surveys over northern Greenland in 1994 and 1999 have been coupled with previously reported data from southern Greenland to analyze the recent mass-balance of the Greenland Ice Sheet. Above 2000 meters elevation, the ice sheet is in balance on average but has some regions of local thickening or thinning. Thinning predominates at lower elevations, with rates exceeding 1 meter per year close to the coast. Interpolation of our results between flight lines indicates a net loss of about 51 cubic kilometers of ice per year from the entire ice sheet, sufficient to raise sea level by 0.13 millimeter per year-approximately 7% of the observed rise.

Abstract A comparative study of raindrop size distribution measurements has been conducted at NASA’s Goddard Space Flight Center where the focus was to evaluate the performance of the upgraded laser-optical OTT Particle Size Velocity (Parsivel2; P2) disdrometer. The experimental setup included a collocated pair of tipping-bucket rain gauges, OTT Parsivel (P1) and P2 disdrometers, and Joss–Waldvogel (JW) disdrometers. Excellent agreement between the two collocated rain gauges enabled their use as a relative reference for event rain totals. A comparison of event total showed that the P2 had a 6% absolute bias with respect to the reference gauges, considerably lower than the P1 and JW disdrometers. Good agreement was also evident between the JW and P2 in hourly raindrop spectra for drop diameters between 0.5 and 4 mm. The P2 drop concentrations mostly increased toward small sizes, and the peak concentrations were mostly observed in the first three measurable size bins. The P1, on the other hand, underestimated small drops and overestimated the large drops, particularly in heavy rain rates. From the analysis performed, it appears that the P2 is an improvement over the P1 model for both drop size and rainfall measurements. P2 mean fall velocities follow accepted terminal fall speed relationships at drop sizes less than 1 mm. As a caveat, the P2 had approximately 1 m s−1 slower mean fall speed with respect to the terminal fall speed near 1 mm, and the difference between the mean measured and terminal fall speeds reduced with increasing drop size. This caveat was recognized as a software bug by the manufacturer and is currently being investigated.

Recent aircraft and satellite laser altimeter surveys of the Amundsen Sea sector of West Antarctica show that local glaciers are discharging about 250 cubic kilometers of ice per year to the ocean, almost 60% more than is accumulated within their catchment basins. This discharge is sufficient to raise sea level by more than 0.2 millimeters per year. Glacier thinning rates near the coast during 2002-2003 are much larger than those observed during the 1990s. Most of these glaciers flow into floating ice shelves over bedrock up to hundreds of meters deeper than previous estimates, providing exit routes for ice from further inland if ice-sheet collapse is under way.

Several new data sets reveal that thinning and speedup of Jakobshavn Isbrae continue, following its recent rapid increase in speed as its floating ice tongue disintegrated. The present speedup rate of ∼5% a −1 over much of the fast‐moving region appears to be a diffusive response to the initial much larger speedup near the front. There is strong seasonality in speed over much of the fast‐flowing main trunk that shows a good inverse correlation with the seasonally varying length of a short (typically ∼6 km) floating ice tongue. This modulation of speed with ice front position supports the hypothesis that the major speedup was caused by loss of the larger floating ice tongue from 1998 to 2003. Analysis of image time series suggests that the transient winter ice tongue is formed when sea ice bonds glacier ice in the fjord to produce a nearly rigid mass that almost entirely suppresses calving. Major calving only resumes in late winter when much of this ice clears from the fjord. The collapse of the ice tongue in the late 1990s followed almost immediately after a sharp decline in winter sea‐ice concentration in Disko Bay. This decline may have extended the length of the calving season for several consecutive years, leading to the ice tongue's collapse.