Mathematical modelling of controlled drug release from polymer micro-spheres : incorporating the effects of swelling, diffusion and dissolution via moving boundary problems

Controlled drug delivery is a key topic in modern pharmacotherapy, where controlled drug delivery devices are required to prolong the period of release, maintain a constant release rate, or release the drug with a predetermined release profile. In the pharmaceutical industry, the development process of a controlled drug delivery device may be facilitated enormously by the mathematical modelling of drug release mechanisms, directly decreasing the number of necessary experiments. Such mathematical modelling is difficult because several mechanisms are involved during the drug release process. The main drug release mechanisms of a controlled release device are based on the device’s physiochemical properties, and include diffusion, swelling and erosion. In this thesis, four controlled drug delivery models are investigated. These four models selectively involve the solvent penetration into the polymeric device, the swelling of the polymer, the polymer erosion and the drug diffusion out of the device but all share two common key features. The first is that the solvent penetration into the polymer causes the transition of the polymer from a glassy state into a rubbery state. The interface between the two states of the polymer is modelled as a moving boundary and the speed of this interface is governed by a kinetic law. The second feature is that drug diffusion only happens in the rubbery region of the polymer, with a nonlinear diffusion coefficient which is dependent on the concentration of solvent. These models are analysed by using both formal asymptotics and numerical computation, where front-fixing methods and the method of lines with finite difference approximations are used to solve these models numerically. This numerical scheme is conservative, accurate and easily implemented to the moving boundary problems and is thoroughly explained in Section 3.2. From the small time asymptotic analysis in Sections 5.3.1, 6.3.1 and 7.2.1, these models exhibit the non-Fickian behaviour referred to as Case II diffusion, and an initial constant rate of drug release which is appealing to the pharmaceutical industry because this indicates zeroorder release. The numerical results of the models qualitatively confirms the experimental behaviour identified in the literature. The knowledge obtained from investigating these models can help to develop more complex multi-layered drug delivery devices in order to achieve sophisticated drug release profiles. A multi-layer matrix tablet, which consists of a number of polymer layers designed to provide sustainable and constant drug release or bimodal drug release, is also discussed in this research. The moving boundary problem describing the solvent penetration into the polymer also arises in melting and freezing problems which have been modelled as the classical onephase Stefan problem. The classical one-phase Stefan problem has unrealistic singularities existed in the problem at the complete melting time. Hence we investigate the effect of including the kinetic undercooling to the melting problem and this problem is called the one-phase Stefan problem with kinetic undercooling. Interestingly we discover the unrealistic singularities existed in the classical one-phase Stefan problem at the complete melting time are regularised and also find out the small time behaviour of the one-phase Stefan problem with kinetic undercooling is different to the classical one-phase Stefan problem from the small time asymptotic analysis in Section 3.3. In the case of melting very small particles, it is known that surface tension effects are important. The effect of including the surface tension to the melting problem for nanoparticles (no kinetic undercooling) has been investigated in the past, however the one-phase Stefan problem with surface tension exhibits finite-time blow-up. Therefore we investigate the effect of including both the surface tension and kinetic undercooling to the melting problem for nanoparticles and find out the the solution continues to exist until complete melting. The investigation of including kinetic undercooling and surface tension to the melting problems reveals more insight into the regularisations of unphysical singularities in the classical one-phase Stefan problem. This investigation gives a better understanding of melting a particle, and contributes to the current body of knowledge related to melting and freezing due to heat conduction.

Enhanced processing of SPOT multispectral satellite imagery for environmental monitoring and modelling

The Taita Hills in southeastern Kenya form the northernmost part of Africa’s Eastern Arc Mountains, which have been identified by Conservation International as one of the top ten biodiversity hotspots on Earth. As with many areas of the developing world, over recent decades the Taita Hills have experienced significant population growth leading to associated major changes in land use and land cover (LULC), as well as escalating land degradation, particularly soil erosion. Multi-temporal medium resolution multispectral optical satellite data, such as imagery from the SPOT HRV, HRVIR, and HRG sensors, provides a valuable source of information for environmental monitoring and modelling at a landscape level at local and regional scales. However, utilization of multi-temporal SPOT data in quantitative remote sensing studies requires the removal of atmospheric effects and the derivation of surface reflectance factor. Furthermore, for areas of rugged terrain, such as the Taita Hills, topographic correction is necessary to derive comparable reflectance throughout a SPOT scene. Reliable monitoring of LULC change over time and modelling of land degradation and human population distribution and abundance are of crucial importance to sustainable development, natural resource management, biodiversity conservation, and understanding and mitigating climate change and its impacts. The main purpose of this thesis was to develop and validate enhanced processing of SPOT satellite imagery for use in environmental monitoring and modelling at a landscape level, in regions of the developing world with limited ancillary data availability. The Taita Hills formed the application study site, whilst the Helsinki metropolitan region was used as a control site for validation and assessment of the applied atmospheric correction techniques, where multiangular reflectance field measurements were taken and where horizontal visibility meteorological data concurrent with image acquisition were available. The proposed historical empirical line method (HELM) for absolute atmospheric correction was found to be the only applied technique that could derive surface reflectance factor within an RMSE of < 0.02 ps in the SPOT visible and near-infrared bands; an accuracy level identified as a benchmark for successful atmospheric correction. A multi-scale segmentation/object relationship modelling (MSS/ORM) approach was applied to map LULC in the Taita Hills from the multi-temporal SPOT imagery. This object-based procedure was shown to derive significant improvements over a uni-scale maximum-likelihood technique. The derived LULC data was used in combination with low cost GIS geospatial layers describing elevation, rainfall and soil type, to model degradation in the Taita Hills in the form of potential soil loss, utilizing the simple universal soil loss equation (USLE). Furthermore, human population distribution and abundance were modelled with satisfactory results using only SPOT and GIS derived data and non-Gaussian predictive modelling techniques. The SPOT derived LULC data was found to be unnecessary as a predictor because the first and second order image texture measurements had greater power to explain variation in dwelling unit occurrence and abundance. The ability of the procedures to be implemented locally in the developing world using low-cost or freely available data and software was considered. The techniques discussed in this thesis are considered equally applicable to other medium- and high-resolution optical satellite imagery, as well the utilized SPOT data.

A GIS-based model of soil erosion and transport

Soil erosion is a natural process that occurs when the force of wind, raindrops or running water on the soil surface exceeds the cohesive forces that bind the soil together. In general, vegetation cover protects the soil from the effects of these erosive forces. However, land management activities such as ploughing, burning or heavy grazing may disturb this protective layer, exposing the underlying soil. The decision making process in rural catchment management is often supported by the predictive modelling of soil erosion and sediment transport processes within the catchment, using established techniques such as the Universal Soil Loss Equation [USLE] and the Agricultural Nonpoint Source pollution model [AGNPS]. In this article, the authors examine the range of erosion models currently available and describe the application of one of these to the Burrishoole catchment on the north-west coast of Ireland, which has suffered heavy erosion of blanket peat in recent years.

Restrospective genetic monitoring of the threatened Yellow marsh saxifrage (Saxifraga hirculus) reveals genetic erosion but provides valuable information for conservation strategies

Aim: Retrospective genetic monitoring, comparing genetic diversity of extant populations with historical samples, can provide valuable and often unique insights into evolutionary processes informing conservation strategies. The Yellow marsh saxifrage (Saxifraga hirculus) is listed as ‘critically endangered’ in Ireland with only two extant populations. We quantified genetic changes over time and identified genotypes in extant populations that could be used as founders for reintroductions to sites where the species is extinct.

Location: Ireland.

Methods: Samples were obtained from both locations where the species is currently found, including the most threatened site at the Garron Plateau, Co. Antrim, which held only 13 individuals during 2011. Herbarium samples covering the period from 1886 to 1957 were obtained including plants from the same area as the most threatened population, as well as three extinct populations. In total, 422 individuals (319 present-day and 103 historical) were genotyped at six microsatellite loci. Species distribution modelling was used to identify areas of potentially suitable habitat for reintroductions.

Results: Level of phenotypic diversity within the most threatened population was significantly lower in the present-day compared with historical samples but levels of observed heterozygosity and number of alleles, whilst reduced, did not differ significantly. However, Bayesian clustering analysis suggested gradual lineage replacement over time. All three measures of genetic diversity were generally lower at the most threatened population compared with the more substantial extant populations in Co. Mayo. Species distribution modelling suggested that habitat at one site where the species is extinct may be suitable for reintroduction.

Main conclusions: The dominant genetic lineage in the most threatened population is rare elsewhere; thus, care needs to be taken when formulating any potential reintroduction programme. Our findings highlight both the need for genetic monitoring of threatened populations, but also for its swift implementation before levels of diversity become critically low.

In Vitro Release from Reverse Poloxamine/alpha-Cyclodextrin Matrices: Modelling and Comparison of Dissolution Profiles

Gels obtained by complexation of octablock star polyethylene oxide/polypropylene oxide copolymers (Tetronic 90R4) with -cyclodextrin (-CD) were evaluated as matrices for drug release. Both molecules are biocompatible so they can be potentially applied to drug delivery systems. Two different types of matrices of Tetronic 90R4 and -CD were evaluated: gels and tablets. These gels are capable to gelifying in situ and show sustained erosion kinetics in aqueous media. Tablets were prepared by freeze-drying and comprising the gels. Using these two different matrices, the release of two model molecules, L-tryptophan (Trp), and a protein, bovine serum albumin (BSA), was evaluated. The release profiles of these molecules from gels and tablets prove that they are suitable for sustained delivery. Mathematical models were applied to the release curves from tablets to elucidate the drug delivery mechanism. Good correlations were found for the fittings of the release curves to different equations. The results point that the release of Trp from different tablets is always governed by Fickian diffusion, whereas the release of BSA is governed by a combination of diffusion and tablet erosion. 

Sensitivity of a hydraulic model to changes in channel erosion during extreme flooding

Recent research into flood modelling has primarily concentrated on the simulation of inundation flow without considering the influences of channel morphology. River channels are often represented by a simplified geometry that is implicitly assumed to remain unchanged during flood simulations. However, field evidence demonstrates that significant morphological changes can occur during floods to mobilise the boundary sediments. Despite this, the effect of channel morphology on model results has been largely unexplored. To address this issue, the impact of channel cross-section geometry and channel long-profile variability on flood dynamics is examined using an ensemble of a 1D-2D hydraulic model (LISFLOOD-FP) of the 1:2102 year recurrence interval floods in Cockermouth, UK, within an uncertainty framework. A series of hypothetical scenarios of channel morphology were constructed based on a simple velocity based model of critical entrainment. A Monte-Carlo simulation framework was used to quantify the effects of channel morphology together with variations in the channel and floodplain roughness coefficients, grain size characteristics, and critical shear stress on measures of flood inundation. The results showed that the bed elevation modifications generated by the simplistic equations reflected a good approximation of the observed patterns of spatial erosion despite its overestimation of erosion depths. The effect of uncertainty on channel long-profile variability only affected the local flood dynamics and did not significantly affect the friction sensitivity and flood inundation mapping. The results imply that hydraulic models generally do not need to account for within event morphodynamic changes of the type and magnitude modelled, as these have a negligible impact that is smaller than other uncertainties, e.g. boundary conditions. Instead morphodynamic change needs to happen over a series of events to become large enough to change the hydrodynamics of floods in supply limited gravel-bed rivers like the one used in this research.

Calibration Strategies Designed For Long Term Sediment Transport Modelling

The presented work deals with the calibration of a 2D numerical model for the simulation of long term bed load transport. A settled basin along an alpine stream was used as a case study. The focus is to parameterise the used multi fractional transport model such that a dynamically balanced behavior regarding erosion and deposition is reached. The used 2D hydrodynamic model utilizes a multi-fraction multi-layer approach to simulate morphological changes and bed load transport. The mass balancing is performed between three layers: a top mixing layer, an intermediate subsurface layer and a bottom layer. Using this approach bears computational limitations in calibration. Due to the high computational demands, the type of calibration strategy is not only crucial for the result, but as well for the time required for calibration. Brute force methods such as Monte Carlo type methods may require a too large number of model runs. All here tested calibration strategies used multiple model runs utilising the parameterization and/or results from previous run. One concept was to reset to initial bed elevations after each run, allowing the resorting process to convert to stable conditions. As an alternative or in combination, the roughness was adapted, based on resulting nodal grading curves, from the previous run. Since the adaptations are a spatial process, the whole model domain is subdivided in homogeneous sections regarding hydraulics and morphological behaviour. For a faster optimization, the adaptation of the parameters is made section wise. Additionally, a systematic variation was done, considering results from previous runs and the interaction between sections. The used approach can be considered as similar to evolutionary type calibration approaches, but using analytical links instead of random parameter changes.

Spatial Assessment of Erosion and its Impact on Soil Fertility in the Tajik Foothills

Efficient planning of soil conservation measures requires, first, to understand the impact of soil erosion on soil fertility with regard to local land cover classes; and second, to identify hot spots of soil erosion and bright spots of soil conservation in a spatially explicit manner. Soil organic carbon (SOC) is an important indicator of soil fertility. The aim of this study was to conduct a spatial assessment of erosion and its impact on SOC for specific land cover classes. Input data consisted of extensive ground truth, a digital elevation model and Landsat 7 imagery from two different seasons. Soil spectral reflectance readings were taken from soil samples in the laboratory and calibrated with results of SOC chemical analysis using regression tree modelling. The resulting model statistics for soil degradation assessments are promising (R2=0.71, RMSEV=0.32). Since the area includes rugged terrain and small agricultural plots, the decision tree models allowed mapping of land cover classes, soil erosion incidence and SOC content classes at an acceptable level of accuracy for preliminary studies. The various datasets were linked in the hot-bright spot matrix, which was developed to combine soil erosion incidence information and SOC content levels (for uniform land cover classes) in a scatter plot. The quarters of the plot show different stages of degradation, from well conserved land to hot spots of soil degradation. The approach helps to gain a better understanding of the impact of soil erosion on soil fertility and to identify hot and bright spots in a spatially explicit manner. The results show distinctly lower SOC content levels on large parts of the test areas, where annual crop cultivation was dominant in the 1990s and where cultivation has now been abandoned. On the other hand, there are strong indications that afforestations and fruit orchards established in the 1980s have been successful in conserving soil resources.

Extensive management of field margins enhances their potential to mitigate soil erosion

Soil erosion is a widespread problem in agricultural landscapes, particularly in regions with strong rainfall events. Vegetated field margins can mitigate negative impacts of soil erosion by trapping eroded material. In this data set, we present data of sediment trapped by 12 field margins during the monsoon season of 2013 in an agricultural landscape in the Haean-myun catchment in South Korea. Prior to the beginning of monsoon season, we equipped a total of 12 sites representing three replicates for each of four different types of field margins ("managed flat", "managed steep", "natural flat" and "natural steep") with Astroturf mats with a size of 34 cm x 25 cm (850 cm**2). The mats (n = 15 / site) were installed at three levels: upslope, immediately before the field margin to quantify the sediments that reach it, in the middle of the field margin to quantify the locally trapped sediments, and after the field margin at the downslope edge to quantify the sediments that leave the field margin to the next field or to the stream. Sediment was collected after each rain event until the end of the monsoon season.

Pollen analysis, particle size analysis, end-member modelling and x-ray fluorescence measurements from sediment cores off Cape Ghir, Morocco

Understanding past human-climate-environment interactions is essential for assessing the vulnerability of landscapes and ecosystems to future climate change. This is particularly important in southern Morocco where the current vegetation is impacted by pastoralism, and the region is highly sensitive to climate variability. Here, we present a 2000-year record of vegetation, sedimentation rate, XRF chemical element intensities, and particle size from two decadal-resolved, marine sediment cores, raised from offshore Cape Ghir, southern Morocco. The results show that between 650 and 850 AD the sedimentation rate increased dramatically from 100 cm/1000 years to 300 cm/1000 years, and the Fe/Ca and pollen flux doubled, together indicating higher inputs of terrestrial sediment. Particle size measurements and end-member modelling suggest increased fluvial transport of the sediment. Beginning at 650 AD pollen levels from Cichorioideae species show a sharp rise from 10% to 20%. Pollen from Atemisia and Plantago, also increase from this time. Deciduous oak pollen percentages show a decline, whereas those of evergreen oak barely change. The abrupt increase in terrestrial/fluvial input from 650 to 850 AD occurs, within the age uncertainty, of the arrival of Islam (Islamisation) in Morocco at around 700 AD. Historical evidence suggests Islamisation led to population increase and development of southern Morocco, including expanded pastoralism, deforestation and agriculture. Livestock pressure may have changed the vegetation structure, accounting for the increase in pollen from Cichorioideae, Plantago, and Artemisia, which include many weedy species. Goats in particular may have played a dominant role as agents of erosion, and intense browsing may have led to the decline in deciduous oak; evergreen oak is more likely to survive as it re-sprouts more vigorously after browsing. From 850 AD to present sedimentation rates, Fe/Ca ratios and fluvial discharge remain stable, whereas pollen results suggest continued degradation. Pollen results from the past 150 years suggest expanded cultivation of olives and the native argan tree, and the introduction of Australian eucalyptus trees. The rapidly increasing population in southern Morocco is causing continued pressure to expand pastoralism and agriculture. The history of land degradation presented here suggests that the vegetation in southern Morocco may have been degraded for a longer period than previously thought and may be particularly sensitive to further land use changes. These results should be included in land management strategies for southern Morocco.