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.     

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.     

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.     

Instrumentation and Numerical Analysis of Cylindrical Diaphragm Wall Movement During Deep Excavation at Coastal Area

The lateral deflection of a cylindrical diaphragm wall and the associated ground movement induced by deep excavation are analyzed by performing site instrumentations and numerical analyses in the coastal area of Korea. Wall lateral deflection, rebar stress, and pore water pressure were measured and analyzed in eight directions. Variations of soil properties with the decrease of confining pressure are compared by performing various in situ tests before ad after excavation. To calculate the wall lateral deflection accurately, the effects of small strain nonlinearity, confining pressure, and the hysteresis loading/unloading loop developed during excavation are considered in the proposed numerical analysis. By comparing numerical results with measured ones, the importances of considering small strain nonlinearity and confining pressure reduction in the nonlinear (FEM) are emphasized. Also, the effects of wall stiffness on the performance of cylindrical diaphragm walls are studied for future similar excavation in the coastal area.     

Spectrophotometric determination of pH in seawater off Taiaroa Head,Otago, New Zealand: Full-spectrum modelling and prediction of pCO2 levels

Accurate measurement of seawater pH has long been sought by marine chemists (for example: [Dickson, A.G. 1993a. The measurement of sea water pH. Marine Chemistry, 44, 131–142, Dickson, A.G. 1993b. pH buffers for sea water media based on the total hydrogen ion concentration scale. Deep-Sea Research, 40, 107–118; Zhang, 1996; Tapp, M., Hunter, K.A., Currie, K. and Macaskill, B. 2000. Apparatus of continuous-flow underway spectrophotometric measurement of surface water pH. Marine Chemistry 72(2–4), 193–202; Friis, K., Koetzinger, A., Wallace, D.W.R. 2004. Spectrophotometric pH measurement in the ocean: Requirements, design and testing of an autonomous charge-coupled device detector system. Limnology and Oceanography: Methods 2, 126–136]. Recently, such attempts have taken on greater significance as anthropogenic carbon dioxide emissions may create rapidly changing oceanic pH. Spectrophotometric techniques have been accepted generally as the best for determination of seawater pH. Here we report a new technique using thymol blue as the indicator dye and fitting the entire spectrum from 400 to 900 nm rather than measuring the absorbance values at only two or three points in the spectrum. This full-spectrum modelling enables a reduction in signal to noise over other techniques. In the laboratory, we find with seawater samples a pH precision increase of five-fold “within” a sample and seven-fold “between” samples when comparing the full spectrum to the three-point method of analysis [Zhang, H., Byrne, R.H. 1996. Spectrophotometric pH measurements of surface seawater at in-situ conditions: absorbance and protonation behaviour of thymol blue. Marine Chemistry 52, 17–25].     

Estimating hydraulic conductivity using grain-size analyses, aquifer tests, and numerical modeling in a riverside alluvial system in South Korea

Hydraulic conductivity (K) for an alluvial system in a riverbank filtration area in Changwon City, South Korea, has been studied using grain-size distribution, pumping and slug tests, and numerical modeling. The alluvial system is composed of layers: upper fine sand, medium sand, lower fine sand, and a highly conductive sand/gravel layer at the base. The geometric mean of K for the sand/gravel layer (9.89?×?10^?4 m s^?1), as determined by grain-size analyses, was 3.33 times greater than the geometric mean obtained from pumping tests (2.97?×?10^?4 m s^?1). The geometric mean of K estimates obtained from slug tests (3.08?×?10^?6 m s^?1) was one to two orders of magnitude lower than that from pumping tests and grain-size analyses. K estimates derived from a numerical model were compared to those derived from the grain-size methods, slug tests and pumping tests in order to determine the degree of deviation from the numerical model. It is considered that the K estimates determined by the slug tests resemble the uppermost part of the alluvial deposit, whereas the K estimates obtained by grain-size analyses and pumping tests are similar to those from the numerical model for the sand/gravel layer of the riverside alluvial system.