Dynamic growth model for Scots pine (Pinus sylvestris L.) plantations in Galicia (north-western Spain)

In this study we developed a dynamic growth model for Scots pine (Pinus sylvestris L.) plantations in Galicia (north-western Spain). The data used to develop the model were obtained from a network of permanent plots, of between 10 and 55-year-old, which the Unidade de Xestión Forestal Sostible (Sustainable Forest Management Unit) of the University of Santiago de Compostela has set up in pure plantations of this species of pine in its area of distribution in Galicia. In this model, the initial stand conditions at any point in time are defined by three state variables (number of trees per hectare, stand basal area and dominant height), and are used to estimate stand volume, classified by commercial classes, for a given projection age. The model uses three transition functions expressed as algebraic difference equations of the three corresponding state variables used to project the stand state at any point in time. In addition, the model incorporates a function for predicting initial stand basal area, which can be used to establish the starting point for the simulation. This alternative should only be used when the stand is not yet established or when no inventory data are available. Once the state variables are known for a specific moment, a distribution function is used to estimate the number of trees in each diameter class, by recovering the parameters of the Weibull function, using the moments of first and second order of the distribution (arithmetic mean diameter and variance, respectively). By using a generalized height–diameter function to estimate the height of the average tree in each diameter class, combined with a taper function that uses the above predicted diameter and height, it is then possible to estimate total or merchantable stand volume.     

Sensitivity of stand dynamics to grazing in mixed Pinus sylvestris and Quercus pubescens forests: A modelling study

Only recently, studies of forest succession have started to include the effects of browsing by wild or domestic ungulates. We aim to contribute to this topic by analysing the influence of goat grazing on the long-term coexistence of Scots pine (Pinus sylvestris L.) and pubescent oak (Quercus pubescens Willd.) in the low-elevation forests of an inner-Alpine dry valley. The forest gap model ForClim was first adapted to these site conditions by examining the site-dependent sensitivity of the model with regard to the species-specific parameterisation of the drought tolerance as well as the light demand of establishing and adult trees. In a second step, the behaviour of the model was investigated with respect to different grazing intensities and species-specific browsing susceptibilities. The last step was the application of a grazing scenario based on forest history, with 150 years of heavy browsing (by goats) at the beginning of the simulated forest succession, followed by less intensive grazing pressure.     

Modelling spatial and temporal variability in a zero-inflated variable: The case of stone pine (Pinus pinea L.) cone production

Modelling masting habit, i.e. the spatial synchronized annual variability in fruit production, is a huge task due to two main circumstances: (1) the identification of main ecological factors controlling fruiting processes, and (2) the common departure of fruit data series from the main basic statistical assumptions of normality and independence. Stone pine (Pinus pinea L.) is one of the main species in the Mediterranean basin that is able to grow under hard limiting conditions (sandy soils and extreme continental climate), and typically defined as a masting species. Considering the high economical value associated with edible nut production, the masting habit of stone pine has been a main concern for the forest management of the species. In the present work we have used annual fruit data series from 740 stone pine trees measured during a 13 years period (1996-2008) in order: (a) to verify our main hypothesis pointing out to the existence of a weather control of the fruiting process in limiting environments, rather than resource depletion or endogenous inherent cycles; (b) to identify those site factors, stand attributes and climate events affecting specific traits involved in fruiting process; and (c) to construct a model for predicting spatial and temporal patterns of variability in stone pine cone production at different spatial extents as region, stand and tree. Given the nature of the data, the model has been formulated as zero-inflated log-normal, incorporating random components to carry out with the observed lack of independence. This model attains efficiencies close to 70-80% in predicting temporal and spatial variability at regional scale. Though efficiencies are reduced according to the spatial extent of the model, it leads to unbiased estimates and efficiencies over 35-50% when predicting annual yields at tree or stand scale, respectively. In this sense, the proposed model is a main tool for facilitating decision making in some management aspects such as the quantification of total amount of cones annually supplied to nut industry, design of cone harvest programs or the optimal application of seedling felling.     

The response of the woodpigeon (Columba palumbus) to relaxation of intraspecific competition: A hybrid modelling approach

The recent rapid growth of the woodpigeon population in the British Isles is a cause for concern for environmental managers. It is unclear what has driven their increase in abundance. Using a mathematical model, we explored two possible mechanisms, reduced intraspecific competition for food and increased reproductive success. We developed an age-structured hybrid model consisting of a system of ordinary differential equations that describes density-dependent mortality and a discrete component, which represents the birth-pulse. We investigated equilibrium population dynamics using our model. The two hypotheses predict contrasting population age profiles at equilibrium. We adapted the model to examine the impacts of control measures. We showed that an annual shooting season that follows the period of density-dependent mortality is the most effective control strategy because it simultaneously removes adult and juvenile woodpigeons. The model is a first step towards understanding the processes that influence the dynamics of woodpigeon populations.     

A simulation model for Japanese sardine (Sardinops melanostictus) migrations in the western North Pacific

A two-dimensional individual-based model coupled with fish bioenergetics was developed to simulate migration and growth of Japanese sardine (Sardinops melanostictus) in the western North Pacific. In the model, fish movement is controlled by feeding and spawning migrations with passive transport by simulated ocean current. Feeding migration was assumed to be governed by search for local optimal habitats, which is estimated by the spatial distribution of net growth rate of a sardine bioenergetics model. The forage density is one of the most important factors which determines the geographical distributions of Japanese sardine during their feeding migrations. Spawning migration was modeled by an artificial neural network (ANN) with an input layer composed of five neurons that receive environmental information (surface temperature, temperature change experienced, current speed, day length and distance from land). Once the weight of the ANN was determined, the fish movement was solved by combining with the feeding migration model. To obtain the weights of the ANN, three experiments were conducted in which (1) the ANN was trained with back propagation (BP) method with optimum training data, (2) genetic algorithm (GA) was used to adjust the weights and (3) the weights of the ANN were decided by the GA with BP, respectively. BP is a supervised learning technique for training ANNs. GA is a search technique used in computing to find approximate solutions, such as optimization of parameters. Condition factor of sardine in the model is used as a factor of optimization in the GA works. The methods using only BP or GA did not work to search the appropriate weights in the ANN for spawning migration. In the third method, which is a combined approach of GA with BP, the model reproduced the most realistic spawning migration of Japanese sardine. The changes in temperature and day length are important factors for the orientation cues of Japanese sardine according to the sensitivity analysis of the weights of the ANN.     

A model of radiocesium cycling in a sand hills—Turkey Oak (Quercus laevis) ecosystem

Digital simulation models of radiocesium cycling in Turkey Oaks were developed from in situ ¹³⁴Cs tagging studies. Predictions of ¹³⁴Cs steady-state distribution for 3-, 4- and 5-compartment, donor-controlled models were compared with the estimated fallout ¹³⁷Cs distribution as a measure of model validation; output from the 5-compartment model compared best. Sensitivity analysis demonstrated that Turkey Oak burden of ¹³⁴Cs was equally sensitive to the output rate from the tree compartment and the availability of ¹³⁴Cs for uptake (i.e., presence in the root zone) but not the rate of uptake by Turkey Oaks. Observed distribution and model predictions indicate that radiocesium is readily bioaccumulated by Turkey Oaks (～13% of the ecosystem burden) from the soil and is cycled within the sand hills—Turkey Oak ecosystem.     

Influence of host migration between woodland and pasture on the population dynamics of the tick Ixodes ricinus: A modelling approach

Ticks act as vectors of pathogens that can be harmful to animals and/or humans. Epidemiological models can be useful tools to investigate the potential effects of control strategies on diseases such as tick-borne diseases. The modelling of tick population dynamics is a prerequisite to simulating tick-borne diseases and the corresponding spread of the pathogen. We have developed a dynamic model to simulate changes in tick density at different stages (egg, larva, nymph and adult) under the influence of temperature. We have focused on the tick Ixodes ricinus, which is widespread in Europe. The main processes governing the biological cycles of ticks were taken into account: egg laying, hatching, development, host (small, mainly rodents, or large, like deer and cattle, mammals) questing, feeding and mortality. This model was first applied to a homogeneous habitat, where simulations showed the ability of the model to reproduce the general patterns of tick population dynamics. We considered thereafter a multi-habitat model, where three different habitats (woodland, ecotone and meadow) were connected through host migration. Based on this second application, it appears that migration from woodland, via the ecotone, is necessary to sustain the presence of ticks in the meadow. Woodland can therefore be considered as a source of ticks for the meadow, which in turn can be regarded as a sink. The influence of woodland on surrounding tick densities increases in line with the area of this habitat before reaching a plateau. A sensitivity analysis to parameter values was carried out and demonstrated that demographic parameters (sex ratio, development, mortality during feeding and questing, host finding) played a crucial role in the determination of questing nymph densities. This type of modelling approach provides insight into the influence of spatial heterogeneity on tick population dynamics.     

Agent-based simulation of effects of stress on forest tent caterpillar (Malacosoma disstria Hübner) population dynamics

The forest tent caterpillar (Malacosoma disstria Hübner) (FTC) has an outbreak cycle of approximately 10 years; however, smaller spatial scale analyses show some regions have longer or more frequent periods of high defoliation. This may be a result of local forest fragmentation, pollution or other sources of stress that may affect FTC directly or indirectly through stress on their hosts or parasitoids. Population dynamics of FTC were examined to investigate how stress may alter the severity and frequency of defoliation. We developed a spatially explicit agent-based model to simulate the host-parasitoid dynamics of FTC. Theoretical and empirically derived parameters were established using past literature and over 50 years of population data of FTC from Ontario, Canada. We find that increasing FTC fecundity, FTC dispersal or parasitoid mortality resulted in more severe outbreaks while a decrease in parasitoid fecundity or searching efficiency resulted in an overall elevation of defoliation. Parasitoid efficiency was the most effective parameter for altering the FTC defoliation. Since plant stress has been shown to alter several of these parameters in nature due to changes in food quality, habitat suitability, and chemical cue interference, our results suggest that forests affected by stressors such as climate change and pollution will have more severe and frequent defoliation from these insects than surrounding unaffected forests. As stressors such as drought and pollution emissions are predicted to increase in frequency or intensity over the next few decades, understanding how they may affect the outbreak cycle of a forest defoliator can aid in planning strategies to reduce the detrimental effects of this insect.     

Population dynamics of splitnose rockfish (Sebastes diploproa) in the Northeast Pacific Ocean

We developed an age-structured population model of splitnose rockfish, Sebastes diploproa, in the Northeast Pacific Ocean. Splitnose rockfish is a bycatch species that co-occurs with several commercially important species that are currently declared overfished. Bycatch species are typically not the focus of stock assessment efforts because of their limited economic importance, but they may suffer the same population declines as species with which they co-occur. To examine the dynamics of splitnose rockfish for the first time, we analyzed data from three groundfish fisheries and four research surveys conducted in the Northeast Pacific Ocean. To develop a model, we used Stock Synthesis, a statistical framework for the construction of a population dynamics models utilizing both fishery-dependent and fishery-independent data. In the model, we reconstructed the total catch of the species back to 1900, estimated the dynamics of the stock spawning output and recruitment and evaluated biomass depletion relative to the stock's unfished state, as well as sources of uncertainty in model outputs. The results indicate that the splitnose rockfish is currently not overfished even though it has experienced several periods of abrupt decline in its biomass. Revisiting age data from earlier years, monitoring fishery discard, and investigating the spatial dynamics of splitnose rockfish is important to further improve the understanding of this species’ population dynamics, and decrease uncertainty in model results.     

Modeling the impacts of hunting on the population dynamics of red howler monkeys (Alouatta seniculus)

Overexploitation of wildlife populations occurs across the humid tropics and is a significant threat to the long-term survival of large-bodied primates. To investigate the impacts of hunting on primates and ways to mitigate them, we developed a spatially explicit, individual-based model for a landscape that included hunted and un-hunted areas. We used the large-bodied neotropical red howler monkey (Alouatta seniculus) as our case study species because its life history characteristics make it vulnerable to hunting. We modeled the influence of different rates of harvest and proportions of landscape dedicated to un-hunted reserves on population persistence, population size, social dynamics, and hunting yields of red howler monkeys. In most scenarios, the un-hunted populations maintained a constant density regardless of hunting pressure elsewhere, and allowed the overall population to persist. Therefore, the overall population was quite resilient to extinction; only in scenarios without any un-hunted areas did the population go extinct. However, the total and hunted populations did experience large declines over 100 years under moderate and high hunting pressure. In addition, when reserve area decreased, population losses and losses per unit area increased disproportionately. Furthermore, hunting disrupted the social structure of troops. The number of male turnovers and infanticides increased in hunted populations, while birth rates decreased and exacerbated population losses due to hunting. Finally, our results indicated that when more than 55% of the landscape was harvested at high (30%) rates, hunting yields, as measured by kilograms of biomass, were less than those obtained from moderate harvest rates. Additionally, hunting yields, expressed as the number of individuals hunted/year/km², increased in proximity to un-hunted areas, and suggested that dispersal from un-hunted areas may have contributed to hunting sustainability. These results indicate that un-hunted areas serve to enhance hunting yields, population size, and population persistence in hunted landscapes. Therefore, spatial regulation of hunting via a reserve system may be an effective management strategy for sustainable hunting, and we recommend it because it may also be more feasible to implement than harvest quotas or restrictions on season length.     

A system dynamics model for the management of the Manila clam, Ruditapes philippinarum (Adams and Reeve, 1850) in the Bay of Arcachon (France)

The Manila clam Ruditapes philippinarum (Adams and Reeve, 1850) is one of the mollusc species that, driven mainly by the shellfish market industry, has extended throughout the world, far beyond the limits of its original habitat. The Manila clam was introduced into France for aquaculture purposes, between 1972 and 1975. In France, this venerid culture became increasingly widespread and, since 1988, this species has colonised most of the embayments along the French Atlantic coast. In 2004, this development resulted in a fishery of ca. 520 t in Arcachon Bay.     

Identifying key life-traits for the dynamics and gene flow in a weedy crop relative: Sensitivity analysis of the GeneSys simulation model for weed beet (Beta vulgaris ssp. vulgaris)

The benefits of genetically modified herbicide-tolerant (GMHT) sugar beet (Beta vulgaris) varieties stem from their presumed ability to improve weed control and reduce its cost, particularly targeting weed beet, a harmful annual weedy form of the genus Beta (i.e. B. vulgaris ssp. vulgaris) frequent in sugar beet fields. As weed beet is totally interfertile with sugar beet, it is thus likely to inherit the herbicide-tolerance transgene through pollen-mediated gene flow. Hence, the foreseeable advent of HT weed beet populations is a serious threat to the sustainability of GM sugar beet cropping systems. For studying and quantifying the long-term effects of cropping system components (crop succession and cultivation techniques) on weed beet population dynamics and gene flow, we developed a biophysical process-based model called GeneSys-Beet in a previous study. In the present paper, the model was employed to identify and rank the weed life-traits as function of their effect on weed beet densities and genotypes, using a global sensitivity analysis to model parameters. Monte Carlo simulations with simultaneous randomization of all life-trait parameters were carried out in three cropping systems contrasting for their risk for infestation by HT weed beets. Simulated weed plants and bolters (i.e. beet plants with flowering and seed-producing stems) were then analysed with regression models as a function of model parameters to rank processes and life-traits and quantify their effects. Key parameters were those determining the timing and success of growth, development, seed maturation and the physiological end of seed production. Timing parameters were usually more important than success parameters, showing for instance that optimal timing of weed management operations is more important than its exact efficacy. The ranking of life-traits though depended on the cropping system and, to a lesser extent, on the target variable (i.e. GM weeds vs. total weed population). For instance, post-emergence parameters were crucial in rotations with frequent sugar beet crops whereas pre-emergence parameters were most important when sugar beet was rare. In the rotations with frequent sugar beet and insufficient weed control, interactions between traits were small, indicating diverse populations with contrasted traits could prosper. Conversely, when sugar beet was rare and weed control optimal, traits had little impact individually, indicating that a small number of optimal combinations of traits would be successful. Based on the analysis of sugar beet parameters and genetic traits, advice for the future selection of sugar beet varieties was also given. In climatic conditions similar to those used here, the priority should be given to limiting the presence of hybrid seeds in seed lots rather than decreasing varietal sensitivity to vernalization.     

Simulating the temporal pattern of waste production in farmed gilthead seabream (Sparus aurata), European seabass (Dicentrarchus labrax) and Atlantic bluefin tuna (Thunnus thynnus)

Most fish farming waste output models provide gross waste rates as a function of stocked or produced biomass for a year or total culture cycle, but without contemplating the temporality of the discharges. This work aims to ascertain the temporal pattern of waste loads by coupling available growth and waste production models and developing simulation under real production rearing conditions, considering the overlapping of batches and management of stocks for three widely cultured species in the Mediterranean Sea: gilthead seabream (Sparus aurata), European seabass (Dicentrarchus labrax) and Atlantic bluefin tuna (Thunnus thynnus). For a similar annual biomass production, the simulations showed that waste output and temporal dumping patterns differ between the three species as a result of the disparities in growth velocity, nutrient digestibility, maintenance metabolic budget and husbandry. The simulations allowed the temporal patterns including the periods of maximum discharge and the dissolved and particulate nitrogen and phosphorus content in the wastes released to be determined, both of which were seen to be species-specific.     

Modeling the population dynamics of cotton bollworm Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae) over a wide area in northern China

The cotton bollworm Helicoverpa armigera (Lepidoptera: Noctuidae) is one of the most serious crop pests in northern China, calling for accurate prediction of pest outbreaks and strategies for pest control. A computer model is developed to simulate the population dynamics of H. armigera over a wide area in northern China. The area considered covers 12 provinces where serious outbreaks of H. armigera have been observed. In this model, pest development is driven by local ambient temperature, and adults migrate long distances between regions and select preferred hosts for oviposition within a region. Six types of host including cotton, wheat, corn, peanut, soybean and a single category composed of all other minor hosts are considered in this model. Survival rates of eggs and larvae are based on life-table data, and simulated as a function of host type, host phenology and temperature. The incidence of diapause depends on temperature and photoperiod experienced during the larval stage. Survival rate of non-diapause pupae is a nonlinear function of rainfall, and overwinter survival rate is a nonlinear function of temperature. Insecticide is applied when population density exceeds the economic threshold on a host crop within a region. Comparisons of model output with light-trap data indicate that our model reflects the pest population dynamics over a wide area, and could potentially be used for testing novel pest control strategies in northern China.     

Application of a temperature-dependent von Bertalanffy growth model to bullhead (Cottus gobio)

The most studied and commonly applied model of fish growth is the von Bertalanffy model. However, this model does not take water temperature into account, which is one of the most important environmental factors affecting the life cycle of fish, as many physiological processes that determine growth, e.g. metabolic rate and oxygen supply, are directly influenced by temperature. In the present study we propose a version of the von Bertalanffy growth model that includes mean annual water temperatures by correlating the growth coefficient, k, explicitly and the asymptotic length, L_∞, implicitly to water temperature. All relationships include parameters with an obvious biological relevance that makes them easier to identify. The model is used to fit growth data of bullhead (Cottus gobio) at different locations in the Bez River network (Drme, France). We show that temperature explains much of the growth variability at the different sampling sites of the network.     

Analysis of annual fluctuations of C. nodosa in the Venice lagoon: Modeling approach

A simple simulation model was developed to describe the growth trends of Cymodocea nodosa (Ucria) Ascherson based on data sets from the Venice lagoon. The model reproduces the seasonal fluctuations in the above and belowground biomass and in shoot density. The modeling results are in good agreement with data on net production, growth rates and chemical–physical parameters of water. It was assumed that light and temperature are the most important factors controlling C. nodosa development, and that the growth was not limited by nutrient availability. The aim was to simulate biomass production as a function of external forcing variables (light, water temperature) and internal control (plant density). A series of simulation experiments were performed with the basic model showing that among the most important phenomena affecting C. nodosa growth are: (1) inhibition of production and recruitment of new shoots by high temperature and (2) light attenuation due to seasonal fluctuation.     

Validation of a multispecies forest dynamics model using 50-year growth from Eucalyptus forests in eastern Australia

One of the key problems confronting ecological forecasting is the validation of computer models. Here we report successful validation of a forest dynamics model Ecosystem Dynamics Simulator (EDS), adapted from the JABOWA-II forest succession model. This model and many variants derived from it have successfully simulated growth dynamics of uneven-aged mixed forests under changing environment with a moderate amount of input data. But rarely are adequate time-series data available for quantitative model validation. This study tested the performance of EDS in projecting the tree density, tree diameter at breast height (dbh), tree height, basal area and aboveground biomass of uneven-aged, mixed species sclerophyll forests in St. Mary state forests of eastern Australia. The test data were collected between 1951 and 2005. Every tree was uniquely numbered, tagged and measured in consecutive re-measurements. Projected growth attributes were compared with those observed in an independent validation dataset. The model produced satisfactory projections of tree density (91.7%), dbh (92.3%), total tree height (82.8%), basal area (89.3%) and aboveground biomass (87.6%) compared to the observed attributes. These results suggest that the EDS model can provide reasonable capability in projecting growth dynamics of uneven-aged, mixed species sclerophyll forests.     

A simulation study of the medium-term effects of field patterns on cross-pollination rates in oilseed rape (Brassica napus L.)

Gene flow between cultivars within a landscape may lead to impurities that reduce harvest value. In OSR, as for most crops, impurity rates are expected to depend on the spatial distribution of crops over the landscape. However, in contrast to other well-studied crops such as maize, OSR crops generate seed banks in European agro-ecosystems. Gene flow is thus a spatio-temporal process which depends on cropping systems. We therefore aimed at identifying spatial variables that have an effect on regional or local harvest impurities, taking account of the time since the introduction of OSR crops in the regions and of cropping system. Gene flow was simulated over 36 field patterns cultivated with either 15% or 30% of OSR fields, among which 10% or 50% were GM, for three contrasted cropping systems, with the GeneSys software already used for EU co-existence studies. Through regression analyses, we determined spatial and agronomic factors that most affected harvest impurity rates of non-GM OSR after one or seven years of OSR cultivation. The cropping system was the main factor explaining regional harvest impurity rates. Its importance increased after six years of OSR cultivation. For a given cropping system, the regional impurity rate after one year increased linearly with the current proportion of GM crop. In contrast, impurity rates after six years largely depended on the proportions of OSR crop (GM or not) in the two preceding years. During the first year of OSR cultivation, local impurity rates were mostly explained by the distance to the closest GM field. After six years, these rates were mostly explained by the density of GM volunteers in the analysed field and, to a lesser degree, to that of volunteers in neighbour non-OSR fields. Cropping systems were most important in determining impurity rates and the way impurity rates related to regional or local factors. Determination of isolation distances to ensure harvest purity should thus consider past history of OSR cultivation in the area and, in particular, how current or future cropping systems will manage volunteers. Regression quantiles were fitted to the simulated data to determine regional rules (i.e. the maximum regional area of GM OSR and isolation distances between GM and non-GM crops) as a function of the risk accepted by the decision-maker (i.e. the % of situations exceeding harvest impurity thresholds), the cropping system and the volunteer infestation.     

A stochastic simulation model of brown shrimp (Penaeus aztecus ives) growth,movement, and survival in Galveston Bay,Texas

A stochastic simulation model of brown shrimp (Penaeus aztecus Ives) population dynamics in Galveston Bay, Texas, is described, validated, and used to evaluate the effects of management alternatives and changing environmental conditions on shrimp dynamics. The model is composed of submodels representing: (1) recruitment, (2) growth, (3) natural mortality, (4) fishing mortality, and (5) emigration of brown shrimp. The model predicts significant changes in total annual harvest from the food shrimp, bait, and recreational fisheries resulting from (1) closure of the bay system to all fishing except during the spring and fall open seasons, (2) two-week postponement of the opening and closing of the open seasons for the food shrimp fishery, (3) a 2.5°C increase and (4) a 2.5°C decrease in mean water temperature, (5) an 80% increase and (6) an 80% decrease in fishing effort. No significant change in the total annual harvest is predicted when the food shrimp fishing season is extended from May 15 through December 15. Sensitivity analysis suggests that field experimentation designed specifically to test the hypothesis of a 60-day time lag between brown shrimp recruitment into the bays and exposure to the fishery should receive high priority. Simulation results are discussed within a management framework.     

Sea level rise and eelgrass (Zostera marina) production: A spatially explicit relative elevation model for Padilla Bay, WA

The dynamics that govern the elevation of a coastal wetland relative to sea level are complex, involving non-linear feedbacks among opposing processes. Changes in the balance between these processes can result in significant alterations to vegetation communities that are adapted to a specific range of water levels. Given that current sedimentation rates in Padilla Bay, Washington are likely less than historical levels and that eustatic sea level rise is accelerating, the extensive Zostera marina (eelgrass) meadows in the bay may be at risk of eventual submergence. We developed a spatially explicit relative elevation model and used it to project changes in the productivity and distribution of eelgrass in Padilla Bay over the next century. The model is mechanistic and incorporates many of the processes and feedbacks that govern coastal wetland elevation change. Accretion estimates made using ²¹⁰Pb dating of sediment cores, sediment characteristics measured within cores, and eelgrass productivity and decomposition data were used to initialize and calibrate the model. Validation was performed using an elevation change rate measured with a network of surface elevation tables. Both the field data and model simulations revealed a net accretion deficit for the bay. Simulations using current rates of sea level rise indicated an overall expansion of eelgrass within Padilla Bay over the next century as it migrates from the center of the bay shoreward.