Modeling of Water Quality, Quantity, and SustainabilityView this Special Issue
Conceptual Model for Simulating the Adjustments of Bankfull Characteristics in the Lower Yellow River, China
We present a conceptual model for simulating the temporal adjustments in the banks of the Lower Yellow River (LYR). Basic conservation equations for mass, friction, and sediment transport capacity and the Exner equation were adopted to simulate the hydrodynamics underlying fluvial processes. The relationship between changing rates in bankfull width and depth, derived from quasiuniversal hydraulic geometries, was used as a closure for the hydrodynamic equations. On inputting the daily flow discharge and sediment load, the conceptual model successfully simulated the 30-year adjustments in the bankfull geometries of typical reaches of the LYR. The square of the correlating coefficient reached 0.74 for Huayuankou Station in the multiple-thread reach and exceeded 0.90 for Lijin Station in the meandering reach. This proposed model allows multiple dependent variables and the input of daily hydrological data for long-term simulations. This links the hydrodynamic and geomorphic processes in a fluvial river and has potential applicability to fluvial rivers undergoing significant adjustments.
The bankfull characteristics of alluvial rivers are basic research topics in fluvial processes [1, 2]. In recent years, various natural and anthropogenic disturbances (such as dam construction, water and soil conservation, enhanced water usage with rapid socioeconomic development, and evaporation and precipitation variation related to climate change) have altered the runoff and sediment yield in many large river basins. Concerns have been raised regarding adjustments to river geometries and their impacts on local flood defence and ecological safety. There is emerging interest in the temporal adjustment of bankfull characteristics for rivers subject to significant climate change or intense human interference .
Two kinds of approaches have been developed to quantify the variation in bankfull characteristics in terms of the timescales at which the channel adjustment is explored: geomorphic and hydrodynamic approaches. The geomorphic approach is usually based on geomorphic laws expressed by power law, hyperbolic, and exponential equations or relatively complex differential equations [4–7]. It is widely recognised that geomorphic systems respond rapidly, immediately following the disruption of the equilibrium state, but thereafter exhibit a declining rate of adjustment as a new equilibrium or relaxed state is approached [5, 8, 9]. Applications with years or decades of geomorphic models can be found for various alluvial river basins, such as the South Island of New Zealand , Lower Hunter Valley in NewSouthWales , sand-bed streams of the American Midwest , Sabie River in South Africa , Karst landscape evolution in the Cave City in Kentucky , Lower Yellow River (LYR) in China [15–18], basins in southeastern Australia , and North Fork Toutle River in Washington . The river channel slope , channel width , bankfull cross-section area and discharge [16, 17], bed level [9, 23], and channel width-depth ratio  have been considered as dependent variables of flow and the sediment volume flux per year.
The hydrodynamic approach reproduces channel response processes by accounting for microscale river dynamics. That is, the Saint Venant equations (continuity and momentum equations) for river flow, the sediment transport equation, and the river bed or bank deformation equations are set up to describe the instantaneous relations between channel adjustment and the incoming flow and sediment conditions. One-, two-, and three-dimensional models for either steady flow [25–29] or unsteady flow [30–32] have been used widely to investigate channel degradation or aggradation during flood events. Recently, much effort has shifted to the long-term simulation of channel responses. The formation and evolution of meandering or wandering rivers has been simulated considering the long-term fluvial actions [33–36]. For example, Sun et al.  reproduced the phenomenon of channel evolution, including the shallowing and widening of the initial channel, the rhythmical sequence of alternate bars and erosional bumps, and the occurrence of meandering abrupt shifts and bank cusps. Asahi et al.  captured the nature process of river meandering in a computational model that considers the effects of bank erosion, the process of land accretion along the inner banks of meander bends, and the formation of channel cut-offs.
In nature, the channel response is achieved by the cumulative effects of fluvial erosion/deposition processes, which relate to both geomorphic processes and hydrodynamic events. The geomorphic approach, particularly the rate law method, can capture the overall behaviour of the channel response to seasonal or annual flow and sediment data, but it cannot account for the effects of flood events. The channel response can be characterised by a set of variables, but the coordination of multiple dependent variables has not been explored. By contrast, the hydrodynamic approach is capable of representing the details of the channel response process at small timescales, such as bank erosion, bar migration, and channel shallowing/widening. When applied on large timescales, the computational limitations and requirements for accurate resistance parameters and boundary conditions may stop it from producing reliable simulation and predictions.
Here we propose a physics-based model linking hydrodynamic and geomorphic scales. It is intended to represent channel response processes subject to both hydrodynamic and geomorphic controls. The model adopts hydrodynamic equations and boundary equations with daily data as input. It enables the characterisation of flood events in large rivers (typically lasting for around 10 days) and the simulation of continuous response behaviour over a long-term period (years or decades). In a scientific sense, linking across hydrodynamic and geomorphic scales is an attempt to overcome the shortcomings of each scale by considering the flood details and obtaining the long-period channel response processes. It may be applied both to simulate fluvial processes over the past few decades and to forecast channel-forming processes with potential flow and sediment series in the future for large alluvial rivers like LYR.
2. Background of This Study
2.1. Overview of the Study Area
The Yellow River is the second longest river in China and supports 12% of the Chinese population (Figure 1). The river is about 5464 km long with a drainage area of 752,000 km2. The Yellow River is famous for its excessive sediment load and deficient flow. The long-term average sediment load at the Sanmenxia Dam is 1.6 billion tons per year, which ranks first in the world.
The LYR stretches from Huayuankou to Lijin, as shown in Figure 2. It is about 520 km long and exhibits variable channel patterns. It is usually divided into three geomorphological distinct reaches . The upper reach from Huayuankou to Gaocun is a typical wandering channel pattern. The lower reach, from Taochengpu to Lijin, is a stable meandering channel pattern. The intervening reach transitions from wandering to meandering. The LYR is a natural alluvial river because of the high suspended sediment load it carries and deposits. According to the observed data, the total deposition volume in the LYR reached about 5.52 billion m3 during the period from 1950 to 1999, of which 60% was deposited in the wandering reach.
Here, the Huayuankou and Lijin Stations were chosen to validate the effects of the model, because these two stations are representative of typical wandering and meandering reaches, respectively. They have different bankfull characteristics as a function of flow and sediment load.
For Lijin Station in the meandering reach, the bankfull characteristics are easy to determine by observing when the water fills the main channel without overtopping the banks of the floodplain (Figure 3). For Huayuankou Station, cross-sections usually exhibit a compound geometric shape (Figure 4). Given nonstationary flow and sediment discharge, the location of the main channel often shifts laterally, which makes it difficult to define the zone of the main channel and bankfull characteristics in the wandering reach. As before, earlier and later measurements of the cross-section and the upstream and downstream cross-sectional profiles were used for reference to determine the bankfull stage .
The deviation of the inner-annual distributions of flow and sediment discharge can cause different channel forms in the LYR. These two chosen reaches had similar annual mean values of flow and sediment discharge, but different annual distributions (Figure 5). The annual mean flow discharge from 1950 to 2002 was 1248 m3/s at Huayuankou Station and 1036 m3/s at Lijin Station, while the annual peak flow discharge (daily peak data) was 6014 and 4713 m3/s, respectively. The annual mean sediment discharge from 1950 to 2002 was 32 m3/s at Huayuankou Station and 26 m3/s at Lijin Station, while the annual peak sediment discharge (daily peak data) was 653 and 330 m3/s, respectively.
Using time series analysis, a decreasing trend was detected for the annual mean and daily peak flow and sediment discharge for the two stations over the past 50 years using the Mann-Kendall method . Two significant times marking mutation of the pattern of the annual mean flow discharge were detected at the stations (1969 and 1986) using the Fisher discriminant criterion . The two years were linked to dam construction on the Yellow River for the Liujiaxia and Longyangxia reservoirs.
2.2. Bankfull Characteristics in Typical Sections
Figure 6 shows the observed bankfull characteristics (bankfull discharge , bankfull area , bankfull width , and bankfull depth ) at Lijin Station from 1964 to 2000. All of the data were obtained from the measured profile data after the flood season each year. Table 1 shows the correlation matrix. From the figure and table, there are linear relations between and , and , and and . However, the data mining indicated that the bankfull width varied from 400 to 600 m. For this, there were no obvious linear relations. Since these bankfull characteristics were all determined by the incoming flow and sediment conditions on a large time scale; there is no coherent geomorphic law to interpret all processes.
The same analysis for Huayuankou in the wandering reach is shown in Figure 7 and Table 2. Nonlinear relations existed more widely here than at Lijin. Except for , none of the other relationships were obviously linear. The channel geometry changed so quickly and markedly that the bankfull width varied from 1800 to 3200 m. This showed greater variation and was more nonlinear than in the meandering reach.
3. Study Methodology
In previous geomorphic models, channel response to disturbances can be described using nonlinear decay functions, for which the rate law has been most widely used to describe relaxation paths and recovery times. The rate law takes the form of a negative exponential equation: Here, is a characteristic parameter of the fluvial system, is the asymptotic value of this parameter, is the decay rate that quantifies how rapidly approaches , and is the time since the disturbance. Many studies have shown that the rate law can simulate and predict the channel response to disturbances with seasonal or annual data [5, 40–42].
As mentioned above, the rate law is used to illustrate the channel forming macroprocess (101-102 years) without any physical mechanism. It smoothes the flood details and cannot account for the effects of the deviation of flood events. This paper focuses on the dynamic processes of the channel geometry of alluvial rivers. Bankfull parameters such as the bankfull discharge , area , width , and depth all vary with time. To capture the major characteristics of the fluvial response, but allow lower computational complexity, a quasiequilibrium state was assumed. The temporal variation of the flow discharge and sediment load was distinguished for typical flood events, but the governing equations of flow and sediment transport took the form for steady flows. This suggests an intermediate timescale between the hydrodynamic and geomorphic scales. Above this, the temporal variation is accounted for, but it is neglected below this. A suitable timescale for a large river like the LYR that has a typical flood event of 10 days is 1 day. This can account for the daily variation in the flow and sediment conditions during the flood event.
The cross-sections of the LYR are mostly wide and shallow with 400–4000 m river top widths and 1–5 m river depths. As a result, the wetted perimeter (m) can nearly be substituted by the top widths (m). The hydraulic radius (m) can be expressed as , where is the flow area (m2) and (m) is the average river depth expressed by . With the assumption of quasiequilibrium [43, 44], the flow continuity, momentum, sediment transport, and bed deformation equations can be expressed as where is the bankfull discharge (m3/s), is the bankfull channel width, is the bankfull channel depth, is the bankfull velocity (m/s), is the bankfull hydraulic radius (m), is the channel slope, and is the Manning coefficient determined by the median grain size of the bed material load and channel depth in the form .
Equation (4) is Zhang Ruijin’s sediment-carrying capacity formula , which is used widely for the LYR, where is the sediment concentration (kg/m3), is the sediment settling velocity (m/s) in muddy water, and and are parameters (at Huayuankou Station, and ; at Lijin Station, and ). is the riverbed porosity, is the bed elevation (m), is the total sediment transport rate in volume per unit stream width (m2/s), and is the streamwise coordinate. Of these equations, (5) (i.e., the Exner equation) is for unsteady flow, whereas the others are for steady flows. This is consistent with the fact that the adjustment of the channel form takes much longer than that of the flow and sediment conditions.
These four governing equations involve two independent variables, and , but five dependent variables in the bankfull state, , , , , and . On a medium timescale of 100-101 years, changes less than and . As is fixed, over long periods of time for a given and , the expected channel cross-section can be obtained in terms of and .
For an alluvial river like the LYR, when . Under these circumstances, the bankfull depth is The left-hand side of (6) is calculated from known data, and the right-hand side of the equation is an increasing function of . A trial method can be used to obtain the value of . Then, and can be obtained as From (6) and (7), and are determined by the incoming steady flow and sediment discharge. Nevertheless, for alluvial rivers like the LYR, the flow discharge and sediment volume are not steady on intra- or interannual timescales. As a result, series of and should be expected instead of a fixed value. In other words, a dynamic equation including and is needed. This is the key problem to be addressed in this paper.
We rewrite (5) in integral form where is the length of the river reach and is the volumetric sediment transport rate (m3/s). With the assumption of a rectangular channel, we have For rivers in a stationary state, the relationship between the geometrical dimension (i.e., and ) and water discharge that passes through this cross-section takes the form [8, 47] where and are coefficients that vary with the characteristics of the river basin. For the statistical results of most alluvial rivers across the world , and are exponents, usually with and .
For a nonstationary state, is a time-dependent variable, and steady values of and do not exist. Nonetheless, the relations between bankfull channel geometry and flow discharge still work as
A plot of the channel response processes is shown in Figure 8, which shows a theoretical situation where the change of water occurs by skipping from one value to another and the channel geometry parameters ( or ) change synchronously. Nevertheless, a process connects the past bankfull values to the current ones. In this process, the value of or is neither the past bankfull value nor the current one. It is also a bankfull value including two impact input processes: the past flow process and current one. Here, we sought a rule that describes the transition between the two stable states.
Equation (11) is differentiated at time leading to The final form is
The computational model in a reach of the LYR is shown in Figure 9. There, inlet and outlet data appear from the entrance and exit stations of the reach. CS0 and are the sediment transport capacity (kg/m3) of the two stations; is the flow discharge of the reach (m3/s), is the grain size of the bed material load in the 65% range (mm), and is the length of the reach (m).
4. Results and Discussion
4.1. Simulation Results
From Jinan Station to the Lijin Station, the model uses the daily water and sediment discharge data from the 1964 to 2000: the bed slope , , and is observed and changed with time every year as shown in Table 3. The simulation results are shown in Figure 10.
Figure 10 shows that two types of year cause nonnegligible errors. The first type is the years with low flow discharge and high sediment load, such as 1973 and 1997. In 1973, the annual sediment load at Lijin Station was 1.2 billion tons per year or 50% more than the mean value in the past 50 years. The annual flow volume was 28.15 billion m3 per year or 10% less than the mean value in the past 50 years. In 1997, the annual sediment load at Lijin Station was 0.016 billion tons per year or 98% less than the mean value in the past 50 years. The annual flow volume was 1.86 billion m3 per year, with 227 days in the zero-flow state. In these years, deposition along the river occurred more heavily at the outlet. This was the main reason why the observed bankfull width was reduced and the calculated value was not. The second type was years with overflood events, such as what occurred in 1967 and 1988. In 1967, the maximum daily discharge was 8510 m3/s and the observed bankfull discharge (after the flood season) was 7500 m3/s. In 1988, the maximum daily discharge was 5220 m3/s and the observed bankfull discharge (after the flood season) was 5000 m3/s. However, overbank flood events were less important than the larger sediment load and smaller flow volume from the error analysis at Lijin Station. One explanation for the smaller impact might be the fact that the overbank flood volume was not much larger than the bankfull flood at the river outlet.
The daily water and sediment discharge data for Xiaolangdi and Huayuankou Stations from 1976 to 1997 was used in the model, together with the daily water discharge data of Heishiguan and Wuzhi Stations. The simulation results are shown in Figure 11.
The calculated errors at Huayuankou station were more obvious than those at Lijin station. Besides the first type, which included 1987 and 1992, the second type, which included , , 1988, and 1996, had more impact on the errors. may be taken as an example: the max daily discharge was 7760, 13,400, and 7580 m3/s, respectively, and the bankfull discharge (after the flood season) was 5320, 6000, and 6800 m3/s. Three years of overbank flood events caused obvious river width expansion, which the methods in this paper cannot reflect.
4.2. Accuracy Impact Factors for Simulations
4.2.1. Lateral Inflow
Two tributaries in the Xiaolangdi-Huayuankou reach and the Yiluo and Qin Rivers affected the flow and sediment discharge in the reach. Given the lack of daily sediment concentration data, a simplified study is required.
According to the statistical data from 1960 to 1996, the annual average sediment discharge of the Yiluo River (Heishiguan station) is , accounting for 0.9% of that in the main stream of the LYR. The annual average sediment discharge of the Qin River (Wuzhi station) is , accounting for 0.4% of that in the main stream of the LYR. Therefore, these two rivers are considered to be clear in the calculation for this section. This simplification might affect the accuracy in some years when the sediment load cannot be ignored.
4.2.2. Overbank Flood
Differences exist in the scour and silting mechanisms between overbank and normal floods. Scouring occurs when the water level rises, while silting occurs when the water level falls. Floodplain silting with channel scouring is common during overbank flood periods in the LYR.
For overbank floods with a low sediment load, for example, the flood that occurred from July 31 to August 8, 1982, the peak discharge at Huayuankou station was 15,300 m3/s and the average sediment concentration was 67 kg/m3. This caused 0.217 billion tons of sediment to be deposited in the floodplain and 0.15 billion tons to be scoured from the main channel in the Huayuankou-Aishan reach.
For overbank floods with high sediment load, for example, the flood that occurred in 1977, the peak discharge at Huayuankou station was 10,800 m3/s and the average sediment concentration was 437 kg/m3. This caused the amount of sediment deposited in the main channel to be three times that deposited in the floodplain in the Huayuankou-Aishan reach.
Consequently, the methods used in this paper are not suitable for overbank situations. The key issue is that the distribution of resistance in the floodplain and main channel is unknown. A simplified statistical method is used to deal with overbank flood events. For flows with a high sediment concentration, sediment is deposited all along the river, and a distributing coefficient reflecting the deposition ratio between the floodplain and main channel is calibrated against all years. Here, is the deposition volume in the main channel (m3) and is in the floodplain.
For flows with a high sediment concentration, sediment is deposited all along the river , while for flows with a low sediment concentration, the main channel is scoured, the floodplain is deposited, and . According to statistical data, the criterion for high sediment loads is 300 kg/m3.
4.2.3. Channel Shape of a Uniform Reach
Here, the chosen river reach is assumed to be uniform so that the depositing or eroding sediment volume can be distributed along the reach evenly. This assumption is much more acceptable in the Jinan-Lijin reach than in the Xiaolangdi-Huayuankou reach, because there is a valley channel beyond Xiaolangdi Station that affects the sediment distribution along the Xiaolangdi-Huayuankou reach. This might be an important impact factor for the accuracy of the model applied to Huayuankou station.
(1)Based on the equations for stable alluvial rivers, namely, the Exner equation, flow continuity, the momentum, and sediment concentration conservation equations and assuming that the riverbed slope is constant on middle timescales (100-101 years), we found stationary geomorphic relations for time-varying fluvial processes. We obtained a time-dependent relation for the channel geometry given nonstationary alluvial processes.(2)The final solutions from the equations for stable alluvial rivers simply reflect the adjusted direction of channel geometry. The adjusted magnitude of the hydraulic geometry each day is determined by the daily input of depositing or eroding sediment volumes. The direction and magnitude give the final channel shape.(3)Using daily flow and sediment data to determine the adjustment direction and magnitude, the variation in the bankfull characteristics for Lijin and Huayuankou Stations in the LYR were calculated. The final results agreed closely with the observed values. This implies that the assumptions are acceptable and the method is applicable for the reservoir and channel regulation in the LYR.(4)Compared with the rate-law and hydrodynamic methods, the proposed method enables multiple dependent variables and the input of detailed hydrological data for long-term simulations. However, because of those assumptions for simplification, the model we provided also has some limitations. Firstly, it cannot be applied to a dramatic time-varying fluvial procedure like an overbank flow or a hyperconcentration flood, only fitting for a quasisteady fluvial process. Secondly, the time-step of one day is still large for a physics-based model which smoothes the deviation of inner-daily flood process. Thirdly, the key equation in this model is derived by the experiential relations between channel geometry and hydraulic relations, which is not verified in theory. These limitations are the main sources of the model error and will be our next work in the future.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
This study was supported financially by the China Ministry of Science and Technology through the Key Project in the National Science & Technology Pillar Program in the Twelfth Five-year Plan Period of China (Grant no. 2012BAB02B02). It was supported by the Open Program Fund of MWR Key Laboratory of Yellow River Sedimentation of the Yellow River Conservancy Commission.
G. Lacey, “Flow in alluvial channels with sandy mobile beds,” Proceedings of the Institute of Civil Engineers, vol. 11, pp. 145–164, 1958.View at: Google Scholar
D. B. Simons and M. L. Albertson, “Uniform water conveyance channels in alluvial materials,” Journal of the Hydraulics Division, vol. 86, pp. 33–71, 1960.View at: Google Scholar
W. L. Graf, “The rate law in fluvial geomorphology,” American Journal of Science, vol. 277, no. 2, pp. 178–191, 1977.View at: Google Scholar
G. P. Williams and M. G. Wolman, Downstream Effects of Dams on Alluvial Rivers, Geological Survey Professional Paper no. 1286, U.S. Government Printing Office, Washington, DC, USA, 1984.
A. Simon, “A model of channel response in disturbed alluvial channels,” Earth Surface Processes & Landforms, vol. 14, no. 1, pp. 11–26, 1989.View at: Google Scholar
C. E. Thorn and M. R. Welford, “The equilibrium concept in geomorphology,” Annals of the Association of American Geographers, vol. 84, no. 4, pp. 666–696, 1994.View at: Google Scholar
D. Knighton, Fluvial Forms and Processes: A New Perspective, Arnold, London, UK, 1998.
M. P. Mosley, “Semi-determinate hydraulic geometry of river channels, South Island, New Zealand,” Earth Surface Processes & Landforms, vol. 6, no. 2, pp. 127–137, 1981.View at: Google Scholar
C. Gippel, “Changes in stream channel morphology at tributary junctions, Lower Hunter Valley, New South Wales,” Australian Geographical Studies, vol. 28, no. 2, pp. 291–307, 1985.View at: Google Scholar
J. E. Pizzuto, “Flow variability and the bankfull depth of sand-bed streams of the American Midwest,” Earth Surface Processes & Landforms, vol. 11, no. 4, pp. 441–450, 1986.View at: Google Scholar
T. J. Pietsch and G. C. Nanson, “Bankfull hydraulic geometry; the role of in-channel vegetation and downstream declining discharges in the anabranching and distributary channels of the Gwydir distributive fluvial system, southeastern Australia,” Geomorphology, vol. 129, no. 1-2, pp. 152–165, 2011.View at: Publisher Site | Google Scholar
A. D. Howard, “Equilibrium and time scales in geomorphology: application to sand-bed alluvial streams,” Earth Surface Processes & Landforms, vol. 7, no. 4, pp. 303–325, 1982.View at: Google Scholar
B. Yu and M. G. Wolman, “Some dynamic aspects of river geometry,” Water Resources Research, vol. 23, no. 3, pp. 501–509, 1987.View at: Google Scholar
E. Pasche and G. Rouvé, “Overbank flow with vegetatively roughened flood plains,” Journal of Hydraulic Engineering, vol. 111, no. 9, pp. 1262–1278, 1985.View at: Google Scholar
S. E. Darby and C. R. Thorne, “Predicting stage-discharge curves in channels with bank vegetation,” Journal of Hydraulic Engineering, vol. 122, no. 10, pp. 583–586, 1996.View at: Google Scholar
E. Mosselman, “Morphological modeling of rivers with erodible banks,” Hydrological Processes, vol. 12, no. 8, pp. 1357–1370, 1998.View at: Google Scholar
Z.-H. Sun, Y.-T. Li, Y. Huang, and K.-C. Gao, “Fluvial process of sandbars and shoals in branching channels of the middle Yangtze River,” Journal of Hydraulic Engineering, vol. 42, no. 12, pp. 1398–1406, 2011.View at: Google Scholar
K. Asahi, Y. Shimizu, J. Nelson, and G. Parker, “Numerical simulation of river meandering with self-evolving banks,” Journal of Geophysical Research: Earth Surface, vol. 118, no. 4, pp. 2208–2229, 2013.View at: Google Scholar
D. S. Wilks, Statistical Methods in the Atmospheric Sciences, Academic Press, Oxford, UK, 3rd edition, 2011.
J. M. Hooke, “River channel adjustment to meander cutoffs on the River Bollin and River Dane, northwest England,” Geomorphology, vol. 14, no. 3, pp. 235–253, 1995.View at: Google Scholar
A. Simon and C. R. Thorne, “Channel adjustment of an unstable coarse-grained stream: opposing trends of boundary and critical shear stress, and the applicability of extremal hypotheses,” Earth Surface Processes & Landforms, vol. 21, no. 2, pp. 155–180, 1996.View at: Google Scholar
G. A. Richard, Quantification and prediction of lateral channel adjustments downstream from Cochiti Dam [Ph.D. thesis], Colorado State University, Rio Grande, NM, USA, 2001.
C. H. Hu, Complex Channel Response and Variation of 13 Flow and Sediment Processes of the Yellow River, Science Press, Beijing, China, 1st edition, 2005.
E. D. Andrews, “Effective and bankfull discharges of streams in the Yampa River basin, Colorado and Wyoming,” Journal of Hydrology, vol. 46, no. 3-4, pp. 311–330, 1980.View at: Google Scholar