Uncertainty analysis of climate change mitigation options in the forestry sector using a generic carbon budget model

Industrialized countries agreed on a reduction of greenhouse gas emissions under the Kyoto Protocol. Many countries elected forest management activities and the resulting net balance of carbon emissions and removals of non-CO₂ greenhouse gases by forest management in their climate change mitigation measures. In this paper a generic dynamic forestry model (FORMICA) is presented. It has an empirical basis. Several modules trace C pools relevant for the Kyoto Protocol and beyond: biomass, litter, deadwood and soil, and harvested wood products. The model also accounts for the substitution of fossil fuels by wood products and bioenergy.     

Modelling carbon and nitrogen dynamics in forest ecosystems of Central Russia under different climate change scenarios and forest management regimes

The individual-based stand-level model EFIMOD was used for large-scale simulations using standard data on forest inventories as model inputs. The model was verified for the case-study of field observations, and possible sources of uncertainties were analysed. The approach developed kept the ability for fine-tuning to account for spatial discontinuity in the simulated area. Several forest management regimes were simulated as well as forest wildfires and climate changes. The greatest carbon and nitrogen accumulations were observed for the regime without cuttings. It was shown that cuttings and wildfires strongly influence the processes of carbon and nitrogen accumulations in both soil and forest vegetation. Modelling also showed that the increase in annual average temperatures resulted in the partial relocation of carbon and nitrogen stocks from soil to plant biomass. However, forest management, particularly harvesting, has a greater effect on the dynamics of forest ecosystems than the prescribed climate change.     

treedyn3 forest simulation model

The treedyn3 forest simulation model is a process model of tree growth, carbon and nitrogen dynamics in a single-species, even-aged forest stand. It is based on the treedyn model. Major changes include the computation of sun angle and radiation as a function of latitude and day of the year, the closed-form integration of canopy production as a function of day and hour, the introduction of tree number, height, and diameter as separate state variables, and different growth strategies, mortalities, and resulting self-thinning as function of crowding competition.The tree/soil system is described by a set of nonlinear ordinary differential equations for the state variables: tree number, base diameter, tree height, wood biomass, nitrogen in wood, leaf mass, fine root mass, fruit biomass, assimilate, carbon and nitrogen in litter, carbon and nitrogen in soil organic matter, and plant-available nitrogen. The model includes explicit formulations of all relevant ecophysiological processes such as: computation of radiation as a function of seasonal time, daytime and cloudiness, light attenuation in the canopy, and canopy photosynthesis as function of latitude, seasonal time, and daytime, respiration of all parts, assimilate allocation, increment formation, nitrogen fixation, mineralization, humification and leaching, forest management (thinning, felling, litter removal, fertilization etc.), temperature effects on respiration and decomposition, and environmental effects (pollution damage to photosynthesis, leaves, and fine roots). Only ecophysiological parameters which can be either directly measured or estimated with reasonable certainty are used. treedyn3 is a generic process model which requires species- and site-specific parametrization. It can be applied to deciduous and coniferous forests under tropical, as well as temperate or boreal conditions.The paper presents a full documentation of the mathematical model as well as representative simulation results for spruce and acacia.     

Optimal resource management control for CO2 emission and reduction of the greenhouse effect

In recent years, the world has witnessed an ever-growing concern towards global warming caused by greenhouse gases, such as carbon dioxide (CO₂). In order to reduce the emissions of CO₂ without limiting economic growth, substantial investments should target the development of clean technology and the expansion of forested areas. Considering the limited availability of resources, investments must be used in the most effective way. The present work proposes a method to efficiently manage these resources by applying the optimal control theory to a new mathematical model that describes the dynamics of the atmospheric CO₂. The contributions of this work are twofold: (1) present a model that describes the dynamic relation of CO₂ emission with investment in reforestation and clean technology and (2) present a method to efficiently manage the available resources by casting an optimal control problem. The mathematical model uses ordinary differential equations to relate the production of CO₂ with forest area and Gross Domestic Product (GDP). The model parameters are adjusted to fit the actual published data. Given an appropriate performance index, the optimal solution is found by numerically solving the Two-Point Boundary Value Problem (TPBVP) that arises from the application of Pontriagyn's Maximum Principle. The sensitivity of the obtained numerical solution is evaluated with respect to the uncertainties in the model parameters. The main objective of this work is to provide a quantitative tool for the efficient allocation of resources to reduce the greenhouse effect caused CO₂ emissions.     

Constructing change: The evolution of land and resource management in Alberta,Canada

Recognizing the limitations common to both centralized and privatized management regimes, institutionalized resource management is beginning to incorporate the knowledge and skills of local resource users, coupled with the enabling policies and legislation of state systems, to arrive at cooperative approaches to resource management. These varying and dynamic approaches to resource management have been compelled largely through the recognition of the limited capabilities of existing management systems to adapt effectively to ecosystem change and the evolving needs of resource users. These cooperative approaches to management should not, however, be considered an institutional end-point, but rather a phase in the perpetual transition of a social system; each unique in character and individually variable depending on the resource being managed, the political climate in which management occurs, as well as the differing strategies employed by resource users to enact institutional change. Drawing from the experiences of the Whitefish Lake First Nation of Alberta, Canada, this paper presents a brief overview of the evolution of resource management theory, grounded in the real-world formation of the Whitefish Lake First Nation — Province of Alberta Cooperative Management Agreement.     

The political ecology of Yukon forestry: managing the forest as if people mattered

Over the past several decades a growing interest has emerged in the relationship between forest dependent communities and their use of non-timber forest resources (NTFR—i.e. berries, roots, barks). Motivated in large part by international concerns over the loss of biological and cultural diversity, efforts are now being made to conserve the world's remaining forests by repositioning NTFR to the forefront of the forest management process. Together with the Little Salmon Carmacks First Nation, we have undertaken an assessment that: 1) addresses the role of NTFR in fulfilling basic human necessities and well-being of community members, and 2) identifies the perceptual differences that exist among First Nation and non-First Nation community members regarding the security and future availability of forest resources. Based on our findings, we argue that if the true value of the forest is to be understood, and thus sustained, those involved in management must consider the full array of products and services the forest provides. Thus the conventional and over-simplified view that the forest is to provide primarily for commercial timber must be challenged in light of the multiplicity of values community members gain from forest use.     

Alternative transportation fuel standards: Welfare effects and climate benefits

This paper develops an integrated model of the fuel and agricultural sectors to analyze the welfare and greenhouse gas emission (GHG) effects of the existing Renewable Fuel Standard (RFS), a Low Carbon Fuel Standard (LCFS) and a carbon price policy. The conceptual framework shows that these policies differ in the incentives they create for the consumption and mix of different types of biofuels and in their effects on food and fuel prices and GHG emissions. We also simulate the welfare and GHG effects of these three policies which are normalized to achieve the same level of US GHG emissions. By promoting greater production of food-crop based biofuels, the RFS is found to lead to a larger reduction in fossil fuel use but also a larger increase in food prices and a smaller reduction in global GHG emissions compared to the LCFS and carbon tax. All three policies increase US social welfare compared to a no-biofuel baseline scenario due to improved terms-of-trade, even when environmental benefits are excluded; global social welfare increases with a carbon tax but decreases with the RFS and LCFS due to the efficiency costs imposed by these policies, even after including the benefits of mitigating GHG emissions.     

Modelling forest management within a global vegetation model—Part 1: Model structure and general behaviour

This article describes a new forest management module (FMM) that explicitly simulates forest stand growth and management within a process-based global vegetation model (GVM) called ORCHIDEE. The net primary productivity simulated by ORCHIDEE is used as an input to the FMM. The FMM then calculates stand and management characteristics such as stand density, tree size distribution, tree growth, the timing and intensity of thinnings and clear-cuts, wood extraction and litter generated after thinning. Some of these variables are then fed back to ORCHIDEE. These computations are made possible with a distribution-based modelling of individual tree size. The model derives natural mortality from the relative density index (rdi), a competition index based on tree size and stand density. Based on the common forestry management principle of avoiding natural mortality, a set of rules is defined to calculate the recurrent intensity and frequency of forestry operations during the stand lifetime. The new-coupled model is called ORCHIDEE-FM (forest management).The general behaviour of ORCHIDEE-FM is analysed for a broadleaf forest in north-eastern France. Flux simulation throughout a forest rotation compare well with the literature values, both in absolute values and dynamics.Results from ORCHIDEE-FM highlight the impact of forest management on ecosystem C-cycling, both in terms of carbon fluxes and stocks. In particular, the average net ecosystem productivity (NEP) of 225 gC m⁻² year⁻¹ is close to the biome average of 311 gC m⁻² year⁻¹. The NEP of the “unmanaged” case is 40% lower, leading us to conclude that management explains 40% of the cumulated carbon sink over 150 years. A sensitivity analysis reveals 4 major avenues for improvement: a better determination of initial conditions, an improved allocation scheme to explain age-related decline in productivity, and an increased specificity of both the self-thinning curve and the biomass-diameter allometry.     

Sanitary felling of Norway spruce due to spruce bark beetles in Slovenia: A model and projections for various climate change scenarios

A model is presented to predict sanitary felling of Norway spruce (Picea abies) due to spruce bark beetles (Ips typographus, Pityogenes chalcographus) in Slovenia according to different climate change scenarios. The model incorporates 21 variables that are directly or indirectly related to the dependent variable, and that can be arranged into five groups: climate, forest, landscape, topography, and soil. The soil properties are represented by 8 variables, 4 variables define the topography, 4 describe the climate, 4 define the landscape, and one additional variable provides the quantity of Norway spruce present in the model cell. The model was developed using the M5′ model tree. The basic spatial unit of the model is 1 km², and the time resolution is 1 year. The model evaluation was performed by three different measures: (1) the correlation coefficient (51.9%), (2) the Theil's inequality coefficient (0.49) and (3) the modelling efficiency (0.32). Validation of the model was carried out by 10-fold cross-validation. The model tree consists of 28 linear models, and model was calculated for three different climate change scenarios extending over a period until 2100, in 10-year intervals. The model is valid for the entire area of Slovenia; however, climate change projections were made only for the Maribor region (596 km²). The model assumes that relationships among the incorporated factors will remain unchanged under climate change, and the influence of humans was not taken into account. The structure of the model reveals the great importance of landscape variables, which proved to be positively correlated with the dependent variable. Variables that describe the water regime in the model cell were also highly correlated with the dependent variable, with evapotranspiration and parent material being of particular importance. The results of the model support the hypothesis that bark beetles do greater damage to Norway spruce artificially planted out of its native range in Slovenia, i.e., lowlands and soils rich in N, P, and K. The model calculation for climate change scenarios in the Maribor region shows an increase in sanitary felling of Norway spruce due to spruce bark beetles, for all scenarios. The model provides a path towards better understanding of the complex ecological interactions involved in bark beetle outbreaks. Potential application of the results in forest management and planning is discussed.     

Inter-annual variability of ecosystem production in boreal jack pine forests (1975-2004) estimated from tree-ring data using CBM-CFS3

We describe and apply a method of using tree-ring data and an ecosystem model to reconstruct past annual rates of ecosystem production. Annual data on merchantable wood volume increment and mortality obtained by dendrochronological stand reconstruction were used as input to the Carbon Budget Model of the Canadian Forest Sector (CBM-CFS3) to estimate net ecosystem production (NEP), net primary production (NPP), and heterotrophic respiration (Rh) annually from 1975 to 2004 at 10 boreal jack pine (Pinus banksiana Lamb.) stands in Saskatchewan and Manitoba, Canada. From 1975 (when sites aged 41-60 years) to 2004 (when they aged 70-89 years), all sites were moderate C sinks except during some warmer than average years where estimated Rh increased. Across all sites and years, estimated annual NEP averaged 57 g Cm⁻² yr⁻¹ (range −31 to 176 g Cm⁻² yr⁻¹), NPP 244 g Cm⁻² yr⁻¹ (147-376 g Cm⁻² yr⁻¹), and Rh 187 g Cm⁻² yr⁻¹ (124-270 g Cm⁻² yr⁻¹). Across all sites, NPP was related to stand age and density, which are proxies for successional changes in leaf area. Regionally, warm spring temperature increased NPP and defoliation by jack pine budworm 1 year previously reduced NPP. Our estimates of NPP, Rh, and NEP were plausible when compared to regional eddy covariance and carbon stock measurements. Inter-annual variability in ecosystem productivity contributes uncertainty to inventory-based assessments of regional forest C budgets that use yield curves predicting averaged growth over time. Our method could expand the spatial and temporal coverage of annual forest productivity estimates, providing additional data for the development of empirical models accounting for factors not presently considered by these models.     

A tree and climate assessment tool for modelling ecosystem response to climate change

Understanding how vulnerable forest ecosystems are to climate change is a key requirement if sustainable forest management is to be achieved. Modelling the response of species in their regeneration niche to phenological and biophysical processes that are directly influenced by climate is one method for achieving this understanding. A model was developed to investigate species resilience and vulnerability to climate change within its fundamental-regeneration niche. The utility of the developed model, tree and climate assessment (TACA), was tested within the interior Douglas-fir ecosystem in south-central British Columbia. TACA modelled the current potential tree species composition of the ecosystem with high accuracy and modelled significant responses amongst tree species to climate change. The response of individual species suggests that the studied ecosystem could transition to a new ecosystem over the next 100 years. TACA showed that it can be an effective tool for identifying species resilience and vulnerability to changes in climate within the most sensitive stage of development, the regeneration phase. The TACA model was able to identify the degree of change in phenological and biophysical variables that control tree establishment, growth and persistence. The response to changes in one or more of these variables resulted in changes in the climatic suitability of the ecosystem for species and enabled a measure of vulnerability to be quantified. TACA could be useful to forest managers as a decision support tool for adaptation actions and by researchers interested in modelling stand dynamics under climate change.     

Managing mangroves in Bangladesh: A strategy analysis

Bangladesh, favoured by a tropical climate, houses the world’s largest stretch of mangroves forests (Sundarbans Reserved Forest) and plantations. Around half of the forests of the country occur in the coastal zone. People extract various goods and services from the mangroves. Nevertheless the mangrove forests are depleting. Although the extent of the Sundarbans forest has not changed much, its decline is of a qualitative nature. Mangrove plantations are increasing in area but they are losing growing stock. To arrest this, Bangladesh has adopted several strategies. The ‘Sustainable Ecosystem Management’ strategy has now been adopted instead of the ‘Sustained Yield Principle’. Biodiversity conservation and enhancement has been taken as a key management goal. A zoning system is being developed for both production and protection purposes. The government facilitates alternative income for the local people by generating activities for the communities which are dependent on the forest. Different non-governmental organizations collaborate with the government in reducing the local people’s dependence on the forest. Coastal plantations are erected to protect people from cyclones and to make the land more suitable for habitation. Through this greening of the coastal belt tree plantation is encouraged in coastal villages. Coastal embankments are being planted and leased to poor settlers in exchange for routine maintenance of the embankments. Plantations on newly accreted mud flats help in stabilizing the land, which can later on be settled by victims of erosion elsewhere. These adopted management measures do not only contribute to forestry resource management but also to the social, environmental and economic wellbeing of the coastal communities. These efforts are at present being integrated into an Integrated Coastal Zone Management (ICZM) project.     

Mathematical analysis of change in forest carbon use efficiency with stand development: A case study on Abies veitchii Lindl.

Changes in carbon use efficiency (CUE), which is defined as the ratio of net primary production (NPP) to gross primary production (GPP), were analyzed for Abies veitchii Lindl. forests with respect to stand development by developing a simple mathematical model incorporating data on physiological variables and leaf mass ratio. A decrease in CUE with stand development was successfully expressed as a function of stand biomass (y) based on the following three assumptions: (1) a power-law relationship between mean respiration and mean individual tree mass, (2) a power-functional relationship between mean gross primary production and mean individual tree mass, and (3) self-thinning relationship between stand biomass and density. Based on this model, a parameter of CUE–y relationship was defined, and it was clarified that CUE decrease with stand development is caused not by the ratio of specific respiration rate to specific gross photosynthetic rate, but by leaf mass ratio. Since CUE is high in young forests, helpful information on selecting woody species when planting seedlings was provided from the viewpoints of reducing CO₂ in the atmosphere and global warming.     

The effects of disturbance based management on the dynamics of Mediterranean vegetation: A hierarchical and spatially explicit modeling approach

The eastern Mediterranean region has been subjected to intensive human disturbance in the past 10,000 years, mainly in the forms of agro-pastoral activities such as grazing, shrub clearing, and prescribed burning. This disturbance history resulted in the formation of highly heterogeneous landscapes, characterized by high biodiversity. Recent changes in human activities have resulted in a decrease of landscape heterogeneity, leading to decreasing biodiversity and increasing fire risk. To conserve heterogeneity, land managers apply disturbance based management practices, using the same activities that created and maintained landscape heterogeneity in the past. However, the long-term and large-scale outcomes of these disturbances are often unknown, due to the complex response of Mediterranean vegetation to disturbance. Here we report on a spatially explicit, hybrid, and spatially hierarchical ecological model developed by us. The model attempts to predict the outcome of various disturbance based management activities on the long-term spatio-temporal dynamics of five common Mediterranean vegetation types. The model uses a spatially explicit state and transition formulation, with continuous transition functions. Model simulations were conducted on a Mediterranean landscape in Northern Israel, incorporating various disturbance practices that are common in the region. Simulation results highlight the potential of disturbance based management as a tool for conserving landscape heterogeneity, as well as the complex interactions between disturbances and the spatial structure of the landscape in Mediterranean regions.     

Dissolved organic carbon concentrations and fluxes in forest catchments and streams: DOC-3 model

Dissolved organic carbon (DOC) concentrations in south-western Nova Scotia streams, sampled at weekly to biweekly intervals, varied across streams from about 3 to 40 mg L⁻¹, being highest mid-summer to fall, and lowest during winter to spring. A 3-parameter model (DOC-3) was proposed to project daily stream DOC concentrations and fluxes from modelled estimates for daily soil temperature and moisture, year-round, and in relation to basin size and wetness. The parameters of this model refer to (i) a basin-specific DOC release parameter “k_DOC^”, related to the wet- and open-water area percentages per basin, (ii) the lag time “τ” between DOC production and subsequent stream DOC emergence, related to the catchment area above the stream sampling location; and (iii) the activation energy “Ea”, to deal with the temperature effect on DOC production. This model was calibrated with the 1988-2006 DOC concentration data from three streams (Pine Marten, Moosepit Brook, and the Mersey River sampled at or near Kejimkujik National Park, or KNP), and was used to interpret the biweekly 1999-2003 DOC concentrations data (stream, ground and lake water, soil lysimeters) of the Pockwock-Bowater Watershed Project near Halifax, Nova Scotia. The data and the model revealed that the DOC concentrations within the streams were not correlated to the DOC concentrations within the soil- and groundwater, but were predictable based on (i) the hydrologically inferred weather-induced changes in soil moisture and temperature next to each stream, and (ii) the topographically inferred basin area and wet- and open-water area percentages associated with each stream (R² = 0.53; RMSE = 3.5 mg L⁻¹). Model-predicted fluxes accounted 74% of the hydrometrically determined DOC exports at KNP.     

Coupling biochemical and biophysical processes at the leaf level: an equilibrium photosynthesis model for leaves of C3 plants

The paper presents a generic computer model for estimating short-term steady-state fluxes of CO₂, water vapor, and heat from broad leaves and needle-leaved coniferous shoots of C₃ plant species. The model explicitly couples all major processes and feedbacks known to impact leaf biochemistry and biophysics including biochemical reactions, stomatal function, and leaf-boundary layer heat- and mass-transport mechanisms. The ability of the model to successfully predict measured photosynthesis and stomatal-conductance data as well as to simulate a variety of observed leaf responses is demonstrated. A model application investigating physiological and environmental regulation of leaf water-use efficiency (WUE) under steady-state conditions is discussed. Simulation results suggest that leaf physiology has a significant control over the environmental sensitivity of leaf WUE. The implementation of a highly efficient solution technique allows the model to be directly incorporated into plant-canopy and terrestrial ecosystem models.     

Response of tree growth to a changing climate in boreal central Canada: A comparison of empirical,process-based,and hybrid modelling approaches

The impact of 2 × CO₂ driven climate change on radial growth of boreal tree species Pinus banksiana Lamb., Populus tremuloides Michx. and Picea mariana (Mill.) BSP growing in the Duck Mountain Provincial Forest of Manitoba (DMPF), Canada, is simulated using empirical and process-based model approaches. First, empirical relationships between growth and climate are developed. Stepwise multiple-regression models are conducted between tree-ring growth increments (TRGI) and monthly drought, precipitation and temperature series. Predictive skills are tested using a calibration–verification scheme. The established relationships are then transferred to climates driven by 1× and 2 × CO₂ scenarios using outputs from the Canadian second-generation coupled global climate model. Second, empirical results are contrasted with process-based projections of net primary productivity allocated to stem development (NPP_s). At the finest scale, a leaf-level model of photosynthesis is used to simulate canopy properties per species and their interaction with the variability in radiation, temperature and vapour pressure deficit. Then, a top-down plot-level model of forest productivity is used to simulate landscape-level productivity by capturing the between-stand variability in forest cover. Results show that the predicted TRGI from the empirical models account for up to 56.3% of the variance in the observed TRGI over the period 1912–1999. Under a 2 × CO₂ scenario, the predicted impact of climate change is a radial growth decline for all three species under study. However, projections obtained from the process-based model suggest that an increasing growing season length in a changing climate could counteract and potentially overwhelm the negative influence of increased drought stress. The divergence between TRGI and NPP_s simulations likely resulted, among others, from assumptions about soil water holding capacity and from calibration of variables affecting gross primary productivity. An attempt was therefore made to bridge the gap between the two modelling approaches by using physiological variables as TRGI predictors. Results obtained in this manner are similar to those obtained using climate variables, and suggest that the positive effect of increasing growing season length would be counteracted by increasing summer temperatures. Notwithstanding uncertainties in these simulations (CO₂ fertilization effect, feedback from disturbance regimes, phenology of species, and uncertainties in future CO₂ emissions), a decrease in forest productivity with climate change should be considered as a plausible scenario in sustainable forest management planning of the DMPF.     

Agency and structure: a grounded theory approach to explain land-use change in the Mindo and western foothills of Pichincha,Ecuador

ABSTRACT

Although life and land decisions are individual, driven by perceptions of reality, they reflect broader social processes. This research aims to understand relevant land-use change processes and the context within which land-use change occurs in the study area. For this, we employ grounded theory techniques and procedures to analyze narratives and life history interviews. Based on these narratives, we re-construct past land-use changes. Additionally, we identify structural conditions that drive change, several dimensions of change, including cultural-cognitive dimensions, and future discourses. The identified structural conditions motivate changes in concepts, actions, and practices, in land-use, and institutions, eventually leading to generational changes. Further research is needed to examine how the so-called structural conditions producing change varies in different settings and contexts. These findings can provide insight into certain patterns and knowledge that may contribute to community planning, policy design, and the conception of sustainable solutions with more grounded knowledge.     

Tradeoffs between forestry resource and conservation values under alternate policy regimes: A spatial analysis of the western Canadian boreal plains

An important element of resource management and conservation is an understanding of the tradeoffs between marketed products, such as timber, and measures of environmental quality, such as biodiversity. In this paper, we develop an integrated economic-ecological spatial optimization model that we then apply to evaluate alternate forest policies on a 560,000 km² study region of managed boreal forest in Alberta and British Columbia, Canada. The integrated model incorporates dynamic forest sector harvesting, current levels of oil and gas sector development, coarse-filter or habitat-based old forest indicators, a set of empirical forest bird abundance models, and statistical models of the natural and current fire regimes. Using our integrated model, economic tradeoff curves, or production possibility frontiers, are developed to illustrate the cost of achieving coarse-filter targets by a set time (50 years) within a 100-year time horizon. We found levels of ecological indicators and economic returns from the timber industry could both be increased if spatial constraints imposed by the current policy environment were relaxed; other factors being equal, this implies current policy should be revised. We explore the production possibility frontier's relationship to the range of natural variation of old forest habitat, and show how this range can be used to guide choices of preferred locations along the frontier. We also show that coarse-filter constraints on the abundance of certain habitat elements are sufficient to satisfy some fine-filter objectives, expressed as the predicted abundances of various species of songbirds.     

Tracking bioaccumulation in aquatic organisms: A dynamic model integrating life history characteristics and environmental change

We have developed a dynamic model to track the evolution of contaminant concentration in an aquatic organism as a function of season and ontogeny throughout its life cycle. We have focused our analysis on the round goby (Apollonia melanostoma), a globally distributed invasive forage fish. By integrating bioenergetics with a bioaccumulation model, we illustrate how life history characteristics interact to influence contaminant accumulation. We use uncertainty and sensitivity analyses to assess how the model output is affected by uncertainty and variability in model parameters. We then demonstrate the influence of important physiological characteristics on contaminant accumulation with two scenarios. First, we probe the influence of sexual dimorphism by comparing gender-specific accumulation of a standard polychlorinated biphenyl congener, PCB153, in male and female round gobies. We hypothesize that lipid loss in female gobies during spawning season leads to a decrease in the PCB body burden compared to male gobies. Second, we compare PCB accumulation in the round goby and in the mottled sculpin (Cottus bairdi), the native forage fish that the round goby displaced in southern Lake Michigan, to determine whether the invasive species has an intrinsically different bioaccumulation potential than the native one. Our non-intuitive findings from these simulations illustrate how the interaction of growth rate with other life history characteristics lead to unexpected bioaccumulation patterns. The model we present is a flexible tool that integrates complex and dynamic interactions among environmental parameters, thus providing a means to better assess the potential for chemical accumulation in human and wildlife populations, and aiding the development of ecological forecasts.