The journal Communications in Numerical Methods in Engineering with Biomedical Applications becomes the International Journal for Numerical Methods in Biomedical Engineering (IJNMBE) from 1st January 2010

Communications in Numerical Methods in Engineering, founded by Roland W. Lewis in 1985, will change its title to the ‘International Journal for Numerical Methods in Biomedical Engineering’ and has a revised Aims and Scope. Copyright © 2009 John Wiley & Sons, Ltd.     

Bifurcation and second‐order work in geomaterials

In this paper, the ability of a material rate‐independent system to evolve toward another mechanical state from an equilibrium configuration, with no change in the control parameters, is investigated. From a mechanical point of view, this means that the system can spontaneously develop kinetic energy with no external disturbance from an equilibrium state, which corresponds to a particular case of bifurcation. The existence of both conjugate incremental strain and stress such that the second‐order work vanishes is established as a necessary and sufficient condition for the appearance of this bifurcation phenomenon. It is proved that this fundamental result is independent of the constitutive relation of the rate‐independent material considered. Then the case of homogeneous loading paths is investigated, and, as an illustration, the subsequent results are applied to interpret the well‐known liquefaction observed under isochoric triaxial loading conditions with loose granular materials. Copyright © 2006 John Wiley & Sons, Ltd.     

Constraining Holocene sea‐surface conditions in the south‐western Black Sea using dinoflagellate cysts

This paper presents two dinoflagellate cyst records from the south‐western shelf of the Black Sea. A new site, MAR05‐13, from the Sakarya shelf is described and placed into context with site MAR02‐45, ～250 km distant on the Thracian shelf. The records provide a centennial resolution of surface water conditions in the Holocene. Analysis of the data suggests that the surface salinity of the south‐western shelf increased in a gradual and progressive manner. In the period ～11 000–9000 cal a BP the assemblages suggest surface‐water salinities between 7–13 psu. The initial arrival of euryhaline species, ～8100 cal a BP, is linked to the reconnection of the Black Sea and Marmara Sea. The suggested surface water changes related to the reconnection took approximately 1000 years. Following this initial change in assemblages, a further increase in the number of euryhaline species is noted between 5000 and 4000 cal a BP. This is linked to the establishment of more saline surface‐water conditions, close to present‐day values. The record for MAR05‐13 highlights the complexity of the changes in cyst assemblages during the mid‐Holocene. Copyright © 2012 John Wiley & Sons, Ltd.     

Formal subdivision of the Holocene Series/Epoch: a Discussion Paper by a Working Group of INTIMATE (Integration of ice‐core,marine and terrestrial records) and the Subcommission on Quaternary Stratigraphy (International Commission on Stratigraphy)

This discussion paper, by a Working Group of INTIMATE (Integration of ice‐core, marine and terrestrial records) and the Subcommision on Quaternary Stratigraphy (SQS) of the International Commission on Stratigraphy (ICS), considers the prospects for a formal subdivision of the Holocene Series/Epoch. Although previous attempts to subdivide the Holocene have proved inconclusive, recent developments in Quaternary stratigraphy, notably the definition of the Pleistocene–Holocene boundary and the emergence of formal subdivisions of the Pleistocene Series/Epoch, mean that it may be timely to revisit this matter. The Quaternary literature reveals a widespread but variable informal usage of a tripartite division of the Holocene (‘early’, ‘middle’ or ‘mid’, and ‘late’), and we argue that this de facto subdivision should now be formalized to ensure consistency in stratigraphic terminology. We propose an Early–Middle Holocene Boundary at 8200 a BP and a Middle–Late Holocene Boundary at 4200 a BP, each of which is linked to a Global Stratotype Section and Point (GSSP). Should the proposal find a broad measure of support from the Quaternary community, a submission will be made to the International Union of Geological Sciences (IUGS), via the SQS and the ICS, for formal ratification of this subdivision of the Holocene Series/Epoch. Copyright © 2012 John Wiley & Sons, Ltd.     

Salinograph trends as indicators of the recession characteristics of stream components

Compared to hydrograph recession analysis, which is widely applied in engineering hydrology, the quantitative assessment of stream salinity with time (i.e. the salinograph) has received significantly less attention. In particular, while in many previous hydrological studies an inverse relationship between hydrograph and salinograph responses is apparent, the concept of salinity accession (the inversely related salinity counterpart to hydrograph recession) has not been introduced nor quantitatively evaluated in previous literature. In this study, we conduct a mathematical analysis of salinograph accession, and determine new quantitative relationships between salinity accession and hydrograph recession parameters. An equation is formulated that reproduces the general trend in salinity accession. A salinity accession parameter k_c is then introduced and is shown to be the ratio of direct runoff to total stream flow recession parameters: k_r/k. The groundwater recession parameter k_g was estimated using a simple and rapid method that uses both salinograph and hydrograph data. Salinity accession type‐curves illustrate that under certain conditions, the relative steepness of individual salinographs is dependent upon the ratio of groundwater salinity to direct runoff salinity: C_g/C_r. The salinity accession algorithms are applied to two contrasting field settings: Scott Creek, South Australia and Sandy Creek, northern Queensland, Australia. It was found that k_g > k during periods of obvious stream flow recession, for the events analysed. Salinograph accession behaviour was fairly similar for both sites, despite contrasting environments. Using assumed end‐member salinities for groundwater and direct runoff based upon field observations, the behaviour of k_c from the Scott Creek site was approximately reproduced by varying the initial groundwater to runoff flow ratio: Q_g0/Q_r0, within reasonable parameter ranges. The use of salinograph information when used in addition to standard hydrograph analyses provided useful information on recession characteristics of stream components. Copyright © 2008 John Wiley & Sons, Ltd.     

Impact of wildfire on source water contributions in Devil Creek,CA: evidence from end‐member mixing analysis

A geochemical and end‐member mixing analysis (EMMA) is undertaken in Devil Canyon catchment, located in southern California, to further understanding of watershed behaviour and source water contributions after an acute and extensive wildfire. Physical and chemical transformations in post‐fire watersheds are known to increase overland flow and decrease infiltration, mainly due to formation of a hydrophobic layer at, or near, the soil surface. However, less is known about subsurface flow response in burned watersheds. The current study incorporates EMMA to evaluate and quantify source water contributions before, and after, a catchment affected by wildfires in southern California during the fall of 2003. Pre‐ and post‐fire stream water data were available at several sampling sites within the catchment, allowing the identification of contributing water sources at varying spatial scales. Proposed end‐member observations (groundwater, overland flow, shallow subsurface flow) were also collected to constrain and develop the catchment mixing model. Post‐fire source water changes are more evident in the smaller and faster responding sub‐basin (interior sampling point). Early post‐fire storm events are dominated by overland flow with no significant soil water or groundwater flow contribution. Inter‐storm streamwater in this smaller basin shows an increase in groundwater and a decrease in soil water. In the larger, baseflow‐dominated system, source water components appear less affected by fire. A slight increase in lateral flow is observed with only a slight decrease in baseflow. Changes in the post‐fire flow regimes affect nutrient loading and chemical response of the basin. Relatively rapid recovery of the chaparral ecosystem is evidenced, with active re‐growth and evapotranspiration evidenced by the fourth post‐fire rainy season. Copyright © 2008 John Wiley & Sons, Ltd.     

Application of integrated back‐propagation network and self organizing map for flood forecasting

Much of the nonlinearity and uncertainty regarding the flood process is because hydrologic data required for estimation are often tremendously difficult to obtain. This study employed a back‐propagation network (BPN) as the main structure in flood forecasting to learn and to demonstrate the sophisticated nonlinear mapping relationship. However, a deterministic BPN model implies high uncertainty and poor consistency for verification work even when the learning performance is satisfactory for flood forecasting. Therefore, a novel procedure was proposed in this investigation which integrates linear transfer function (LTF) and self‐organizing map (SOM) to efficiently determine the intervals of weights and biases of a flood forecasting neural network to avoid the above problems. A SOM network with classification ability was applied to the solutions and parameters of the BPN model in the learning stage, to classify the network parameter rules and to obtain the winning parameters. The outcomes from the previous stage were then used as the ranges of the parameters in the recall stage. Finally, a case study was carried out in Wu‐Shi basin to demonstrate the effectiveness of the proposal. Copyright © 2009 John Wiley & Sons, Ltd.     

A fuzzy dynamic flood routing model for natural channels

A fuzzy dynamic flood routing model (FDFRM) for natural channels is presented, wherein the flood wave can be approximated to a monoclinal wave. This study is based on modification of an earlier published work by the same authors, where the nature of the wave was of gravity type. Momentum equation of the dynamic wave model is replaced by a fuzzy rule based model, while retaining the continuity equation in its complete form. Hence, the FDFRM gets rid of the assumptions associated with the momentum equation. Also, it overcomes the necessity of calculating friction slope (S_f) in flood routing and hence the associated uncertainties are eliminated. The fuzzy rule based model is developed on an equation for wave velocity, which is obtained in terms of discontinuities in the gradient of flow parameters. The channel reach is divided into a number of approximately uniform sub‐reaches. Training set required for development of the fuzzy rule based model for each sub‐reach is obtained from discharge‐area relationship at its mean section. For highly heterogeneous sub‐reaches, optimized fuzzy rule based models are obtained by means of a neuro‐fuzzy algorithm. For demonstration, the FDFRM is applied to flood routing problems in a fictitious channel with single uniform reach, in a fictitious channel with two uniform sub‐reaches and also in a natural channel with a number of approximately uniform sub‐reaches. It is observed that in cases of the fictitious channels, the FDFRM outputs match well with those of an implicit numerical model (INM), which solves the dynamic wave equations using an implicit numerical scheme. For the natural channel, the FDFRM outputs are comparable to those of the HEC‐RAS model. Copyright © 2009 John Wiley & Sons, Ltd.     

Spectrum of Fluctuations of Solar Noise Storms

The calculation of the Fourier transform of noise storm (NS) fluctuations showed that the power spectrum was adequately described by the expression G(F)∼1/F. Our results rule out the possibility that NS radiation is formed from random, short-term bursts (so-called type I bursts), since the spectrum of the sum of random short fluctuations is flat, but the real NS has a hyperbolic spectrum. This spectrum is monotonic and does not contain any components that exceed the level of the statistical fluctuations (i.e., the results of observations do not reveal the presence of periodic or resonant properties of the emission source). The hyperbolic shape of the spectrum shows that the main energy of a NS is contained in the slower temporal fluctuations.     

Effects of land‐use changes on hydrological processes in the middle basin of the Heihe River,northwest China

The effects of land‐use changes on the runoff process in the midstream plain of this arid inland river basin are a key factor in the rational allocation of water resources to the middle and lower reaches. The question is whether and by how much increasingly heavy land use impacts the hydrological processes in such an arid inland river basin. The catchment of the Heihe River, one of the largest inland rivers in the arid region of northwest China, was chosen to investigate the hydrological responses to land‐use change. Flow duration curves were used to detect trends and variations in runoff between the upper and lower reaches. Relationships among precipitation, upstream runoff, and hydrological variables were identified to distinguish the effects of climatic changes and upstream runoff changes on middle and downstream runoff processes. The quantitative relation between midstream cultivated land use and various parameters of downstream runoff processes were analysed using the four periods of land‐use data since 1956. The Volterra numerical function relation of the hydrological non‐linear system response was utilized to develop a multifactor hydrological response simulation model based on the three factors of precipitation, upstream runoff, and cultivated land area. The results showed that, since 1967, the medium‐ and high‐coverage natural grassland area in the midstream region has decreased by 80·1%, and the downstream runoff has declined by 27·32% due to the continuous expansion of the cultivated land area. The contribution of cultivated land expansion to the impact on the annual total runoff is 14–31%, on the annual, spring and winter base flow it is 44–75%, and on spring and winter discharge it is 23–64%. Once the water conservation plan dominated by land‐use structural adjustments is implemented over the next 5 years, the mean annual discharge in the lower reach could increase by 8·98% and the spring discharge by 26·28%. This will significantly alleviate the imbalance between water supply and demand in both its quantity and temporal distribution in the middle and lower reaches. Copyright © 2006 John Wiley & Sons, Ltd.