Mathematical modeling for improved greenhouse gas balances, agro-ecosystems, and policy development: Lessons from the Australian experience

If the land sector is to make significant contributions to mitigating anthropogenic greenhouse gas (GHG) emissions in coming decades, it must do so while concurrently expanding production of food and fiber. In our view, mathematical modeling will be required to provide scientific guidance to meet this challenge. In order to be useful in GHG mitigation policy measures, models must simultaneously meet scientific, software engineering, and human capacity requirements. They can be used to understand GHG fluxes, to evaluate proposed GHG mitigation actions, and to predict and monitor the effects of specific actions; the latter applications require a change in mindset that has parallels with the shift from research modeling to decision support. We compare and contrast 6 agro-ecosystem models (FullCAM, DayCent, DNDC, APSIM, WNMM, and AgMod), chosen because they are used in Australian agriculture and forestry. Underlying structural similarities in the representations of carbon flows though plants and soils in these models are complemented by a diverse range of emphases and approaches to the subprocesses within the agro-ecosystem. None of these agro-ecosystem models handles all land sector GHG fluxes, and considerable model-based uncertainty exists for soil C fluxes and enteric methane emissions. The models also show diverse approaches to the initialisation of model simulations, software implementation, distribution, licensing, and software quality assurance; each of these will differentially affect their usefulness for policy-driven GHG mitigation prediction and monitoring. Specific requirements imposed on the use of models by Australian mitigation policy settings are discussed, and areas for further scientific development of agro-ecosystem models for use in GHG mitigation policy are proposed.

In vitro modeling of the prostate cancer microenvironment

Prostate cancer is the most commonly diagnosed malignancy in men and advanced disease is incurable. Model systems are a fundamental tool for research and many in vitro models of prostate cancer use cancer cell lines in monoculture. Although these have yielded significant insight they are inherently limited by virtue of their two-dimensional (2D) growth and inability to include the influence of tumour microenvironment. These major limitations can be overcome with the development of newer systems that more faithfully recreate and mimic the complex in vivo multi-cellular, three-dimensional (3D) microenvironment. This article presents the current state of in vitro models for prostate cancer, with particular emphasis on 3D systems and the challenges that remain before their potential to advance our understanding of prostate disease and aid in the development and testing of new therapeutic agents can be realised.

Dynamic, invivo, real-time detection of retinal oxidative status in a model of elevated intraocular pressure using a novel, reversibly responsive, profluorescent nitroxide probe

Changes to the redox status of biological systems have been implicated in the pathogenesis of a wide variety of disorders including cancer, Ischemia-reperfusion (I/R) injury and neurodegeneration. In times of metabolic stress e.g. ischaemia/reperfusion, reactive oxygen species (ROS) production overwhelms the intrinsic antioxidant capacity of the cell, damaging vital cellular components. The ability to quantify ROS changes in vivo, is therefore essential to understanding their biological role. Here we evaluate the suitability of a novel reversible profluorescent probe containing a redox-sensitive nitroxide moiety (methyl ester tetraethylrhodamine nitroxide, ME-TRN), as an in vivo, real-time reporter of retinal oxidative status. The reversible nature of the probe's response offers the unique advantage of being able to monitor redox changes in both oxidizing and reducing directions in real time. After intravitreal administration of the ME-TRN probe, we induced ROS production in rat retina using an established model of complete, acute retinal ischaemia followed by reperfusion. After restoration of blood flow, retinas were imaged using a Micron III rodent fundus fluorescence imaging system, to quantify the redox-response of the probe. Fluorescent intensity declined during the first 60 min of reperfusion. The ROS-induced change in probe fluorescence was ameliorated with the retinal antioxidant, lutein. Fluorescence intensity in non-Ischemia eyes did not change significantly. This new probe and imaging technology provide a reversible and real-time response to oxidative changes and may allow the in vivo testing of antioxidant therapies of potential benefit to a range of diseases linked to oxidative stress

Data driven modeling for power transformer lifespan evaluation

Large sized power transformers are important parts of the power supply chain. These very critical networks of engineering assets are an essential base of a nation’s energy resource infrastructure. This research identifies the key factors influencing transformer normal operating conditions and predicts the asset management lifespan. Engineering asset research has developed few lifespan forecasting methods combining real-time monitoring solutions for transformer maintenance and replacement. Utilizing the rich data source from a remote terminal unit (RTU) system for sensor-data driven analysis, this research develops an innovative real-time lifespan forecasting approach applying logistic regression based on the Weibull distribution. The methodology and the implementation prototype are verified using a data series from 161 kV transformers to evaluate the efficiency and accuracy for energy sector applications. The asset stakeholders and suppliers significantly benefit from the real-time power transformer lifespan evaluation for maintenance and replacement decision support.

Dynamic pattern of suicide in Australia, 1986-2005: A descriptive-analytic study

Objective This study explores the spatiotemporal variations of suicide across Australia from 1986 to 2005, discusses the reasons for dynamic changes, and considers future suicide research and prevention strategies. Design Suicide (1986–2005) and population data were obtained from the Australian Bureau of Statistics. A series of analyses were conducted to examine the suicide pattern by sex, method and age group over time and geography. Results Differences in suicide rates across sex, age groups and suicide methods were found across geographical areas. Male suicides were mainly completed by hanging, firearms, gases and self-poisoning. Female suicides were primarily completed by hanging and self-poisoning. Suicide rates were higher in rural areas than in urban areas (capital cities and regional centres). Suicide rates by firearms were higher in rural areas than in urban areas, while the pattern for self-poisoning showed the reverse trend. Suicide rates had relatively stable trend for the total population and those aged between 15 and 54, while suicide decreased among 55 years and over during the study period. There was a decrease in suicides by firearms during the study period especially after 1996 when a new firearm control law was implemented, while suicide by hanging continued to increase. Areas with a high proportion of indigenous population (eg, northwest of Queensland and top north of the Northern Territory) had shown a substantial increase in suicide incidence after 1995. Conclusions Suicide rates varied over time and space and across sexes, age groups and suicide methods. This study provides detailed patterns of suicide to inform suicide control and prevention strategies for specific subgroups and areas of high and increased risk.

Development and validation of a multilevel model for predicting workload under routine and nonroutine conditions in an air traffic management center

Objective: The aim of this study was to develop a model capable of predicting variability in the mental workload experienced by frontline operators under routine and nonroutine conditions. Background: Excess workload is a risk that needs to be managed in safety-critical industries. Predictive models are needed to manage this risk effectively yet are difficult to develop. Much of the difficulty stems from the fact that workload prediction is a multilevel problem. Method: A multilevel workload model was developed in Study 1 with data collected from an en route air traffic management center. Dynamic density metrics were used to predict variability in workload within and between work units while controlling for variability among raters. The model was cross-validated in Studies 2 and 3 with the use of a high-fidelity simulator. Results: Reported workload generally remained within the bounds of the 90% prediction interval in Studies 2 and 3. Workload crossed the upper bound of the prediction interval only under nonroutine conditions. Qualitative analyses suggest that nonroutine events caused workload to cross the upper bound of the prediction interval because the controllers could not manage their workload strategically. Conclusion: The model performed well under both routine and nonroutine conditions and over different patterns of workload variation. Application: Workload prediction models can be used to support both strategic and tactical workload management. Strategic uses include the analysis of historical and projected workflows and the assessment of staffing needs. Tactical uses include the dynamic reallocation of resources to meet changes in demand.

Bayesian predictive modeling and comparison of oil samples

Statistical comparison of oil samples is an integral part of oil spill identification, which deals with the process of linking an oil spill with its source of origin. In current practice, a frequentist hypothesis test is often used to evaluate evidence in support of a match between a spill and a source sample. As frequentist tests are only able to evaluate evidence against a hypothesis but not in support of it, we argue that this leads to unsound statistical reasoning. Moreover, currently only verbal conclusions on a very coarse scale can be made about the match between two samples, whereas a finer quantitative assessment would often be preferred. To address these issues, we propose a Bayesian predictive approach for evaluating the similarity between the chemical compositions of two oil samples. We derive the underlying statistical model from some basic assumptions on modeling assays in analytical chemistry, and to further facilitate and improve numerical evaluations, we develop analytical expressions for the key elements of Bayesian inference for this model. The approach is illustrated with both simulated and real data and is shown to have appealing properties in comparison with both standard frequentist and Bayesian approaches

Spatially explicit structural equation modeling

Structural equation modeling (SEM) is a powerful statistical approach for the testing of networks of direct and indirect theoretical causal relationships in complex data sets with intercorrelated dependent and independent variables. SEM is commonly applied in ecology, but the spatial information commonly found in ecological data remains difficult to model in a SEM framework. Here we propose a simple method for spatially explicit SEM (SE-SEM) based on the analysis of variance/covariance matrices calculated across a range of lag distances. This method provides readily interpretable plots of the change in path coefficients across scale and can be implemented using any standard SEM software package. We demonstrate the application of this method using three studies examining the relationships between environmental factors, plant community structure, nitrogen fixation, and plant competition. By design, these data sets had a spatial component, but were previously analyzed using standard SEM models. Using these data sets, we demonstrate the application of SE-SEM to regularly spaced, irregularly spaced, and ad hoc spatial sampling designs and discuss the increased inferential capability of this approach compared with standard SEM. We provide an R package, sesem, to easily implement spatial structural equation modeling.

Modeling meat yield based on measurements of body traits in genetically improved giant freshwater prawn (GFP) Macrobrachium rosenbergii

A single-generation dataset consisting of 1,730 records from a selection program for high growth rate in giant freshwater prawn (GFP, Macrobrachium rosenbergii) was used to derive prediction equations for meat weight and meat yield. Models were based on body traits [body weight, total length and abdominal width (AW)] and carcass measurements (tail weight and exoskeleton-off weight). Lengths and width were adjusted for the systematic effects of selection line, male morphotypes and female reproductive status, and for the covariables of age at slaughter within sex and body weight. Body and meat weights adjusted for the same effects (except body weight) were used to calculate meat yield (expressed as percentage of tail weight/body weight and exoskeleton-off weight/body weight). The edible meat weight and yield in this GFP population ranged from 12 to 15 g and 37 to 45 %, respectively. The simple (Pearson) correlation coefficients between body traits (body weight, total length and AW) and meat weight were moderate to very high and positive (0.75–0.94), but the correlations between body traits and meat yield were negative (−0.47 to −0.74). There were strong linear positive relationships between measurements of body traits and meat weight, whereas relationships of body traits with meat yield were moderate and negative. Step-wise multiple regression analysis showed that the best model to predict meat weight included all body traits, with a coefficient of determination (R 2) of 0.99 and a correlation between observed and predicted values of meat weight of 0.99. The corresponding figures for meat yield were 0.91 and 0.95, respectively. Body weight or length was the best predictor of meat weight, explaining 91–94 % of observed variance when it was fitted alone in the model. By contrast, tail width explained a lower proportion (69–82 %) of total variance in the single trait models. It is concluded that in practical breeding programs, improvement of meat weight can be easily made through indirect selection for body trait combinations. The improvement of meat yield, albeit being more difficult, is possible by genetic means, with 91 % of the variation in the trait explained by the body and carcass traits examined in this study.

Analysis of parameter uncertainty and sensitivity in PCPF-1 modeling for predicting concentrations of rice herbicides

This paper demonstrates the procedures for probabilistic assessment of a pesticide fate and transport model, PCPF-1, to elucidate the modeling uncertainty using the Monte Carlo technique. Sensitivity analyses are performed to investigate the influence of herbicide characteristics and related soil properties on model outputs using four popular rice herbicides: mefenacet, pretilachlor, bensulfuron-methyl and imazosulfuron. Uncertainty quantification showed that the simulated concentrations in paddy water varied more than those of paddy soil. This tendency decreased as the simulation proceeded to a later period but remained important for herbicides having either high solubility or a high 1st-order dissolution rate. The sensitivity analysis indicated that PCPF-1 parameters requiring careful determination are primarily those involve with herbicide adsorption (the organic carbon content, the bulk density and the volumetric saturated water content), secondary parameters related with herbicide mass distribution between paddy water and soil (1st-order desorption and dissolution rates) and lastly, those involving herbicide degradations. © Pesticide Science Society of Japan.

Modeling approaches for pesticide exposure assessment in rice paddies

Pesticide use in paddy rice production may contribute to adverse ecological effects in surface waters. Risk assessments conducted for regulatory purposes depend on the use of simulation models to determine predicted environment concentrations (PEC) of pesticides. Often tiered approaches are used, in which assessments at lower tiers are based on relatively simple models with conservative scenarios, while those at higher tiers have more realistic representations of physical and biochemical processes. This chapter reviews models commonly used for predicting the environmental fate of pesticides in rice paddies. Theoretical considerations, unique features, and applications are discussed. This review is expected to provide information to guide model selection for pesticide registration, regulation, and mitigation in rice production areas.

Modeling of Channel Potential and Subthreshold Slope of Symmetric Double-Gate Transistor

A new physically based classical continuous potential distribution model, particularly considering the channel center, is proposed for a short-channel undoped body symmetrical double-gate transistor. It involves a novel technique for solving the 2-D nonlinear Poisson's equation in a rectangular coordinate system, which makes the model valid from weak to strong inversion regimes and from the channel center to the surface. We demonstrated, using the proposed model, that the channel potential versus gate voltage characteristics for the devices having equal channel lengths but different thicknesses pass through a single common point (termed ``crossover point''). Based on the potential model, a new compact model for the subthreshold swing is formulated. It is shown that for the devices having very high short-channel effects (SCE), the effective subthreshold slope factor is mainly dictated by the potential close to the channel center rather than the surface. SCEs and drain-induced barrier lowering are also assessed using the proposed model and validated against a professional numerical device simulator.

Complex crystallization dynamics in amorphous germanium observed with dynamic transmission electron microscopy

Crystallization of amorphous germanium (a-Ge) by laser or electron beam heating is a remarkably complex process that involves several distinct modes of crystal growth and the development of intricate microstructural patterns on the nanosecond to ten microsecond time scales. Here we use dynamic transmission electron microscopy (DTEM) to study the fast, complex crystallization dynamics with 10 nm spatial and 15 ns temporal resolution. We have obtained time-resolved real-space images of nanosecond laser-induced crystallization in a-Ge with unprecedentedly high spatial resolution. Direct visualization of the crystallization front allows for time-resolved snapshots of the initiation and roughening of the dendrites on submicrosecond time scales. This growth is followed by a rapid transition to a ledgelike growth mechanism that produces a layered microstructure on a time scale of several microseconds. This study provides insights into the mechanisms governing this complex crystallization process and is a dramatic demonstration of the power of DTEM for studying time-dependent material processes far from equilibrium.

In situ investigation of explosive crystallization in a-Ge: Experimental determination of the interface response function using dynamic transmission electron microscopy

The crystallization of amorphous semiconductors is a strongly exothermic process. Once initiated the release of latent heat can be sufficient to drive a self-sustaining crystallization front through the material in a manner that has been described as explosive. Here, we perform a quantitative in situ study of explosive crystallization in amorphous germanium using dynamic transmission electron microscopy. Direct observations of the speed of the explosive crystallization front as it evolves along a laser-imprinted temperature gradient are used to experimentally determine the complete interface response function (i.e., the temperature-dependent front propagation speed) for this process, which reaches a peak of 16 m/s. Fitting to the Frenkel-Wilson kinetic law demonstrates that the diffusivity of the material locally/immediately in advance of the explosive crystallization front is inconsistent with those of a liquid phase. This result suggests a modification to the liquid-mediated mechanism commonly used to describe this process that replaces the phase change at the leading amorphous-liquid interface with a change in bonding character (from covalent to metallic) occurring in the hot amorphous material.