Somatic condition of the fish, Notopterus notopterus (Pallas) during different phases of the reproductive cycle

The study on the mathematical relationship between length and weight of fishes has considerable importance in fishery research especially for the study of fish population dynamics and patterns of growth in fish stocks. This aspect is useful in rationally managing the population of fishes, and is an index of relative robustness of the fish). Variations from expected weight or length of an individual fish or group of fish indicate the general well being and state of gonadal development. The condition of fish can be influenced by the onset of maturity. The results indicate that the overall condition of the female and male fish seems to get improved towards the approach of maturing and mature stages of the gonads during pre-spawning phases. The reduction in the condition factor (K) and somatic condition factor (Ks) during spawning phase may be because of the fish becoming exhausted due to spawning activity. Although improvement in the condition of the fish was observed on the approach of breeding period, the overall condition seems to be better in the aquatic body from which the fish was collected.     

Estimating the functional form for the density dependence from life history data

Two contrasting approaches to the analysis of population dynamics are currently popular: demographic approaches where the associations between demographic rates and statistics summarizing the population dynamics are identified; and time series approaches where the associations between population dynamics, population density, and environmental covariates are investigated. In this paper, we develop an approach to combine these methods and apply it to detailed data from Soay sheep (Ovis aries). We examine how density dependence and climate contribute to fluctuations in population size via age- and sex-specific demographic rates, and how fluctuations in demographic structure influence population dynamics. Density dependence contributes most, followed by climatic variation, age structure fluctuations and interactions between density and climate. We then simplify the density-dependent, stochastic, age-structured demographic model and derive a new phenomenological time series which captures the dynamics better than previously selected functions. The simple method we develop has potential to provide substantial insight into the relative contributions of population and individual-level processes to the dynamics of populations in stochastic environments.     

Biological responses to environmental forcing: the linear tracking window hypothesis

Determining the relative contributions of intrinsic and extrinsic processes to the regulation of biological populations has been a recurrent ecological issue. Recent discussions concerning ecosystem "regime shifts" again raise the question of whether population fluctuations are mainly controlled by external forcing. Results of nonlinear time series analyses indicate that pelagic populations typically do not passively track stochastic environmental variables. Rather, population dynamics are better described as nonlinear amplification of physical forcing by biological interactions. However, we illustrate that in some cases populations do show linear tracking of the physical environment. To explain why population dynamics can sometimes be linear, we propose the linear tracking window hypothesis: populations are most likely to track the stochastic environmental forcing when their generation time matches the characteristic time scale of the environmental signal. While our observations follow this hypothesis well, our results indicate that the linear tracking window is a necessary but not a sufficient condition.     

Flow and transport characteristics at an Epinephelus guttatus (red hind grouper) spawning aggregation site in St. Thomas (US Virgin Islands)

This study characterizes the flow field at a spawning site located at the shelf break of a Caribbean island for the Epinephelus guttatus (red hind grouper) in relation to this species spawning events. In order to understand the oceanographic dynamics targeted by the fishes, current measurements were profiled throughout the water column for almost a year at the spawning site. The characteristics of the flow field and its evolution after spawning were investigated by using a numerical ocean model that resolved the observed tide and simulated the island scale flow where passive, neutrally buoyant virtual particles were released for 10 days to trace the flow pathways.Observed currents during the spawning period revealed that the flow was vertically sheared, to the south and weakest at the bottom, and to the west or east at the surface. The tidal analysis revealed that the flow at the time of spawning was directed across and on-shelf, although weaker close to the bottom. The model showed that the initial on-shelf transport was counteracted by the bottom flow directed to the shelf break, where virtual particles were entrained by the downwelling flow. A significant percent of particles resided less than two hundred meters deep, in the vicinity of the chlorophyll maximum and returned to the shelf break, close to the release location within 8-10 days. This journey was largely controlled by the timing between downwelling at the spawning site and upwelling further east at the shelf break, which was driven by the coupling between wind and tide induced vertical movements at the shelf break and deeper. The release location, vertical rotation of its flow field, and its transport properties were shown to be relatively resilient to the passage of transient sub-mesoscale eddies as well as to acute mesoscale flow reversals, suggesting that physical retention is maximized in the area surrounding the spawning site.     

Modelling the effect of constant and fluctuating food supply on egg production rates of Acartia grani

A mathematical model of the individual budget of a spawning female of the copepod Acartia grani (Sars) has been used to simulate the time-scale of egg production over various external forcings (or inputs) of food fluctuation conditions. The budget matter in the body of the copepod females is distributed through four compartments: the whole digestive tract (globally named as gut), the hemolymph (which include the body fluid with available nutriments for the organs), the structural body weight, and the gonad. This small calanoid species does not carry lipid reserves but cumulate some labile reserves in its body, according to food availability. The model results show how the continuous spawning varies with food fluctuations, and suggest the mechanisms inducing the delay of response to starvation by using the metabolic reserves. Three different patterns in egg production response are observed: food fluctuations with frequencies below 12 h have no effect on egg production; food fluctuations of 12 h to 5 days induce synchronous egg production fluctuations; beyond 5 days the strong physiological changes induced by long starvation durations create delays in the responses to food replenishment. The available data of cultivated cohorts under laboratory conditions are used to validate the model. The properties underlined by the model, in particular its weak capacity to respond to starvation, allow explaining A. grani distribution in specific habitats. Different experimental protocols for complementary experiments are proposed to complete the model validation in other forcing conditions.     

The pattern of species turnover resulting from stochastic population dynamics: The model and field data

The model of random population dynamics in a sampling site returns geometric distribution of longevities of continuous presence (=persistence) and Poisson distribution of the presence–absence transitions. This discrete-time stochastic process describes the presence–absence pattern observed in the beetles surveyed 6 years on Mount Carmel, Israel. Homogeneous pools of species mostly on the Families rank, exhibit the predicted by the model patterns. Conformity to an ergodic hypothesis is the criterion of ecological homogeneity. This criterion assumes the equivalence of short-term behavior of entire pool and long-term behavior of any species from this pool. The pool of all 801 species of Order Coleoptera does not match the model. Thus a taxon of an arbitrary rank may not be considered a priory as a unit of ecological study. Determined from field data parameters of the model are biased and magnitude of the bias depends on longevity of the survey. Parameter of distribution depends also on species tolerance, which is the level adaptation of given species to given environment in given time interval. Random process of species turnover may be considered as a game of species to gain their presence against deteriorative fluctuations of environmental conditions.     

Meso-scale hydrodynamic and reproductive asynchrony affects the source–sink metapopulation structure of the coastal polychaete Pectinaria koreni

The polychaete Pectinaria koreni exhibits a complex life cycle characterized by non-overlapping generations and widespread larval dispersal. To explore how “local” metapopulation genetic structure varies spatially and temporally during population turnover, we combined observations on demography, larval dispersal through hydrodynamic modelling and population genetics of successive age cohorts in the Baie de Seine (eastern English Channel, France). Mature adults (March), newly settled (July) and later-stage juveniles (September) were sampled in 2003 on the edge and in the main demes of the metapopulation. Demes displayed an asynchronous dynamics due to variations in habitat quality affecting reproductive timing (e.g. three distinct spawning events observed) and in local larval supply linked to temporal fluctuations of hydrodynamism. Two-source populations were identified among dense areas with the greatest larval retention and self-recruitment rates: one with a single recruitment event, stable temporal genetic variation and a strong spatial genetic re-homogenization during turnover, and the other with two recruitment events and significant allele frequency changes over time. Sink habitats displayed single recruitment event and experienced strong inter-generational (juveniles vs adults) genetic changes due to genetic drift associated with strong winter mortalities. Altogether, results suggested that adult spawning asynchrony and sweepstakes reproductive success, together with genetic drift, played a greater role than hydrodynamics itself in determining effective recruitment rates at some sites and generating genetic patchiness.     

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.     

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.     

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.     

Analysis of lipid components for determining the condition of anchovy larvae,Engraulis mordax

It has been presumed that intertidal spawning by Limulus polyphemus minimizes the loss of egges to subtidal predators; however, this strategy involves considerable risks. Massive beach strandings of adults accompany seasonal spawning migrations of crabs along Cape May in Delaware Bay, (USA). At least 190000 horseshoe crabs, approximating 10% of the adult population, died from beach stranding along the New Jersey shore of Delaware Bay during the 1986 (May to June) spawning season. Abnormalities of the telson (which is used in righting behavior) were significantly more common among stranded crabs than among individuals actively spawning on the intertidal beach. The number of stranded crabs per day was not correlated with tidal height or environmental variables (wind speed, wave height) which characterized the conditions at spawning. A complex suite of factors, including the size of the available spawning population, tidal and weather conditions, and beach slope, influence the number stranded during the breeding season. Horseshoe crab stranding results in a large loss of gravid females from the population, and may represent a major input of organic matter to intertidal sandy beaches in certain regions of Delaware Bay.     

Individual based model of slug population and spatial dynamics

The slug, Deroceras reticulatum, is one of the most important pests of agricultural and horticultural crops in UK and Europe. In this paper, a spatially explicit individual based model (IbM) is developed to study the dynamics of a population of D. reticulatum. The IbM establishes a virtual field within which slug spatial dynamics and changes in abundance were simulated. The strong dependence of slug behaviour on environmental conditions is built into the model, which is based upon previous work on the environmental dependence of slug population dynamics. The simulation results show that the IbM described well changes in the slug population. The IbM proved capable of describing slug populations over 3.5 years, including the presence, magnitude and duration of D. reticulatum population crashes within this period. Moreover, the model was capable of reproducing slug population dynamics at two sites, with distinct weather and some 100 km apart, with minor changes in initialisation values but no change in model structure and parameter values. A study of field heterogeneity, which might simulate various field designs, indicated the importance of spatial structuring to slug population dynamics and the utility of the IbM for simulating a range of potential spatial management treatments for slug control to maximise crop yield. This IbM system performs well and is currently being used as part of an integrated approach to predict slug population dynamics and control in the UK.     

Potential latitudinal variation in egg size and number of a geographically widespread reef fish,revealed by common-environment experiments

Variation in maternal reproductive traits was examined in field and reared populations of a geographically widespread reef fish, Pomacentrus coelestis (Pomacentridae), drawn from three different latitudes in Japan. Size-specific clutch size and clutch weight of wild fish increased with increasing latitude. Conversely, latitudinal variation in egg size of wild fish was obscure in same-season comparisons, probably because of the temperature effect on egg size. Common-environment experiments conducted at three temperatures showed that egg size decreased with increasing temperature in all populations. In the experiments, egg size, clutch size and clutch weight differed among populations at all temperatures, showing clear latitudinal clines. Females from low latitude spawned larger eggs at every experimental temperature. Size-specific clutch size and weight were greater in females from high latitude. Thus, the northern fish had a larger reproductive output per spawning and a larger number of smaller eggs in a spawning. Such interpopualtion variation in this fish is likely to be partially genetically based, although environmental effects on the variation cannot be entirely ruled out. This study provides evidence of potential latitudinal variation in the egg size and number in a coastal fish, by common-environment experiments. The close correspondence between latitudes and these maternal reproductive traits may be a consequence of local adaptation.     

Spatial synchrony in coral reef fish populations and the influence of climate

We investigated spatial patterns of synchrony among coral reef fish populations and environmental variables over an eight-year period on the Great Barrier Reef, Australia. Our aims were to determine the spatial scale of intra- and interspecific synchrony of fluctuations in abundance of nine damselfish species (genus Pomacentrus) and assess whether environmental factors could have influenced population synchrony. All species showed intraspecific synchrony among populations on reefs separated by < or =100 km, and interspecific synchrony was also common at this scale. At greater spatial scales, only four species showed intraspecific synchrony, over distances ranging from 100-300 km to 500-800 km, and no cases of interspecific synchrony were recorded. The two mechanisms most likely to cause population synchrony are dispersal and environmental forcing through regionally correlated climate (the Moran effect). Dispersal may have influenced population synchrony over distances up to 100 km as this is the expected spatial range for ecologically significant reef fish dispersal. Environmental factors are also likely to have synchronized population fluctuations via the Moran effect for three reasons: (1) dispersal could not have caused interspecific synchrony that was common over distances < or =100 km because dispersal cannot link populations of different species, (2) variations in both sea surface temperature and wind speed were synchronized over greater spatial scales (>800 km) than fluctuations in damselfish abundance (< or =800 km) and were correlated with an index of global climate variability, the El Ni?o-Southern Oscillation (ENSO), and (3) synchronous population fluctuations of most damselfish species were correlated with ENSO; large population increases often followed ENSO events. We recorded regional variations in the strength of population synchrony that we suspect are due to spatial differences in geophysical, oceanographic, and population characteristics, which act to dilute or enhance the effects of synchronizing mechanisms. We conclude that synchrony is common among Pomacentrus populations separated by tens of kilometers but less prevalent at greater spatial scales, and that environmental variation linked to global climate is likely to be a driving force behind damselfish population synchrony at all spatial scales on the Great Barrier Reef.     

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.     

Can synchronous spawning be predicted from environmental parameters? A case study of the lugworm Arenicola marina

The annual epidemic spawning period of a Scottish population of Arenicola marina (L.) has been recorded over a period of 13 yr. This population spawns between mid-October and mid-November in a discrete spawning event over a period of 4 to 5 d. Endocrine manipulation experiments showed that spawning is induced in females only if sufficient titres of PMH (prostomial maturation hormone) are present in the prostomia. These levels are attained during the 2 to 3 wk prior to the natural spawning date. The East Sands, St. Andrews population always spawns during periods of spring tides regardless of tidal amplitude or whether they are full- or new-moon tides. Meteorological data, including sea-temperature data were collected for each year, and correlation of the environmental data with spawning time was attempted. Correlation of spawning times with weather patterns showed that mean daily air pressures were significantly higher during the spawning period than from September to November as a whole. Evidence also suggests that a reduction in sea temperature is required prior to spawning. A significant moderate negative correlation was found between May to July air temperatures and spawning date, suggesting that higher May to July temperatures may induce early spawning. Daily rainfall and wind speed were also lower during the spawning period, but not significantly so. These results indicate that air pressure (or changes therein) may act as a final spawning cue, and the advantages of this are discussed in relation to fertilization success. A model of the interplay between environmental parameters and the endocrine mechanisms controlling the induction of spawning is proposed. Higher than average summer temperatures may advance gametogenesis to bring the population into a state of maturity (full-size oocytes, well-developed sperm morulae), and may also advance spawning time. Once the population has completed gametogenesis, a drop in sea temperature is then required to trigger an increase in endocrine titres within the prostomium, without which spawning cannot be induced by prostomial injection. The population spawns on spring tides; however a lack of clement weather coinciding with the spring tide will result in population-wide spawning being aborted, as in 1996. Clement weather (high pressure, low rainfall and wind speed) in conjunction with spring tides permits spawning to proceed to completion. Received: 21 June 1999 / Accepted: 25 January 2000     

A hydrologically driven model of swamp water mosquito population dynamics

We develop a swamp water mosquito population model that is forced solely by environmental variability. Measured temperature and land surface wetness conditions are used to simulate Anopheles walkeri population dynamics in a northern New Jersey habitat. Land surface wetness conditions, which represent oviposition habitat availability, are derived from simulations using a dynamic hydrology model. Using only these two density-independent effects, population model simulations of biting Anoph. walkeri correlate significantly with light trap collections. These results suggest that prediction of mosquito populations and the diseases they transmit could be better constrained by inclusion of environmental variability.     

Persistent Unequal Sex Ratio in a Population of Grayling (Salmonidae) and Possible Role of Temperature Increase

In some fishes, water chemistry or temperature affects sex determination or creates sex‐specific selection pressures. The resulting population sex ratios are hard to predict from laboratory studies if the environmental triggers interact with other factors, whereas in field studies, singular observations of unusual sex ratios may be particularly prone to selective reporting. Long‐term monitoring largely avoids these problems. We studied a population of grayling (Thymallus thymallus) in Lake Thun, Switzerland, that has been monitored since 1948. Samples of spawning fish have been caught about 3 times/week around spawning season, and water temperature at the spawning site has been continuously recorded since 1970. We used scale samples collected in different years to determine the average age of spawners (for life‐stage specific analyses) and to identify the cohort born in 2003 (an extraordinarily warm year). Recent tissue samples were genotyped on microsatellite markers to test for genetic bottlenecks in the past and to estimate the genetically effective population size (N_e). Operational sex ratios changed from approximately 65% males before 1993 to approximately 85% males from 1993 to 2011. Sex ratios correlated with the water temperatures the fish experienced in their first year of life. Sex ratios were best explained by the average temperature juvenile fish experienced during their first summer. Grayling abundance is declining, but we found no evidence of a strong genetic bottleneck that would explain the apparent lack of evolutionary response to the unequal sex ratio. Results of other studies show no evidence of endocrine disruptors in the study area. Our findings suggest temperature affects population sex ratio and thereby contributes to population decline. Persistencia de Proporción de Sexos Desigual en una Población de Tímalos (Salmonidae) y el Posible Papel del Incremento de la Temperatura     

Role of the nutritive value of sediment in regulating population dynamics of the deposit-feeding polychaete Streblospio shrubsolii

The importance of changes in nutritive value of the sediment in regulating the dynamics of the small deposit-feeding polychaete Streblospio shrubsolii (Buchanam, 1890) was examined. First, fluctuations in the population dynamics of S. shrubsolii and in the nutritive content of the sediment were studied on a seasonal scale. Subsequently, a manipulative field experiment was used to test the hypothesis that, in summer, availability of food regulates aspects of the population dynamics of S. shrubsolii. In summer 1997, sharp increases in protein concentrations and in the number of individuals of S. shrubsolii were found compared to all the other seasons. Densities of S. shrubsolii also showed an increase in the subsequent season. Newly recruited juveniles (19-31 setigers) were abundant in summer and autumn 1997, whereas the number of juveniles (32-42 setigers) increased in autumn 1997. Large numbers of adults were found in summer 1997. Significant negative correlations occurred in summer 1997 between the concentrations of proteins and the total number of individuals, the number of newly recruited juveniles and the density of adults. No differences in densities of juveniles, adults or total number of individuals were found among plots before starting the experiment. After 20 days of experimental enrichment, proteins showed a tendency to increase. After 40 days of enrichment, the total number of individuals was significantly greater in the enriched plots. Large numbers of newly recruited larvae and postlarval adults were found in the enriched plots. The experiment corroborated the model that the growth of population was limited by the availability of food in summer. Thus, in this season, patches of food could have a fundamental role in determining the population dynamics of small deposit-feeders under natural conditions and over short time scales.     

Longleaf pine (Pinus palustris) and hardwood dynamics in a fire-maintained ecosystem: A simulation approach

Longleaf pine (Pinus palustris) savannas of the southeastern U.S. represent an archetype of a fire dependent ecosystem. They are known to have very short fire return intervals (∼1-3 years) that perpetuate understory plant diversity (up to 50 species m⁻²), support pine recruitment, and suppress fire sensitive hardwoods. Understanding the relationships that regulate longleaf and southern hardwoods is especially critical. With decreased fire frequency, insufficient intensity, or lack of underground competition, a woody mid-story rapidly develops, dominated by fire sensitive trees and shrubs that in-turn suppress more fire dependent species (including pine seedlings). This may occur in forest gaps, where pine-needle abundance is diminished, reducing fire spread potential. The interactions between longleaf pine, hardwoods, forest fuels, and fire frequency are complex and difficult to understand spatially. The objective of this study was to develop a spatially explicit longleaf pine-hardwood stochastic simulation model (LLM), incorporating tree demography, plant competition, and fuel and fire characteristics. Data from two longleaf pine study sites were used to develop and evaluate the model with the goal to incorporate simple site-specific calibration parameters for model versatility. Specific model components included pine seed masting, hardwood clonal sprouting, response to fire (re-sprouting, mortality), and tree density driven competition effects. LLM spatial outputs were consistent with observed forest gap dynamics associated with pine seedling establishment and hardwood encroachment. Changes in fire frequency (i.e., fire probability = 0.35-0.05) illustrated a shift in community structure from longleaf pine dominated to a hardwood dominated community. This approach to assessing model response may be useful in characterizing longleaf ecosystem resilience, especially at intermediate fire frequencies (e.g., 0.15) where the community may be sensitive to small changes in the fire regime. Height distributions and population densities were similar to in situ findings (field and LIDAR data) for both study sites. Height distributions output by the LLM illustrated fluctuations in population structure. The LLM was especially useful in determining knowledge gaps associated with fuel and fire heterogeneity, plant-plant interactions, population structure and its temporal fluctuations, and hardwood demography. This is the first known modeling work to simulate interactions between longleaf pine and hardwoods and provides a foundation for further studies on fire and forest management, especially in relation to ecological forestry practices, restoration, and site-specific applications.