\documentclass[12pt]{article}

\usepackage{graphicx}
\usepackage{latexsym}
\usepackage{caption2}
\usepackage{flafter}
\usepackage{pdflatex}

\oddsidemargin 0.3cm

\textwidth 15.5cm

\textheight 23cm

\pagestyle{myheadings}
\markboth{C. J. Meli\'an, and J. Bascompte}{C. J. Meli\'an, and J. Bascompte}

\title{Food web cohesion}
\author{Carlos J. Meli\'an,$^{} \footnote{Corresponding author;
phone:+34 954 232340; fax:+34 954 621125; e-mail:
cmelian@ebd.csic.es}$\hspace{0.05 in} and  Jordi Bascompte
\\
\\Integrative Ecology Group
\\Estaci\'on Biol\'ogica de Do\~nana, CSIC
\\Apdo. 1056, E-41080 Sevilla, Spain.}
\\
\\

\begin{document}



\date{}
\maketitle
\baselineskip=8.75 mm

\vspace{1 in}
\centerline{\large Running head: Food web cohesion.}

  
\newpage


\centerline{ABSTRACT}

Both dynamic and topologic approaches in food webs have shown
how structure alters conditions for stability. However, while most studies
concerning the structure of food webs have shown a non-random pattern, it
still remains unclear how this structure is related to compartmentalization
and to responses to perturbations.
Here we build a bridge between connectance, food web structure, and compartmentalization by studying how links are distributed within and
between subwebs. A $k$-subweb is defined as a subset of species which are connected to at least $k$ species from the same subset.
We study the subweb $k$-frequency distribution (i.e. the number of
$k$-subwebs in each food web). This distribution is highly skewed, decaying
in all cases as a power law. The most dense subweb has the most interactions
despite containing a small number of species, and shows connectivity values
independent of species richness. The removal of the most dense subweb implies
multiple fragmentation. Our results show a cohesive organization, that is, a
high number of small subwebs highly connected among them through the most dense subweb. We discuss the implications of this organization in relation to different types of disturbances.

\vspace{0.2 in}

\centerline{Key words: Connectance. Food web structure. Compartmentalization.}
\centerline{Subweb. Null model. Cohesion}

\end{abstract}

\vspace{0.3 in}


\newpage


\section{ Introduction}

The structure of food webs is an important
property for understanding dynamic (May 1972, DeAngelis 1975, Pimm 1979,
Lawlor 1980) and topologic (Pimm 1982) stability. Both theoretical and empirical approximations have shown food web structure
represented by guilds (Root 1967), blocks and modules (May 1972), cliques and
dominant cliques (Cohen 1978, Yodzis 1982), compartments (Pimm 1979), subwebs
(Paine 1980), block submatrices (Critchlow and Stearns 1982), and simplicial
complex (Sugihara 1983). From these studies it is well known that food webs are
not randomly assembled. However, it still
remains unclear how the non-random structure of food webs is related to
compartmentalization and its topologic and dynamic implications for
stability following perturbations (Pimm and Lawton 1980, Polis
1991, Strong 1992, Raffaelli and Hall 1992, Solow et al. 1999). This is specially relevant after
studies showing a much larger complexity of food webs than previously expected
(Polis 1991, Strong 1992, Hall and Raffaelli 1993, Polis and Strong 1996).

Current studies show groups of species more connected
than they are with other groups of species (Solow and Beet 1998, Montoya and
Sol\'e 2002). However, these studies do not
make explicit reference to the number of modules and their heterogeneity (see
Ravasz et al. 2002). Here, we build
 a bridge between connectance, food web structure, and
compartmentalization by studying how links are distributed within and between
subsets of species in twelve highly resolved food webs. 

Specifically, we address the following questions:
(1) How are subwebs structured within highly resolved food webs?
(2) What is the relation between food web structure and compartmentalization?
(3) What are the implications of subweb structure for responses to perturbations? 
In order to answer these questions we develop an operative
definition of subweb.

\newpage
\begin{document}
\includepdf{meliancoh.pdf}
\end{document}


\section{Measures of food web structural organization}

\subsection{$k$-Subweb}

A $k$-subweb is defined here as a subset of species
which are connected to at least $k$ prey and/or predators from the same
subset. A $k$-subweb has the following features: 

\begin{enumerate}
\item Subwebs are defined using only information on the presence/absence of
interactions. 

\item Each species belongs only to one subweb, the subset where
each species has the highest value of $k$. 

\item Each subweb contains species from different trophic levels.
\end{enumerate}

Fig. 1 makes explicit this concept. As noted, different subwebs with
the same $k$-value are disjointed in the
web. The sum of the total number of disjointed subwebs with at least $k$
interactions represents the frequency of $k$-subwebs. If we denote by $S_T$ and
$S_k$ the total number of subwebs and the number of $k$-subwebs respectively,
the subweb $k$-frequency distribution is thus $p(S_{k})=S_k/S_T$ (represented
as the cumulative distribution, $P(S_{k})$ throughout this paper).

\subsection{The most dense subweb}

The most dense subweb is the subset of connected species with
the largest number of interactions per species (white circles in Fig. 1 and
Fig. 2). In order to get a measure of cohesion we calculate
and compare connectance for the twelve food webs studied
here (see Table 1). If real food webs are cohesive, we will find a value of connectance of the
most dense subweb significantly larger than both global connectance and the
connectance of the most dense subweb for a series of food webs models. Global connectance is defined as:

\begin{displaymath}
C = \frac{L}{S^2},\eqno(1)
\end{displaymath}
where $L$ is the number of links in the web and $S^2$ is the maximum
number of possible links, including cannibalism and mutual
predation (Martinez 1991). Similarly, we can define the connectance of the
densest subweb ($C_d$) as:

\begin{displaymath}
C_d = \frac{L_d}{S_d^2},\eqno(2)
\end{displaymath}
where $L_d$ is the number of interactions within the most dense subweb,
and $S_d^2$ is the maximum number of possible interactions within
the most dense subweb.

\subsection{Null models of food web structure}

 Can the most dense subweb observed in food webs be reproduced by
 models with different levels of complexity? To answer this question, five food web models were
  tested. We generated $50$ replicates of each model with the same number of
  species and global connectance as the real food webs. Our statistic ($P$) is
  the probability that a random replicate has a $C_d$ value equal
  or higher than the observed value (Manly 1998).   

 In the first model, the basic null model, any link among species occurs with
 the same probability, equal to the global connectance ($C$) of the
 empirical web (Cohen 1978). The
 second model (Cohen et al.'s 1990 cascade model), assigns each species a random value drawn uniformly from the interval [0,1] and each species has the probability $P=2CS/(S-1)$ of
 consuming only species with values less than its own. The third model is the niche
 model by Williams and Martinez (2000). This model assigns a randomly drawn
``niche value" to each species, similarly to the cascade model. Species
 are then constrained to consume all prey species within one range of values
 whose randomly chosen center is less than the consumer's niche
 value. In the preferential attachment model (Barab\'asi and
 Albert 1999), the probability that a new species will be connected to a
 previous species is proportional to the connectivity of the later (both for resources and predators ($j$) of each new
 species), so that $P(k_j)= \frac{k_j}{\sum_i(k_i)}$. Finally, the local rewiring algorithm randomizes the empirical data yet strictly conserves ingoing
 and outgoing links (Connor and Simberloff
 1979; Gotelli 2001). In this algorithm, a pair of directed links $A-B$ and
 $C-D$ are randomly selected. They are rewired in such a way that $A$ becomes
 connected to $D$, and $C$ to $B$, provided that none of these links already
 existed in the network, in which case the rewiring stops, and a new
 pair of edges is selected (Maslov and Sneppen 2002). A library of codes in
 $Matlab$ to generate these matrices is available from the authors upon
 request.

\section{Results}

For the five largest food webs, we calculated the subweb $k$-frequency
distribution. The distribution was found to be strongly skewed with the best
fit following a power law in all webs (see cumulative distribution in Fig. 2).
The average $\pm SD$ of the exponent ($\gamma$) for the five food webs was
$-1.34 \pm 0.57$. This means that subwebs show an extreme heterogeneity, with
most subwebs with a small number of interactions per species and a unique most dense subweb. 

In Silwood Park (Fig. 2a), species belonging to the most dense part ($9\%$ of species in the web) embody $70\%$ of the 
interactions ($26\%$ among the species of the most dense subweb and $44\%$
among the subweb-species and the rest of the web). In Ythan Estuary
(Fig. 2b), the most dense subweb ($21\%$ of species in the web)
holds $74\%$ of all the links in the web ($30\%$ among the species of the most
dense subweb and $44\%$ among the species of the most dense subweb and the
rest of the web). The fraction of interactions in the most dense subweb of
El Verde (Fig. 2c), Little Lake Rock (Fig. 2d),
and Caribbean Coral Reef (Fig. 2e), (with $27\%$, $22\%$, and $31\%$ of
species in the web, respectively) represents $78\%$, $77\%$, and $89\%$ of the
total interactions, respectively ($35\%$, $24\%$, and $33\%$ among the species
of the most dense subweb and $43\%$, $53\%$, and $56\%$ among the most dense
subweb-species and the rest of the web, respectively).

 The average $\pm SD$
percentage of species in the most dense subweb is $22\%\pm 8\%$, and the
average $\pm SD$ percentage of interactions within the
most dense part is $78\%\pm 6\%$. This means that a small number of
species contain the most interactions. The average $\pm SD$ percentage of
species in the most dense subweb in the five null models tested is $86\%\pm
5\%$ for the basic model, $84\%\pm6\%$ for the cascade model, $43\%\pm10\%$ for
the niche model, $37\%\pm15\%$ for the preferential attachment model, and
$28\%\pm13\%$ for the local rewiring algorithm model. 

Table 1 shows global connectance ($C$), the connectance of the most dense
subweb for real data ($C_d$) and the average for each one of the null models
tested (the basic, $\bar{C}_d_b$; cascade, $\bar{C}_d_c$; niche,
$\bar{C}_d_n$; local rewiring algorithm, $\bar{C}_d_l_r_a$ and preferential attachment, $\bar{C}_d_p_a$).
The values of $C_d$ are significantly higher ($P<0.01$) in the twelve food webs for
the basic and cascade model (see Table 1), with the exception of St. Martin in
the cascade model ($0.05<P<0.1$). For the Niche model, three of the most 
resolved food webs (Silwood Park, El Verde and Little Rock Lake), departed
significantly ($P<0.01$) and the rest of the most resolved food webs departed
marginally ($0.05<P<0.1$) (with the exception of Ythan, $P=0.18$). In the
local rewiring algorithm, two
of the most resolved food webs, El Verde and Little Rock Lake, departed
significantly ($P<0.05$ and $P<0.01$, respectively), and the rest of the most resolved food webs
departed marginally ($0.05<P<0.1$), with the exception of the Caribbean food
web ($P>0.1$). Finally, in the preferential attachment model the most resolved
food webs departed marginally ($0.05<P<0.1$) (see Table 1), with the exception
of the Caribean food web ($P>0.1$).  

While $C$, $\bar{C}_d_b$, $\bar{C}_d_c$, $\bar{C}_d_n$ and $\bar{C}_d_p_a$
decay as a power-law as the number of species increases
($r^2=0.53, P< 0.01$; $r^2=0.56, P< 0.01$; $r^2=0.6, P< 0.01$; $r^2=0.73, P<
0.01$; $r^2=0.8, P< 0.01$ respectively), $C_d$ is independent of species
richness ($r^2\leq0.16, P\geq0.47$), which  suggests a scale-invariant
property in the structure of food webs (similarly to the empirical data, the
average value of the $C_d$ in the local rewiring algorithm, $\bar{C}_d_l_r_a$ is independent
of species richness, $r^2\leq0.23, P\geq0.24$).

To further confirm the potential cohesion of the most dense subweb, we removed
it and checked whether the remaining web is fragmented, and if so, in how many
pieces. The web becomes fragmented in $54$ parts in Silwood Park, $37$ parts
in Ythan Estuary, $29$ parts in the Caribbean Coral Reef, $7$ parts in El
Verde, and it does not get fragmented in Little Rock Lake. This multiple
fragmentation shows the role of the most dense subweb.

\section{Discussion}

It is well known that (i) connectance has a very narrow range of values
(Warren 1990, 1994, Martinez 1992), and (ii) food webs are not randomly
assembled (Lawlor 1978, Cohen 1978, Pimm 1980, Ulanowicz and Wolff 1991, Solow et al. 1999). However, little is known about how different subweb frequency distributions are compatible with a specific connectance value and its implications for dynamic and topologic stability.

In this paper we have studied the statistical
properties of the structure in subwebs ($k$-frequency distribution) and its
heterogeneous pattern.
 If this pattern were homogeneous, a single macroscopic description such as connectance would
adequately characterize the organization of food webs. But this is not the
case. There is a need to
move beyond descriptions based on mean field properties such as mean
connectance (Cohen 1978, Pimm 1980, Critchlow and Stearns 1982, Yodzis
1982, Sugihara 1983) to consider these other variables characterizing the
structural organization of food webs.

Our results indicate both a high level of structure (with well-defined
$k$-subwebs) and a cohesive organization (the most dense subweb). While connectance is a scale-variant
property (May 1974, Rejm\'anek and Stary 1979, Yodzis 1980, Jordano 1987,
Sugihara et al. 1989, Bersier et al. 1999, Winemiller et al. 2001), the connectance within the most dense
subweb in the twelve food webs studied is not correlated to species
richness. This is in striking contrast to the null models explored with the
exception of the local rewiring algorithm. Although the degree of connectance (see Table 1), the types of
historical and current human disturbances (Baird and Ulanowicz 1993, Raffaelli 1999), as well as other
ecological and geographic factors were different in the food webs explored, a
similar structural organization was found. This confers a remarkable level of
generality to our results.

What type of mechanisms are underlying this cohesive pattern?
As we have shown, food web models with increasing heterogeneity
in links' distribution don't capture (niche model with the exception of
Ythan and Caribbean) or marginally capture (local rewiring algorithm
and preferential attachment with the exception of the Caribbean) the internal structure of the most resolved food
webs. The biological mechanisms explaining the pattern here reported could be
elucidated by comparing the identity and attributes of the species forming the
most dense subweb across different food webs. If the species composing the most dense subweb in
each food web are taxonomically and phylogenetically different, an ecological
explanation should be suggested (Schoener 1989). However, if the species
forming the most dense subweb are phylogenetically related, evolutionary
mechanisms should be proposed (Williams and Martinez 2000). An intermediate
case would be that in which there are phylogenetic
differences but there is correlation with any biological attribute such as
body size (Cohen et al. 2003) or other physiological and behavioral feature
(Kondoh 2003). In this case, intermediate mechanisms should be suggested.

 These results have implications relative to the previously proposed
 hypothesis about the propagation of perturbations (Pimm and Lawton 1980). The presence of a
high number of small subwebs highly connected among them through the most
dense subweb supports a structured view of the reticulate hypothesis. How are these highly structured
and reticulated webs responding to disturbances? On one hand, the significantly larger probability of
interactions between highly connected intermediate species 
may favor the propagation of disturbances (i.e. a contaminant) through the web (Meli\'an and
Bascompte 2002, Williams et al. 2002). On the other hand, this cohesive
structure may decrease the probability of
network fragmentation when species are removed (Albert et al. 2000, Sol\'e and
Montoya 2001, Dunne et al. 2002). Also, the results presented here may be 
relevant to studies trying to address whether the pattern of subweb structure
may affect the likelihood of trophic cascades (Polis 1991, Strong 1992, Berlow
1999, Pace et al. 1999, Yodzis 2000, Shurin et al. 2002). 

\section{Acknowledgments}

We thank Sandy Liebhold, Pedro Jordano, George Sugihara, Louis-F\'elix
Bersier and Miguel A. Fortuna for useful comments on a previous draft, Enric Sala for bringing to
our attention Opitz's work, and Enrique Collado for computer
assistance. Funding was provided by a Grant to J.B. from the Spanish Ministry
of Science and Technology (BOS2000-1366-C02-02) and a Ph.D. Fellowship to C.J.M. (FP2000-6137).


\newpage

\section{Literature Cited}

Albert, R., Jeong, H., and Barab\'asi, A-L. (2000). Error and attack tolerance
of complex networks. Nature (London) {\bf 406}:378-382.

Almunia, J., Basterretxea, G., Aristegui, J., and Ulanowicz,
R. E. (1999). Benthic-pelagic switching in a coastal subtropical lagoon.
Estuarine Coastal and Shelf Science {\bf 49}:363-384.

Barab\'asi, A-L. \& Albert, R. (1999). Emergence of scaling in random
networks. Science {\bf 286}:509-512.

Baird, D. and Ulanowicz, R. E. (1989). The seasonal dynamics of the Chesapeake
Bay ecosystem. Ecological Monographs {\bf 59}:329-364.

Baird, D. and Ulanowicz, R. E. (1993). Comparative study on the trophic
structure, cycling and ecosystem properties of four tidal estuaries. Marine
Ecology Progress Series {\bf 99}:221-237.

Berlow, E. L. (1999). Strong effects of weak interactions in ecological
communities. Nature (London) {\bf 398}:330-334.

Bersier, L. F., Dixon, P., and Sugihara, G. (1999). Scale-Invariant or scale
dependent behaviour of the link density property in food webs: a matter of
sampling effort? American Naturalist {\bf 153}:676-682.

Connor, E. F. \& Simberloff, D. (1979). The assembly of species
communities: chance or competition? Ecology {\bf 60}:1132-1140.

Cohen, J. E. (1978). Food Webs and Niche Space. Princeton University Press.

Cohen, J. E. Briand, F. \& Newman, C. M. (1990). Community Food Webs: Data and
Theory. Springer-Verlag. Berlin. Germany.

Cohen, J. E., Jonsson, T., \& Carpenter, S. R. (2003). Ecological community
description using the food web, species abundance, and body size.
 Proceedings of the National Academy of Sciences (USA) {\bf 100}:1781-1786. 

Critchlow, R. E. and Stearns, S. C. (1982). The structure of food
webs. American Naturalist {\bf 120}:478-499. 

DeAngelis, D. L. (1975). Stability and connectance in food webs
models. Ecology {\bf 56}:238-243.

Dunne, J. A., Williams, R. J. and Martinez, N. D. (2002). Food web structure
and network theory: The role of connectance and size. 
Proceedings of the National Academy of Sciences (USA) {\bf 99}:12917-12922.

Goldwasser, Ll., and Roughgarden, J. (1993). Construction and analysis
of a large caribbean food web. Ecology {\bf 74}:1216-1233.

Gotelli, N. J. (2001). Research frontiers in null model analysis. Global
Ecology \& Biogeography {\bf 10}:337-343.

Hall, S. J., and Raffaelli, D. G. (1993). Food webs: Theory and
reality. Advances in Ecological Research {\bf 24}:187-239.

Huxham, M., Beaney, S., and Raffaelli, D. (1996). Do parasites reduce the
change of triangulation in a real food web? Oikos {\bf 76}:284-300.

Jordano, P. (1987). Patterns of mutualistic interactions in pollination and
seed dispersal: connectance, dependence asymmetries and coevolution. American
Naturalist {\bf 129}:657-677.

Kondoh, M. (2003). Foraging adaptation and the relationship between food web
complexity and stability. Science {\bf 299}:1388-1391.

Lawlor, L. R. (1978). A comment on randomly constructed model
ecosystems. American Naturalist {\bf 112}:445-447.

Lawlor, L. R. (1980). Structure and stability in natural and randomly
constructed competitive communities. American Naturalist {\bf 116}:394-408. 

Manly, B. F. J. (1998). Randomization, Bootstrap and Montecarlo Methods in
Biology. Chapman and Hall. $2^{nd}$ edition. London.

May, R. M. (1972). Will a large complex system be stable? Nature
(London) {\bf 238}:413-414.

May, R. M. (1974). Stability and Complexity in Model Ecosystems. $2^{nd}$
edition. Princeton University Press, Princeton, New Yersey, USA.

Martinez, N. D. (1991). Artifacts or attributes? Effects of resolution
on the Little Rock Lake food web. Ecological Monographs {\bf 61}:367-392.

Martinez, N. D. and Lawton, J. H. (1992). Constant connectance in community
food webs. American Naturalist {\bf 139}:1208-1218.

Martinez, N. D., Hawkins, B. A., Dawah, H. A., and Feifarek, B. P. (1999).
Effects of sampling effort on characterization of food-web structure. Ecology
{\bf 80}:1044-1055.

Maslov, S., and Sneppen, K. (2002). Specificity and stability in topology of
protein networks. Science {\bf 296}:910-913

Meli\'an, C. J., and Bascompte, J. (2002). Complex networks: two ways to be
robust? Ecology Letters {\bf 5}:705-708.

Memmott, J., Martinez, N. D., and Cohen, J. E. (2000) . Predators,
parasitoids and pathogens: species richness, trophic generality and
body sizes in a natural food web. Journal of Animal Ecology {\bf 69}:1-15.

Montoya, J. M., and Sol\' e, R. V. (2002). Small world pattern in food
webs. Journal of Theoretical Biology {\bf 214}:405-412.

Opitz, S. (1996). Trophic interactions in Caribbean coral reefs.
ICLARM Tech Ed.

Pace, M. L., Cole, J. J., Carpenter, S. R., and Kitchell,
J. F. (1999). Trophic cascades revealed in diverse ecosystems. Trends in
Ecology and Evolution {\bf 14}:483-488.

Paine, R. T. (1980). Food webs: linkage, interaction strength and
community infrastructure. Journal of Animal Ecology {\bf 49}:667-685.

Pimm, S. L. (1979). The Structure of food webs. Theoretical Population Biology
{\bf 16}:144-158.

Pimm, S. L. (1980). Properties of food web. Ecology {\bf 61}:219-225.

Pimm, S. L., and Lawton, J. H. (1980). Are food webs divided into compartments?
Journal of Animal Ecology {\bf 49}:879-898.

Pimm, S. L. (1982). Food Webs. Chapman \& Hall. London.

Polis, G. A. (1991). Complex trophic interactions in deserts: an empirical
critique of food web theory. American Naturalist {\bf 138}:123-155.

Polis, G. A. and Strong, D. R. (1996). Food web complexity and community
dynamics. American Naturalist {\bf 147}:813-846.

Ravasz, E., Somera, A. L., Mongru, D. A., Oltvai, Z. N., and Barab\'asi,
A.-L. (2002). Hierarchical organization of modularity in metabolic
networks. Science {\bf 297}:1551-1555.

Raffaelli, D. and Hall, S. (1992). Compartments and predation in an estuarine
food web. Journal of Animal Ecology {\bf 61}:551-560.

Raffaelli, D. (1999). Nutrient enrichment and trophic organisation in an
estuarine food web. Acta Oecologica {\bf 20}:449-461.

Reagan, D. P. and Waide, R. B. (1996). The Food Web of a Tropical Rain
Forest. University of Chicago Press.

Rejm\'anek, M. and Stary, P. (1979). Connectance in real biotic communities and
critical values for stability of model ecosystems. Nature (London) {\bf 280}:311-315.

Root, R. B. (1967). The niche exploitation pattern of the blue-gray
gnatcatcher. Ecological Monographs {\bf 37}:317-350.

Schoener, T. W. (1989). Food webs from the small to the large. Ecology {\bf 70}:1559-1589.

Shurin, J. B., Borer, E. T., Seabloom, E. W., Anderson, K., Blanchette, C. A.,
Broitman, B., et al. (2002). A cross-ecosystem comparison of the strength of
trophic cascades. Ecology Letters {\bf 5}:785-791.

Solow, A. R. and Beet, A. R. (1998). On lumping species in food webs.
Ecology {\bf 79}:2013-2018.

Solow, A. R., Costello, C., and Beet, A. (1999). On an early result on
stability and complexity. American Naturalist {\bf 154}:587-588.

Sol\'e, R. V. and Montoya, J. M. (2001). Complexity and fragility in
ecological networks. Proceedings of the Royal Society of London B {\bf 268}:2039-2045.

Sugihara, G. (1983). Holes in niche space: a derived assembly rule and its
relation to intervality. In current trends in food web theory, D.L. DeAngelis,
W. Post, and G. Sugihara, eds. (Oak Ridge, TN:Report 5983, Oak Ridge National
Laboratory) , pp. 25.35.

Sugihara, G., Schoenly, K., and Trombla, A. (1989). Scale invariance in food
web properties. Science {\bf 245}:48-52.

Strong, D. R. (1992). Are trophic cascades all wet? differentiation and
donor-control in speciose ecosystems. Ecology {\bf 73}:747-754.

Ulanowicz, R. E., and Wolff, W. F. (1991). Ecosystem flow networks: loaded
dice. Mathematical Biosciences {\bf 103}:45-68.

Warren, P. H. (1989). Spatial and temporal variation in the structure of a
freshwater food web. Oikos {\bf 55}:299-311.

Warren, P. H. (1990). Variation in food-web structure: the determinants of
connectance. American Naturalist {\bf 136}:689-700.

Warren, P. H. (1994). Making connections in food webs. Trends in Ecology and
Evolution {\bf 9}:136-141.

Williams, R. J., and Martinez, N. D. (2000). Simple rules yield complex food
webs. Nature (London) {\bf 404}:180-183.

Williams, R. J., Berlow, E. L., Dunne, J. A., Barab\'asi, A-L., and Martinez,
N. D. (2002). Two degrees of separation in complex food webs. Proceedings of
the National Academy of Sciences (USA) {\bf 99}:12913-12916.

Winemiller, K., Pianka, E. R., Vitt, L. J., and Joern, A. (2001). Food web laws
or niche theory? Six independent empirical tests. American Naturalist
{\bf 158}:193-199.

Yodzis, P. (1980). The connectance of real ecosystems. Nature (London)
{\bf 284}:544-545.

Yodzis, P. (1982). The compartmentation of real and assembled
ecosystems. American Naturalist {\bf 120}:551-570.

Yodzis, P. (1998). Local trophodynamics and the interaction of marine mammals
and fisheries in the Benguela ecosystem. Journal of Animal Ecology
{\bf 67}:635-658.

Yozdis, P. (2000). Diffuse effects in food webs. Ecology {\bf 81}:261-266.

\newpage

\section{Table} 

\centerline{\caption{Table 1: Food webs studied and their statistical properties}}
\vspace{0.1 in}
\hspace{-0.6 in}
\begin{tabular}{|c|c|c|c|c|c|c|c|c|c|c|}\hline
{\small$FW$ & \small$S$ & \small\bar{k}\pm SD & \small$C$ & \small$C_d$ &\small$\bar{C}_d_b$ & \small$\bar{C}_d_c$ & \small$\bar{C}_d_n$ & \small$\bar{C}_d_l_r_a$ & \small$\bar{C}_d_p_a$ & \small$Distribution$\\\hline\hline
MAS & 23 & $6 \pm 3$ & 0.13 & 0.26 & 0.16$^{***}$ & 0.17$^{***}$ & 0.32$^{NS}$ & 0.24$^{NS}$ & 0.39$^{NS}$ & $-$\\
\hline
BEN & 29 & $14 \pm 6$ & 0.24 & 0.34 & 0.26$^{***}$ & 0.265$^{***}$ & 0.38$^{NS}$ & 0.32$^{NS}$ & 0.36$^{NS}$ & $-$\\
\hline
COA & 30 & $19 \pm 8$ & 0.32 & 0.47 & 0.34$^{***}$ & 0.35$^{***}$ & 0.47$^{NS}$ & 0.44$^{**}$ & 0.42$^{NS}$ & $-$\\
\hline
CHE & 36 & $5 \pm 3$ & 0.06 & 0.14 & 0.08$^{***}$ & 0.09$^{***}$ & 0.18$^{NS}$ & 0.15$^{NS}$ & 0.33$^{NS}$ & $-$\\
\hline
SKI & 37 & $21 \pm 9$ & 0.27 & 0.51 & 0.285$^{***}$ & 0.29$^{***}$ & 0.4$^{*}$ & 0.43$^{**}$ & 0.39$^{NS}$ & $-$\\
\hline
STM & 44 & $10 \pm 6$ & 0.11 & 0.16 & 0.13$^{***}$ & 0.135$^{*}$ & 0.22$^{NS}$ & 0.17$^{NS}$ & 0.29$^{NS}$ & $-$\\
\hline 
UKG & 75 & $3 \pm 3$ & 0.02 & 0.26 & 0.033$^{***}$ & 0.04$^{***}$ & 0.14$^{*}$ & 0.14$^{*}$ & 0.16$^{*}$ & $-$\\
\hline 
YE & 134 & $9 \pm 10$ & 0.033 & 0.23 & 0.04$^{***}$ & 0.04$^{***}$ & 0.12$^{NS}$ & 0.19$^{*}$ & 0.12$^{*}$ & \small$PL$ $({\gamma}= -1.87)$\\
\hline 
SP & 154 & $5 \pm 7$ & 0.02 & 0.38 & 0.02$^{***}$ & 0.025$^{***}$ & 0.12$^{***}$ & 0.32$^{*}$ & 0.21$^{*}$ & \small$PL$ $({\gamma}= -1.98)$\\
\hline
EV & 156 & $19 \pm 18$ & 0.06 & 0.3 & 0.067$^{***}$ & 0.07$^{***}$ & 0.14$^{***}$ & 0.26$^{**}$ & 0.17$^{*}$ & \small$PL$ $({\gamma}= -1.22)$\\ 
\hline
LRL & 182 & $26 \pm 22$ & 0.07 & 0.36 & 0.07$^{***}$ & 0.08$^{***}$ & 0.16$^{***}$ & 0.19$^{***}$ & 0.17$^{*}$ & \small$PL$ $({\gamma}= -0.97)$\\ 
\hline 
CAR & 237 & $26 \pm 34$ & 0.05 & 0.19 & 0.057$^{***}$ & 0.06$^{***}$ & 0.12$^{*}$ & 0.2$^{NS}$ & 0.15$^{NS}$ & \small$PL$ $({\gamma}= -0.65)$\\ 
\hline
\end{tabular}
\hspace{0.2 in}
\vspace{0.1 in}

$\bullet$ [Table 1] Maspalomas (MAS), Almunia et al. 1999; Benguela (BEN), Yodzis 1998;
Coachella (COA), Polis 1991; Chesapeake Bay (CHE), Baird \& Ulanowicz 1989;
Skipwith Pond (SKI), Warren 1989; St. Martin (STM), Goldwasser \& Roughgarden
1993; United Kingdom Grassland (UKG), Martinez et al. 1999;
Ythan Estuary (YE), Huxam et al. 1996; Silwood Park (SP), Memmott et al.
2000; El Verde (EV), Reagan \& Waide 1996; Little Rock Lake (LRL), Martinez
1991; and Caribbean Coral Reef (CAR), Opitz 1996.
$S$, number of species; $\bar{k}\pm SD$, mean and standard deviation of the
number of links per species. $C$, connectance, $C_d$, connectance of the most dense subweb for the empirical webs.
$\bar{C}_d_b$, $\bar{C}_d_c$, $\bar{C}_d_n$, $\bar{C}_d_l_r_a$ and
$\bar{C}_d_p_a$, mean connectance of the densest subweb for
$50$ replicates of the basic, cascade, niche, local rewiring algorithm, and
preferential attachment, respectively. Level of significance, $P<0.01$ ($***$), $P<0.05$ ($**$), $0.05<P<0.1$ ($*$), and not significant ($NS$). $Distribution$ refers to the best fit of the
subweb $k$-frequency distribution (cumulative distribution calculated only for the five largest food webs), $PL$, power-law, and the number refers to the scaling exponent (slope).


\section{Figure Legends}

$\bullet$ [Fig. 1] A hypothetical food web graph. A subset of vertices
is called a $k$-subweb if every species within the subset is connected to at
least $k$ prey and/or predators from the same subset. We can observe the
following subwebs:
four separate 0-subwebs (i.e., species have no
links with other species within the same subset, but have one or more links
with other $k$-subwebs of higher degree; black nodes); one 2-subweb; one 3-subweb;
and one 5-subweb, the most dense subweb (white nodes). The links within the most dense
subweb are represented as clearer lines. The density of such interactions
represents the connectance of the most dense subweb ($C_d$). Broken lines represent
the interactions between the densest subweb and the rest of the web. The
density of such interactions represents the intersubweb connectance between the
most dense subweb and the rest of the web.
Note that the web becomes fragmented in five parts if we eliminate the
densest subweb.

$\bullet$ [Fig. 2] Food web structure for Silwood Park (a),
Ythan Estuary (b), El Verde (c), Little Rock Lake (d) and Caribbean
(e). Grey level and line type code as in Fig. 1. 
The Subweb $k$-frequency distribution (represented as the
cumulative distribution $P(S k)$) for Silwood Park (f), Ythan Estuary (g), El
Verde (h), Little Rock Lake (i) and Caribbean (j). As noted, the subweb $k$-frequency distribution is highly skewed, decaying in all cases as a power law with an exponent $\gamma$= $-1.34 \pm 0.57$. The network visualization was done using the Pajek program for large network
analysis:$<http://vlado.fmf.uni-lj.si/pub/networks/pajek/pajekman.htm>$.

\newpage

\section{Figures}
\vspace{0.4 in}

\begin{center}
\includegraphics[width=10cm] {ksubweb.eps}
\vspace{0.1 in}
\centerline{\caption{Fig. 1}}
\end{center}

\newpage

 
\begin{center}
\includegraphics[width=12cm] {Fig2Cohesion1.eps}
\end{center}

\end{document}

\begin{figure}
\hspace{-1in}\begin{minipage} [th] {0.55\textwidth}
\includegraphics[width=4in] {SilSwKFD2.eps}
\end{minipage}
\begin{minipage} [th] {0.7\textwidth}
\includegraphics[width=5in] {KSubElVerde.eps}
\end{minipage}

\hspace{-1.5in}
\begin{minipage}{0.7\textwidth}
\includegraphics[width=4.5in] {YtSubweb.eps}
\end{minipage}
\begin{minipage}{0.7\textwidth}
\includegraphics[width=4.5in] {LRLsubweb1.eps}
\end{minipage}
\renewcommand{\figurename}
\caption{Fig. 2}
\end{figure}

\newpage

\begin{center}
\includegraphics[width=16cm] {kSubFreq1.eps}
\vspace{0.1 in}
\centerline{\caption{Fig. 3}}
\end{center}

\end{document}


\newpage

\begin{center}
\includegraphics[width=13cm] {Cd.eps}
\vspace{0.1 in}
\centerline{\caption{Fig. 4}}
\end{center}








\newpage
\vspace{4cm in}
\begin{center}
\includegraphics[width=13cm] {Cohesion3.eps}
\vspace{0.l in}
\centerline{\caption{Fig. 5}}
\end{center}

\newpage
\vspace{8cm in}
\begin{center}
\includegraphics[width=10cm] {DatRea12.eps}
\vspace{0.l in}
\centerline{\caption{Fig. 6}}
\end{centerJ











