Simulating Carbon Stocks and Fluxes of an African Tropical Montane Forest with an Individual-Based Forest Model

Tropical forests are carbon-dense and highly productive ecosystems. Consequently, they play an important role in the global carbon cycle. In the present study we used an individual-based forest model (FORMIND) to analyze the carbon balances of a tropical forest. The main processes of this model are tree growth, mortality, regeneration, and competition. Model parameters were calibrated using forest inventory data from a tropical forest at Mt. Kilimanjaro. The simulation results showed that the model successfully reproduces important characteristics of tropical forests (aboveground biomass, stem size distribution and leaf area index). The estimated aboveground biomass (385 t/ha) is comparable to biomass values in the Amazon and other tropical forests in Africa. The simulated forest reveals a gross primary production of 24 tcha-1yr-1. Modeling above- and belowground carbon stocks, we analyzed the carbon balance of the investigated tropical forest. The simulated carbon balance of this old-growth forest is zero on average. This study provides an example of how forest models can be used in combination with forest inventory data to investigate forest structure and local carbon balances.

Tuberculosis in Cape Town: An age-structured transmission model.

BACKGROUND Tuberculosis (TB) is the leading cause of death in South Africa. The burden of disease varies by age, with peaks in TB notification rates in the HIV-negative population at ages 0-5, 20-24, and 45-49 years. There is little variation between age groups in the rates in the HIV-positive population. The drivers of this age pattern remain unknown. METHODS We developed an age-structured simulation model of Mycobacterium tuberculosis (Mtb) transmission in Cape Town, South Africa. We considered five states of TB progression: susceptible, infected (latent TB), active TB, treated TB, and treatment default. Latently infected individuals could be re-infected; a previous Mtb infection slowed progression to active disease. We further considered three states of HIV progression: HIV negative, HIV positive, on antiretroviral therapy. To parameterize the model, we analysed treatment outcomes from the Cape Town electronic TB register, social mixing patterns from a Cape Town community and used literature estimates for other parameters. To investigate the main drivers behind the age patterns, we conducted sensitivity analyses on all parameters related to the age structure. RESULTS The model replicated the age patterns in HIV-negative TB notification rates of Cape Town in 2009. Simulated TB notification rate in HIV-negative patients was 1000/100,000 person-years (pyrs) in children aged <5 years and decreased to 51/100,000 in children 5-15 years. The peak in early adulthood occurred at 25-29 years (463/100,000 pyrs). After a subsequent decline, simulated TB notification rates gradually increased from the age of 30 years. Sensitivity analyses showed that the dip after the early adult peak was due to the protective effect of latent TB and that retreatment TB was mainly responsible for the rise in TB notification rates from the age of 30 years. CONCLUSION The protective effect of a first latent infection on subsequent infections and the faster progression in previously treated patients are the key determinants of the age-structure of TB notification rates in Cape Town.

An ordinal logistic regression model with misclassification of the outcome variable and categorical covariate

Ordinal outcomes are frequently employed in diagnosis and clinical trials. Clinical trials of Alzheimer's disease (AD) treatments are a case in point using the status of mild, moderate or severe disease as outcome measures. As in many other outcome oriented studies, the disease status may be misclassified. This study estimates the extent of misclassification in an ordinal outcome such as disease status. Also, this study estimates the extent of misclassification of a predictor variable such as genotype status. An ordinal logistic regression model is commonly used to model the relationship between disease status, the effect of treatment, and other predictive factors. A simulation study was done. First, data based on a set of hypothetical parameters and hypothetical rates of misclassification was created. Next, the maximum likelihood method was employed to generate likelihood equations accounting for misclassification. The Nelder-Mead Simplex method was used to solve for the misclassification and model parameters. Finally, this method was applied to an AD dataset to detect the amount of misclassification present. The estimates of the ordinal regression model parameters were close to the hypothetical parameters. β1 was hypothesized at 0.50 and the mean estimate was 0.488, β2 was hypothesized at 0.04 and the mean of the estimates was 0.04. Although the estimates for the rates of misclassification of X1 were not as close as β1 and β2, they validate this method. X 1 0-1 misclassification was hypothesized as 2.98% and the mean of the simulated estimates was 1.54% and, in the best case, the misclassification of k from high to medium was hypothesized at 4.87% and had a sample mean of 3.62%. In the AD dataset, the estimate for the odds ratio of X 1 of having both copies of the APOE 4 allele changed from an estimate of 1.377 to an estimate 1.418, demonstrating that the estimates of the odds ratio changed when the analysis includes adjustment for misclassification. ^

Medem Channel profiles based on stand-alone model simulation and corrections with nudging

State-of-the-art process-based models have shown to be applicable to the simulation and prediction of coastal morphodynamics. On annual to decadal temporal scales, these models may show limitations in reproducing complex natural morphological evolution patterns, such as the movement of bars and tidal channels, e.g. the observed decadal migration of the Medem Channel in the Elbe Estuary, German Bight. Here a morphodynamic model is shown to simulate the hydrodynamics and sediment budgets of the domain to some extent, but fails to adequately reproduce the pronounced channel migration, due to the insufficient implementation of bank erosion processes. In order to allow for long-term simulations of the domain, a nudging method has been introduced to update the model-predicted bathymetries with observations. The model-predicted bathymetry is nudged towards true states in annual time steps. Sensitivity analysis of a user-defined correlation length scale, for the definition of the background error covariance matrix during the nudging procedure, suggests that the optimal error correlation length is similar to the grid cell size, here 80-90 m. Additionally, spatially heterogeneous correlation lengths produce more realistic channel depths than do spatially homogeneous correlation lengths. Consecutive application of the nudging method compensates for the (stand-alone) model prediction errors and corrects the channel migration pattern, with a Brier skill score of 0.78. The proposed nudging method in this study serves as an analytical approach to update model predictions towards a predefined 'true' state for the spatiotemporal interpolation of incomplete morphological data in long-term simulations.

Regional ocean model for the Mediterranean, links to model files

A regional ocean general circulation model of the Mediterranean is used to study the climate of the Last Glacial Maximum. The atmospheric forcing for these simulations has been derived from simulations with an atmospheric general circulation model, which in turn was forced with surface conditions from a coarse resolution earth system model. The model is successful in reproducing the general patterns of reconstructed sea surface temperature anomalies with the strongest cooling in summer in the northwestern Mediterranean and weak cooling in the Levantine, although the model underestimates the extent of the summer cooling in the western Mediterranean. However, there is a strong vertical gradient associated with this pattern of summer cooling, which makes the comparison with reconstructions complicated. The exchange with the Atlantic is decreased to roughly one half of its present value, which can be explained by the shallower Strait of Gibraltar as a consequence of lower global sea level. This reduced exchange causes a strong increase of salinity in the Mediterranean in spite of reduced net evaporation.

Efficient parallelization of a regional ocean model for the western Mediterranean Sea

This paper focuses on the parallelization of an ocean model applying current multicore processor-based cluster architectures to an irregular computational mesh. The aim is to maximize the efficiency of the computational resources used. To make the best use of the resources offered by these architectures, this parallelization has been addressed at all the hardware levels of modern supercomputers: firstly, exploiting the internal parallelism of the CPU through vectorization; secondly, taking advantage of the multiple cores of each node using OpenMP; and finally, using the cluster nodes to distribute the computational mesh, using MPI for communication within the nodes. The speedup obtained with each parallelization technique as well as the combined overall speedup have been measured for the western Mediterranean Sea for different cluster configurations, achieving a speedup factor of 73.3 using 256 processors. The results also show the efficiency achieved in the different cluster nodes and the advantages obtained by combining OpenMP and MPI versus using only OpenMP or MPI. Finally, the scalability of the model has been analysed by examining computation and communication times as well as the communication and synchronization overhead due to parallelization.

An Agent-based Classification Model

The major function of this model is to access the UCI Wisconsin Breast Cancer data-set[1] and classify the data items into two categories, which are normal and anomalous. This kind of classification can be referred as anomaly detection, which discriminates anomalous behaviour from normal behaviour in computer systems. One popular solution for anomaly detection is Artificial Immune Systems (AIS). AIS are adaptive systems inspired by theoretical immunology and observed immune functions, principles and models which are applied to problem solving. The Dendritic Cell Algorithm (DCA)[2] is an AIS algorithm that is developed specifically for anomaly detection. It has been successfully applied to intrusion detection in computer security. It is believed that agent-based modelling is an ideal approach for implementing AIS, as intelligent agents could be the perfect representations of immune entities in AIS. This model evaluates the feasibility of re-implementing the DCA in an agent-based simulation environment called AnyLogic, where the immune entities in the DCA are represented by intelligent agents. If this model can be successfully implemented, it makes it possible to implement more complicated and adaptive AIS models in the agent-based simulation environment.

Experimental and diagnostic protocol for the physical component of the CMIP6 Ocean Model Intercomparison Project (OMIP)

The Ocean Model Intercomparison Project (OMIP) aims to provide a framework for evaluating, understanding, and improving the ocean and sea-ice components of global climate and earth system models contributing to the Coupled Model Intercomparison Project Phase 6 (CMIP6). OMIP addresses these aims in two complementary manners: (A) by providing an experimental protocol for global ocean/sea-ice models run with a prescribed atmospheric forcing, (B) by providing a protocol for ocean diagnostics to be saved as part of CMIP6. We focus here on the physical component of OMIP, with a companion paper (Orr et al., 2016) offering details for the inert chemistry and interactive biogeochemistry. The physical portion of the OMIP experimental protocol follows that of the interannual Coordinated Ocean-ice Reference Experiments (CORE-II). Since 2009, CORE-I (Normal Year Forcing) and CORE-II have become the standard method to evaluate global ocean/sea-ice simulations and to examine mechanisms for forced ocean climate variability. The OMIP diagnostic protocol is relevant for any ocean model component of CMIP6, including the DECK (Diagnostic, Evaluation and Characterization of Klima experiments), historical simulations, FAFMIP (Flux Anomaly Forced MIP), C4MIP (Coupled Carbon Cycle Climate MIP), DAMIP (Detection and Attribution MIP), DCPP (Decadal Climate Prediction Project), ScenarioMIP (Scenario MIP), as well as the ocean-sea ice OMIP simulations. The bulk of this paper offers scientific rationale for saving these diagnostics.

OMIP contribution to CMIP6: experimental and diagnostic protocol for the physical component of the Ocean Model Intercomparison Project

The Ocean Model Intercomparison Project (OMIP) is an endorsed project in the Coupled Model Intercomparison Project Phase 6 (CMIP6). OMIP addresses CMIP6 science questions, investigating the origins and consequences of systematic model biases. It does so by providing a framework for evaluating (including assessment of systematic biases), understanding, and improving ocean, sea-ice, tracer, and biogeochemical components of climate and earth system models contributing to CMIP6. Among the WCRP Grand Challenges in climate science (GCs), OMIP primarily contributes to the regional sea level change and near-term (climate/decadal) prediction GCs. OMIP provides (a) an experimental protocol for global ocean/sea-ice models run with a prescribed atmospheric forcing; and (b) a protocol for ocean diagnostics to be saved as part of CMIP6. We focus here on the physical component of OMIP, with a companion paper (Orr et al., 2016) detailing methods for the inert chemistry and interactive biogeochemistry. The physical portion of the OMIP experimental protocol follows the interannual Coordinated Ocean-ice Reference Experiments (CORE-II). Since 2009, CORE-I (Normal Year Forcing) and CORE-II (Interannual Forcing) have become the standard methods to evaluate global ocean/sea-ice simulations and to examine mechanisms for forced ocean climate variability. The OMIP diagnostic protocol is relevant for any ocean model component of CMIP6, including the DECK (Diagnostic, Evaluation and Characterization of Klima experiments), historical simulations, FAFMIP (Flux Anomaly Forced MIP), C4MIP (Coupled Carbon Cycle Climate MIP), DAMIP (Detection and Attribution MIP), DCPP (Decadal Climate Prediction Project), ScenarioMIP, HighResMIP (High Resolution MIP), as well as the ocean/sea-ice OMIP simulations.

The routes of the Mediterranean Sea circulation

The Mediterranean Sea is a semi-enclosed basin connected to the Atlantic Ocean through the narrow and shallow Strait of Gibraltar and further subdivided in two different sub-basins, the Eastern Mediterranean and the Western Mediterranean, connected through the Stait of Sicily. On annual basis, a net heat budget of −7 W/m2, combined with exceeding evaporation over precipation and runoff together with wind stress, is responsible for the antiestuarine character of the zonal thermoaline circulation. The outflow at Gibraltar Strait is mainly composed of Levantine Intermediate Water (LIW) and deep water masses formed in the Western Mediterranean Sea. The aim of this thesis is to validate and quantitatively assess the main routes of water masses composing the ouflow at Gibraltar Strait, using for the first time in the Mediterranean Sea a lagrangian interpretation of the eulerian velocity field produced from an eddy-resolving reanalysis dataset, spanning from 2000 to 2012. A lagrangian model named Ariane is used to map out three-dimensional trajectories in order to describe the pathways of water mass transport from the Strait of Sicily, the Gulf of Lyon and the Northern Tyrrhenian Sea to the Gibraltar Strait. Numerical experiments were carried out by seeding millions of particles in the Strait of Gibraltar and following them backwards in time to track the origins of water masses and transport exchanged between the different sections of the Mediterranean. Finally, the main routes of the intermediate and deep water masses are reconstructed from virtual particles trajectories, which highlight the role of the Western Mediterranean Deep Water (WMDW) as the main contributor to the Gibraltar Strait outflow. For the first time, the quantitative description of the flow of water masses coming from the Eastern Mediterranean towards the Gibraltar Strait is provided and a new route that directly links the Northern Tyrrhenian Sea to Gibralatr Strait has been detected.

Evaluation of an innovative critical power model in intermittent vertical jump

The aim of this study was to test if the critical power model can be used to determine the critical rest interval (CRI) between vertical jumps. Ten males performed intermittent countermovement jumps on a force platform with different resting periods (4.1 +/- 0.3 s, 5.0 +/- 0.4 s, 5.9 +/- 0.6 s). Jump trials were interrupted when participants could no longer maintain 95% of their maximal jump height. After interruption, number of jumps, total exercise duration and total external work were computed. Time to exhaustion (s) and total external work (J) were used to solve the equation Work = a + b . time. The CRI (corresponding to the shortest resting interval that allowed jump height to be maintained for a long time without fatigue) was determined dividing the average external work needed to jump at a fixed height (J) by b parameter (J/s). in the final session, participants jumped at their calculated CRI. A high coefficient of determination (0.995 +/- 0.007) and the CRI (7.5 +/- 1.6 s) were obtained. In addition, the longer the resting period, the greater the number of jumps (44 13, 71 28, 105 30, 169 53 jumps; p<0.0001), time to exhaustion (179 +/- 50, 351 +/- 120, 610 +/- 141, 1,282 +/- 417 s; p<0.0001) and total external work (28.0 +/- 8.3, 45.0 +/- 16.6, 67.6 +/- 17.8, 111.9 +/- 34.6 kJ; p<0.0001). Therefore, the critical power model may be an alternative approach to determine the CRI during intermittent vertical jumps.

An electromechanical finite element model for piezoelectric energy harvester plates

Vibration-based energy harvesting has been investigated by several researchers over the last decade. The goal in this research field is to power small electronic components by converting the waste vibration energy available in their environment into electrical energy. Recent literature shows that piezoelectric transduction has received the most attention for vibration-to-electricity conversion. In practice, cantilevered beams and plates with piezoceramic layers are employed as piezoelectric energy harvesters. The existing piezoelectric energy harvester models are beam-type lumped parameter, approximate distributed parameter and analytical distributed parameter solutions. However, aspect ratios of piezoelectric energy harvesters in several cases are plate-like and predicting the power output to general (symmetric and asymmetric) excitations requires a plate-type formulation which has not been covered in the energy harvesting literature. In this paper. an electromechanically coupled finite element (FE) plate model is presented for predicting the electrical power output of piezoelectric energy harvester plates. Generalized Hamilton`s principle for electroelastic bodies is reviewed and the FE model is derived based on the Kirchhoff plate assumptions as typical piezoelectric energy harvesters are thin structures. Presence of conductive electrodes is taken into account in the FE model. The predictions of the FE model are verified against the analytical solution for a unimorph cantilever and then against the experimental and analytical results of a bimorph cantilever with a tip mass reported in the literature. Finally, an optimization problem is solved where the aluminum wing spar of an unmanned air vehicle (UAV) is modified to obtain a generator spar by embedding piezoceramics for the maximum electrical power without exceeding a prescribed mass addition limit. (C) 2009 Elsevier Ltd. All rights reserved.

An open-boundary integrable model of three coupled XY spin chains

The integrable open-boundary conditions for the model of three coupled one-dimensional XY spin chains are considered in the framework of the quantum inverse scattering method. The diagonal boundary K-matrices are found and a class of integrable boundary terms is determined. The boundary model Hamiltonian is solved by using the coordinate space Bethe ansatz technique and Bethe ansatz equations are derived. (C) 1998 Elsevier Science B.V.

Determining forest structural attributes using an inverted geometric-optical model in mixed eucalypt forests, Southeast Queensland, Australia

The Montreal Process indicators are intended to provide a common framework for assessing and reviewing progress toward sustainable forest management. The potential of a combined geometrical-optical/spectral mixture analysis model was assessed for mapping the Montreal Process age class and successional age indicators at a regional scale using Landsat Thematic data. The project location is an area of eucalyptus forest in Emu Creek State Forest, Southeast Queensland, Australia. A quantitative model relating the spectral reflectance of a forest to the illumination geometry, slope, and aspect of the terrain surface and the size, shape, and density, and canopy size. Inversion of this model necessitated the use of spectral mixture analysis to recover subpixel information on the fractional extent of ground scene elements (such as sunlit canopy, shaded canopy, sunlit background, and shaded background). Results obtained fron a sensitivity analysis allowed improved allocation of resources to maximize the predictive accuracy of the model. It was found that modeled estimates of crown cover projection, canopy size, and tree densities had significant agreement with field and air photo-interpreted estimates. However, the accuracy of the successional stage classification was limited. The results obtained highlight the potential for future integration of high and moderate spatial resolution-imaging sensors for monitoring forest structure and condition. (C) Elsevier Science Inc., 2000.

A two-dimensional analytical solution of groundwater responses to tidal loading in an estuary and ocean

Previous studies on tidal dynamics of coastal aquifers have focussed on the inland propagation of oceanic tides in the cross-shore direction, a configuration that is essentially one-dimensional. Aquifers at natural coasts can also be influenced by tidal waves in nearby estuaries, resulting in a more complex behaviour of head fluctuations in the aquifers. We present an analytical solution to the two-dimensional depth-averaged groundwater flow equation for a semi-infinite aquifer subject to oscillating head conditions at the boundaries. The solution describes the tidal dynamics of a coastal aquifer that is adjacent to a cross-shore estuary. Both the effects of oceanic and estuarine tides on the aquifer are included in the solution. The analytical prediction of the head fluctuations is verified by comparison with numerical solutions computed using a standard finite-difference method. An essential feature of the present analytical solution is the interaction between the cross- and along-shore tidal waves in the aquifer area near the estuary's entry. As the distance from the estuary or coastline increases, the wave interaction is weakened and the aquifer response is reduced, respectively, to the one-dimensional solution for oceanic tides or the solution of Sun (Sun H. A two-dimensional analytical solution of groundwater response to tidal loading in an estuary, Water Resour Res 1997;33:1429-35) for two-dimensional non-interacting tidal waves. (C) 2000 Elsevier Science Ltd. All rights reserved.