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.     

Modeling the carbon dynamics of the La Grande hydroelectric complex in northern Quebec

A mechanistic semi-empirical carbon cycle model of the La Grande reservoir complex in northern Quebec, Canada was conceived in order to investigate the climate impact of such a large alteration of the continental water cycle. The model includes inputs from the drainage basin, organic matter release from flooded soils, CO₂ emissions across the water-atmosphere interface and sedimentation. Most input data stems from previous research by our group on those ecosystems. The model includes the seven reservoirs of the La Grande complex and was run for periods of 50 and 100 years. Terrigeneous dissolved, particulate and suspended soil carbon fluxes and concentrations were computed. Over 100 years, 31.3 × 10¹² g C are released from flooded soils, equivalent to 28-29% of inputs from the drainage basin. 40-74% of dissolved organic carbon is mineralized. CO₂ fluxes over 100 years are 50.5-79.8 × 10¹² g C, 46.4-67.9 × 10¹² g C more than in the absence of reservoirs. The increase in mineralization of organic matter and in CO₂ emissions is a result of the increase in cumulated water residence time due to the creation of the reservoirs. Changes in other carbon sinks and sources likely offset a part of this additional carbon flux to the atmosphere. In the first years following flooding of the reservoir, organic carbon release from flooded soils exceeds CO₂ emissions, implying the downstream export of large quantities of eroded soil organic carbon. After this initial period, CO₂ emissions are fuelled by organic carbon originating from the drainage basin.     

Individual-based modeling of flooding and salinity effects on a coastal swamp forest

Coastal swamps are among the rapidly vanishing wetland habitats in Louisiana. Increased flooding, nutrient and sediment deprivation, and salt-water intrusion have been implicated as probable causes of the decline of coastal swamps. We developed a two-species individual-based forest succession model to compare the growth and composition of a cypress-tupelo swamp under various combinations of flooding intensity and salinity levels, using historical time-series of stage and salinity data as inputs. Our model simulates forest succession over 500 years by representing the growth, mortality, and reproduction of individual Taxodium distichum (baldcypress) and Nyssa aquatica (water tupelo) trees in a 1-km² spatial grid of 10 m × 10 m cells that vary in water levels and salinity through differences in elevation. We independently adjusted the elevations of each cell to obtain different grid-wide mean elevations and standard deviations of elevation; this affected the temporal and spatial pattern of flooding. We calibrated the model by adjusting selected parameters until averaged basal area, stem density and wood production rates under two different mean elevations (partially versus highly flooded) were qualitatively similar to comparable values reported for swamps in the literature. Corroboration involved comparing model predictions to four well-monitored contrasting habitat sites within the Maurepas Basin, Louisiana, USA. Model predictions of both species combined showed the same patterns among sites as the data, but the model overestimated wood production and the dominance of T. distichum. Exploratory simulations predicted that increased flooding leads to swamps with reduced basal areas and stem densities, while increased salinity resulted in lower basal areas at low salinity concentration (∼1-3 psu) and complete tree mortality at higher salinity concentrations (∼2-6 psu). Our model can provide insight into the succession dynamics of coastal swamps and information for the effective design of restoration actions.     

Calculating ecological carrying capacity of shellfish aquaculture using mass-balance modeling: Narragansett Bay, Rhode Island

Increasing growth in the aquaculture industry demands ecosystem-based techniques for management if that growth is to be ecologically sustainable and promote equity among users of the ecosystems in which it occurs. Models of carrying capacity can be used to responsibly limit the growth of aquaculture in increasingly crowded coastal areas. Narragansett Bay, Rhode Island, USA is one such crowded coastal region experiencing a rapid increase in bivalve aquaculture. An ecosystem mass-balance model was used to calculate the ecological carrying capacity of bivalve aquaculture. Cultured oyster biomass is currently at 0.47 t km⁻² and could be increased 625 times without exceeding the ecological carrying capacity of 297 t km⁻². This translates to approximately 38,950 t of harvested cultured oysters annually which is 4 times the total estimated annual harvest of finfish. This potential for growth is due to the high primary productivity and large energy throughput to detritus of this ecosystem. Shellfish aquaculture has potential for continued growth and is unlikely to become food limited due, in part, to the large detritus pool.     

Interpreting spatial heterogeneity of crop yield with a process model and remote sensing

A process-based crop growth model (Vegetation Interface Processes (VIP) model) is used to estimate crop yield with remote sensing over the North China Plain. Spatial pattern of the key parameter—maximum catalytic capacity of Rubisco (V_cmax) for assimilation is retrieved from Normalized Difference of Vegetation Index (NDVI) from Terra-MODIS and statistical yield records. The regional simulation shows that the agreements between the simulated winter wheat yields and census data at county-level are quite well with R² being 0.41-0.50 during 2001-2005. Spatial variability of photosynthetic capacity and yield in irrigated regions depend greatly on nitrogen input. Due to the heavy soil salinity, the photosynthetic capacity and yield in coastal region is less than 50 μmol C m⁻² s⁻¹ and 3000 kg ha⁻¹, respectively, which are much lower than that in non-salinized region, 84.5 μmol C m⁻² s⁻¹ and 5700 kg ha⁻¹. The predicted yield for irrigated wheat ranges from 4000 to 7800 kg ha⁻¹, which is significantly larger than that of rainfed, 1500-3000 kg ha⁻¹. According to the path coefficient analysis, nitrogen significantly affects yield, by which water exerts noticeably indirect influences on yield. The effect of water on yield is regulated, to a certain extent, by crop photosynthetic capacity and nitrogen application. It is believed that photosynthetic parameters retrieved from remote sensing are reliable for regional production prediction with a process-based model.     

Modeling net primary production of a fast-growing forest using a light use efficiency model

As interest grows in the quantification of global carbon cycles, Light Use Efficiency (LUE) model predictions of the forest net primary production (NPP) are being developed at an accelerating rate. Such models can provide useful predictions at large scales, but evaluating their performance has been difficult. In this study, a remote sensing-based LUE model was established to estimate forest NPP. Using the forest inventory data (FID) from the regional forest inventory survey in China and established allometric biomass equations, we calculated the biomass, the biomass increment, and the NPP of Eucalyptus urophylla (E. urophylla) plantation plots in the forestry jurisdiction of the Leizhou Forestry Bureau, Southern China. The FID-based NPP and the NPP from LUE model predictions were then compared to each other. Results show that the NPP from model predictions at a spatial resolution of 30 m × 30 m varied from 0 to 265 gC/(m² month) and showed regional differences. In addition, the stand age had variable effects on the average individual biomass of the E. urophylla plantation plots. The average individual biomass of the young and mid-age forests increased exponentially and logarithmically with the stand age (R² = 0.9178 and R² = 0.8683), respectively. For young and mid-age E. urophylla plantation plots, the LUE model-predicted NPP was fairly consistent with the FID-based NPP, but the model predictions of the NPP were higher than the estimates from FID. Through the analysis of the causes of uncertainty and the possible reasons for the discrepancy between the model-based NPP and FID-based NPP, the FID-derived estimates provided a foundation for model evaluation.     

Application of a Random Forest algorithm to predict spatial distribution of the potential yield of Ruditapes philippinarum in the Venice lagoon, Italy

We present a modelling framework that combines machine learning techniques and Geographic Information Systems to support the management of an important aquaculture species, Manila clam (Ruditapes philippinarum). We use the Venice lagoon (Italy), the first site in Europe for the production of R. philippinarum, to illustrate the potential of this modelling approach. To investigate the relationship between the yield of R. philippinarum and a set of environmental factors, we used a Random Forest (RF) algorithm. The RF model was tuned with a large data set (n = 1698) and validated by an independent data set (n = 841). Overall, the model provided good predictions of site-specific yields and the analysis of marginal effect of predictors showed substantial agreement among the modelled responses and available ecological knowledge for R. philippinarum. The most influent environmental factors for yield estimation were percentage of sand in the sediment, salinity, and water depth. Our results agree with findings from other North Adriatic lagoons. The application of the fitted RF model to continuous maps of all the environmental variables allowed estimates of the potential yield for the whole basin. Such a spatial representation enabled site-specific estimates of yield in different farming areas within the lagoon. We present a possible management application of our model by estimating the potential yield under the current farming distribution and comparing it to a proposed re-organization of the farming areas. Our analysis suggests a reduction of total yield is likely to result from the proposed re-organization.     

Calculation of the unit emergy value of water in an Italian watershed

Emergy is a thermodynamics-based entity that enables the implementation of a holistic environmental accounting system. It contributes to identify and measure all the inputs (energy and matter) supporting a given system, expressed in a common unit, namely solar emergy joule (sej). The emergy per unit product (called unit emergy value, UEV), is a measure of the environmental cost of a given resource. It is specific of the system/process and gives information on the dynamics, components and functioning of it. This paper presents the emergy evaluation of water resources within the watershed of the river Sieve, located in the Province of Florence (Italy). Along the river, an artificial basin has been created by means of a dam to preserve water quantity and quality, and to protect the Florentine area from dangerous floods and inundations. Different UEVs of water can be identified along the course of the river, especially upstream and downstream of the dam. These values quantify both the environmental and human efforts made to ensure and regulate the presence of water at different points of the river. The UEV of water flowing in the river increases from 1.35 × 10⁵ sej/g upstream, to 5.80 × 10⁵ sej/g downstream of the dam, depending mainly on man-made infrastructure. Along the watershed, three different systems of extraction, purification and distribution of water have been chosen on the basis of their dimension, type and location. UEVs of water distributed and the emergy investment necessary to implement different water management strategies are presented. The value of water purified and distributed decreases from 2.00 × 10⁶ sej/g for the smallest plant in the mountainous area, to 1.72 × 10⁶ sej/g for the largest plant, in the city of Florence, depending on production efficiency.     

Increasing accuracy of dispersal kernels in grid-based population models

Dispersal kernels in grid-based population models specify the proportion, distance and direction of movements within the model landscape. Spatial errors in dispersal kernels can have large compounding effects on model accuracy. Circular Gaussian and Laplacian dispersal kernels at a range of spatial resolutions were investigated, and methods for minimizing errors caused by the discretizing process were explored. Kernels of progressively smaller sizes relative to the landscape grid size were calculated using cell-integration and cell-center methods. These kernels were convolved repeatedly, and the final distribution was compared with a reference analytical solution. For large Gaussian kernels (σ > 10 cells), the total kernel error was <10⁻¹¹ compared to analytical results. Using an invasion model that tracked the time a population took to reach a defined goal, the discrete model results were comparable to the analytical reference. With Gaussian kernels that had σ ≤ 0.12 using the cell integration method, or σ ≤ 0.22 using the cell center method, the kernel error was greater than 10%, which resulted in invasion times that were orders of magnitude different than theoretical results. A goal-seeking routine was developed to adjust the kernels to minimize overall error. With this, corrections for small kernels were found that decreased overall kernel error to <10⁻¹¹ and invasion time error to <5%.     

Modeling amphibian energetics, habitat suitability, and movements of western toads, Anaxyrus (=Bufo) boreas, across present and future landscapes

Effective conservation of amphibian populations requires the prediction of how amphibians use and move through a landscape. Amphibians are closely coupled to their physical environment. Thus an approach that uses the physiological attributes of amphibians, together with knowledge of their natural history, should be helpful. We used Niche Mapper™ to model the known movements and habitat use patterns of a population of Western toads (Anaxyrus (=Bufo) boreas) occupying forested habitats in southeastern Idaho. Niche Mapper uses first principles of environmental biophysics to combine features of topography, climate, land cover, and animal features to model microclimates and animal physiology and behavior across landscapes. Niche Mapper reproduced core body temperatures (T_c) and evaporation rates of live toads with average errors of 1.6 ± 0.4 °C and 0.8 ± 0.2 g/h, respectively. For four different habitat types, it reproduced similar mid-summer daily temperature patterns as those measured in the field and calculated evaporation rates (g/h) with an average error rate of 7.2 ± 5.5%. Sensitivity analyses indicate these errors do not significantly affect estimates of food consumption or activity. Using Niche Mapper we predicted the daily habitats used by free-ranging toads; our accuracy for female toads was greater than for male toads (74.2 ± 6.8% and 53.6 ± 15.8%, respectively), reflecting the stronger patterns of habitat selection among females. Using these changing to construct a cost surface, we also reconstructed movement paths that were consistent with field observations. The effect of climate warming on toads depends on the interaction of temperature and atmospheric moisture. If climate change occurs as predicted, results from Niche Mapper suggests that climate warming will increase the physiological cost of landscapes thereby limiting the activity for toads in different habitats.     

Multi-metric evaluation of the models WARM,CropSyst, and WOFOST for rice

WARM (Water Accounting Rice Model) simulates paddy rice (Oryza sativa L.), based on temperature-driven development and radiation-driven crop growth. It also simulates: biomass partitioning, floodwater effect on temperature, spikelet sterility, floodwater and chemicals management, and soil hydrology. Biomass estimates from WARM were evaluated and compared with the ones from two generic crop models (CropSyst, WOFOST). The test-area was the Po Valley (Italy). Data collected at six sites from 1989 to 2004 from rice crops grown under flooded and non-limiting conditions were split into a calibration (to estimate some model parameters) and a validation set. For model evaluation, a fuzzy-logic based multiple-metrics indicator (MQI) was used: 0 (best) ≤ MQI ≤ 1 (worst). WARM estimates compared well with the actual data (mean MQI = 0.037 against 0.167 and 0.173 with CropSyst and WOFOST, respectively). On an average, the three models performed similarly for individual validation metrics such as modelling efficiency (EF > 0.90) and correlation coefficient (R > 0.98). WARM performed best in a weighed measure of the Akaike Information Criterion: (worst) 0<wk<1$0<w_k<1$ (best), considering estimation accuracy and number of parameters required to achieve it (mean wk=0.983$w_k=0.983$ against 0.007 and ∼0.000 with CropSyst and WOFOST, respectively). WARM results were sensitive to 30% of the model parameters (ratio being lower with both CropSyst, <10%, and WOFOST, <20%), but appeared the easiest model to use because of the lowest number of crop parameters required (10 against 15 and 34 with CropSyst and WOFOST, respectively). This study provides a concrete example of the possibilities offered using a range of assessment metrics to evaluate model estimates, predictive capabilities, and complexity.     

Structural changes in lake functioning induced from nutrient loading and climate variability

Climate variability is increasingly recognized as an important regulatory factor, capable of influencing the structural properties of aquatic ecosystems. Lakes appear to be particularly sensitive to the ecological impacts of climate variability, and several long time series have shown a close coupling between climate, lake thermal properties and individual organism physiology, population abundance, community structure, and food web dynamics. Thus, understanding the complex interplay among meteorological forcing, hydrological variability, and ecosystem functioning is essential for improving the credibility of model-based water resources/fisheries management. Our objective herein is to examine the relative importance of the ecological mechanisms underlying plankton seasonal variability in Lake Washington, Washington State (USA), over a 35-year period (1964–1998). Our analysis is founded upon an intermediate complexity plankton model that is used to reproduce the limiting nutrient (phosphate)–phytoplankton–zooplankton–detritus (particulate phosphorus) dynamics in the lake. Model parameterization is based on a Bayesian calibration scheme that offers insights into the degree of information the data contain about model inputs and allows obtaining predictions along with uncertainty bounds for modeled output variables. The model accurately reproduces the key seasonal planktonic patterns in Lake Washington and provides realistic estimates of predictive uncertainty for water quality variables of environmental management interest. A principal component analysis of the annual estimates of the underlying ecological processes highlighted the significant role of the phosphorus recycling stemming from the zooplankton excretion on the planktonic food web variability. We also identified a moderately significant signature of the local climatic conditions (air temperature) on phytoplankton growth (r = 0.41), herbivorous grazing (r = 0.38), and detritus mineralization (r = 0.39). Our study seeks linkages with the conceptual food web model proposed by Hampton et al. [Hampton, S.E., Scheuerell, M.D., Schindler, D.E., 2006b. Coalescence in the Lake Washington story: interaction strengths in a planktonic food web. Limnol. Oceanogr. 51, 2042–2051.] to emphasize the “bottom-up” control of the Lake Washington plankton phenology. The posterior predictive distributions of the plankton model are also used to assess the exceedance frequency and confidence of compliance with total phosphorus (15 μg L⁻¹) and chlorophyll a (4 μg L⁻¹) threshold levels during the summer-stratified period in Lake Washington. Finally, we conclude by underscoring the importance of explicitly acknowledging the uncertainty in ecological forecasts to the management of freshwater ecosystems under a changing global environment.     

A discrete model for estimating the development time from egg to infecting larva of Ostertagia ostertagi parametrized using a fuzzy rule-based system

Ostertagia ostertagi is a nematode, predominantly affecting cattle in the Pampean region of Argentina. A mathematical model parametrized using fuzzy rule-based systems of the Takagi-Sugeno-Kant type (FTSK) for estimating the development time from egg to infecting larval stage L3 of the gastrointestinal parasite O. ostertagi is here proposed. The estimation of development time of O. ostertagi is essential for the generation of appropriate control mechanisms, since this provides information about the time when parasites are ready to migrate to pastures. For the purpose of reflecting the natural environmental conditions, the mean daily temperature is taken as the main and only regulator of the development time. Humidity conditions are considered to be sufficient for the normal development of the larvae. Hence the individual's daily growth is a function of its length and the mean temperature recorded on the previous day. It is expressed in terms of a difference equation with fuzzy parameters, which are defined using laboratory data. Model outputs are tested against results of field experiments. Simulation results are very satisfactory, yielding a mean estimation error (MEE) of 0.64 weeks, with variance 0.34, and a determination coefficient R² = 0.74. The model clearly exhibits an inverse relationship between development time and temperature both in controlled and in field conditions. It also exhibits a very sensitive response both to the order in which the temperature sequence occurs, - reproducing the differences observed between spring and autumn - and to the amplitude of the temperature range.     

Mechanisms of moisture stress in a mid-latitude temperate forest: Implications for feedforward and feedback controls from an irrigation experiment

While it is well established that stomata close during moisture stress, strong correlations among environmental (e.g., vapor pressure deficit, soil moisture, air temperature, radiation) and internal (e.g., leaf water potential, sap flow, root-shoot signaling) variables obscure the identification of causal mechanisms from field experiments. Models of stomatal control fitted to field data therefore suffer from ambiguous parameter identification, with multiple acceptable (i.e., nearly optimal) model structures emphasizing different moisture status indicators and different processes. In an effort to minimize these correlations and improve parameter and process identification, we conducted an irrigation experiment on red maples (Acer rubrum L.) at Harvard Forest (summers of 2005 and 2006). Control and irrigated trees experienced similar radiative and boundary layer forcings, but different soil moisture status, and thus presumably different diurnal cycles of internal leaf water potential. Measured soil moisture and atmospheric forcing were used to drive a transient tree hydraulic model that incorporated a Jarvis-type leaf conductance in a Penman–Monteith framework with a Cowan-type (resistance and capacitance) tree hydraulic representation. The leaf conductance model included dependence on both leaf matric potential, Ψ_L (so-called feedback control) and on vapor pressure deficit, D (so-called feedforward control). Model parameters were estimated by minimizing the error between predicted and measured sap flow. The whole-tree irrigation treatment had the effect of elevating measured transpiration during summer dry-downs, demonstrating the limiting effect that subsurface resistance may have on transpiration during these times of moisture stress. From the best fitted model, we infer that during dry downs, moisture stress manifests itself in an increase of soil resistance with a resulting decrease in Ψ_L, leading to both feedforward and feedback controls in the control trees, but only feedforward control for the irrigated set. Increases in the sum-of-squares error when individual model components were disabled allow us to reject the following three null hypotheses: (1) the f(D) stress is statistically insignificant (p = 0.01); (2) the f(Ψ_L) stress is statistically insignificant (p = 0.07); and (3) plant storage capacitance is independent of moisture status (p = 0.07).     

On the use of growth rate parameters for projecting population sizes: Application to aphids

This paper extends the application of the cumulative size based mechanistic model, which has previously been shown to describe diverse aphid population size data well. The mechanistic model is reviewed with a focus on the explanatory role of the birth and death rate formulation. An analysis of two data sets, one on the mustard aphid and the other on the pecan aphid, indicates that multiple linear regression equations based on the estimated birth and death rate parameters alone account for nearly all (R² > 0.95) of the variability in two key population attributes, namely the peak count and the cumulative density. This indicates that population size variables may be projected directly from the growth rate parameters using linear equations. Such linear relationships based on the birth and death rate parameters are shown to hold also for certain generalized mechanistic models for which the analytical solution is not available. The birth and death rate coefficients, therefore, constitute a new succinct set of variables that could be included in the predictive modeling of aphid populations, as well as other insect and animal populations with local collapse which follow similar growth dynamics.     

Estimating Plecoglossus altivelis altivelis migration using a mass balance model expressed by hydrological distribution parameters in a major limpid river basin in Japan

We examined the influence of several hydrological and meteorological parameters on the migratory movements of ayu Plecoglossus altivelis altivelis in central Japan. When comprehensively evaluating rivers and ayu behaviour on a catchment scale, the subjects of analysis typically include human activities and hydrological and meteorological phenomena. However, limiting analyses to such factors may be too restrictive when human activities are being conducted. Accordingly, we incorporated a biological viewpoint into the evaluation method, analysing hydrological data (river discharge, river water temperature, sea water temperature) to determine watershed characteristics and examining the relationship between these characteristics and the habitat conditions of ayu. Then we constructed a numerical model for ayu migratory runs that incorporated ayu ecology and watershed characteristics. Analyses of ayu movements from a lower estuarine dam demonstrated that downstream displacements were associated with high water flows of more than 200 m³ s⁻¹ at the beginning of summer. We conclude that it is important to consider the effects of environmental parameters on the movements of different fish species to understand the causes of spatial variation in fish distribution in lowland rivers.     

Formulation and testing of a novel river nitrification model

The nitrification process in many river water quality models has been approximated by a simple first order dependency on the water column ammonia concentration, while the benthic contribution has routinely been neglected. In this study a mathematical framework was developed for sediment bed nitrification based on mass transfer theory and Monod bacterial growth kinetics. The model describes ammonia transport across the boundary layer and consumption within the biofilm to quantify the overall nitrification flux. Model results suggest that nitrification is usually controlled by the boundary layer thickness, and we estimated a nitrification velocity range between 0.14 and 0.97 m d⁻¹, assuming typical boundary thicknesses of 0.1–1.0 mm. These ranges compared favorably with reported literature values, including our own measurements. The model was applied to several river systems of different depths where nitrification rates and river depths were available. Assuming that nitrification is exclusively a benthic process, the average velocity of all the rivers evaluated was 0.85 m d⁻¹ (r² = 0.72).     

Carrying capacity of an uninhabited island off the southwestern coast of Korea

An emergy evaluation was carried out to assess the carrying capacity of a small, uninhabited island (Woosedo) off the southwestern coast of Korea. The sea area within 1 km from the high tide level of the island was included in the evaluation. The total environmental emergy input to Woosedo was 1.66E19 sej/yr, with the most emergy contribution from the tidal energy. The land and marine ecosystems of Woosedo contributed 4.97 million Em$ (7600 Em$/ha/yr) to the Korean economy annually. If Woosedo was developed to the national average at the emergy investment ratio of 2.86, its carrying capacity was estimated at 1034 people at the current living standard of Korea. With this population, the island system would not be sustainable with a very low emergy sustainability index of 0.36. At the same living standard used in the developed scenario, the carrying capacity of the island would be 370 people for a sustainable development scenario and 270 people if the renewable emergy were the only source to support the population. The emergy contribution of the marine ecosystem of the island was the major source of support in determining the level of carrying capacity of the island.     

Modelling potential habitat for cougars in midwestern North America

Cougars (Puma concolor) are of considerable interest to wildlife biologists and the general public in midwestern North America, yet no researchers have modelled potential habitat in the region. We created a model of potential cougar habitat in 9 midwestern states using geospatial data, expert-opinion surveys, the analytical hierarchy process, and a GIS. About 8% of the study region contained highly favorable habitat (with favorability scores ≥75%) for cougars; the states of Arkansas (19%) and Missouri (16%) contained the highest proportions of potentially favorable habitat. We identified 6 large (≥2500 km² in size), contiguous areas of highly favorable habitat for cougars. Model testing procedures indicated a valid model when compared to an independent set of cougar locations, a null dataset, and similar studies. Our model is useful as a planning tool to proactively address future human-cougar conflicts should cougars re-colonize the Midwest via subadult dispersal.