To examine how endogenous reserves may influence avian life history, it is often necessary to quantify carcass composition. However, proximate analyses are expensive, time-consuming and difficult to perform under field conditions. Consequently, carcass composition is often estimated from easily measured data. We evaluate methods of estimating carcass composition of the common eider duck Somateria mollissima. We measured, dissected and completed proximate analyses of 92 eiders. Predictive models were derived using multiple regressions of 70 birds, while the remaining 22 were used as an independent test of the models. Each model's accuracy was evaluated by comparing estimates against known values of protein and lipids, using root mean square error (RMSE). Abdominal and leg fat pad mass were highly correlated with total lipid (r = 0.92), and body mass was highly correlated with total protein (r = 0.80). Models that used body mass, fat depots and/or muscle group data were the most accurate (lipids adjusted R2 = 0.93, RMSE = 14.60; protein adjusted R2 = 0.74, RMSE = 11.14). By using these equations it is possible, using dissection data, to accurately estimate carcass composition of eiders. If dissection data are not available, one can still estimate carcass composition using equations that require only morphometrics although in our lipid analysis such equations had relatively low accuracy (lipids adjusted R2 = 0.54, RMSE = 32.74).
Knowledge of body condition can provide insight into avian life history. Birds commonly store fat to ensure that their future energetic needs are met (Griminger 1986), and protein can also be catabolised in times of energy or nutrient shortage (King & Murphy 1985). The amount of reserves that a bird stores is a trade-off between starvation risk and predation risk (Lima 1986). Reserve levels are also influenced by food availability (Joyner et al. 1984, Oosterhuis & van Dijk 2002), environmental conditions (Whyte & Bolen 1984a, Lovvorn 1994), and life history stage (Korschgen 1977, McLandress & Raveling 1981, Parker & Holm 1990).
Body condition is a general term that refers to a bird's ability to meet current and future energy needs (Owen & Cook 1977). It is often inferred from body mass, condition indices or mass of muscles and/or lipid depots (e.g. breast muscle, abdominal fat). To measure carcass composition, it is standard practice to extract lipids from dry tissue using a Soxhlet apparatus, with a solvent such as petroleum ether, and then to combust lean dry tissue to determine the protein content of a carcass (Reynolds & Kunz 2001). Such proximate analyses are expensive, time-consuming and difficult to perform under most field conditions. They also require that birds be sacrificed, and this prevents repeated measurements of the same bird over time and often limits sample size. Consequently, many researchers have attempted to estimate carcass composition using more easily obtained data.
A common technique for indexing size of nutrient reserves is to derive predictive models from regressions of carcass lipids and protein against measurements of body mass, abdominal fat pad mass, or breast muscle mass (Thomas et al. 1983, Piersma 1984, Miller 1989, Boos et al. 2000). Sometimes ratios of body mass over a measure of body size are used as predictor variables in regressions or as indices of body condition on their own (e.g. Wishart 1979, Whyte & Bolen 1984b). However, using ratios as a means of expressing body condition should be done only when the two variables are related isometrically, and they rarely are (Packard & Boardman 1999).
Although several studies have investigated ways to estimate carcass composition of waterfowl, there are no published studies examining any of the sea ducks (Mergini tribe). Sea ducks may have unique physiologies because they spend most of their lives in the marine environment, with several species inhabiting polar environments. The body mass of common eiders Somateria mollissima is known to fluctuate greatly throughout the year (Gorman & Milne 1971, Milne 1976), and it is suspected that eiders in poor condition during winter and spring refrain from breeding in summer (Coulson 1984, Oosterhuis & van Dijk 2002). However, little is known about their body condition in winter (Gorman & Milne 1971, Korschgen 1977). In this study, we evaluate several methods for estimating lipid and protein content of northern common eiders S. m. borealis collected during the non-breeding season in Greenland. Specifically, we evaluate models taken from the literature and models we derived from various external morphological measurements and dissection data. Finally, we determine the minimum analysis necessary to rigorously quantify endogenous reserve levels of sea ducks.
Material and methods
Inuit hunters and fishermen collected 748 common eiders as part of a subsistence harvest from the waters of the southwest coast of Greenland. The birds were either shot or retrieved drowned from fishnets where they were caught unintentionally. From this sample, 92 birds were selected for complete carcass analysis of lipid and protein. They were selected to ensure that a wide range of endogenous reserve levels, and sex and age classes were analysed (18 first-year females; 28 after-first-year females; 15 first-year males; 31 after-first-year males). We also ensured that these 92 birds had not sustained any damage during collection that would prevent us from gathering all the data.
Ducks were aged by plumage characteristics (Baker 1993) and length of the bursa of Fabricius (Mather & Esler 1999), and sexed by syrinx morphology (Beer 1963). The total length of the head and bill and the length of the tarso-metatarsus bone were measured to the nearest 0.1 mm using calipers. Flattened wing length was measured with a wing board to the nearest mm. Birds were weighed to the nearest 1 g with an electronic scale.
Eiders were plucked, except for the head region, and reweighed. Breast muscles (both pectoralis major and minor), leg muscles (all muscles that originate or insert in the femur or tibiotarsus bones), and leg fat pad were dissected from the right side of each bird and weighed. The gizzard, heart and abdominal fat pad were extracted and weighed. The digestive system was removed, emptied of its contents and reweighed. For exact methods of dissection see Jamieson et al. (submitted). After dissection, all excised tissues were returned to the carcass and refrozen.
Dissected birds were shipped frozen to the Avian Energetics Laboratory at the University of Western Ontario, London, ON. There, each bird was thawed and cut into ~ 2-cm3 pieces and subsequently dried to a constant weight in an oven at 80°C. Constant weight was achieved if the bird lost < 1 g per day in the oven (~ 1 week). After dry carcass mass was measured, each bird was homogenised using a hand meat grinder. From each homogenate, a 10-g sample was placed in a single thickness cellulose thimble and placed in a Soxhlet extractor for 16 hours. Petroleum ether was used as the solvent because it extracts mostly storage lipids rather than structural lipids (e.g. phospholipids) or non-lipid compounds (Dobush et al. 1985). After 16 hours of extraction, each sample was dried at 80°C in an oven, weighed and returned to the extractor for another four hours. This process was repeated until a weight change of < 0.01 g was obtained. The lean dry homogenate was placed in a crucible and heated in a muffle furnace for 16 hours at 620°C, to burn off the non-mineral content of the sample. Dry ash weight was recorded. Total carcass water, lipid and protein were calculated as:
The assumption of normality was checked using the Kolmogorov-Smirnov test (Zar 1999). Homoscedasticity was examined using Levene's test (SYSTAT Software Inc. 2002). Linearity was inspected using residual plots of total carcass lipids and protein, and of body mass (Zar 1999). Total body mass was used rather than body mass minus digestive system contents because it was most relevant to field studies where the entire bird is weighed.
Approximately 80% of the birds were randomly selected from each age/sex class and this source group was used to derive models that estimated total lipids and protein (Tabachnick & Fidell 2001). The remaining 20% (test group) provided an independent means of evaluating the performance of predictive models through cross-validation. We ran both unadjusted and Bonferroni adjusted t-tests to ensure that there were no significant differences between carcasses of test and source groups. We compared their carcass composition, morphometrics and masses of muscles groups, organs and lipid depots.
We also considered using a principle component analysis to index body size and ran some preliminary analyses with PC1 as an independent variable. However, we were concerned that this approach might obscure some morphological variation that might otherwise be significant in the regression analyses. Furthermore, initial analyses using PC1 to index body size did not increase the accuracy of the resulting predictive equations. Therefore, we felt it was appropriate to use individual morphometrics.
To estimate total carcass lipid and protein, we applied stepwise backwards multiple regression using combinations of 12 independent variables (Table 1). At each step the variable with the highest P-value was removed from the analysis and the regression was rerun. This process was repeated until only those variables that contributed significantly to explaining variation in the dependent variable remained (P < 0.05). Regressions were run using: a) all independent variables, b) only independent variables that could be measured on live birds, and c) independent variables that were strongly correlated with the dependent variable (Pearson correlation coefficients > 0.60). These variables were run both individually and together as a group.
Independent variables included in various combinations in stepwise backwards multiple regressions used to estimate total carcass lipid and total carcass protein of non-breeding common eiders.
There are several predictive models involving water content that have been applied to birds (Child & Marshall 1970, Campbell & Leatherland 1980, Briggs & Thornton 1988). Percent water content is often correlated negatively with lipid content of a carcass (e.g. Briggs & Thornton 1988, Miller 1989). Therefore, we also ran a regression that included percent water to estimate total lipid content. We applied two previous published methods, the Child-Marshall model (1970) to estimate lipid content and the Campbell-Leatherland model (1980) to estimate lipid and protein content. The Child-Marshall model uses the ratio of carcass water to fat-free weight, whereas the Campbell-Leatherland model uses: 1) the ratio of carcass water to carcass protein, and 2) the ratio of total carcass lipid plus total carcass protein to dry carcass weight (Child & Marshall 1970, Campbell & Leather-land 1980).
All regressions originally included dummy variables for age (1 = first-year; 0 = after-first-year) and sex (0 = female; 1 = male) because during preliminary analyses we found that there were significant differences between the slopes and intercepts of the age and sex classes when body mass or masses of individual fat pads were regressed against carcass lipid.
Traditional condition indices (e.g. body mass divided by a measure of body size) were not used because the data were not related isometrically; that is, they did not produce straight lines through the origin when plotted against each other (Packard & Boardman 1999). We also refrained from using residuals as our data did not meet all the necessary assumptions (Green 2001).
Derived models were applied to the morphological and dissection data of the test group (N = 22), whose lipid and protein content were known, and total carcass lipid and protein were estimated. The accuracy of each model was evaluated by comparing estimates against known values of protein and lipids from each carcass. This was done by calculating the root mean square error of prediction (RMSE; Olden & Jackson 2000) with a small value indicating a model of high accuracy. Each model was ranked according to its RMSE, which was calculated as follows:
All statistical analyses were carried out using SYSTAT 10.2 (SYSTAT Inc. 2002).
Data for body mass, total carcass lipid and total carcass protein met assumptions of normality (Kolmogorov-Smirnov test: P > 0.05), homoscedasticity (Levene's test: P > 0.05), and linearity. Strong multicollinearity was not found among the independent variables (variance inflation factors ranged within 1.1–5.0).
There were no statistically significant differences between carcass components of the source and test groups (Table 2).
Carcass components of northern common eiders wintering in Greenland. Ducks were randomly divided into a source (N = 70) or test (N = 22) group and compared using both unadjusted and Bonferroni adjusted t-tests. No significant differences were found between the two groups.
The variables of abdominal fat pad mass, leg fat pad mass, and body mass were strongly correlated to total carcass lipids (Pearson correlation coefficients > 0.60; Table 3). The locomotory muscles (breast and leg) and body mass were strongly correlated with total carcass protein (Pearson correlation coefficients > 0.60; see Table 3). Due to the high correlations with the dependent variables, these variables were run in their own regression analyses. Lipid and protein contents were not strongly correlated to any measure of body size.
Pearson correlation coefficients of the relationships between various carcass components and total carcass lipid and protein of non-breeding common eiders.
Seven different models were produced using backwards stepwise multiple regressions to estimate total carcass lipids (Table 4). Model 1, which included body mass, leg fat pad mass, abdominal fat pad mass, and flattened wing length, most accurately estimated total carcass lipids (see Table 4). Previously published methods (Marshall-Child Method, Campbell-Leatherland Method, and percent water content) all scored lower in comparison to models containing fat pad weights, as did models derived using only external measures from live birds (see Table 4).
Models derived using backward stepwise multiple regressions to estimate total carcass lipid (TCL) of common eiders during the non-breeding season. Models were ranked according to their root mean square error of prediction (RMSE; see ‘Model evaluation’). All models originally included the variables sex and age.
Six different models were derived to estimate total carcass protein (Table 5). Model 8 had the highest accuracy for predicting total carcass protein, which included body mass, breast muscle mass, gizzard mass and abdominal fat pad mass (see Table 5). Models using only variables easily measured on living birds, i.e. models 9 and 13, scored higher than models that included only muscle groups (see Table 5). The Campbell-Leatherland Method ranked low when compared to derived models (see Table 5).
Models derived using backward stepwise multiple regressions to estimate total carcass protein (TCP) of common eiders during the non-breeding season. Models were ranked according to their root mean square error of prediction (RMSE; see ‘Model evaluation’). All models originally included the variables sex and age.
The primary goal of our study was to derive models using relatively easily measured data to estimate carcass composition of common eiders. However, we found poor fit between simple external measurements of eiders and actual carcass composition, in particular lipid content. We did, however, find that it is possible to accurately estimate the carcass composition of northern common eiders through dissection using models derived from proximate carcass analyses. Among this sample of common eider ducks, Model 1 most accurately estimated total carcass lipid (Adjusted R2 = 0.93, RMSE = 14.60). Variables in Model 1 that significantly contributed to the equation included body mass, leg fat pad mass, abdominal fat pad mass, and flattened wing length; all parameters easily quantified through dissection. Total carcass protein was most accurately estimated by Model 8 which included the independent variables of body mass, breast muscles mass, gizzard mass and abdominal fat pad mass (Adjusted R2 = 0.74, RMSE = 11.14).
Both models 1 and 8 were derived by entering all 12 measured variables into the regression analysis and removing those that did not contribute significantly to explaining the variance observed in carcass composition. Recall that the goal of our study was to develop models to estimate carcass composition and not to make any inferences about biological relationships between variables. Therefore, it was appropriate to include variables that did not appear obviously related to the dependent variable because they may add significantly to the model for unknown reasons (e.g. abdominal fat pad mass explaining some of the observed variation in total carcass protein).
Despite their strong predictive power, we caution that the specific models developed here should be used only to estimate the composition of common eiders that fall within the range of the data used to derive the equation (Zar 1999). Models derived for one population may not be appropriate for others of the same or similar species (Castro & Myers 1990, Sparling et al. 1992). Therefore, the models developed in our study should be applied to northern common eiders that weigh 1,542–2,602 g (the range of body masses included in this study) and only to non-breeding birds. Eiders undergo large physiological changes during the breeding season, which could influence the accuracy of the models. However, the results of our study have several general implications concerning the methods used to study endogenous reserves of other bird species, for example the use of fat pad mass as an index of body condition.
The abdominal fat pad and leg fat pads from dissected eiders were good indices of body condition for several reasons. First, fat pads consist primarily of lipids that are available as an energy source (i.e. triglycerides) rather than structural lipids (e.g. phospholipids). Second, they can be easily removed without causing much damage to the carcass, and this might increase the participation of hunters offering their take for research purposes. Finally, both abdominal and leg fat pad were highly correlated with total carcass lipid (see Table 3). Abdominal fat pad has been found to be highly related to total carcass lipids in other species of waterfowl (red-billed teal Anas erythrorhyncha: Woodall 1978; Canada geese Branta canadensis and lesser snow geese Chen caerulescens caerulescens: Thomas et al. 1983; ring-necked ducks Aythya collaris: Hohman & Taylor 1986; northern pintails Anas acuta: Miller 1989). However, leg fat pad is often overlooked as a possible index of body condition, even though it was also strongly related to total carcass lipids in our study (see Table 3). The leg fat pad is more easily removed than abdominal fat because it is well-defined and not entwined in the internal mesentery or organs.
Models that incorporated external structural measures and body mass predicted carcass protein more accurately than they did carcass lipid content (RMSE for protein ranged within 13.20–13.27 and adjusted R2 ranged within 0.67–0.71, while RMSE for lipids ranged within 32.74–38.93 and adjusted R2 ranged within 0.54–0.61). This has also been found in other waterfowl studies (maned ducks Chenonetta jubata: Briggs 1989; northern pintails: Miller 1989; mallards Anas platyrhynchos: Boos et al. 2000).
We found that it was not possible to accurately predict total carcass lipids from percent carcass water (Adjusted R2 = 0.48, RMSE = 33.10). Unlike Miller (1989) and Johnson et al. (1985), we found that the Child-Marshall and Campbell-Leatherland models did not accurately predict carcass composition, nor did percent water content. Many of the eiders examined in our study had drowned in fishing nets (N = 52). Even though we could account for water held in feathers, we were not able to account for any water contained in their lungs and air sacs. We speculated that this additional water contributed to the inaccuracy of the models predicting carcass lipids from carcass water among our sample of drowned eiders. However, this does not appear to be the case. Independently, we ran the Child-Marshall and Campbell-Leatherland models on a sample of birds that had been shot, and the new RMSE values were even greater than those calculated from drowned and shot birds combined (58.18 compared to 44.57 for the Child-Marshall Model; 45.67 compared to 31.67 for the Campbell-Leatherland Model).
In conclusion, methods using only measures of structural size and body mass were insufficient to predict lipid content in northern eiders. In contrast, the models derived from measures of fat depots and muscle groups attained through simple dissection gave accurate estimates of both lipid and protein content which would be sufficient for detailed studies of avian energetics. We suggest that these methods should be considered to avoid the time and significant costs associated with proximate carcass analyses (i.e. protein, lipid and ash) of all birds within a sample.
Funding for this project was provided by the Canadian Wildlife Service, the Greenland Institute of Natural Resources, Danish Environmental Protection Agency, and the Atlantic Cooperative Wildlife Ecology Research Network. This is ACWERN Publication No. UNB-39. We thank C. Davison Ankney for his many helpful comments and suggestions, and for the use of his Avian Energetic Laboratory at the University of Western Ontario. The comments of Dan Esler, Graham Forbes, Hannu Pöysä, Greg Robertson and two anonymous reviewers greatly improved earlier versions of this manuscript. Our gratitude is also extended to Joël Bêty, Mark Dionne, Diana Hamilton and Gary Packard for many valuable discussions.
- K. Baker 1993. Identification guide to European non-passerines. British Trust for Ornithology, BTO guide 24, Thetford. 332. pp. Google Scholar
- J. V. Beer 1963. The rapid sexing of downy Anatidae by the structure of the syrinx. Wildfowl 14:160–162. Google Scholar
- M. Boos, T. Zorn, A. Koch, Y. Le Maho, and J-P. Robin . 2000. Determining body fuels of wintering mallards. Comptes Rendus de l' Academie des Sciences Serie III- Sciences de la Vie 323:183–193. Google Scholar
- S. V. Briggs 1989. Morphological prediction of body condition in maned ducks. Australian Wildlife Research 16:605–609. Google Scholar
- S. V. Briggs and S. A. Thornton . 1988. Abdominal fat and percent water as predictors of body fat in adult maned duck, Chenonetta jubata. Australian Wildlife Research 15:231–234. Google Scholar
- R. R. Campbell and J. F. Leatherland . 1980. Estimating body protein and fat from water content in lesser snow geese. Journal of Wildlife Management 44:438–446. Google Scholar
- G. Castro and J. P. Myers . 1990. Validity of predictive equations for total body fat in sanderlings form different non-breeding areas. Condor 92:205–209. Google Scholar
- G. I. Child and S. G. Marshall . 1970. A method of estimating carcass fat and fat-free weights in migrant birds from water content of specimens. Condor 72:116–119. Google Scholar
- J. C. Coulson 1984. The population dynamics of the eider duck Somateria mollissima and evidence of extensive non-breeding by adult ducks. Ibis 126:525–543. Google Scholar
- G. R. Dobush, C. D. Ankney, and D. G. Krementz . 1985. The effect of apparatus, extraction time, and solvent type on lipid extraction of snow geese. Canadian Journal of Zoology 63:1917–1920. Google Scholar
- M. L. Gorman and H. Milne . 1971. Seasonal changes in the adrenal steroid tissue of the common eider Somateria mollissima and its relation to organic metabolism in normal and oil-polluted birds. Ibis 113:218–228. Google Scholar
- A. J. Green 2001. Mass/length residuals: measures of body condition or generators of spurious results? Ecology 82:1473–1483. Google Scholar
- P. Griminger 1986. Lipid Metabolism. In: P. D. Sturkie , editor. (Ed.). Avian physiology. Springer-Verlag. New York. pp. 345–358. Google Scholar
- W. L. Hohman and T. S. Taylor . 1986. Indices of fat and protein for ring-necked ducks. Journal of Wildlife Management 50:209–211. Google Scholar
- S. E. Jamieson, G. J. Robertson, K. Falk, F. R. Merkel, and H. G. Gilchrist . submitted. Waterfowl dissection manual. Canadian Wildlife Service Report. Google Scholar
- D. H. Johnson, G. L. Krapu, K. J. Reinecke, and D. G. Jorde . 1985. An evaluation of condition indices for birds. Journal of Wildlife Management 49:569–575. Google Scholar
- D. E. Joyner, R. D. Arthur, and B. N. Jacobson . 1984. Winter weight dynamics, grain consumption and reproductive potential in Canada geese. Condor 86:275–280. Google Scholar
- J. R. King and M. E. Murphy . 1985. Periods of nutritional stress in the annual cycle of endotherms: fact or fiction? American Zoologist 25:955–964. Google Scholar
- C. E. Korschgen 1977. Breeding stress of female eiders in Maine. Journal of Wildlife Management 41:360–373. Google Scholar
- S. L. Lima 1986. Predation risk and unpredictable feeding conditions: determinants of body mass in birds. Ecology 67:377–385. Google Scholar
- J. R. Lovvorn 1994. Nutrient reserves, probability of cold spells and the question of reserve regulation in wintering canvasbacks. Journal of Animal Ecology 63:11–23. Google Scholar
- D. D. Mather and D. Esler . 1999. Evaluation of bursal depth as an indicator of age class of harlequin ducks. Journal of Field Ornithology 70:200–205. Google Scholar
- M. R. McLandress and D. G. Raveling . 1981. Changes in diet and body composition of Canada geese before spring migration. Auk 98:65–79. Google Scholar
- M. R. Miller 1989. Estimating carcass fat and protein in northern pintails during the nonbreeding season. Journal of Wildlife Management 53:123–129. Google Scholar
- H. Milne 1976. Body weights and carcass composition of the common eider. Wildfowl 27:115–122. Google Scholar
- J. D. Olden and D. A. Jackson . 2000. Torturing data for the sake of generality: how valid are our regression models? Ecoscience 7:501–510. Google Scholar
- R. Oosterhuis and K. van Dijk . 2002. Effect of food shortage on the reproductive output of common eiders Somateria mollissima breeding at Griend (Wadden Sea). Atlantic Sea-birds 4:29–37. Google Scholar
- M. Owen and W. A. Cook . 1977. Variation in body weight, wing length and condition of mallards Anas platyrhynchos and their relationship to environmental changes. Journal of Zoology (London) 183:377–395. Google Scholar
- G. C. Packard and T. J. Boardman . 1999. The use of percentages and size-specific indices to normalize physiological data for variation in body size: wasted time, wasted effort? Comparative Biochemistry and Physiology 122:37–44. Google Scholar
- H. Parker and H. Holm . 1990. Patterns of nutrient and energy expenditure in female common eiders nesting in the high Arctic. Auk 107:660–668. Google Scholar
- T. Piersma 1984. Estimating energy reserves of great crested grebes Podiceps cristatus on the basis of body dimensions. Ardea 72:119–126. Google Scholar
- D. S. Reynolds and T. H. Kunz . 2001. Standard methods for destructive body composition analysis. In: J. R. Speakman , editor. (Ed.). Body composition analysis of animals: a handbook of non-destructive methods. Cambridge University Press. U.K. pp. 39–55. Google Scholar
- D. W. Sparling, J. A. Barzen, J. R. Lovvorn, and J. R. Serie . 1992. An evaluation of regression methods to estimate nutritional condition of canvasbacks and other water birds. U.S. Department of the Interior Fish and Wildlife Service. Washington, D.C., USA. 11. pp. Google Scholar
- SYSTAT Software Inc. 2002. SYSTAT 10.2 Statistics I. SYSTAT Software Inc. Richmond California. 686. pp. Google Scholar
- B. G. Tabachnick and L. S. Fidell . 2001. Using Multivariate Statistics. Allyn and Bacon. Toronto, Ontario. 966. pp. Google Scholar
- V. G. Thomas, S. K. Mainguy, and J. P. Prevett . 1983. Predicting fat content of geese from abdominal fat weight. Journal of Wildlife Management 47:1115–1119. Google Scholar
- R. J. Whyte and E. G. Bolen . 1984a. Impacts of winter stress on mallard body composition. Condor 86:477–482. Google Scholar
- R. J. Whyte and E. G. Bolen . 1984b. Variation in winter fat depots and condition indices of mallards. Journal of Wildlife Management 48:1370–1373. Google Scholar
- R. A. Wishart 1979. Indices of structural size and condition of American wigeon (Anas americana). Canadian Journal of Zoology 57:2369–2374. Google Scholar
- P. F. Woodall 1978. Omental fat: a condition index for red-billed teal. Journal of Wildlife Management 42:188–190. Google Scholar
- J. H. Zar 1999. Biostatistical Analysis. Prentice Hall. Upper Saddle River, New Jersey. 663. pp. Google Scholar