Extending process and understanding for the development of complex ecosystem models, with application to the Chatham Rise Atlantis model

<p>The Chatham Rise is a highly productive deep-sea ecosystem that supports numerous substantial commercial fisheries, and is therefore a likely candidate for an ecosystem based approach to fisheries management in New Zealand. This thesis describes model construction, calibration and validation, for the first end-to-end ecosystem model of the Chatham Rise, New Zealand. The work extends beyond what has previously been done for validating such models, and explores uncertainty analyses through bootstrapping the oceanographic variables, perturbing the model's initial conditions, and analysing species interaction effects, with the results further analysed with respect to known data gaps. This enables the inclusion of uncertainty in simulated scenarios using the Chatham Rise Atlantis model, thus providing an envelope of results with which to analyse and understand the likely responses of the Chatham Rise ecosystem. The model was designed with 24 dynamic polygons, 5 water column depth bins, 55 species functional groups, and used 12-hour timesteps. The transfer of energy was tracked throughout the system using nitrogen as the model's main currency. The model simulated the system from 1900–2015, preceded by a 35 year burn-in period. The model produced very similar biomass trajectories in response to historical fishing to corresponding fisheries stock assessment models for key fisheries species. Population dynamics and system interactions were considered realistic with respect to growth rates, mortality rates, diets and species group interactions. The model was found to be generally stable under perturbations to the initial conditions, with lower trophic level species groups having the most variability. The specification of the Spawning Stock Recruitment curve was explored, as it relates to the multi-species and ecosystem models within which it is now applied. Close attention needs to be given to population dynamics specific to multi-species interactions such as predation-release, in particular the Spawning Stock Recruitment curve. Potentially misleading dynamics under predation-release were identified, and the simple solution of applying a cap to recruitment when biomass exceeds virgin levels was explored. The population dynamics of myctophids under fishing induced predation release were analysed with and without limiting recruitment to virgin levels. The effects were evident in several ecosystem indicators, suggesting unintentional mis-specification could lead to erroneous model results. It raises several questions around the specification of the Spawning Stock Recruitment relationship for multispecies models, and more generally, whether the concept of ‘virgin’ (or ‘unfished’) biomass should be reconsidered to reflect dynamic natural mortality and potentially changing unfished states. The model components that had knowledge gaps and were found to most likely to influence model results were the initial conditions, oceanographic variables, and the aggregate species groups ‘seabird’ and ‘cetacean other’. It is recommended that applications of the model, such as forecasting biomasses under various fishing regimes, should include alternatives that vary these components, and present appropriate levels of uncertainty in results. Initial conditions should be perturbed, with greater variability applied to species groups modelled as biomass-pools, and age-structured species groups that have little data available from the literature.</p>

Extending process and understanding for the development of complex ecosystem models, with application to the Chatham Rise Atlantis model

<p>The Chatham Rise is a highly productive deep-sea ecosystem that supports numerous substantial commercial fisheries, and is therefore a likely candidate for an ecosystem based approach to fisheries management in New Zealand. This thesis describes model construction, calibration and validation, for the first end-to-end ecosystem model of the Chatham Rise, New Zealand. The work extends beyond what has previously been done for validating such models, and explores uncertainty analyses through bootstrapping the oceanographic variables, perturbing the model's initial conditions, and analysing species interaction effects, with the results further analysed with respect to known data gaps. This enables the inclusion of uncertainty in simulated scenarios using the Chatham Rise Atlantis model, thus providing an envelope of results with which to analyse and understand the likely responses of the Chatham Rise ecosystem. The model was designed with 24 dynamic polygons, 5 water column depth bins, 55 species functional groups, and used 12-hour timesteps. The transfer of energy was tracked throughout the system using nitrogen as the model's main currency. The model simulated the system from 1900–2015, preceded by a 35 year burn-in period. The model produced very similar biomass trajectories in response to historical fishing to corresponding fisheries stock assessment models for key fisheries species. Population dynamics and system interactions were considered realistic with respect to growth rates, mortality rates, diets and species group interactions. The model was found to be generally stable under perturbations to the initial conditions, with lower trophic level species groups having the most variability. The specification of the Spawning Stock Recruitment curve was explored, as it relates to the multi-species and ecosystem models within which it is now applied. Close attention needs to be given to population dynamics specific to multi-species interactions such as predation-release, in particular the Spawning Stock Recruitment curve. Potentially misleading dynamics under predation-release were identified, and the simple solution of applying a cap to recruitment when biomass exceeds virgin levels was explored. The population dynamics of myctophids under fishing induced predation release were analysed with and without limiting recruitment to virgin levels. The effects were evident in several ecosystem indicators, suggesting unintentional mis-specification could lead to erroneous model results. It raises several questions around the specification of the Spawning Stock Recruitment relationship for multispecies models, and more generally, whether the concept of ‘virgin’ (or ‘unfished’) biomass should be reconsidered to reflect dynamic natural mortality and potentially changing unfished states. The model components that had knowledge gaps and were found to most likely to influence model results were the initial conditions, oceanographic variables, and the aggregate species groups ‘seabird’ and ‘cetacean other’. It is recommended that applications of the model, such as forecasting biomasses under various fishing regimes, should include alternatives that vary these components, and present appropriate levels of uncertainty in results. Initial conditions should be perturbed, with greater variability applied to species groups modelled as biomass-pools, and age-structured species groups that have little data available from the literature.</p>

Endurance Trained Athletes Do Not per se Have Higher Hoffmann Reflexes Than Recreationally Active Controls

The impact of endurance training on spinal neural circuitries remains largely unknown. Some studies have reported higher H-reflexes in endurance trained athletes and therefore, adaptations within the Ia afferent pathways after long term endurance training have been suggested. In the present study we tested the hypothesis that cyclists (n = 12) demonstrate higher Hoffmann reflexes (H-reflexes) compared to recreationally active controls (n = 10). Notwithstanding, highly significant differences in endurance performance (VO2peak: 60.6 for cyclists vs. 46.3 ml/min/kg for controls (p &lt; 0.001) there was no difference in the size of the SOL H-reflex between cyclists and controls (Hmax/Mmax ratio 61.3 vs. 60.0%, respectively (p = 0.840). Further analyses of the H and M recruitment curves for SOL revealed a significant steeper slope of the M recruitment curve in the group of cyclists (76.2 ± 3.8° vs. 72.0 ± 4.4°, p = 0.046) without a difference in the H-recruitment curve (84.6 ± 3.0° vs. 85.0 ± 2.8°, p = 0.784) compared to the control group. Cycling is classified as an endurance sport and thus the findings of the present study do not further support the assumption that long-term aerobic training leads to a general increase of the H-reflex. Amongst methodological differences in assessing the H-reflex, the training-specific sensorimotor control of the endurance sport itself might differently affect the responsiveness of spinal motoneurons on Ia-afferent inputs.

Influence of different transcranial magnetic stimulation current directions on the corticomotor control of lumbar erector spinae muscles during a static task

Different directions of transcranial magnetic stimulation (TMS) can activate different neuronal circuits. While posteroanterior current (PA-TMS) depolarizes mainly interneurons in primary motor cortex (M1), an anteroposterior current (AP-TMS) has been suggested to activate different M1 circuits and perhaps axons from the premotor regions. Although M1 is also involved in the control of axial muscles, no study has explored if different current directions activate different M1 circuits that may have distinct functional role. The aim of the study was to compare the effect of different current directions (PA- and AP-TMS) on the corticomotor control and spatial cortical organisation of the lumbar erector spinae muscle (LES). Thirthy-four healthy participants were recruited for two independent experiments and LES motor-evoked potentials (MEP) were recorded. In experiment 1 (n=17), active motor threshold (AMT), MEP latencies, recruitment curve (90 to 160% AMT), excitatory and inhibitory intracortical mechanisms using paired-pulse TMS (80% followed by 120% AMT stimuli at 2-3-10 and 15ms inter-stimulus intervals) were tested using a double cone (n=12) and a figure-of-eight (n=5) coils. In experiment 2 (n=17), LES cortical representations were tested using PA- and AP-TMS. AMT was higher for AP- compared to PA-TMS (p=0.002). Longer latencies with AP-TMS were compared to PA-TMS (p=0.017). AP-TMS produced more inhibition compared to PA-TMS at 2ms and 3ms (p=0.010), but no difference was observed for longer intervals. No difference was found for recruitment curve and mapping. Those findings suggest that each PA- and AP-TMS may activate different cortical circuits controlling low back muscles as proposed for hand muscles.

Transcutaneous Electrical Neuromodulation of the Cervical Spinal Cord Depends Both on the Stimulation Intensity and the Degree of Voluntary Activity for Training. A Pilot Study

Electrical enabling motor control (eEmc) through transcutaneous spinal cord stimulation offers promise in improving hand function. However, it is still unknown which stimulus intensity or which muscle force level could be better for this improvement. Nine healthy individuals received the following interventions: (i) eEmc intensities at 80%, 90% and 110% of abductor pollicis brevis motor threshold combined with hand training consisting in 100% handgrip strength; (ii) hand training consisting in 100% and 50% of maximal handgrip strength combined with 90% eEmc intensity. The evaluations included box and blocks test (BBT), maximal voluntary contraction (MVC), F wave persistency, F/M ratio, spinal and cortical motor evoked potentials (MEP), recruitment curves of spinal MEP and cortical MEP and short-interval intracortical inhibition. The results showed that: (i) 90% eEmc intensity increased BBT, MVC, F wave persistency, F/M ratio and cortical MEP recruitment curve; 110% eEmc intensity increased BBT, F wave persistency and cortical MEP and recruitment curve of cortical MEP; (ii) 100% handgrip strength training significantly modulated MVC, F wave persistency, F/M wave and cortical MEP recruitment curve in comparison to 50% handgrip strength. In conclusion, eEmc intensity and muscle strength during training both influence the results for neuromodulation at the cervical level.

Assessing head acceleration to identify a motor threshold to galvanic vestibular stimulation

Galvanic vestibular stimulation (GVS) is used to assess vestibular function, but vestibular responses can exhibit variability depending on protocols or intensities used. We measured head acceleration in healthy subjects to identify an objective motor threshold on which to base GVS intensity when assessing postural responses. Thirteen healthy right-handed subjects stood on a force platform, eyes closed, head facing forward. An accelerometer was placed on the vertex to detect head acceleration, and electromyography activity of the right soleus was recorded. GVS (200 ms; current steps 0.5;1-4mA) was applied in a binaural and bipolar configuration. 1) GVS induced a biphasic accelerometer response at a latency of 15 ms. Based on response amplitude, we constructed a recruitment curve for all participants and determined the motor threshold. In parallel, the method of limits was used to devise a more rapid approach to determine motor threshold. 2) We observed significant differences between motor threshold based on therecruitment curve and perceptual thresholds (sensation/perception of movement). No significant difference was observed between the motor threshold based on the method of limits and perceptual thresholds . 3) Using orthogonal polynomial contrasts, we observed a linear progression between multiples of the objective motor threshold (0.5, 0.75, 1, 1.5x motor threshold) and the 95% confidence ellipse area, the first peak of center of pressure velocity, and the short and medium latency responses in the soleus. Hence, an objective motor threshold and a recruitment curve for GVS were determined based on head acceleration, which could increase understanding of the vestibular system.

Methodological aspects of the stock-recruitment curve analysis

Choosing of mathematical model for describing of the pacific salmon progeny dependence on the parents abundance is discussed. Results of different approximations are interpreted. The analysis could be useful for forecasting of the pacific salmon stocks by specialists without deep skills in numerical modeling (like the author).

Therapeutic Elastic Tapes Applied in Different Directions Over the Triceps Surae Do Not Modulate Reflex Excitability of the Soleus Muscle

Context: Elastic taping has been widely used for either to facilitate or to inhibit muscle contraction. The efficacy of elastic taping is allegedly ascribed to physiological mechanisms related to subcutaneous tissue and muscle stimulation as a result of tape tension and direction. However, the underlying mechanisms that support the use of elastic taping are still unclear. Objective: To investigate changes in electrophysiological responses after 48 hours of tape application in different directions on the calf muscles of healthy individuals. Design: Within-subjects design. Setting: Research laboratory. Participants: Twenty-seven physically active males (age 18.0 [4.2] y, height 1.65 [0.07] m, body mass 62.3 [10.3] kg) participated. Interventions: Soleus H-reflex responses were evoked through stimulation of the tibial posterior nerve with 2- to 4-second interval between stimuli (32 sweeps) for each condition (baseline: without tape; facilitation: tape applied from muscle origin to insertion; inhibition: tape applied from muscle insertion to origin). Main Outcome Measures: The H-reflex amplitude values were normalized by the maximal direct response (Mmax). Parameters were estimated from a sigmoidal fit of the H-reflex recruitment curve (ascending limb). Results: No significant differences were found for the parameters derived from the recruitment curve of the H-reflex among the conditions (P > .05). Conclusions: The authors’ findings showed that, irrespective of the direction of tape application, the elastic tape applied over the triceps surae does not generate any significant alteration on the excitability of the reflex pathway for different subpopulations of motor units. The authors therefore suggest a re-examination of the current recommendations on taping direction in clinical and sports activities.

Integration of Convergent Sensorimotor Inputs Within Spinal Reflex Circuits in Healthy Adults

The output from motor neuron pools is influenced by the integration of synaptic inputs originating from descending corticomotor and spinal reflex pathways. In this study, using paired non-invasive brain and peripheral nerve stimulation, we investigated how descending corticomotor pathways influence the physiologic recruitment order of the soleus Hoffmann (H-) reflex. Eleven neurologically unimpaired adults (9 females; mean age 25 ± 3 years) completed an assessment of transcranial magnetic stimulation (TMS)-conditioning of the soleus H-reflex over a range of peripheral nerve stimulation (PNS) intensities. Unconditioned H-reflex recruitment curves were obtained by delivering PNS pulses to the posterior tibial nerve. Subsequently, TMS-conditioned H-reflex recruitment curves were obtained by pairing PNS with subthreshold TMS at short (−1.5 ms) and long (+10 ms) intervals. We evaluated unconditioned and TMS-conditioned H-reflex amplitudes along the ascending limb, peak, and descending limb of the H-reflex recruitment curve. Our results revealed that, for long-interval facilitation, TMS-conditioned H-reflex amplitudes were significantly larger than unconditioned H-reflex amplitudes along the ascending limb and peak of the H-reflex recruitment curve. Additionally, significantly lower PNS intensities were needed to elicit peak H-reflex amplitude (Hmax) for long-interval facilitation compared to unconditioned. These findings suggest that the influence of descending corticomotor pathways, particularly those mediating long-interval facilitation, contribute to changing the recruitment gain of the motor neuron pool, and can inform future methodological protocols for TMS-conditioning of H-reflexes. By characterizing and inducing short-term plasticity in circuitry mediating short- and long-interval TMS-conditioning of H-reflex amplitudes, future studies can investigate supraspinal and spinal circuit contributions to abnormal motor control, as well as develop novel therapeutic targets for neuromodulation.