Ice‐cover variability on shallow lakes at high latitudes: model simulations and observations

A one‐dimensional thermodynamic model for simulating lake‐ice phenology is presented and evaluated. The model can be driven with observed daily or hourly atmospheric forcing of air temperature, relative humidity, wind speed, cloud amount and snowfall. In addition to computing the energy balance components, key model output includes the temperature profile at an arbitrary number of levels within the ice/snow (or the water temperature if there is no ice) and ice thickness (clear ice and snow‐ice) on a daily basis, as well as freeze‐up and break‐up dates. The lake‐ice model is used to simulate ice‐growth processes on shallow lakes in arctic, sub‐arctic, and high‐boreal forest environments. Model output is compared with field and remote sensing observations gathered over several ice seasons. Simulated ice thickness, including snow‐ice formation, compares favourably with field measurements. Ice‐on and ice‐off dates are also well simulated when compared with field and satellite observations, with a mean absolute difference of 2 days. Model simulations and observations illustrate the key role that snow cover plays on the seasonal evolution of ice thickness and the timing of spring break‐up. It is also shown that lake morphometry, depth in particular, is a determinant of ice‐off dates for shallow lakes at high latitudes. Copyright © 2003 John Wiley & Sons, Ltd.     

Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.     

CO<Subscript>2</Subscript> Flux Measurements over the Sea Surface by Eddy Correlation and Aerodynamic Techniques

The measurements of the vertical transport of CO₂ were carried out over the Sea of Japan using the specially designed pier of Kyoto University on September 20 to 22, 2000. CO₂ fluxes were measured by the eddy correlation and aerodynamic techniques. Both techniques showed comparable CO₂ fluxes during sea breeze conditions: −0.001 to −0.08 mg m⁻²s⁻¹ with the mean of −0.05 mg m⁻²s⁻¹. This means that the measuring site satisfies the fetch requirement for meteorological observations under sea breeze conditions. Moreover, the eddy diffusivity coefficient used in the aerodynamic technique is found to be consistent with the coefficient used in the eddy correlation technique. The present result leads us to conclude that the aerodynamic technique may be applicable to underway CO₂ flux measurements over the ocean and may be used in place of the bulk technique. The important point is the need to maintain a measuring accuracy of CO₂ concentration difference of the order of 0.1 ppmv on the research vessels or the buoys.     

Freely floating-body simulation by a 2D fully nonlinear numerical wave tank

Nonlinear interactions between large waves and freely floating bodies are investigated by a 2D fully nonlinear numerical wave tank (NWT). The fully nonlinear 2D NWT is developed based on the potential theory, MEL/material-node time-marching approach, and boundary element method (BEM). A robust and stable 4th-order Runge–Kutta fully updated time-integration scheme is used with regriding (every time step) and smoothing (every five steps). A special φ_n-η type numerical beach on the free surface is developed to minimize wave reflection from end-wall and wave maker. The acceleration-potential formulation and direct mode-decomposition method are used for calculating the time derivative of velocity potential. The indirect mode-decomposition method is also independently developed for cross-checking. The present fully nonlinear simulations for a 2D freely floating barge are compared with the corresponding linear results, Nojiri and Murayama’s (Trans. West-Jpn. Soc. Nav. Archit. 51 (1975)) experimental results, and Tanizawa and Minami’s (Abstract for the 6th Symposium on Nonlinear and Free-surface Flow, 1998) fully nonlinear simulation results. It is shown that the fully nonlinear results converge to the corresponding linear results as incident wave heights decrease. A noticeable discrepancy between linear and fully nonlinear simulations is observed near the resonance area, where the second and third harmonic sway forces are even bigger than the first harmonic component causing highly nonlinear features in sway time series. The surprisingly large second harmonic heave forces in short waves are also successfully reproduced. The fully updated time-marching scheme is found to be much more robust than the frozen-coefficient method in fully nonlinear simulations with floating bodies. To compare the role of free-surface and body-surface nonlinearities, the body-nonlinear-only case with linearized free-surface condition was separately developed and simulated.     

On the use of poorly soluble hygroscopic substances for modification of warm clouds and fogs

From analysis of the relationships between the equilibrium state of condensation nuclei and the relative humidity, a conclusion is made concerning the preferred use of poorly soluble substances for the generation of artificial condensation nuclei (ACN) designed for modification of warm clouds and fogs with the purpose of their dissipation and precipitation formation. The advantage of poorly soluble substances over commonly used soluble salts is that the finely dispersed part of the spectrum of the poorly soluble ACN does not deliquesce and so is not involved in the formation of cloud droplets. For experimental testing of the conclusions, preference was given to cement, whose main soluble substance is calcium oxide with a solubility of about 1 g/l. The spectrum of dispersed cement particles was measured and compared with the spectra of pyrotechnic flares widely used for modification at present. The process of formation of the cloud droplet spectra was simulated in the aerosol chamber of the Institute of Experimental Meteorology by decreasing the preliminarily generated excess pressure. It was found in these experiments that, compared to the spectra of particles formed on background condensation nuclei, the introduction of dispersed cement leads to the broadening of spectra and to a decrease in the concentration of droplets. Even at the early stage of condensation, droplets with radii of ∼20 μm appear. In this case, no “overseeding” phenomenon is observed, which, for soluble substances, manifests itself in an increase in the concentration of cloud droplets with a large ACN concentration. These effects indicate that, according to the existing concepts about the mechanism of warm-cloud modification with hygroscopic substances, the introduction of poorly soluble ACN (in particular, dispersed cement) below the base of cumulus clouds should stimulate coagulation processes and accelerate rain-formation processes. Considerations are given that the introduction of poorly soluble ACN into the already existing cumulus or stratocumulus clouds or fogs should also result in the acceleration of precipitation-formation processes or cloud dispersal. Comparison of cement powder with the well-known means of warm-cloud and fog modification is carried out. Original Russian Text ? N.P. Romanov, A.S. Drofa, N.S. Kim, A.V. Savchenko, G.F. Yaskevich, 2006, published in Izvestiya AN. Fizika Atmosfery i Okeana, 2006, Vol. 42, No. 1, pp. 80–91.     

Scaling analysis of deposition from turbidity currents

Many oil-bearing sedimentary deposits are formed by the settling of particles from turbidity currents. Modeling sedimentary processes that form these turbidites enables the calculation of properties such as extent, depth, porosity and permeability of hydrocarbon-bearing reservoirs. This paper estimates the extent and thickness of turbidites from the initial conditions of the turbidity flow. This is achieved by the application of scaling analysis of the partial differential equations that govern the dynamics of and deposition from turbidity currents. We apply the results of scaling analysis to five modern submarine fans. The predicted and actual values of the dimensions of the fan deposits match well. We then compare the derived results against tabulated sizes of ancient turbidites. The comparisons are good as long as we correctly identify the flow regimes in which the deposition took place. The good agreements observed in the two cases show that the estimates obtained using scaling analysis can provide useful first-guess values for the dimensions of the deposits.     

A comparison of several wave spectra for the random wave diffraction by a semi-infinite breakwater

A comparison of the diffraction of multidirectional random waves using several selected wave spectrum models is presented in this paper. Six wave spectrum models, Bretschneider, Pierson–Moskowitz, ISSC, ITTC, Mitsuyasu, and JONSWAP spectrum, are considered. A discrete form for each of the given spectrum models is used to specify the incident wave conditions. Analytical solutions based on both the Fresnel integrals and polynomial approximations of the Fresnel integrals and numerical solutions of a boundary integral approach have been used to obtain the two-dimensional wave diffraction by a semi-infinite breakwater at uniform water depth. The diffraction of random waves is based on the cumulative superposition of linear diffraction solution. The results of predicted random wave diffraction for each of the given spectrum models are compared with those of the published physical model presented by Briggs et al. [1995. Wave diffraction around breakwater. Journal of Waterway, Port, Coastal and Ocean Engineering—ASCE 121(1), 23–35]. Reasonable agreement is obtained in all cases. The effect of the directional spreading function is also examined from the results of the random wave diffraction. Based on these comparisons, the present model for the analysis of various wave spectra is found to be an accurate and efficient tool for predicting the random wave field around a semi-infinite breakwater or inside a harbor of arbitrary geometry in practical applications.     

Sedimentology of a subtidal,tide-dominated sand body in the Yellow Sea,southwest Korea

Odanam Satoe, a subtidal, tide-dominated sand body in the Yellow Sea, Korea, is linear in plan and asymmetrical in cross-section. It consists of fine- to medium-grained, well-sorted subangular sand. Bedforms consist of high-amplitude (1–2 m) sandwaves on the lower flanks of the gentler-sloping bar surface, and medium-amplitude (0.5-1 m) sandwaves on the sand body trough adjoining the steeper face, the bar crest and shallower parts of the gently sloping bar surface. Bedforms are absent on the relatively steeper bar surface, which is characterized by 2° slopes. Bedform orientation on the gentler slope is oblique by 30° to the bar crest, parallel to the sand-body crest on the crest itself, and opposite to the steeper sand-body face in the trough below the steeper slope of the bar.Bottom current velocity data show that tidal currents are semi-rotary with a flood time—velocity asymmetry over the gentler slope, and ebb time—velocity asymmetry over the steeper slope during most of the tidal cycle. Tidal-current flow parallels bar elongation over the steeper slope, whereas over the gentler slope, tidal-current flow is directed at 30° to the bar crest and changes to normal to the crest one hour prior to low tide. Bedform orientation mapped with side-scan sonar shows agreement with these flow directions.Sand dispersal around the sand body is controlled by time—velocity asymmetry and partial rotary flow directions of tidal currents. This circulation causes not only a trapezoidal mode of grain dispersal, but also westerly migration of the sand body documented from comparative bathymetric surveys in 1964 and 1980.     

Fully nonlinear numerical wave tank (NWT) simulations and wave run-up prediction around 3-D structures

A finite-difference scheme and a modified marker-and-cell (MAC) algorithm have been developed to investigate the interactions of fully nonlinear waves with two- or three-dimensional structures of arbitrary shape. The Navier–Stokes (NS) and continuity equations are solved in the computational domain and the boundary values are updated at each time step by the finite-difference time-marching scheme in the framework of a rectangular coordinate system. The fully nonlinear kinematic free-surface condition is implemented by the marker-density function (MDF) technique developed for two fluid layers.To demonstrate the capability and accuracy of the present method, the numerical simulation of backstep flows with free-surface, and the numerical tests of the MDF technique with limit functions are conducted. The 3D program was then applied to nonlinear wave interactions with conical gravity platforms of circular and octagonal cross-sections. The numerical prediction of maximum wave run-up on arctic structures is compared with the prediction of the Shore Protection Manual (SPM) method and those of linear and second-order diffraction analyses based on potential theory and boundary element method (BEM). Through this comparison, the effects of non-linearity and viscosity on wave loading and run-up are discussed.