Pixel‐scale evaluation of SSM/I sea‐ice algorithms in the marginal ice zone during early fall freeze‐up

Observed reduction in recent sea ice areal extent and thickness has focused attention on the fact that the Arctic marine system appears to be responding to global‐scale climate variability and change. Passive microwave remote‐sensing data are the primary source underpinning these reports, yet problems remain in geophysical inversion of information on ice type and concentration. Uncertainty in sea‐ice concentration (SIC) retrievals is highest in the summer and fall, when water occurs in liquid phase within the snow–sea‐ice system. Of particular scientific interest is the timing and rate of new ice formation due to the control that this form of sea ice has on mass, energy and gas fluxes across the ocean–sea‐ice–atmosphere interface. In this paper we examine the critical fall freeze‐up period using in situ data from a ship‐based and aerial survey programme known as the Canadian Arctic Shelf Exchange study combined with microwave and optical Earth observations data. Results show that: (1) the overall physical conditions observed from aerial survey photography were well matched with coincident moderate‐resolution imaging spectroradiometer data and Radarsat ScanSAR imagery; (2) the shortwave albedo was linearly related to old ice concentration derived from survey photography; (3) the three SSM/I SIC algorithms (NASA Team (NT), NASA Team 2 (NT2), and Bootstrap (BT)) showed considerable discrepancies in pixel‐scale comparison with the Radarsat ScanSAR SICs well calibrated by the aerial survey data. The major causes of the discrepancies are attributed to (1) the inherent inability to detect the new thin ice in the NT and BT algorithms, (2) mismatches of the thin‐ice tie point of the NT2 algorithm, and (3) sub‐pixel ambiguity between the thin ice and the mixture of open water and sea ice. These results suggest the need for finer resolution of passive microwave sensors, such as AMSR‐E, to improve the precision of the SSM/I SIC algorithms in the marginal ice zone during early fall freeze‐up. Copyright © 2006 John Wiley & Sons, Ltd.     

Bulk organic δ13C and C/N ratios as palaeosalinity indicators within a Scottish isolation basin

Microfossils in isolation basin sediments are frequently used to reconstruct sea‐level change, but preservation problems and non‐analogue situations can limit their usefulness. Here we investigate the potential of stable carbon isotopes (δ¹³C) and C/N ratios from bulk organic matter, as an alternative proxy of salinity within isolation basin sediments from a basin in northwest Scotland. Within the Holocene sediment δ¹³C and C/N are determined largely by the mean weighted values of the predominant source of the organic material. Analysis of modern materials and comparison with the diatom record shows that the marine parts of the sequence are dominated by high δ¹³C and variable C/N. In the fresh water sequences the organic material is a mixture of both freshwater aquatic and terrestrial plant input that have relatively low δ¹³C and high C/N. The application of δ¹³C and C/N ratios in the studied basin in general follow the environmental change recorded by the diatoms and shows the potential of bulk organic matter in the investigation of salinity change in isolation basins. Copyright © 2005 John Wiley & Sons, Ltd.     

A non‐coaxial critical state soil model and its application to simple shear simulations

The yield vertex non‐coaxial theory is implemented into a critical state soil model, CASM (Int. J. Numer. Anal. Meth. Geomech. 1998; 22 :621–653) to investigate the non‐coaxial influences on the stress–strain simulations of real soil behaviour in the presence of principal stress rotations. The CASM is a unified clay and sand model, developed based on the soil critical state concept and the state parameter concept. Without loss of simplicity, it is capable of simulating the behaviour of sands and clays within a wide range of densities. The non‐coaxial CASM is employed to simulate the simple shear responses of Erksak sand and Weald clay under different densities and initial stress states. Dependence of the soil behaviour on the Lode angle and different plastic flow rules in the deviatoric plane are also considered in the study of non‐coaxial influences. All the predictions indicate that the use of the non‐coaxial model makes the orientations of the principal stress and the principal strain rate different during the early stage of shearing, and they approach the same ultimate values with an increase in loading. These ultimate orientations are dependent on the density of soils, and independent of their initial stress states. The use of the non‐coaxial model also softens the shear stress evolutions, compared with the coaxial model. It is also found that the ultimate shear strengths by using the coaxial and non‐coaxial models are dependent on the plastic flow rules in the deviatoric plane. Copyright © 2006 John Wiley & Sons, Ltd.     

Visualization of material stiffness in geomechanics analysis

This paper presents novel visualization techniques to simplify representation of the fourth‐order material stiffness tensor as a set of three‐dimensional geometric objects. Stiffness visualization aids in understanding the complex stiffness characteristics of highly non‐linear constitutive models including modelled material anisotropy and loading path dependent stiffness variation. Stiffness visualization is relevant for understanding the relationship of material stiffness to global behaviour in the analysis of a boundary value problem. The spherical pulse stiffness visualization method, developed in the acoustics field, is extended to visualize stiffness of geomaterials using three three‐dimensional objects. This method is limited to relatively simple constitutive models with symmetric stiffness matrices insensitive to loading magnitude and direction. A strain dependent stiffness visualization method is developed that allows the examination of material stiffness for a range of loading directions and is suitable for highly non‐linear and path dependent material models. The proposed stiffness visualization can be represented as 3‐D, 2‐D and 1‐D objects. The visualization technique is used to represent material stiffness and its evolution during simulated soil laboratory tests and deep excavation construction. Copyright © 2005 John Wiley & Sons, Ltd.     

Hydrometeorological behaviour of pine and larch forests in eastern Siberia

Seasonal changes in the water and energy exchanges over a pine forest in eastern Siberia were investigated and compared with published data from a nearby larch forest. Continuous observations (April to August 2000) were made of the eddy‐correlation sensible heat flux and latent heat flux above the canopy. The energy balance was almost closed, although the sum of the turbulent fluxes sometimes exceeded the available energy flux (R_n ? G) when the latent heat flux was large; this was related to the wind direction. We examined the seasonal variation in energy balance components at this site. The seasonal variation and magnitude of the sensible heat flux (H) was similar to that of the latent heat flux (λE), with maximum values occurring in mid‐June. Consequently, the Bowen ratio was around 1·0 on many days during the study period. On some clear days just after rainfall, λE was very large and the sum of H and λE exceeded R_n ? G. The evapotranspiration rate above the dry canopy from May to August was 2·2 mm day^?1. The contributions of understory evapotranspiration (E_u) and overstory transpiration (E_o) to the evapotranspiration of the entire ecosystem (E_t) were both from 25 to 50% throughout the period analysed. These results suggest that E_u plays a very important role in the water cycle at this site. From snowmelt through the tree growth season (23 April to 19 August 2000), the total incoming water, comprised of the sum of precipitation and the water equivalent of the snow at the beginning of the melt season, was 228 mm. Total evapotranspiration from the forest, including interception loss and evaporation from the soil when the canopy was wet, was 208–254 mm. The difference between the incoming and outgoing amounts in the water balance was from +20 to ?26 mm. The water and energy exchanges of the pine and larch forest differed in that λE and H increased slowly in the pine forest, whereas λE increased rapidly in the larch forest and H decreased sharply after the melting season. Consequently, the shape of the Bowen ratio curves at the two sites differed over the period analysed, as a result of the differences in the species in each forest and in soil thawing. Copyright © 2003 John Wiley & Sons, Ltd.     

Dynamics and radiation of young type-Ia supernova remnants: Important physical processes

Astronomy Letters - We analyze the physical processes that should be taken into account when modeling young type-Ia supernova remnants (SNRs) with ages of several hundred years in which forward...     

Comparison of infiltration models

Fourteen popular, representative infiltration models, some physically based, some semi‐empirical and some empirical, were selected for a comparative evaluation. Using the Nash and Sutcliffe efficiency criterion, the models were evaluated and compared for 243 sets of infiltration data collected from field and laboratory tests conducted in India and the USA on soils ranging from coarse sand to fine clay. Based on a relative grading scale, the semi‐empirical Singh–Yu general model, Holtan model and Horton model were graded respectively as 6·52, 5·57 and 5·48 out of 10. The empirical Huggins and Monke model, modified Kostiakov and Kostiakov model were graded as 5·57, 5·30 and 5·22, respectively. The physically based non‐linear and linear models of Smith–Parlange were graded as 5·48 and 5·22, respectively. Other models were ranked lower than these models. All the models generally performed poorly in field tests on Georgia's sandy soils, except the Robertsdale loamy sand. Copyright © 2003 John Wiley & Sons, Ltd.     

Seasonal evolution of meltwater generation,storage and discharge at a non‐temperate glacier in Svalbard

In glacierized catchments, meteorological inputs driving surface melting are translated into runoff outputs mediated by the glacier hydrological system: analysis of the relationship between meteorology and diurnal and seasonal patterns of runoff should reflect the functioning of that system, with the role of meltwater storage likely to be of particular importance. Daily meltwater storage is determined for a glacier at 78 °N in the Svalbard archipelago, by comparing inputs calculated from a surface energy balance model with measured outputs (proglacial discharge). Solar radiation, air temperature, wind speed and proglacial discharge are then analysed by regression and time‐series methods, in order to assess the meteorology–discharge relationship and its variation at diurnal and seasonal time‐scales. The recorded discharge time‐series can be divided into two contrasting intervals: up to early August, proglacial discharge was high and variable, mean hydrographs showed little indication of diurnal cycling, ARIMA models of discharge indicated a non‐seasonal, moving‐average generating process, and there was a net loss of meltwater from storage; from early August, proglacial discharge was low and relatively invariable, but with clearer diurnal cycles, regression models of discharge showed substantially improved correlations with air temperature and solar radiation, ARIMA models indicated a non‐seasonal, autoregressive generating process, and eventually a seasonal component, and there was a net gain in meltwater storage. The transition between the two periods is brief compared with the duration of the melt season. The runoff response to meteorology therefore lacks the strongly progressive element previously identified in mid‐latitude glacierized catchments. In particular, the glacier hydrological system only appears responsive to diurnal forcing following the depletion of the seasonal snowpack meltwater store. Copyright © 2001 John Wiley & Sons, Ltd.     

Performance of Jacobi preconditioning in Krylov subspace solution of finite element equations

This paper examines the performance of the Jacobi preconditioner when used with two Krylov subspace iterative methods. The number of iterations needed for convergence was shown to be different for drained, undrained and consolidation problems, even for similar condition number. The differences were due to differences in the eigenvalue distribution, which cannot be completely described by the condition number alone. For drained problems involving large stiffness ratios between different material zones, ill‐conditioning is caused by these large stiffness ratios. Since Jacobi preconditioning operates on degrees‐of‐freedom, it effectively homogenizes the different spatial sub‐domains. The undrained problem, modelled as a nearly incompressible problem, is much more resistant to Jacobi preconditioning, because its ill‐conditioning arises from the large stiffness ratios between volumetric and distortional deformational modes, many of which involve the similar spatial domains or sub‐domains. The consolidation problem has two sets of degrees‐of‐freedom, namely displacement and pore pressure. Some of the eigenvalues are displacement dominated whereas others are excess pore pressure dominated. Jacobi preconditioning compresses the displacement‐dominated eigenvalues in a similar manner as the drained problem, but pore‐pressure‐dominated eigenvalues are often over‐scaled. Convergence can be accelerated if this over‐scaling is recognized and corrected for. Copyright © 2002 John Wiley & Sons, Ltd.     

Eagle Nebula Pillars: From Models to Observations

Over the past few years, our group has been developing hydrodynamic models to simulate formation of the Eagle Nebula pillars. The true test of any model is, of course, how well it can reproduce the observations. Here, we discuss how we go about testing our models against observations. We describe the process by which we “observe” the model data to create synthetic maps. We show an example of this technique using one of our model runs and compare the resultant synthetic map to the real one.