Preliminary comparison of natural versus model-predicted recovery of vessel-generated seagrass injuries in Florida Keys National Marine Sanctuary

Each year, more than 500 motorized vessel groundings cause widespread damage to seagrasses in Florida Keys National Marine Sanctuary (FKNMS). Under Section 312 of the National Marine Sanctuaries Act (NMSA), any party responsible for the loss, injury, or destruction of any Sanctuary resource, including seagrass, is liable to the United States for response costs and resulting damages. As part of the damage assessment process, a cellular automata model is utilized to forecast seagrass recovery rates. Field validation of these forecasts was accomplished by comparing model-predicted percent recovery to that which was observed to be occurring naturally for 30 documented vessel grounding sites. Model recovery forecasts for both Thalassia testudinum and Syringodium filiforme exceeded natural recovery estimates for 93.1% and 89.5% of the sites, respectively. For Halodule wrightii, the number of over- and under-predictions by the model was similar. However, where under-estimation occurred, it was often severe, reflecting the well-known extraordinary growth potential of this opportunistic species. These preliminary findings indicate that the recovery model is consistently generous to Responsible Parties in that the model forecasts a much faster recovery than was observed to occur naturally, particularly for T. testudinum, the dominant seagrass species in the region and the species most often affected. Environmental setting (i.e., location, wave exposure) influences local seagrass landscape pattern and may also play a role in the recovery dynamics for a particular injury site. An examination of the relationship between selected environmental factors and injury recovery dynamics is currently underway. (PDF file contains 20 pages.)

Vascular endothelial growth factor regulates melanoma cell adhesion and growth in the bone marrow microenvironment via tumor cyclooxygenase-2

Background: Human melanoma frequently colonizes bone marrow (BM) since its earliest stage of systemic dissemination, prior to clinical metastasis occurrence. However, how melanoma cell adhesion and proliferation mechanisms are regulated within bone marrow stromal cell (BMSC) microenvironment remain unclear. Consistent with the prometastatic role of inflammatory and angiogenic factors, several studies have reported elevated levels of cyclooxygenase-2 (COX-2) in melanoma although its pathogenic role in bone marrow melanoma metastasis is unknown. Methods: Herein we analyzed the effect of cyclooxygenase-2 (COX-2) inhibitor celecoxib in a model of generalized BM dissemination of left cardiac ventricle-injected B16 melanoma (B16M) cells into healthy and bacterial endotoxin lipopolysaccharide (LPS)-pretreated mice to induce inflammation. In addition, B16M and human A375 melanoma (A375M) cells were exposed to conditioned media from basal and LPS-treated primary cultured murine and human BMSCs, and the contribution of COX-2 to the adhesion and proliferation of melanoma cells was also studied. Results: Mice given one single intravenous injection of LPS 6 hour prior to cancer cells significantly increased B16M metastasis in BM compared to untreated mice; however, administration of oral celecoxib reduced BM metastasis incidence and volume in healthy mice, and almost completely abrogated LPS-dependent melanoma metastases. In vitro, untreated and LPS-treated murine and human BMSC-conditioned medium (CM) increased VCAM-1-dependent BMSC adherence and proliferation of B16M and A375M cells, respectively, as compared to basal medium-treated melanoma cells. Addition of celecoxib to both B16M and A375M cells abolished adhesion and proliferation increments induced by BMSC-CM. TNF alpha and VEGF secretion increased in the supernatant of LPS-treated BMSCs; however, anti-VEGF neutralizing antibodies added to B16M and A375M cells prior to LPS-treated BMSC-CM resulted in a complete abrogation of both adhesion-and proliferation-stimulating effect of BMSC on melanoma cells. Conversely, recombinant VEGF increased adherence to BMSC and proliferation of both B16M and A375M cells, compared to basal medium-treated cells, while addition of celecoxib neutralized VEGF effects on melanoma. Recombinant TNFa induced B16M production of VEGF via COX-2-dependent mechanism. Moreover, exogenous PGE2 also increased B16M cell adhesion to immobilized recombinant VCAM-1. Conclusions: We demonstrate the contribution of VEGF-induced tumor COX-2 to the regulation of adhesion-and proliferation-stimulating effects of TNFa, from endotoxin-activated bone marrow stromal cells, on VLA-4-expressing

Growth, movement, and attrition of northern bluefin tuna, Thunnus thynnus, in the Pacific Ocean, as determined by tagging

ENGLISH: The growth of northern bluefin tuna is described by a two-stanza model. For fish between 191 and 564 mm in length the Gompertz curve, with values of 581 mm and 4.32 for Loo and K (annual), respectively, is used. The fish between 564 and 1530 mm grow linearly, at the rate of 0.709 mm per day. Age-O fish tagged and released in the western Pacific Ocean have been recaptured in the western, central, and eastern Pacific. The minimum time between release in the western Pacific and recapture in the eastern Pacific is 215 days. Older fish, mostly Land 2-year olds, tagged and released in the eastern Pacific have been recaptured in the eastern and western Pacific. The minimum time between release in eastern Pacific and recapture in the western Pacific is 674 days. The coefficient of natural mortality is estimated from data on growth and ambient temperature to be 0.276 on an annual basis, with 90-percent confidence limits of 0.161 and 0.47L Spawning of northern bluefin takes place only in the western Pacific. Some of the juveniles migrate to the eastern Pacific, where they reside for several months to several years before returning to the western Pacific. The portion of fish which migrate to the eastern Pacific varies among years, and this appears to be an important cause of the annual variation in the catches in the eastern Pacific Ocean. SPANISH: El crecimiento del atún aleta azul del norte es descrito por un modelo de dos estadios. Para los peces de entre 191 y 564 mm de talla se usa la curva de Gompertz, con valores de 581 mm y 4.32 para Loo y K (anual), respectivamente. Los peces de entre 564 y 1530 mm crecen de forma lineal, a 0.709 mm por día. Peces de edad Omarcados y liberados en el Pacífico occidental han sido recapturados en el Pacífico occidental, central, y oriental. La demora mínima entre la liberación en el Pacífico occidental y la recaptura en el Pacífico oriental es de 215 días. Peces mayores, principalmente de 1 ó 2 años de edad, marcados y liberados en el Pacífico oriental han sido re capturados en el Pacífico occidental y oriental. La demora mínima entre la liberación en el Pacífico oriental y la recaptura en el Pacífico occidental es de 674 días. Se estima el coeficiente de mortalidad natural a partir de los datos de crecimiento y temperatura ambiental en un 0.276 anual, con límites de confianza al 90% de 0.161 y 0.471. El aleta azul del norte desova únicamente en el Pacífico occidental. Algunos de los juveniles migran al Pacífico oriental, donde permanecen entre varios meses y varios años antes de regresar al Pacífico occidental. La porción de los peces que migran al Pacífico oriental varía entre años, y ésto parece ser una causa importante de la variación anual en las capturas en el Océano Pacífico oriental. (PDF contains 94 pages.)

Earthquakes of the Nepal Himalaya : towards a physical model of the seismic cycle

Home to hundreds of millions of souls and land of excessiveness, the Himalaya is also the locus of a unique seismicity whose scope and peculiarities still remain to this day somewhat mysterious. Having claimed the lives of kings, or turned ancient timeworn cities into heaps of rubbles and ruins, earthquakes eerily inhabit Nepalese folk tales with the fatalistic message that nothing lasts forever. From a scientific point of view as much as from a human perspective, solving the mysteries of Himalayan seismicity thus represents a challenge of prime importance. Documenting geodetic strain across the Nepal Himalaya with various GPS and leveling data, we show that unlike other subduction zones that exhibit a heterogeneous and patchy coupling pattern along strike, the last hundred kilometers of the Main Himalayan Thrust fault, or MHT, appear to be uniformly locked, devoid of any of the “creeping barriers” that traditionally ward off the propagation of large events. The approximately 20 mm/yr of reckoned convergence across the Himalaya matching previously established estimates of the secular deformation at the front of the arc, the slip accumulated at depth has to somehow elastically propagate all the way to the surface at some point. And yet, neither large events from the past nor currently recorded microseismicity nearly compensate for the massive moment deficit that quietly builds up under the giant mountains. Along with this large unbalanced moment deficit, the uncommonly homogeneous coupling pattern on the MHT raises the question of whether or not the locked portion of the MHT can rupture all at once in a giant earthquake. Univocally answering this question appears contingent on the still elusive estimate of the magnitude of the largest possible earthquake in the Himalaya, and requires tight constraints on local fault properties. What makes the Himalaya enigmatic also makes it the potential source of an incredible wealth of information, and we exploit some of the oddities of Himalayan seismicity in an effort to improve the understanding of earthquake physics and cipher out the properties of the MHT. Thanks to the Himalaya, the Indo-Gangetic plain is deluged each year under a tremendous amount of water during the annual summer monsoon that collects and bears down on the Indian plate enough to pull it away from the Eurasian plate slightly, temporarily relieving a small portion of the stress mounting on the MHT. As the rainwater evaporates in the dry winter season, the plate rebounds and tension is increased back on the fault. Interestingly, the mild waggle of stress induced by the monsoon rains is about the same size as that from solid-Earth tides which gently tug at the planets solid layers, but whereas changes in earthquake frequency correspond with the annually occurring monsoon, there is no such correlation with Earth tides, which oscillate back-and-forth twice a day. We therefore investigate the general response of the creeping and seismogenic parts of MHT to periodic stresses in order to link these observations to physical parameters. First, the response of the creeping part of the MHT is analyzed with a simple spring-and-slider system bearing rate-strengthening rheology, and we show that at the transition with the locked zone, where the friction becomes near velocity neutral, the response of the slip rate may be amplified at some periods, which values are analytically related to the physical parameters of the problem. Such predictions therefore hold the potential of constraining fault properties on the MHT, but still await observational counterparts to be applied, as nothing indicates that the variations of seismicity rate on the locked part of the MHT are the direct expressions of variations of the slip rate on its creeping part, and no variations of the slip rate have been singled out from the GPS measurements to this day. When shifting to the locked seismogenic part of the MHT, spring-and-slider models with rate-weakening rheology are insufficient to explain the contrasted responses of the seismicity to the periodic loads that tides and monsoon both place on the MHT. Instead, we resort to numerical simulations using the Boundary Integral CYCLes of Earthquakes algorithm and examine the response of a 2D finite fault embedded with a rate-weakening patch to harmonic stress perturbations of various periods. We show that such simulations are able to reproduce results consistent with a gradual amplification of sensitivity as the perturbing period get larger, up to a critical period corresponding to the characteristic time of evolution of the seismicity in response to a step-like perturbation of stress. This increase of sensitivity was not reproduced by simple 1D-spring-slider systems, probably because of the complexity of the nucleation process, reproduced only by 2D-fault models. When the nucleation zone is close to its critical unstable size, its growth becomes highly sensitive to any external perturbations and the timings of produced events may therefore find themselves highly affected. A fully analytical framework has yet to be developed and further work is needed to fully describe the behavior of the fault in terms of physical parameters, which will likely provide the keys to deduce constitutive properties of the MHT from seismological observations.

Low Climate Stabilisation under Diverse Growth and Convergence Scenarios

4 p.

Energy Supply and the Sustainability of Endogenous Growth

29 p.

The DICER Model: Methodological Issues and Initial Results

31 p.

Low Climate Stabilisation under Diverse Growth and Convergence Scenarios

9 p.

The Uniform World Model: A Methodology for Predicting the Health Impacts of Air Pollution

47 p.

A simple approach for the estimation of food consumption from growth rates at different environmental conditions and its application to juvenile cod (Gadus morhua L.) of a fjordic sea loch on the west coast of Scotland

Abstract Growth and condition of fish are functions of available food and environmental conditions. This led to the idea of using fish as a “consumption sensor” for the measurement of food intake over a defined period of time. A bio-physical model for the estimation of food consumption was developed based on the von Bertalanffy model. Whereas some of the input variables of the model, the initial and final lengths and masses of a fish and the temperature within the time period considered can easily be measured, internal characteristics of the species have to be determined indirectly. Three internal parameters are used in the model: the maintenance consumption at 0°C, the temperature dependence of this consumption and the food efficiency, the percentage of the ingested food utilized. Estimates of the parameters for a given species can be determined by feeding experiments. Here, data from published feeding experiments on juvenile cod, Gadus morhua L., were used to validate the model. The average of the relative error for the food intake predicted by the model for individual fish was about 24 %, indicating that fish used the food with different efficiencies. However, grouping the fish according to size classes and temperature lowered the relative error of the predicted food intake for the group to typically 5 %. For a group containing all fish of the feeding experiment the relative prediction error was about 2 %. Zusammenfassung Wachstum und Kondition der Fische sind von der verfügbaren Nahrung und von Umweltbedingungen abhängig. Dies führte zur Idee, Fisch als „Konsum-Sensor“ für die Messung der Nahrungsaufnahme über einen definierten Zeitraum zu verwenden. Auf Grundlage des von Bertalanffy-Modells wurde ein bio-physikalisches Modell zur Schätzung der Futteraufnahme entwickelt. Während einige der Eingangsgrößen des Modells leicht gemessen werden können (Anfangs- und Endlänge und -körpermasse der Fische und die Temperatur innerhalb des betrachteten Zeitraum), können interne Parameter der betrachteten Art nur indirekt bestimmt werden. Drei interne Parameter werden in dem Modell verwendet: Die Erhaltungskonsumtion bei 0° C, die Temperaturabhängigkeit dieser Rate und der Wirkungsgrad der Nahrung (der Anteil der Nahrung ,der aufgenommen und verwendet und nicht ungenutzt wieder ausgeschieden wird). Die Modellparameter für eine bestimmte Art können durch Fütterungsversuche bestimmt werden. Um das Modell zu validieren wurden Daten aus veröffentlichten Fütterungsversuchen mit juvenilen Kabeljau (Gadus morhua L.) verwendet. Modell und Wirklichkeit weichen in der Regel voneinander ab. Der durchschnittliche relative Fehler der durch das Modell vorhergesagten Nahrungsaufnahme betrug für Einzelfische etwa 24%, was darauf hinweist, dass einzelne Fisch die Nahrung mit unterschiedlichen Wirkungsgraden verwerten. Allerdings senkte die Gruppierung der Fische nach Größenklassen und Temperatur den relativen Vorhersagefehler für die Nahrungsaufnahme der Gruppe auf etwa 5%. Für alle Fische im Fütterungsversuch ist der relative Vorhersagefehler etwa 2%.

Stability of Hypervelocity Boundary Layers

The early stage of laminar-turbulent transition in a hypervelocity boundary layer is studied using a combination of modal linear stability analysis, transient growth analysis, and direct numerical simulation. Modal stability analysis is used to clarify the behavior of first and second mode instabilities on flat plates and sharp cones for a wide range of high enthalpy flow conditions relevant to experiments in impulse facilities. Vibrational nonequilibrium is included in this analysis, its influence on the stability properties is investigated, and simple models for predicting when it is important are described.

Transient growth analysis is used to determine the optimal initial conditions that lead to the largest possible energy amplification within the flow. Such analysis is performed for both spatially and temporally evolving disturbances. The analysis again targets flows that have large stagnation enthalpy, such as those found in shock tunnels, expansion tubes, and atmospheric flight at high Mach numbers, and clarifies the effects of Mach number and wall temperature on the amplification achieved. Direct comparisons between modal and non-modal growth are made to determine the relative importance of these mechanisms under different flow regimes.

Conventional stability analysis employs the assumption that disturbances evolve with either a fixed frequency (spatial analysis) or a fixed wavenumber (temporal analysis). Direct numerical simulations are employed to relax these assumptions and investigate the downstream propagation of wave packets that are localized in space and time, and hence contain a distribution of frequencies and wavenumbers. Such wave packets are commonly observed in experiments and hence their amplification is highly relevant to boundary layer transition prediction. It is demonstrated that such localized wave packets experience much less growth than is predicted by spatial stability analysis, and therefore it is essential that the bandwidth of localized noise sources that excite the instability be taken into account in making transition estimates. A simple model based on linear stability theory is also developed which yields comparable results with an enormous reduction in computational expense. This enables the amplification of finite-width wave packets to be taken into account in transition prediction.

Growth parameters for Lates niloticus (L.), Bagrus docmak (Forsskal), Oreochromis niloticus (L.), Clarias gariepinus (Burchell) and Synodontis species derived from tag returns

In a tagging experiment carried out in the Kenyan waters of Lake Victoria, an annual growth increment of 29 cm yr was obtained for Lates niloticus (L.). Growth parameters obtained using the von Bertalanffy model on the growth curve fitted by eye were L (inf.) = 122 cm yr and k = 0.26 yr. Data for other species tagged were inadequate to obtain meaningful results.

The impact of physical processes on algal growth

Mixing and transport processes in surface waters strongly influence the structure of aquatic ecosystems. The impact of mixing on algal growth is species-dependent, affecting the competition among species and acting as a selective factor for the composition of the biocoenose. Were it not for the ever-changing ”aquatic weather”, the composition of pelagic ecosystems would be relatively simple. Probably just a few optimally adapted algal species would survive in a given water-body. In contrast to terrestrial ecosystems, in which the spatial heterogeneity is primarily responsible for the abundance of niches, in aquatic systems (especially in the pelagic zone) the niches are provided by the temporal structure of physical processes. The latter are discussed in terms of the relative sizes of physical versus biological time-scales. The relevant time-scales of mixing and transport cover the range between seconds and years. Correspondingly, their influence on growth of algae is based on different mechanisms: rapid changes are relevant for the fast biological processes such as nutrient uptake and photosynthesis, and the slower changes are relevant for the less dynamic processes such as growth, respiration, mineralization, and settling of algal cells. Mixing time-scales are combined with a dynamic model of photosynthesis to demonstrate their influence on algal growth.

Tsunamis - a model of their generation and propagation

A general solution is presented for water waves generated by an arbitrary movement of the bed (in space and time) in a two-dimensional fluid domain with a uniform depth. The integral solution which is developed is based on a linearized approximation to the complete (nonlinear) set of governing equations. The general solution is evaluated for the specific case of a uniform upthrust or downthrow of a block section of the bed; two time-displacement histories of the bed movement are considered.

An integral solution (based on a linear theory) is also developed for a three-dimensional fluid domain of uniform depth for a class of bed movements which are axially symmetric. The integral solution is evaluated for the specific case of a block upthrust or downthrow of a section of the bed, circular in planform, with a time-displacement history identical to one of the motions used in the two-dimensional model.

Since the linear solutions are developed from a linearized approximation of the complete nonlinear description of wave behavior, the applicability of these solutions is investigated. Two types of non-linear effects are found which limit the applicability of the linear theory: (1) large nonlinear effects which occur in the region of generation during the bed movement, and (2) the gradual growth of nonlinear effects during wave propagation.

A model of wave behavior, which includes, in an approximate manner, both linear and nonlinear effects is presented for computing wave profiles after the linear theory has become invalid due to the growth of nonlinearities during wave propagation.

An experimental program has been conducted to confirm both the linear model for the two-dimensional fluid domain and the strategy suggested for determining wave profiles during propagation after the linear theory becomes invalid. The effect of a more general time-displacement history of the moving bed than those employed in the theoretical models is also investigated experimentally.

The linear theory is found to accurately approximate the wave behavior in the region of generation whenever the total displacement of the bed is much less than the water depth. Curves are developed and confirmed by the experiments which predict gross features of the lead wave propagating from the region of generation once the values of certain nondimensional parameters (which characterize the generation process) are known. For example, the maximum amplitude of the lead wave propagating from the region of generation has been found to never exceed approximately one-half of the total bed displacement. The gross features of the tsunami resulting from the Alaskan earthquake of 27 March 1964 can be estimated from the results of this study.

Estimation of the growth curve parameters in Macrobrachium rosenbergii

Growth is one of the most important characteristics of cultured species. The objective of this study was to determine the fitness of linear, log linear, polynomial, exponential and Logistic functions to the growth curves of Macrobrachium rosenbergii obtained by using weekly records of live weight, total length, head length, claw length, and last segment length from 20 to 192 days of age. The models were evaluated according to the coefficient of determination (R2), and error sum off square (ESS) and helps in formulating breeders in selective breeding programs. Twenty full-sib families consisting 400 PLs each were stocked in 20 different hapas and reared till 8 weeks after which a total of 1200 animals were transferred to earthen ponds and reared up to 192 days. The R2 values of the models ranged from 56 – 96 in case of overall body weight with logistic model being the highest. The R2 value for total length ranged from 62 to 90 with logistic model being the highest. In case of head length, the R2 value ranged between 55 and 95 with logistic model being the highest. The R2 value for claw length ranged from 44 to 94 with logistic model being the highest. For last segment length, R2 value ranged from 55 – 80 with polynomial model being the highest. However, the log linear model registered low ESS value followed by linear model for overall body weight while exponential model showed low ESS value followed by log linear model in case of head length. For total length the low ESS value was given by log linear model followed by logistic model and for claw length exponential model showed low ESS value followed by log linear model. In case of last segment length, linear model showed lowest ESS value followed by log linear model. Since, the model that shows highest R2 value with low ESS value is generally considered as the best fit model. Among the five models tested, logistic model, log linear model and linear models were found to be the best models for overall body weight, total length and head length respectively. For claw length and last segment length, log linear model was found to be the best model. These models can be used to predict growth rates in M. rosenbergii. However, further studies need to be conducted with more growth traits taken into consideration