Engineering topics
Publications and source records attributed to Freilich, M. H..
The European Center for Medium-Range Weather Forecasts (ECMWF) and Institut Francais Pour la Recherche et 1'Exploitation de la Mer (IFREMER) European Remote-Sensing Satellite (ERS-1), named IFR2, surface wind velocity data products are compared during July 1995 over the Arabian Sea.
Spatial variations of the east-west and north-south components of surface wind stress are critical in studies of ocean circulation and biological-physical interactions because surface wind stress curl produces a vertical velocity in the upper ocean at the bottom of the Ekman Layer.The ERS-1 scatterometer provides reasonable coverage and direct measurements of vector of winds. Three schemes are evaluated relative to high-quality moored-bouy wind observations recorded in the central Arabian Sea, where high surface waves and high atmospheric water content during the southeast monsoon adversely affect the estimation of satellite-derived winds.
A statistical technique is developed for determining fully empirical model functions relating altimeter backscatter (sigma(sub 0)) measurements to near-surface neutral stability wind speed. By assuming that sigma(sub 0) varies monotonically and uniquely with wind speed, the method requires knowledge only of the separate, rather than joint distribution functions of sigma(sub 0) and wind speed. Analytic simplifications result from using a Weibull distribution to approximate the global ocean wind speed distribution; several different wind data sets are used to demonstrate the validity of the Weibull approximation. The technique has been applied to 1 year of Geosat data. Validation of the new and historical model functions using an independent buoy data set demonstrates that the present model function not only has small overall bias and root mean square (RMS) errors, but yields smaller systematic error trends with wind speed and pseudowave age than previously published models. The present analysis suggests that generally accuracte altimeter model functions can be derived without the use of colocated measurements, nor is additional significant wave height information measured by the altimeter necessary.
A 10 m height wind vector data set, named CMODFD, from ERS-1 Active Microwave Instrument (AMI) measurements, was compared with moored buoy wind observations. The CMODFD wind vectors compared favorably with moored buoy wind observations at about 60 sites. Satisfactory results were obtained from intercomparisons between CMODFD and moored buoy estimates of horizontal wind divergence and wind stress curl. The CMODFD positions of the Intertropical Convergence Zone (ITCZ) at 28 deg west were nearly equal to those determined from satellite images of clouds. The CMODFD data were employed to compute the Ekman vertical velocity in the region of the Pacific North Equatorial Countercurrent (NECC), and a correlative relation was found between the annual cycles and longitudinal variations of the NECC and Ekman vertical motion.
Backscatter cross section measurements from the ERS-1 Active Microwave Instrument (AMI) scatterometer were reprocessed to vector winds using the Freilich-Dunbar model function and a meteorologically aided ambiguity removal scheme. This consistent data set was used to examine the wind field over the Southern Ocean from 20 to 60 deg South. The large number of ERS-1 measurements allows relatively accurate calculation of annual mean wind, stress, and curl fields as well as overall statistics of the winds at mid to high southern latitudes. The long duration of the data time series allows preliminary examination of low frequency (semi annual) wind variability.
Shoaling wave fields generated in laboratory experiments were analyzed to determine the sensitivity of nonlinear interactions to the directional distributions of incident waves. Peaks in the directional spectra observed in shallow water were consistent with near resonating, quadratic interactions between two primary waves transferring energy to a third wave with the sum frequency and vector sum wavenumber of the primary waves. Directionally colinear waves forced a higher-frequency wave propagating in the same direction as the primary waves, while directionally spread (i.e., noncolinear) primary waves forced a higher-frequency wave that propagated in a direction between those of the interacting primary waves. Deepwater wave fields with similar frequency spectra but different directional spectra evolved to different shallow-water directional spectra, yet their shallow-water frequency spectra were remarkably similar. This result suggests that the shape of the directional spectrum of the incident wave field has only a small effect on the magnitudes of nonlinear energy transfer during shoaling. The principal effect of directionality in the incident wave field is on the directions, not the amplitudes, of the nonlinearly generated waves. The laboratory data demonstrate clearly the importance of triad interactions between noncolinear and colinear shoaling waves.
Freilich and Dunbar (1993) produced a 10-m height wind vector data set, named CMODFD, for the European Space Agency (ESA) first European Remonte Sensing (ERS-1) Active Microwave Instrument (AMI) measurement.The CMODFD wind vectors compared favorably with moored-buoy wind observations at about 60 sites.
Backscatter cross section (sigma(sub 0)) measurements from the ERS-1 scatterometer are collocated with surface wind velocity data from operational Numerical Weather Prediction (NWP) analyses to yield a fully empirical C band model function. The empirical model does not obey a power law at speeds less than 5 m/s, and upwind/crosswind ratios are small for all incidence angles at low wind speeds and for small incidence angles at all wind speeds. Downwind sigma(sub 0) exceeds upwind values for incidence angles below approximately 28 degrees. The full ERS-1 scatterometer data set has been reprocessed using the NWP model function.
A Jan. - Jun. 1992 10 m height wind vector field was computed with the CMOD2 model function. Freilich and Dunbar used ESA measurements of AMI (Active Microwave Instrument) sigma-naught to derive a refined model function and to produce an alternate Jan. - Jun. 1992 10 m height wind vector field, named FD data. The ESA and FD model functions and algorithms were substantially different. The FD wind vectors were in better agreement with the climatological mean wind field and with moored buoy wind measurements. The average correlation coefficients between monthly mean and daily averaged FD and moored buoy matchups of the Cartesian wind components were 0.84 and 0.74, respectively. The average root mean square differences between the monthly mean and daily averaged FD and moored buoy matchups of the Cartesian wind components were 1.25 and 3.0 m/s, respectively.
Scatterometry and scatterometer design issues are reviewed. The design of the NASA Scatterometer (NSCAT) to be flown on the Japanese ADEOS mission is presented. Building on Seasat experience, the NSCAT system includes several enhancements, such as three antenna azimuths in each of two swaths, and an onboard digital Doppler processor to allow backscatter measurements to be colocated everywhere within the orbit. These enhancements will greatly increase the quality of the NSCAT wind data. The ground processing of data is discussed, and scatterometers of the next decade are briefly described.
The predictions of linear and nonlinear (Boussinesq) shoaling wave models for nonbreaking unidirectional surface gravity waves are compared to field observations, with particular emphasis on quantities that may be important for cross-shore sediment transport. The extensive data sets were obtained on two natural beaches, span water depths between 1 and 10 m, and include incident wave power spectra with narrow, broad, and bimodal shapes. Significant wave heights varied between approximately 30 and 100 cm, and peak periods between approximately 8 and 18 s. The evolution of total variances of sea surface elevation, cross-shore velocity, and horizontal acceleration is modeled at least qualitatively well by both linear and nonlinear theories. Only the nonlinear theory predicts the increasingly asymmetric sea surface elevations and horizontal velocities (pitched-forward wave shapes) and the weaker variation of skewness (difference between crest and trough profiles) which are observed to occur during shoaling. The nonlinear theory also models qualitatively well the large skewed accelerations which occur during the passage of asymmetric waves.
The rapid spatial recurrence of weakly nonlinear and weakly dispersive progressive shallow-water waves is examined using a numerical integration technique on the discretized and truncated form of the Boussinesq equations. This study primarily examines recurrence in wave fields with Ursell number O(1) and characterizes the sensitivity of recurrence to initial spectral shape and number of allowed frequency modes. It is shown that the rapid spatial recurrence is not an inherent property of the considered Boussinesq systems for evolution distances of 10-50 wavelengths. The main result of the study is that highly truncated Boussinesq models of resonant shallow-water ocean surface gravity waves predict rapid multiple recurrence cycles, but that this is an artifact dependent on the number of allowed modes. For initial conditions consisting of essentially all energy concentrated in a single mode, damping of the recurrence cycles increases as the number of low-power background modes increases. When more than 32 modes are allowed, the recurrence behavior is relatively insensitive to the number of allowed modes.
The spatial evolution of a directionally spread wave field on a near-planar natural beach is examined using data from longshore arrays of pressure sensors and wave staffs at depths 10.3 and 4.1 m. High-resolution frequency-directional spectra from the deeper array are used to initialize a linear refraction model, and the resulting model predictions are compared with frequency-directional measurements at the shallow array. Linear theory inaccurately predicts both the shapes of directional spectra in shallow water and the total variances in some frequency bands. The discrepancies are largest for frequencies associated with maxima in the bicoherence spectrum, suggesting the importance of nonlinear effects.
Laboratory and field measurements of suspended sediment in the nearshore suggest that fluid accelerations are an important factor in sediment transport by oscillatory waves. Here, Eulerian accelerations of the cross-shore velocity are calculated from measurements of velocity obtained by an array of bottom-mounted electromagnetic flow meters spanning a natural surf zone. Large shoreward accelerations of brief duration are associated with the steep front faces of both near-breaking and breaking waves. Weaker offshore accelerations of longer duration occur during passage of the more gently sloped rear faces. The acceleration field is thus strongly skewed in the shoreward direction. Power spectra and bispectra indicate, as expected, that statistics of the acceleration field are significantly influenced by high-frequency motions but are rather insensitive to surf beat.
The statistics of Sxy estimates derived from orthogonal-component measurements are examined. Based on results of Goodman (1957), the probability density function (pdf) for Sxy(f) estimates is derived, and a closed-form solution for arbitrary moments of the distribution is obtained. Characteristic functions are used to derive the exact pdf of Sxy(tot). In practice, a simple Gaussian approximation is found to be highly accurate even for relatively few degrees of freedom. Implications for experiment design are discussed, and a maximum-likelihood estimator for a posterior estimation is outlined.
Vector winds measured by the Seasat-A Satellite Scatterometer (SASS) are analyzed to determine the spatial structure of oceanic surface winds over wavelengths from 200 to 2200 km. The analysis is performed in four latitudinal bands in the Pacific Ocean. The salient features of the results are summarized as follows: (1) for each of the four geographic regions, the spectra of meridional and zonal wind components and of kinetic energy are consistent with a power-law dependence on wavenumber; for midlatitude regions in both the Northern and Southern hemispheres the wave-number dependence of kinetic energy is k exp -2.2, while for tropical regions in both hemispheres it is k exp-1.9, (2) for each individual region, the spectral dependence on wavenumber is nearly the same for both velocity components and for kinetic energy, (3) comparisons of zonal and meridional component spectra indicate that midlatitude winds may be isotropic, while tropical winds may be significantly anisotropic, and (4) the coherence between wind components is small everywhere.
The National Aeronautics and Space Administration scatterometer (NSCAT) is to be flown as part of the Navy Remote Ocean Sensing System (N-ROSS) scheduled for launch in 1989. The NSCAT will provide frequent accurate and high-resolution measurements of vector winds over the global oceans. NSCAT data will be applicable to a wide range of studies in oceanography, meteorology, and instrument science. The N-ROSS mission, is outlined, are described. The capabilities of the NSCAT flight instrument and an associated NASA research ground data-processing and distribution system, and representative oceanographic meteorological, and instrument science studies that may benefit from NSCAT data are surveyed.