Engineering PapersSearch

Engineering topics

Zukor, Dorothy

Publications and source records attributed to Zukor, Dorothy.

Software Engineering Support of the Third Round of Scientific Grand Challenge Investigations: Earth System Modeling Software Framework Survey

One of the most significant challenges in large-scale climate modeling, as well as in high-performance computing in other scientific fields, is that of effectively integrating many software models from multiple contributors. A software framework facilitates the integration task, both in the development and runtime stages of the simulation. Effective software frameworks reduce the programming burden for the investigators, freeing them to focus more on the science and less on the parallel communication implementation. while maintaining high performance across numerous supercomputer and workstation architectures. This document surveys numerous software frameworks for potential use in Earth science modeling. Several frameworks are evaluated in depth, including Parallel Object-Oriented Methods and Applications (POOMA), Cactus (from (he relativistic physics community), Overture, Goddard Earth Modeling System (GEMS), the National Center for Atmospheric Research Flux Coupler, and UCLA/UCB Distributed Data Broker (DDB). Frameworks evaluated in less detail include ROOT, Parallel Application Workspace (PAWS), and Advanced Large-Scale Integrated Computational Environment (ALICE). A host of other frameworks and related tools are referenced in this context. The frameworks are evaluated individually and also compared with each other.

Talbot, Bryan

Robust Matching of Wavelet Features for Sub-Pixel Registration of Landsat Data

For many Earth and Space Science applications, automatic geo-registration at sub-pixel accuracy has become a necessity. In this work, we are focusing on building an operational system, which will provide a sub-pixel accuracy registration of Landsat-5 and Landsat-7 data. The input to our registration method consists of scenes that have been geometrically and radiometrically corrected. Such pre-processed scenes are then geo-registered relative to a database of Landsat chips. The method assumes a transformation composed of a rotation and a translation, and utilizes rotation- and translation-invariant wavelets to extract image features that are matched using statistically robust feature matching and a generalized Hausdorff distance metric. The registration process is described and results on four Landsat input scenes of the Washington, D.C. area are presented.

LeMoigne, Jacqueline

Semiclassical Models for Virtual Antiparticle Pairs, the Unit of Charge e, and the QCD Coupling alpha(sub s)

New semiclassical models of virtual antiparticle pairs are used to compute the pair lifetimes, and good agreement with the Heisenberg lifetimes from quantum field theory (QFT) is found. The modeling method applies to both the electromagnetic and color forces. Evaluation of the action integral of potential field fluctuation for each interaction potential yields approximately Planck's constant/2 for both electromagnetic and color fluctuations, in agreement with QFT. Thus each model is a quantized semiclassical representation for such virtual antiparticle pairs, to good approximation. When the results of the new models and QFT are combined, formulae for e and alpha(sub s)(q) are derived in terms of only Planck's constant and c.

Batchelor, David

Improved Humidity Sensing with the Chilled Mirror: Really?

Testing of the chilled mirror sensor was initiated at Wallops Island in 1997. The chilled mirror dew point system is integrated with the Sippican, Inc., MK2 radiosonde providing a relatively inexpensive instrument. Early tests suggested that better stratospheric humidity measurements might be available. But, recent tests with different configurations of the mirror's cooler indicated that the systems capability to cool to the very low dew point temperatures required in the stratosphere were not being met with the present system. Nonetheless, the present mirror technology, while still undergoing development gives better humidity information between 400 and 100 hPa than the current routine radiosonde sensor. Comparisons are given of the chilled mirror and the typical operational humidity sensors.

Schmidlin, Francis J.

ECC Ozonesonde Calibration and Observations: Satellite Validation

The reliability of the Electrochemical Concentration Cell (ECC) ozonesonde depends on the care exercised in preparing the instrument for use. Although the ECC can be quickly prepared and flown, generally within less then one day if necessary, it is best to prepare the instrument at least one week prior to use, and as our tests have confirmed even 2-3 weeks prior to use may actually be better. There are a number of factors that must be considered when preparing an ECC. These basically are the pump efficiency, volumetric flow rate, temperature of the air entering the pump, and the background current. Also of importance is the concentration of the potassium iodide solution. Tests conducted at Wallops Island (38 N) has enabled us to identify potential problem areas and ways to avoid them. The calibration and pre-flight preparation methods will be discussed. The method of calibrating the ECC also is used at Ascension Island (8 S) and Natal, Brazil (5 S). Comparisons between vertical profiles of the ECC instrument and satellites will be reviewed as well as comparison with ground based instruments, such as, the Dobson Spectrophotometer and hand held Microtops photometers.

Schmidlin, Francis J.

Semi-Classical Models for Virtual Antiparticle Pairs

Virtual particle-antiparticle pairs of massive elementary particle& are predicted in Quantum Field Theory (QFT) to appear from the vacuum and annihilate each other again within their Heisenberg lifetimes h/4mc(exp 2). In this work, semiclassical models of this process - for the cases of massive leptons, quarks, and the massive weak bosons W and Z - are constructed. It is shown that the dynamical lifetime of the particle- antiparticle system in each case equals the Heisenberg lifetime to good approximation, and obeys appropriate quantization conditions on the field fluctuation action. In other words, the dynamical lifetime of the semiclassical model agrees with QED and QCD to good approximation. But the formula for the dynamical lifetime in each model includes the force strength coupling constant (e in the lepton case, alpha(sup s) (q(exp 2)) in the quark cases), while the Heisenberg lifetime formula does not. Observing the agreement of the Heisenberg and dynamical lifetimes, we may derive the QED and QCD coupling constants in terms of h, c, and numerical factors only.

Batchelor, David

Comparison of Implicit Collocation Methods for the Heat Equation

We combine a high-order compact finite difference scheme to approximate spatial derivatives arid collocation techniques for the time component to numerically solve the two dimensional heat equation. We use two approaches to implement the collocation methods. The first one is based on an explicit computation of the coefficients of polynomials and the second one relies on differential quadrature. We compare them by studying their merits and analyzing their numerical performance. All our computations, based on parallel algorithms, are carried out on the CRAY SV1.

Kouatchou, Jules

Enhancement of Tropical Land Cover Mapping with Wavelet-Based Fusion and Unsupervised Clustering of SAR and Landsat Image Data

The characterization and the mapping of land cover/land use of forest areas, such as the Central African rainforest, is a very complex task. This complexity is mainly due to the extent of such areas and, as a consequence, to the lack of full and continuous cloud-free coverage of those large regions by one single remote sensing instrument, In order to provide improved vegetation maps of Central Africa and to develop forest monitoring techniques for applications at the local and regional scales, we propose to utilize multi-sensor remote sensing observations coupled with in-situ data. Fusion and clustering of multi-sensor data are the first steps towards the development of such a forest monitoring system. In this paper, we will describe some preliminary experiments involving the fusion of SAR and Landsat image data of the Lope Reserve in Gabon. Similarly to previous fusion studies, our fusion method is wavelet-based. The fusion provides a new image data set which contains more detailed texture features and preserves the large homogeneous regions that are observed by the Thematic Mapper sensor. The fusion step is followed by unsupervised clustering and provides a vegetation map of the area.

LeMoigne, Jacqueline

Multi-Sensor Registration of Earth Remotely Sensed Imagery

Assuming that approximate registration is given within a few pixels by a systematic correction system, we develop automatic image registration methods for multi-sensor data with the goal of achieving sub-pixel accuracy. Automatic image registration is usually defined by three steps; feature extraction, feature matching, and data resampling or fusion. Our previous work focused on image correlation methods based on the use of different features. In this paper, we study different feature matching techniques and present five algorithms where the features are either original gray levels or wavelet-like features, and the feature matching is based on gradient descent optimization, statistical robust matching, and mutual information. These algorithms are tested and compared on several multi-sensor datasets covering one of the EOS Core Sites, the Konza Prairie in Kansas, from four different sensors: IKONOS (4m), Landsat-7/ETM+ (30m), MODIS (500m), and SeaWIFS (1000m).

LeMoigne, Jacqueline

Performance of the Wavelet Decomposition on Massively Parallel Architectures

Traditionally, Fourier Transforms have been utilized for performing signal analysis and representation. But although it is straightforward to reconstruct a signal from its Fourier transform, no local description of the signal is included in its Fourier representation. To alleviate this problem, Windowed Fourier transforms and then wavelet transforms have been introduced, and it has been proven that wavelets give a better localization than traditional Fourier transforms, as well as a better division of the time- or space-frequency plane than Windowed Fourier transforms. Because of these properties and after the development of several fast algorithms for computing the wavelet representation of any signal, in particular the Multi-Resolution Analysis (MRA) developed by Mallat, wavelet transforms have increasingly been applied to signal analysis problems, especially real-life problems, in which speed is critical. In this paper we present and compare efficient wavelet decomposition algorithms on different parallel architectures. We report and analyze experimental measurements, using NASA remotely sensed images. Results show that our algorithms achieve significant performance gains on current high performance parallel systems, and meet scientific applications and multimedia requirements. The extensive performance measurements collected over a number of high-performance computer systems have revealed important architectural characteristics of these systems, in relation to the processing demands of the wavelet decomposition of digital images.

El-Ghazawi, Tarek A.

Satellite Observed Variability in Antarctic and Arctic Surface Temperatures and Their Correlation to Open Water Areas

Recent studies using meterological station data have indicated that global surface air temperature has been increasing at a rate of 0.05 K/decade. Using the same set of data but for stations in the Antarctic and Arctic regions (>50 N) only, the increases in temperature were 0.08, and 0.22 K/decade, when record lengths of 100 and 50 years, respectively, were used. To gain insights into the increasing rate of warming, satellite infrared and passive microwave observations over the Arctic region during the last 20 years were processed and analyzed. The results show that during this period, the ice extent in the Antarctic has been increasing at the rate of 1.2% per decade while the surface temperature has been decreasing at about 0.08 K per decade. Conversely, in the Northern Hemisphere, the ice extent has been decreasing at a rate of 2.8% per decade, while the surface temperatures have been increasing at the rate of 0.38 K per decade. In the Antarctic, it is surprising that there is a short term trend of cooling during a global period of warming. Very large anomalies in open water areas in the Arctic were observed especially in the western region, that includes the Beaufort Sea, where the observed open water area was about 1x10(exp 6) sq km, about twice the average for the region, during the summer of 1998. In the eastern region, that includes the Laptev Sea, the area of open water was also abnormally large in the summer of 1995. Note that globally, the warmest and second warmest years in this century, were 1998 and 1995, respectively. The data, however, show large spatial variability with the open water area distribution showing a cyclic periodicity of about ten years, which is akin to the North Atlantic and Arctic Oscillations. This was observed in both western and eastern regions but with the phase of one lagging the other by about two years. This makes it difficult to interpret what the trends really mean. But although the record length of satellite data is still relatively short and the climate trend difficult to establish, the immediate impact of a continued warming trend may be very profound.

Comiso, Josefino C.

Sea Ice Motion from Wavelet Analysis of Satellite Data

Wavelet analysis of NASA scatterometer (NSCAT) backscatter and Defense Meteorological Satellite Program (DMSP) Special Sensor Microwave/Imager (SSM/I) radiance data can be used to obtain daily sea ice drift information for the Arctic region. This technique provides improved spatial coverage over the existing array of Arctic Ocean buoys and better temporal resolution over techniques utilizing data from satellite synthetic aperture radars. Comparisons with ice motion derived from ocean buoys give good quantitative agreement. Both comparison results from NSCAT and SSM/I are compatible, and the results from NSCAT can definitely complement that from SSM/I when there are cloud or surface effects. Then three sea-ice drift daily results from NSCAT, SSM/I, and buoy data can be merged as a composite map by some data fusion techniques. The ice flow streamlines are highly correlated with surface air pressure contours. Examples of derived ice-drift maps in December 1996 illustrate large-scale circulation reversals over a period of four days. A method for deriving divergence and shear at the large-scale has been developed and comparison between buoys and satellite results shows a good agreement. These calibrated/validated results indicate that NSCAT, SSM/I merged daily ice motion are suitably accurate to identify and closely locate sea ice processes, and to improve our current knowledge of sea ice drift and related processes through the data assimilation of ocean-ice numerical model. For demonstration purpose, the ice velocities derived from satellite data are compared with the ice velocities derived from a coupled ice-ocean interaction model. The comparison reveals that the general circulation patterns of the two are quite similar but the ice velocity differences between the two are quite significant. In order to quantify the wind effects on ice motion, empirical orthogonal functions (EOF) are used in the principal component analysis for both ice motion and pressure field. Some preliminary results of sea-ice motion from QuikScat will also be presented.

Liu, Antony K.

Internal Wave Study in the South China Sea Using SAR

Recently, the internal wave distribution maps in the China Seas have been compiled from hundreds of ERS-1/2, RADARSAT, and Space Shuttle SAR (Synthetic Aperture Radar) images from 1993 to 1999. Based on internal wave distribution map, most of internal waves in the northeast part of South China Sea were propagating westward. The wave crest can be as long as 200 km with amplitude of 100 m due to strong current from the Kuroshio branching out into the South China Sea. Based on the observations from drilling rigs near DongSha Island by Amoco Production Co., the solitons may be generated in a 4 km wide channel between Batan and Sabtang islands in Luzon Strait. The proposed generation mechanism is similar to the lee wave formation from a shallow topography. Both depression and elevation internal waves have been observed in the same RADARSAT ScanSAR image on May 4, 1998 near DongSha Island. Furthermore, depression and elevation internal waves have also been observed by SAR at the same location on the shelf in April and June, 1993 (in different seasons) respectively. Numerical models have been used to interpret their generation mechanism and evolution processes. Based on the SAR images, near DongSha Island, the westward propagating huge internal solitons are often encountered and diffracted/broken by the coral reefs on the shelf. After passing the island, the diffracted waves will re-merge or interact with each other. It has been observed that after the nonlinear wave-wave interaction, the phase of wave packet is shifted and wavelength is also changed. Examples of mesoscale features observed in SAR images, such as fronts, raincells, bathymetry, ship wakes, and oil spills will be presented. Recent mooring measurements in April 1999 near Dongsha Island, future field test ASIAEX (Asian Seas International Acoustics Experiment) planned for April 2001, and some pretest survey data will be discussed in this paper.

Liu, Antony K.

The Dependence of Ocean Surface Emissivity on Wind Vector as Measured with TMR

Global TMR brightness temperature observations at 18, 21, and 37 GHz have been colocated with near-simultaneous SeaWinds wind vector data as well as with a monthly SST climatological product. The combined data allow us to study the dependence of ocean surface emissivity, at each frequency, upon both wind speed and direction. Results show a clear two-branch wind speed dependence; weak and linear below 6 m/s with an abrupt increase in sensitivity above that point. Our analysis also shows that the nadir-view ocean surface emissivity depends on the angle between the wind direction and TMR's antenna polarization orientation.

Tran, N.

An Altimeter Wind Speed Model Using Both Radar Backscatter and Significant Wave Height

The development and validation of a new altimeter wind speed model will be presented. This algorithm provides a direct mapping of TOPEX-measured backscatter and significant wave height to 10 m wind speed. A large scatterometer/altimeter crossover data set was assembled to develop the routine and several large ancillary data sets have been assembled for validation purposes. Validation results suggest that this two input routine provides marginal, yet measurable improvements over the standard single-parameter MCW algorithm.

Vandemark, D.

Sea Ice Variability in the Sea of Okhotsk from Passive Microwave Satellite Observations

The Sea of Okhotsk, located between 50 and 60 N, is bounded by the Kamchatka Peninsula, Siberia, Sakhalin Island, and the Kuril Island chain and is the largest midlatitude seasonal sea ice zone in the Northern Hemisphere. The winter sea ice cover begins to form in November and expands to cover most of the sea by March. Over the following three months, the ice retreats with only small ice-covered areas remaining by the beginning of June. The sea is ice free or nearly ice free on average for six months of the year, from June through November. The recent compilation of a consistent, long-term record of Northern Hemisphere sea ice extents based on passive microwave satellite observations from the Nimbus 7 Scanning Multichannel Microwave Radiometer and from four Defense Meteorological Satellite Program Special Sensor Microwave Imagers provides the basis for assessing long-term sea ice extent variability in the Sea of Okhotsk. Analysis of this 20-year data record (1979-1998) shows that based on yearly averages the overall extent of the Sea of Okhotsk ice cover is decreasing at the rate of -8.1+/-2.1x10(exp 3) sq km/yr (-17.2%/decade), in contrast to the rate of decrease of -33.3+/-0.7x10(exp 3) sq km/yr (-2.7%/decade) for the Northern Hemisphere as a whole. There is large regional sea ice extent variability of the Arctic ice cover. Two of the nine Arctic regions analyzed, the Bering Sea and the Gulf of St. Lawrence, show increases of 0.8+/-1.4xl0(exp 3) sq km/yr (2.7%/decade) and 1.2+/-0.5xl0(exp 3) sq km/yr (17.1%/decade), respectively. Interestingly, the Sea of Okhotsk and the Gulf of St. Lawrence show about equal percentage changes, but of opposite sign. The Sea of Okhotsk exhibits its greatest percent decrease (-24.3%/decade) during spring (April-June). The year of maximum winter sea ice extent for the Sea of Okhotsk was 1979, whereas the minimum winter sea ice extent occurred in 1984.

Cavalieri, Donald J.