Estimation of snow water equivalent using microwave radiometry over Arctic first‐year sea ice

The magnitude and spatial distribution of snow on sea ice are both integral components of the ocean–sea‐ice–atmosphere system. Although there exists a number of algorithms to estimate the snow water equivalent (SWE) on terrestrial surfaces, to date there is no precise method to estimate SWE on sea ice. Physical snow properties and in situ microwave radiometry at 19, 37 and 85 GHz, V and H polarization were collected for a 10‐day period over 20 first‐year sea ice sites. We present and compare the in situ physical, electrical and microwave emission properties of snow over smooth Arctic first‐year sea ice for 19 of the 20 sites sampled. Physical processes creating the observed vertical patterns in the physical and electrical properties are discussed. An algorithm is then developed from the relationship between the SWE and the brightness temperature measured at 37 GHz (55°) H polarization and the air temperature. The multiple regression between these variables is able to account for over 90% of the variability in the measured SWE. This algorithm is validated with a small in situ data set collected during the 1999 field experiment. We then compare our data against the NASA snow thickness algorithm, designed as part of the NASA Earth Enterprise Program. The results indicated a lack of agreement between the NASA algorithm and the algorithm developed here. This lack of agreement is attributed to differences in scale between the Special Sensor Microwave/Imager and surface radiometers and to differences in the Antarctic versus Arctic snow physical and electrical properties. Copyright © 2003 John Wiley & Sons, Ltd.     

An Instantaneous-Frequency Filtering Technique for High-Frequency Radar Oceanography

High-frequency (HF) radar systems are remote sensing tools that can be used to measure oceanographic parameters. Problems can occur when using the conventional periodogram (PG) method for computing power spectral estimates from backscattered radar signals. Temporal and spatial inhomogeneities within the radar measurement region can cause distortion in the spectra. This paper describes an instantaneous-frequency (IF) filtering technique that has been developed to measure the first-order modulation contained within the radar signal. Successful removal of this modulation is shown to yield an increased quality and quantity of ocean measurements     

Quaternary stratigraphy of the buried valleys of the lower Red Deer River,Alberta, Canada

The Quaternary fills of the buried valleys of southern Alberta and Saskatchewan have provided a wealth of information for the reconstruction of the glacial-interglacial record of the western plains of Canada, and this paper reports on the previously unstudied stratigraphy of the buried Calgary Valley and its former tributaries in the lower Red Deer River area. We attempt to differentiate Empress Group sediments, which potentially relate to pre-glacial, interglacial/ interstadial and post-glacial lake and river deposition, using sedimentology, stratigraphy and palaeoecology. Twenty-nine stratigraphical logs indicate that Empress Group sediments have infilled a considerably large area of badlands and tributary coulees that once drained into the Calgary Valley, located 15 km to the north of Dinosaur Provincial Park. Radiocarbon dates of 52.4 ka, 27.4 ka and > 42.4 ka and glacially modified quartz grains suggest that at least some of the valley fills date to interglacial or interstadial periods and may be mid-Wisconsinan in age. However, outcrops of an older till overlying other valley fills suggest that the buried valleys were only partially excavated during interglacials/interstadials and that older (even pre-glacial) sediments could have survived. Subglacial channels, recognisable on air photographs, largely coincide with buried valley positions due to the preferential excavation of the Quaternary sediment by meltwater and are filled with post-glacial lake sediment from which a radiocarbon date of 16 ka BP was obtained. Pre-glacial and glacial/post-glacial Empress Group sediments are lithologically indistinct but cover a large time span in southern Alberta.     

Timing and tectonic setting of Grenvillian metamorphism—constraints from a transect along Georgian Bay,Ontario*

Systematic mapping of a transect along the well-exposed shores of Georgian Bay, Ontario, combined with the preliminary results of structural analysis, geochronology and metamorphic petrology, places some constraints on the geological setting of high-grade metamorphism in this part of the Central Gneiss Belt. Correlations within and between map units (gneiss associations) have allowed us to recognize five tectonic units that differ in various aspects of their lithology, metamorphic and plutonic history, and structural style. The lowest unit, which forms the footwall to a regional decollement, locally preserves relic pre-Grenvillian granulite facies assemblages reworked under amphibolite facies conditions during the Grenvillian orogeny. Tectonic units above the decollement apparently lack the early granulite facies metamorphism; out-of-sequence thrusting in the south produced a duplex-like structure. Two distinct stages of Grenvillian metamorphism are apparent. The earlier stage (c. 1160–1120 Ma) produced granulite facies assemblages in the Parry Sound domain and upper amphibolite facies assemblages in the Parry Island thrust sheet. The later stage (c. 1040–1020 Ma) involved widespread, dominantly upper amphibolite facies metamorphism within and beneath the duplex. Deformation and metamorphism recently reported from south and east of the Parry Sound domain at c. 1100–1040 Ma have not yet been documented along the Georgian Bay transect. The data suggest that early convergence was followed by a period of crustal thickening in the orogenic core south-east of the transect area, with further advance to the north-west during and after the waning stages of this deformation.     

Modelling the April 15 spectrum of Supernova 1993 J

We present the modeling of the ultraviolet and optical spectra obtained simultaneously on 1993 April 15 with the HST and at Lick Observatory. A Monte Carlo code is employed in the modeling and a comparison is made between models reported by different groups. With an atmosphere similar to the Sun in chemical composition, the observed spectral lines are well reproduced by a power law density structure of index around 20 except the strong H and HeI λ5876 lines which have peculiar absorption profiles. The photospheric velocity is found to be 9500 km/s and the blackbody temperature of the spectrum is 7990 K. For H and HeI λ5876, we suggest a two-component density structure which has a smoother layer located immediately outside the steeply decreasing inner envelope. The power law indices are most probably 20 and 3, respectively, with the transition point at about 13 000 km/s. In addition, this outer smooth layer serves to flatten the far UV spectrum as observed.     

Obtaining basic simulation data for a heterogeneous field with stratified marine soils

Thinly stratified sedimentary deposits in a heterogeneous field were investigated to obtain basic physical data for the simulation of water flow. A procedure is described which translates a thinly stratified soil profile into a number of functional layers using functional hydrological properties. A functional layer is defined as a combination of one or more soil horizons and should (i) be recognizable during a soil survey using an auger and (ii) show significantly different functional hydrological properties when compared with another functional layer. This procedure gave three easily recognizable functional layers. Sets of hydrological characteristics of these three functional layers were obtained by physical measurements of the soil and by estimation, using textural data for classification into a standard Dutch series. The performance of several combinations of these sets was tested by comparing simulated and measured soil matric potentials for seven plots during one year. The best simulation results were obtained if measured soil hydraulic characteristics were used for relatively homogeneous functional layers and if the soil hydraulic characteristics were estimated at each location for the most heterogeneous layer.     

Water transport from the North-east Irish Sea to western Scottish coastal waters: Further observations from time-trend matching of Sellafield radiocaesium

Radiocaesium isotopes, discharged into the North-east Irish Sea from the Sellafield (formerly Windscale) nuclear fuel reprocessing plant in Cumbria, have been employed as flow monitors to update and extend the record of coastal water movement from the Irish Sea to the Clyde Sea area and, further north, to Loch Etive. The temporal trends in radiocaesium levels have been used to determine the extent of water mixing en route and to define mean advection rates. Flow conditions from the Irish Sea have changed considerably since the mid-1970s, the residence time of northern Irish Sea waters being ～12 months during 1978–1980 inclusive. Average transport times of four and six months are estimated for the Sellafield to Clyde and Sellafield to Etive transects respectively. Sellafield ¹³⁷Cs levels in seawater were diluted by factors of 27 and 50 respectively during current movement to the Clyde and Etive areas. The decrease in salinity-corrected ¹³⁷Cs concentrations between the Clyde and Etive suggests that dilution by Atlantic water occurs, the latter mainly entering the Firth of Lorne from the west. The majority (～94%) of the radiocaesium supply to Loch Etive enters the Firth of Lorne via the portion of the coastal current circulating west of Islay, only ～6% arriving via the Sound of Jura.     

Stable nonlinear methods for sensor array processing;

Most nonlinear high-resolution bearing estimators are unstable in the presence of correlated noise, system phase errors, and multipath arrivals because they inadvertently place too much emphasis on unstable eigenvectors of the cross-sensor correlation matrix. For moderately correlated noise there will be sufficiently many stable eigenvectors to resolve and localize discrete sources. A method is given (the "stable nonlinear method" or SNLM) whereby reweighting of the eigenvectors is achieved implicitly, without actual calculation of the eigenvectors. This SNLM is compared with Capon's maximum likelihood method (MLM) in simulations of correlated noise, partially correlated signals, and phase errors, and is shown to provide good stability in the cases considered.