Figure 1.
Mass fraction evolution across the membrane layer for different values of fKn for Kn = 0.1.
Citation: Alozie Alaoma, Nikolaos Voulvoulis. Mineral resource active regions: The need for systems thinking in management[J]. AIMS Environmental Science, 2018, 5(2): 78-95. doi: 10.3934/environsci.2018.2.78
| [1] | Tianyuan Xu, Shanming Ji, Chunhua Jin, Ming Mei, Jingxue Yin . EARLY AND LATE STAGE PROFILES FOR A CHEMOTAXIS MODEL WITH DENSITY-DEPENDENT JUMP PROBABILITY. Mathematical Biosciences and Engineering, 2018, 15(6): 1345-1385. doi: 10.3934/mbe.2018062 |
| [2] | Wenjie Zhang, Lu Xu, Qiao Xin . Global boundedness of a higher-dimensional chemotaxis system on alopecia areata. Mathematical Biosciences and Engineering, 2023, 20(5): 7922-7942. doi: 10.3934/mbe.2023343 |
| [3] | Sunwoo Hwang, Seongwon Lee, Hyung Ju Hwang . Neural network approach to data-driven estimation of chemotactic sensitivity in the Keller-Segel model. Mathematical Biosciences and Engineering, 2021, 18(6): 8524-8534. doi: 10.3934/mbe.2021421 |
| [4] | Qianhong Zhang, Fubiao Lin, Xiaoying Zhong . On discrete time Beverton-Holt population model with fuzzy environment. Mathematical Biosciences and Engineering, 2019, 16(3): 1471-1488. doi: 10.3934/mbe.2019071 |
| [5] | Chichia Chiu, Jui-Ling Yu . An optimal adaptive time-stepping scheme for solving reaction-diffusion-chemotaxis systems. Mathematical Biosciences and Engineering, 2007, 4(2): 187-203. doi: 10.3934/mbe.2007.4.187 |
| [6] | Xu Song, Jingyu Li . Asymptotic stability of spiky steady states for a singular chemotaxis model with signal-suppressed motility. Mathematical Biosciences and Engineering, 2022, 19(12): 13988-14028. doi: 10.3934/mbe.2022652 |
| [7] | Tingting Yu, Sanling Yuan . Dynamics of a stochastic turbidostat model with sampled and delayed measurements. Mathematical Biosciences and Engineering, 2023, 20(4): 6215-6236. doi: 10.3934/mbe.2023268 |
| [8] | Lin Zhang, Yongbin Ge, Zhi Wang . Positivity-preserving high-order compact difference method for the Keller-Segel chemotaxis model. Mathematical Biosciences and Engineering, 2022, 19(7): 6764-6794. doi: 10.3934/mbe.2022319 |
| [9] | Changwook Yoon, Sewoong Kim, Hyung Ju Hwang . Global well-posedness and pattern formations of the immune system induced by chemotaxis. Mathematical Biosciences and Engineering, 2020, 17(4): 3426-3449. doi: 10.3934/mbe.2020194 |
| [10] | Marcin Choiński, Mariusz Bodzioch, Urszula Foryś . A non-standard discretized SIS model of epidemics. Mathematical Biosciences and Engineering, 2022, 19(1): 115-133. doi: 10.3934/mbe.2022006 |
Nanomedicine knows an increasing interest due to its estimated importance in efficient medication delivery [1]. Crucial medical topics deal with its application in several fields, such as overcoming the blood-brain barrier [2], anti-diabetic medication [3] and uptake in endothelial cells [4]. Nanomedicine deals with the novel chemical, physical or biological properties of medication at nanometric scale, as well as the application of nanotechnology aiming at making a medical diagnosis, treating or preventing the cause of illness. Therefore, nanomedicines have the ability, due to their particular size or surface functionality, to demonstrate physicochemical and/or biological properties that make them interesting for medication delivery, controlled medication release and/or improved medication transport through biological barriers, such as cell membranes, which would otherwise be much more difficult [5]. One of the key parameters for successful medication delivery is controlling or improving the permeation of medication through the destined place (e.g. vascular system, cell) in the human body [6,7,8]. When speaking about cell permeation and the body uptake efficiency of drugs into the body that it entails, many works take the apparent cell permeability, Papp, through monolayers of human epithelial cells as a widely accepted standard [9,10,11]. Values of Papp (often presented in medicine as a volumetric flux per unity of area) of such cells are also often used for the estimation of the blood-brain barrier [12,13]. It is useful to understand the apparent cell permeation as a series of parallel and serial resistances (from layer to another). This vision is based on the concept of solubility diffusion [14]. The permeability is seen here, therefore, not as an imposed force, but is rather passive in a way that only spontaneous diffusion governs the permeation process. This diffusion process is determined by a concentration gradient. The concentration gradient is classically described by Fick's law or, in case of sizes (within the membranes separating the monolayers) at nanoscale, by a more general evolution equation, describing the extended diffusion law, as will be discussed later. Active transport, as could be caused by transporter proteins, is not considered here. Although, in reality, chemicals/drugs could be transported by proteins (and thus have an effect on Papp), it is not taken into account here. It is well admitted that generally three pathways are possible through a cellular monolayer, i.e. in between the cells (the paracellular way), through the cell (considered here) and along the membrane solely (the lateral way) [15,16]. In this case, parallel resistances of the three pathways would apply. Here, we only consider the pathway through the cell. Following the pH partition hypothesis [17], it can be assumed that only the neutral form of the drug chemicals can pass the membrane resistances. Therefore, chemicals that are ionized should undergo an acid-base reaction in between the water layer and the membrane before they are able to pass the membrane as a neutral species. According to the aforementioned pH partition hypothesis, such a reaction can be considered to occur instantaneously, meaning that no kinetic hinderness (and therefore no additional resistance) needs to be considered.
It is often suggested that the diffusion through membranes are of non-Fickian nature. We propose that this is either because of non-local effects (the size of the membranes are of the same order of magnitude as the drug molecules or not much bigger) or because the relaxation time becomes important. As controlling and enhancing drug delivery is an important issue in the field of nanomedicine, next to in vitro experiments, modelling is crucial to understand the underlying mechanisms. In nanomedicine, the inertia of the diffusing particles is relevant in addition to the small size of the channels through which the medication passes, being in some cases of the same order of magnitude as the mean free path of the particles. Furthermore, purely microscopic models are often difficult to implement directly to real-life applications. Therefore, further progress in this field can be enhanced by models that are easily applicable, whilst keeping a qualitative efficient and understandable way to link to experiments in a multi-disciplinary way. Therefore, we propose here a model based on thermodynamic principles, where we extend them by incorporating effects, in a rigourous way, that account for the relaxation and size of the drug molecules with respect to the channels through which they diffuse. With the purpose to demonstrate the applicability of our new theory, we treat two cases: drug diffusion through a cell membrane and drug release out of a scaffold. For the first case, let us imagine the existence of two compartments of well-mixed water that are separated by a barrier, i.e. a membrane: the first containing the drug molecules as the donor and the second represents the inner cell or the receptor. The purpose is to follow the evolution of the mass fraction profile through the membrane at the beginning of the diffusion process, where we assume that the drug molecules (in that beginning process) concentration is that high that the mass fraction at the inlet of the membrane is constant. For the second case, we discuss drug release out of a scaffold and follow the concentration profile to assess the rapidity of the drug release. To account for the importance of the relaxation time of the drug molecules and the non-local effects, we propose the formulation of mass transport using Extended Non-Equilibrium Thermodynamics [18,19], whose main characteristic is to upgrade the flux and higher order fluxes to the rank of independent variables.
At micro and nanoscales, mass transport is mostly influenced by non-local effects and relaxation times. The classical Fick law
| J=−ρD∇c, | (1) |
relating the mass flux vector J to the mass fraction gradient ∇c, with ρ and D denoting the density and mass diffusion coefficient, respectively, is not applicable at short diffusional time (with respect to the relaxation time, i.e. large relaxation times) and small spatial scales. In order to account for large relaxation times, Fick's law can be generalized by a Cattaneo-like equation (in analogy to the one used for heat transport [20]) under the form
| τr∂tJ+J=−ρD∇c | (2) |
where τr designates the relaxation time of the mass flux and ∂t the partial time derivative. To account furthermore for non-local effects, Extended Non-Equilibrium Thermodynamics (see for underlying developments in previous work [18]) proposes a further generalization
| τr∂tJ+J=−ρD∇c+ℓ2∇2J | (3) |
It is convenient to reformulate (3) in dimensionless form. Therefore, the mass flux J is rescaled with ρDδ, the space coordinate with a characteristic length δ (standing also for the membrane thickness), the time t with the characteristic diffusion time τd=δ2D. Introducing the Knudsen number Kn≡ℓδ, a dimensional number M=τrτd and using the same notation for the non-dimensional quantities, Eq (3) becomes
| M∂tJ+J=−∇c+Kn2∇2J | (4) |
It should be noted here that by choosing τd as time scale, one implies that any change in M means a change in the relaxation time, which evokes non-local effects. This motivates proposing an adaptation for the expression for M. Indeed, kinetically speaking, one may define a relation D≡13ℓ2τr, in analogy to the thermal conductivity [1]. Using the definition of τd and Kn, one finds a relation M=13Kn2. However, in suspensions, these kinetic considerations are not always valid and differences could occur. Acknowledging that some partial proportionality is admitted between M and the Kn number (the relaxation time and the mean free path are related in that they both can cause non-local effects), this motivates to define rather M=fKn13Kn2. Note that using this definition does not change anything for the analysis to be performed (indeed, in principle, any number could be chosen for M in the framework of a dimensionless analysis conditioned by correct interpretation of the results). Moreover, fKn is a factor that accounts for the observation that in nanofluids the relaxation time is not necessarily fully proportional to the Kn number. This factor also allows investigating separately the effects of the relaxation time and the mean free path on drug diffusion.
We can see clearly differences between (3) and the conventional constitutional Fick's law for mass diffusion (1), but also some similarities. When we deal with systems where the characteristic size is found to be much larger than the mean free path of the molecules l, Eq (3) reduces to the Cattaneo-like equation (2). This would also correspond to neglect any higher order of mass fluxes. Now, when the changes in the mass flux occur in a time span that is much larger than the relaxation time, we obtain Fick's law. Let us stress that the latter issue depends on the physicochemical relation the medication molecule can have with the carrier fluid or the host environment. Other well-known methods that are used for solving models that deal with nanoscale systems are the phonon-Boltzmann equation, Monte Carlo simulations and Molecular Dynamics. The Monte Carlo simulations and Molecular Dynamics method are numerical methods that generally do not require higher order variables, but instead demand high computing times, whilst being less flexible for modelling various scales at a time, although some efforts have been performed concerning multi-scale molecular dynamics [21]. Nonetheless, the flexibility in scale remains less than the present approach. The phonon-Boltzmann equation also demands high computational time for solving a model. The aforementioned aspects are generally less of a problem when it comes to the model that is proposed in this work. It should be noted that the present development could be extended by introducing non-linear mass fluxes or higher-order (higher than the second order) fluxes. However, implementing such an extension, albeit theoretically not an impossible task (see [18]), would introduce numerical difficulties and higher computing time, which is an absent aspect of the present approach. Furthermore, it would also loose the ability to propose analytical solutions for some particular cases, offering a pedagogical presentation allowing easier interpretations. Non-linear contributions of the mass flux (at first and higher orders) could be interesting to consider when the length scale becomes much smaller than the mean free path of the molecules. For the purposes of nanomedicine applications, this would in many cases lead to unnecessary complications. Indeed, for many cases a linear mass flux limited to second-order mass fluxes (which leads to the term ∇2J in (3)) is sufficient. In comparison, it already appeared to be sufficient for heat transport applications [22] and for predicting nanoflow properties [23].
From typical clinical observations that lead to data in the literature, it appears that one of the most important issues to be studied for improving drug delivery are through enhancing the permeability. In porous media, the concentration gradient can be negligible with respect to the flux, which is often the case for nanoporous flow that is pressure-driven by some external pressure. There are some cases in the human body that are comparable from a conceptual point of view, such as oxygen and carbon dioxide transfer from the lungs to the veins within the alveoles. Nonetheless, porous flow within cell membrane pores presents often the difference that, due to the absence of externally imposed pressure, the concentration gradient can be quite important with respect to the flux. In case of passive permeability, one speaks often of diffusional processes. As the diffusion process through the membrane is often the limiting pathway, the flux across the membrane will undergo non-local effects. Flux variation in the perpendicular direction to the flux across the membrane can be neglected, assuming a well-dispersed medication solution at the entrance of the whole membrane. One would then use (3) to describe such a process, reduced to a one-dimensional equation, replacing the ∇ operator by the derivative ∂y. However, the concentration gradient is not simply a gradient of the concentrations from both sides of the membrane, since concentration jumps can exist at the entry or exit of those membranes due to lipid bilayer interaction or wettability effects. The gradient that counts is the one within the membrane. Therefore, it is important to consider the fact that flux by diffusion across a membrane (even if it were uniform) depends on the solubility. As there are different solubilities at both sides of the membrane and within the membrane as well, we can define partition coefficients ϑin and ϑout for the entry and exit of the membrane, respectively:
| ϑin=Nmb,inNin | (5) |
| ϑout=Nmb,outNout | (6) |
where Nmb, in and Nmb, out are the concentrations within the membrane from the entry and exit sides, respectively. Also, Nin and Nout stand for the same concentrations, but at the medium sides (outside the membrane), respectively. It is more convenient to write (5) and (6) in mass fractions, c, where the subscripts mean the same thing as for the concentrations. Assuming ϑin = ϑout≡ϑ, we have
| ϑ=cmb,in(1−cin)cin(1−cmb,in) | (7) |
| ϑ=cmb,out(1−cout)cout(1−cmb,out) | (8) |
where ϑ is to be determined experimentally as a material property. We can readily see that for a higher partition coefficient, the concentration gradient within the membrane is then larger than on the outer sides of the membrane, which results into a higher membrane solubility and thus a higher permeation. To model the whole process, next to the one-dimensional version of Eq (3), the mass fraction evolution is given by the classical species balance equation, which in dimensionless form gives
| ∂tc=−∇⋅J | (9) |
Since the thickness of a cell membrane is much smaller than its size in the other directions, we consider a one-dimensional model through the membrane with thickness δ, assuming that in all directions other than across the membrane the variables remain constant. The entrance of the membrane (at the drug donor side) is denoted by the dimensionless coordinate y = 0 and the exit (at the drug receptor side) by y = 1. The transport of the drug molecules through a cell membrane will then be modeled by the one-dimensional versions of Eqs. (4) and (9)
| M∂tJ+J=−∂yc+Kn2∂2yJ | (10) |
| ∂tc=−∂yJ | (11) |
Here, we omitted the subscript mb for the membrane and will do so in the following. The drug at the donor side of, but outside, the membrane is given by cin. The inner-side drug mass fraction of the membrane at y = 0 is then, via the partition theory, given by
| c|y=0=ϑcin1+(ϑ−−1)cin | (12) |
At the cell-side of the membrane, y = 1, we assume a continuity of fluxes so that
| −D∂c∂y=−Dce∂cce∂y | (13) |
where the subscript ce denotes the inner cell. Assuming that the inner cell is an unstirred well-mixed environment, the drug uptake can be considered to occur under a constant rising concentration. Therefore, the gradient is nearly zero, so that we can say that ∂cce∂y≈0. This allows to state approximately that at y = 1
| ∂c∂y|y=1=0 | (14) |
The initial mass fraction within the membrane is given by
| c|t=0=c0 | (15) |
A supplementary interest is the continuous delivery of drugs out of a scaffold. In this case, a scaffold, containing the drugs, is introduced into the area of interest (the manner in which this is done is not the point of discussion here). This scaffold is a porous structure, where we are interested to provide a tool for analysing the diffusional flow through this medium. One side of the scaffold contains the drugs, maintained at a certain concentration or mass fraction, cin. The other side is where the drugs enter the area of interest, where we assume that the uptake is immediate so that the mass fraction is zero there, i.e. c0 = 0. Let us approximate such a scaffold as cylindrical pores. Contrary to the previous subsection, the dimension of the flux direction (axial direction of the pores) is larger than the radius of such a cylinder (radial direction of the pores). In that case, we can say (and also motivated by previous work in porous systems [23]) that the flux in the axial direction of the pores (the direction of drug release) will be approximately constant. However, the flux will be variable over the radial direction of the flux. Indeed, due to the small pore size with respect to the axial dimension of the pores, the wall effects will influence the flux and a flux profile exists in the radial direction. As for the mass fraction gradient (being a slow variable), its radial gradient is negligible with respect to the axial one (influenced by the drug release). As such, Eq (4) gives in dimensionless form (using ρDR, R and L as scaling for the mass flux, the radial spatial coordinate and the axial spatial coordinate, respectively, where R is the pore radius and L the length of the pore)
| J=−α∂c∂x+Kn2r1r∂∂r(r∂J∂r), | (16) |
with Knr≡ℓR, r the radial spatial coordinate, x the axial spatial coordinate in the pores and α≡RL. Assuming maximum flux in the middle of the pore (r = 0), i.e. ∂J∂r|r=0=0, and first-order slip (defining Bs≡ℓsR, where ℓs stands for the slip length, related to the wall material) at the wall (r = 1), i.e.
| J=−Bs∂J∂r|r=1, | (17) |
the flux can be found as a function of r. Taking subsequently an average value over a cross-section of the cylinder gives the following mass flux
| J=−α∂c∂x(1−2Kn2rB[1,1Knr]KnrB[0,1Knr]+BsB[1,1Knr]) | (18) |
with B the Bessel-I function. It should be noted that Eq (18) depends on how slip is observed and how the mass flux is defined. Other effects can influence the mass flux, but here it is the purpose to show how physical models can be used in the field of nanomedicine. Such an effective flux (taking into account the material of the scaffold (via the slip factor) and non-local effects via the Knudsen number) is now filled in the dimensionless (time scale is taken as L2D) species balance in the flow direction
| ∂tc=−1α∂xJ | (19) |
The evolution of the mass fraction is then described as
| ∂tc=∂2c∂x2(1−2Kn2B[1,1Kn]KnB[0,1Kn]+BsB[1,1Kn]) | (20) |
The boundary condition on one side of the scaffold is given by
| ∂c∂x|x=0=0 | (21) |
at x = 0, which stands for assuming that the drug only goes on one direction and
| c|x=1=c0e−t | (22) |
at x = 1, which allows for a drug uptake relaxation time of the same order of magnitude as the characteristic diffusion time (which is reasonable, considering that transport in cells is mainly of diffusive character), describing the mass fraction going from c0 to 0 (complete drug uptake). This depends of course of the system in question and the area in the body. Nevertheless, studying the effect of non-locality and slip factor on drug uptake, the choice of a relaxation time is not of importance here and does not influence our analysis. Eq (22) also shows the initial condition by evaluating it at t = 0.
The model is solved using a semi-implicit difference method for the spatial discretization and the forward Euler method for the temporal discretization for each time step t. At the beginning of the simulation, the initial conditions are imposed. At a certain time t, the equations for the mass flux and mass fractions are written in matrix form A*V = B, where A is a matrix of dimension (n-2)×(n-2) (n being the number of nodes in the discretized spatial domain and two nodes being dedicated to the boundary conditions, the bulk equations are thus discretized on (n-2) nodes) containing the coefficients in the equations for the mass fractions and mass fluxes at time t+1, V a vector of dimension (n-2) containing the unknown mass fraction and mass flux values at the spatial nodes at time t+1 and B a vector of dimension (n-2) containing the known mass fraction and mass flulx values at time t. Due to the spatial gradients in the bulk equations, the values at nodes 2 and n-1 depend on the boundary values at nodes 1 and n, respectively. The vector V is calculated as V = A-1B. The boundary condition for at the membrane entrance is imposed and held at that value at each time step. The boundary condition at the exit of the membrane at time t+1 is then obtained through the bulk values at time t+1. This procedure is repeated until the relative errors of the bulk and boundary values of both the mass fraction and mass flux between the previous and present loop become smaller than 10-4. The obtained values at time t+1 are stocked in the output matrix and used as the known values for the next time step t+2.
We use our model, given by Eqs. (10)–(12), (14)–(15), and solve it by means of a semi-implicite finite difference algorithm, explained in section 4.1. We take cin = 0.8, c0 = 0 and ϑ = 1 (we neglect the effect of the partition coefficient for the moment). The results are presented in Figures 1-5. Figures 1-3 show the evolution of the mass fraction across the membrane for Kn numbers Kn = 0.1, Kn = 1 and Kn = 10, respectively. The different Kn numbers stand for different sizes of the membrane with respect to the molecular mean free path of the drug molecules in the fluid present in the membrane (which can be taken as a watery solution). Furthermore, for each Kn number different fKn values are used: fKn = 0.1, fKn = 1 and fKn = 10.
We can see that for a given Kn number a higher fKn number (higher relaxation time with respect to a fixed diffusion characteristic time), the evolution of the mass flux is significantly slower, with a large gradient at the membrane entrance. This implies that the mass flux is very high at the membrane entrance but very low at the membrane exit. This effect is precisely referred to as non-locality. It is important here that by choosing the diffusion characteristic time as the scale time, any change in fKn implies a change in the relaxation time. This means that one compares systems with different relaxation times, which illustrates the effect of medication diffusion as a function of different medication-carrier fluid systems. A higher Kn number has the same non-local effect, albeit on the geometrical level (for the same system). This latter effect is not so clearly visible in Figures 1–3, but can be more appreciated looking to the evolution of the mass fraction at the membrane exit y = 1 in Figure 4.
Effectively, Figure 4 shows that a higher Kn number, for all fKn values, results into a lower mass fraction, which is representative for a lower drug diffusion/solubility in the membrane. The non-local effect at high Kn numbers is responsible for a non-Fickian mass fraction profile, which is indicated by a mass fraction at the exit of the membrane that is less than the initial mass fraction at steady state (this in comparison with Fickian behaviour, where for a zero-flux boundary condition, the profile should become flat). This means that for drugs that have large relaxation times or mean free paths, the drug uptake could be troublesome. It is interesting to note that the values shown in Figure 4 are exactly at the membrane exit, but that a mass fraction jump exists just before that exit. Not being well discernible in Figures 1-3, a zoom of some of those curves has been plotted in Figure 5.
This mass fraction jump can be seen clearly in Figure 5 and stems from the non-local behaviour from both the mass flux time derivative (∂tJ) and mass flux diffusion term (∇2J). Now, if the partition coefficient was not at unity, we would have not only a jump due to a non-locality at the membrane side of the interface, but also at the cell side, due the partition coefficient. Figure 6 shows the influence of the partition coefficient on the mass fraction evolution.
It can be seen that a higher partition coefficient results into a higher mass fraction and a higher drug uptake (higher mass fraction at the membrane exit). Since a zero-flux boundary condition (standing for a quick drug diffusion in the cell and not for a wall-like boundary, as is discussed under Eq (13)) is maintained at the membrane exit, the partition coefficient does not affect the mass fraction here. However, it does effect the quantitative value of the mass fraction just outside the membrane (in the cell). The higher the partition coefficient, the lower the mass fraction will be here. Also, as Figure 6 seems to show, a higher partition coefficient does not exhibit a different behaviour of the mass fraction at the membrane exit. That is, a higher partition coefficient gives a higher mass fraction of the drug at the membrane entrance, which translates itself into a higher mass fraction at the membrane exit. Nonetheless, this increase is cancelled by the partition coefficient itself at the membrane exit. This motivates using drugs that have a higher partition coefficient at the membrane entrance, but a lower one at the membrane exit. This would, of course, not lead to a higher total uptake (mass conservation implies that the total quantity of the uptake is equal to the administered quantity), but may lead to a higher efficiency if a higher mass fraction is needed for a better healing in a shorter time period. Having discussed drug uptake by means of drug diffusion through a cell membrane, it is also interesting to provide tools to drug release by scaffolds, applicable in the field of nanomedicine. This is treated in the next subsection.
We use our model, given by Eqs (20)–(22), and solve it in the same way as for the non-local drug diffusion model. We take c0 = 0.8 and present the results in Figures 7 and 8 for different Kn numbers and slip factors Bs. When we take Bs = 0, this stands for a zero slip length of the mass flux. This means that the mass flux at the wall of the pores in the scaffold is equal to zero. However, when Bs→∞, this means that there is full slip at the pore walls and that the mass flux is everywhere equal in the radial direction. We take a value Bs = 10 to be quite large, standing for "large slip" and compare its results to "no-slip". The slip behaviour of the mass flux is not clearly defined in the literature, but is rather often used in the framework of fluid flow, related to its velocity. The very definition of a mass flux in a porous medium being [23]
| J=ϕρf(vf−v) | (23) |
where ϕ is the porosity of the porous medium, ρf the density of the fluid flowing through it, vf the fluid velocity and v the barycentric velocity, we see that the mass flux is directly related to the fluid velocity. If we note that in a two-component system the barycentric velocity is expressed as ρv = ϕρfvf+(1-ϕ)ρsvs, wherein the total mass density ρ is ρ = ϕρf+(1-ϕ)ρs with ρs and vs the density and velocity of the solid component, respectively, and that in case of a scaffold vs = 0, we can quickly see that in (23) the only variable the mass flux is related directly to is the fluid velocity vf. As such, the discussion on the slip factor becomes more physical and understandable.
Figures 7 and 8 show that the Kn number tends to delay the drug release. This suggests that for pores of sizes of the order of magnitude of the mean free path of the drug molecules the drug release is slowed down. On the other hand, in some applications, it is necessary to have pores of such a size. Figure 8 shows that the slower drug release can be compensated by a higher slip factor, which accelerates the release of medication, illustrated by mass fraction profiles that evolve faster with respect to Figure 7. Figure 9 shows the mass fraction evolution for different Kn numbers for "no-slip" (Bs = 0) and "high-slip" (Bs = 10) at the membrane exit (y = 1). It can be seen that indeed, when Bs is increased, the reduction of the mass fraction due to an increase in the Kn number (from Kn = 0.1 to Kn = 1) is compensated, resulting into a similar profile. This confirms that by increasing the slip factor, which amounts to say altering the material properties of the scaffold pores to properties that reduce the interaction with the drug molecule solutions, compensates the reduction of drug release in nanopores.
This work presents new developments in Extended Non-Equilibrium Thermodynamics, proposing a model for studying drug uptake through cell membranes and drug release out of scaffolds. The theory is based on elevating fluxes to the same independent status of the usual variables. In this particular case, the mass flux has the same status as the mass fraction or concentration of the medication solution. This means that the mass flux follows its own evolution next to that of the mass fraction, being connected on an equal level via a species balance. By taking into account the relative importance of the relaxation time with respect to the diffusion characteristic time, and of the mean free path with respect to the system characteristic size, non-local effects are included in this way.
As for the diffusion through a cell membrane, it appeared that for sizes and relaxation times of the order of magnitude of or smaller than the mean free path (Kn = O(1)) and diffusion characteristic time (M = O(1)), respectively, leads to a lower drug uptake. This is due to the non-local effect, mentioned before. Non-Fickian profiles result from this, where it can be seen that at high non-locality a large mass fraction gradient is present close to the membrane entrance, but a nearly horizontal one near the exit. The latter contributes to the reduction of drug uptake. Furthermore, it has also been shown that at the membrane exit, a mass fraction jump (within the membrane) occurs that is not due to the boundary condition, but also due the non-locality effect. This is an effect that has not only been observed for heat transport [25], where it is already an established phenomenon, but also for mass transport [26]. The effect of the partition coefficient on the mass fraction is also studied, since they are typically introduced to account for the typical mass fraction jump (outside the membrane) observed in membrane diffusion studies, though still presenting experimental uncertainty to its precise value [26,27,28]. No particular effect has been observed other than that a higher partition coefficient causes a higher mass fraction jump at the membrane boundaries. If the medication could be chosen in such a way that the partition coefficient is high at the entrance and low at the exit, the efficiency of the medication could be increased if a higher concentration would lead to a better healing. Of course, mass conservation implies that this would not influence the total quantity of the uptake, but rather the speed of uptake, which could influence the kinetics of healing.
As for the drug release from the scaffold, the non-locality effect also showed to have a reducing impact on the release of drugs from the scaffold. It appeared that such a reduction could be countered by choosing the material of the scaffold in such a way that (besides biocompatibility, which is not a complicated task on its own) the interaction between the scaffold pore walls and the medication solution is weakened. This results into a higher slip factor enhancing drug release.
Future developments of this work include experimental demonstration, using the proposed model in this work. As quantitative agreement between experiments and models are difficult to find in the literature, comparison with experimental data should first be performed by nanomaterials of which material properties are well known and documented. As example, we can mention the delivery of nucleic acids with DNA nanostructures as new means of drug delivery by nanocarriers [29,30]. In the case of the nanocarriers, which can also be seen as scaffolds, the diffusion of the carriers through the cell membranes could be described by the set of Eqs (10)–(12) and (14)–(15). Subsequently, the delivery of the drugs out of the carriers into the cells could be described as release from a scaffold by Eqs (20)–(22). The difficulty lies in determining the interaction the medication can have with the nanocarrier, such as the DNA nanostructures [29]. In other words, it is important to know the slip length or slipping coefficient in Eq (17). Enhancement of drug release does not depend only on the release of the drug out of the carrier, but also the solubility of the carrier itself through the cell membranes. It is reported that DNA nanostructures can overcome the usually large circulation times of carriers [29,30]. The aforementioned discussion implies a necessity of measuring some material properties in order to use the presented model in this paper. In some cases, information about material properties could be obtained by in vitro experiments with a controlled experimental setup in a similar way performed by [31] for high-precision regulation of cell behaviour. This brings the intermediate step between modelling and full experimental demonstration to be identified by bio-mimicking cell membranes or scaffolds. Finally, such information can complete an experimental demonstration of our model.
Financial support from BelSPo is acknowledged.
The author declares there is no conflict of interest.
| [1] | Asafu-Adjaye J, Environmental Economics for Non-Economists. Techniques and Policies for Sustainable Development. World Scientific, 2005. |
| [2] |
Bridge G (2004) Contested Terrain: Mining and the Environment. Annu Rev Environ Resour 29: 205–259. doi: 10.1146/annurev.energy.28.011503.163434
|
| [3] | Roseland M (2000) Sustainable community development : integrating environmental , economic, and social objectives. Prog Plann 54: 73–132. |
| [4] |
Voulvoulis N, Skolout JWF, Oates CJ, et al. (2013) From chemical risk assessment to environmental resources management: the challenge for mining. Environ Sci Pollut Res 20: 7815–7826. doi: 10.1007/s11356-013-1785-8
|
| [5] |
Bryan BA, Grandgirard A, Ward JR (2010) Quantifying and exploring strategic regional priorities for managing natural capital and ecosystem services given multiple stakeholder perspectives. Ecosystems 13: 539–555. doi: 10.1007/s10021-010-9339-0
|
| [6] |
De Groot RS, Wilson MA, Boumans RMJ (2002) A typology for the classification, description and valuation of ecosystem functions, goods and services. Ecol Econ 41: 393–408. doi: 10.1016/S0921-8009(02)00089-7
|
| [7] | Boyd J, Banzhaf S (2007) What are ecosystem services? The need for standardized environmental accounting units. Ecol Econ 63: 616–626. |
| [8] |
Osuji LC, Erondu ES, Ogali RE (2010) Upstream petroleum degradation of mangroves and intertidal shores: The Niger Delta experience. Chem Biodivers 7: 116–128. doi: 10.1002/cbdv.200900203
|
| [9] |
Solomon F, Katz E, Lovel R (2008) Social dimensions of mining: Research, policy and practice challenges for the minerals industry in Australia. Resour Policy 33: 142–149. doi: 10.1016/j.resourpol.2008.01.005
|
| [10] | Humphreys M, Sachs J, Stiglitz J (2007) Escaping the Resource Curse. Columbia University Press. |
| [11] |
Bhattacharyya S, Hodler R (2010) Natural resources, democracy and corruption. Eur Econ Rev 54: 608–621. doi: 10.1016/j.euroecorev.2009.10.004
|
| [12] | Ross M (2006) A closer look at oil, diamonds, and civil wars. Annu Rev Polit Sci 8: 265–300. |
| [13] | Leite C, Weidmann J (1999) Does Mother Nature Corrupt? Natural Resources, Corruption, and Economic Growth. International Monetary Fund. |
| [14] | Ekins P (2002) Economic and Environmental Sustainability: The prospects for Green Growth. Routledge, London, 374p. |
| [15] | O'Lear S (2012) Oil and energy, In: Environmental Politics Scale and Power, S. O'Lear, Ed. Cambridge: Cambridge University Press, 55–86. |
| [16] |
Mehlum H, Moene KO, Torvik R (2002) Plunder & Protection Inc. J Peace Res 39: 447–459. doi: 10.1177/0022343302039004005
|
| [17] | Commission of the European Communities (2000) Promoting Sustainable development in the EU non-energy extractive industry. Communication from the commission.COM (2000) 265. |
| [18] |
Franks DM, Brereton D, Moran CJ (2013) The cumulative dimensions of impact in resource regions. Resour Policy 38: 640–647. doi: 10.1016/j.resourpol.2013.07.002
|
| [19] | Corvalan C, Hales S, McMichael AJ (2005) Ecosystems and Human Well-being: Health Synthesis. World health organization, 2005. |
| [20] | Laurenti R (2013) Applications of Systems Thinking within the Sustainability Domain : Product Design, Product Systems and Stakeholder Perspectives. Ph.D. dissertation, School of Industrial Engineering and Management KTH Royal Institute of Technology Stockholm, Sweden. |
| [21] |
Young N, Matthews R (2007) Resource economies and neoliberal experimentation: the reform of industry and community in rural British Columbia. AREA 39: 176–185. doi: 10.1111/j.1475-4762.2007.00739.x
|
| [22] |
Doloreux D, Amara N, Landry R (2008) Mapping regional and sectoral characteristics of knowledge-intensive business services: Evidence from the Province of Quebec (Canada). Growth Change 39: 464–496. doi: 10.1111/j.1468-2257.2008.00434.x
|
| [23] | Pan J, Oates CJ, Ihlenfeld C, et al. (2010) Screening and prioritisation of chemical risks from metal mining operations, identifying exposure media of concern. Environ Monit Assess 163 555–571. |
| [24] |
Moran CJ, Kunz NC (2014) Sustainability as it pertains to minerals and energy supply and demand: a new interpretative perspective for assessing progress. J Clean Prod 84: 16–26. doi: 10.1016/j.jclepro.2014.09.008
|
| [25] | MMSD (Mining, Minerals and Sustainable Development) (2002) The Report of the Mining, Minerals and Sustainable Development Project. International Institute for Environment and Development and World Business Council for Sustainable Development, TD195.M5 B74 2002. |
| [26] |
Peck P, Sinding K (2003) Environmental and social disclosure and data richness in the mining industry. Bus Strateg Environ 12: 131–146. doi: 10.1002/bse.358
|
| [27] |
Kitula AGN (2006) The environmental and socio-economic impacts of mining on local livelihoods in Tanzania: A case study of Geita District. J Clean Prod 14: 405–414. doi: 10.1016/j.jclepro.2004.01.012
|
| [28] | Kindzierski WB (1999) Importance of human environmental exposure to hazardous air pollutants from gas flares. Environ Rev 8: 41–62. |
| [29] | Rourke DO, Connolly S (2003) Just Oil? The Distribution of Environmental and Social Impacts of oil Production and consumption. Energy 28: 587–617. |
| [30] |
Azapagic A (2004) Developing a framework for sustainable development indicators for the mining and minerals industry. J Clean Prod 12: 639–662. doi: 10.1016/S0959-6526(03)00075-1
|
| [31] |
Kuemmerle T, Baskin L, Leitão PJ, et al. (2014) Potential impacts of oil and gas development and climate change on migratory reindeer calving grounds across the Russian Arctic. Divers Distrib 20: 416–429. doi: 10.1111/ddi.12167
|
| [32] |
Ite AE, Ibok UJ (2013) Gas Flaring and Venting Associated with Petroleum Exploration and Production in the Nigeria's Niger Delta. Am J Environ Prot 1: 70–77. doi: 10.12691/env-1-4-1
|
| [33] |
Esteves AM (2008) Mining and social development: Refocusing community investment using multi-criteria decision analysis. Resour Policy 33: 39–47. doi: 10.1016/j.resourpol.2008.01.002
|
| [34] |
Kohl B, Farthing L (2012) Material constraints to popular imaginaries: The extractive economy and resource nationalism in Bolivia. Polit Geogr 31: 225–235. doi: 10.1016/j.polgeo.2012.03.002
|
| [35] | Khoday K, Perch L (2012) Development from below: Social Accountability in Natural Resource Management. Working Paper, International Policy Centre for Inclusive Growth, No. 91. |
| [36] |
Isham J, Woolcock M, Pritchett L, et al. (2005) The varieties of resource experience: Natural resource export structures and the political economy of economic growth. World Bank Econ Rev 19: 141–174. doi: 10.1093/wber/lhi010
|
| [37] |
Auty RM (1995) Industrial Policy, Sectoral Maturation, and Postwar Economic Growth in Brazil: The Resource Curse Thesis. Econo Geography 71: 257–272. doi: 10.2307/144311
|
| [38] | Sachs JD,Warner AM (1997) Natural Resource Abundance and Economic Growth. Center for International Development and Harvard Institute for International Development. |
| [39] |
Sachs JD,Warner AM (2001) The curse of natural resources. Eur Econ Rev 45: 827–838. doi: 10.1016/S0014-2921(01)00125-8
|
| [40] | Australian Academy of Science (1995) Notes of the Fenner Conference on the Environment: Risk & Uncertainty in Environmental Management, 13–16 November, Becker House, Canberra. |
| [41] |
Bebbington DH (2013) Extraction, inequality and indigenous peoples: Insights from Bolivia. Environ Sci Policy 33: 438–446. doi: 10.1016/j.envsci.2012.07.027
|
| [42] |
Bulte EH, Damania R, Deacon RT (2005) Resource intensity, institutions, and development. World Dev 33: 1029–1044. doi: 10.1016/j.worlddev.2005.04.004
|
| [43] |
Akpabio EM, Akpan NS (2010) Governance and Oil Politics in Nigeria's Niger Delta: The Question of Distributive Equity. J Hum Ecol 30: 111–121. doi: 10.1080/09709274.2010.11906280
|
| [44] |
Jaskoski M (2014) Environmental Licensing and Conflict in Peru's Mining Sector: A Path-Dependent Analysis. World Dev 64: 873–883. doi: 10.1016/j.worlddev.2014.07.010
|
| [45] |
Tiainen H, Sairinen R, Novikov V (2014) Mining in the Chatkal Valley in Kyrgyzstan-Challenge of social sustainability. Resour Policy 39: 80–87. doi: 10.1016/j.resourpol.2013.11.005
|
| [46] | Gylfason T (2001) Natural Resources, Education And Economic Develpoment. Europ Econ Rev 45: 847–859. |
| [47] |
Cockx L, Francken N (2016) Natural resources: A curse on education spending? Energy Policy 92: 394–408. doi: 10.1016/j.enpol.2016.02.027
|
| [48] | Vásquez JA, Vega JMA, Matsuhiro B, et al. (1999) The ecological effects of mining discharges on subtidal habitats dominated by macroalgae in northern Chile: Population and community level studies. Sixteenth Internatinal Seaweed Symposium Springer, Dordrecht, 217–229. |
| [49] |
Urkidi L (2010) A glocal environmental movement against gold mining: Pascua-Lama in Chile. Ecol Econ 70: 219–227. doi: 10.1016/j.ecolecon.2010.05.004
|
| [50] |
Mendoza-Cantú A, Heydrich SC, Cervantes IS, et al. (2011) Identification of environmentally vulnerable areas with priority for prevention and management of pipeline crude oil spills. J Environ Manage 92: 1706–1713. doi: 10.1016/j.jenvman.2011.02.008
|
| [51] |
LeBillon P (2001) Angola's Political Economy of War: The Role of Oil and Diamonds , 1975–2000. Afr Aff (Lond) 100: 55–80. doi: 10.1093/afraf/100.398.55
|
| [52] | Reed K (2009) Crude Existence. Environment and the politics of oil in northern angola. University of California Press, 340p. |
| [53] | Mantovani ET (2017) Indigenous Warao (Orinoco Delta) contaminated by Corporacion Venezolana de Guayana and illegal mining, Venezuela. Ejolt - Environmental Justice, 2017. Available from: https://ejatlas.org/conflict/indigenas-warao-en-el-bajo-delta-del-orinoco-contaminados-por-desechos-de-la-corporacion-venezolana-de-guayana-y-de-la-mineria-ilegal. |
| [54] | Watts M (2004) Resource curse? governmentality, oil and power in the Niger Delta, Nigeria. Geopolitics 9: 50–80. |
| [55] |
Odeyemi O, Ogunseitan OA (1985) Petroleum industry and its pollution potential in Nigeria. Oil Petrochemical Pollut 2: 223–229. doi: 10.1016/S0143-7127(85)90218-2
|
| [56] |
Nduka JKC, Orisakwe OE, Ezenweke LO, et al. (2008) Acid rain phenomenon in niger delta region of Nigeria: economic, biodiversity, and public health concern. Scientific World J 8: 811–818. doi: 10.1100/tsw.2008.47
|
| [57] |
Orisakwe OE, Asomugha R, Obi E, et al. (2001) Ecotoxicological study of the Niger-Delta area of the River Niger. Bull Environ Contam Toxicol 66: 548–552. doi: 10.1007/s00128-001-0042-x
|
| [58] | Ana GREE, Sridhar MKC, Asuzu MC (2010) Environmental risk factors and hospital-based cancers in two Nigerian cities. J Pub Health Epidemiol 2: 216–223.. |
| [59] |
Ekpoh IJ, Obia AE (2010) The role of gas flaring in the rapid corrosion of zinc roofs in the Niger Delta Region of Nigeria. Environmentalist 30: 347–352. doi: 10.1007/s10669-010-9292-7
|
| [60] | Donwa P (2011) Environmental accounting and host community agitation in Nigeria: The petroleum industry experience. Int Rev Bus Res Pap 7: 98–108. |
| [61] |
Zaidi S (1994) Human health effects of oil development in the Ecuadorian Amazon: A challenge to legal thinking. Environ Impact Assess Rev 14: 337–348. doi: 10.1016/0195-9255(94)90005-1
|
| [62] |
Hurtig AK, San Sebastián M (2002) Geographical differences in cancer incidence in the Amazon basin of Ecuador in relation to residence near oil fields. Int J Epidemiol 31: 1021–1027. doi: 10.1093/ije/31.5.1021
|
| [63] | San Sebastián M, Hurtig AK (2005) Oil development and health in the Amazon basin of Ecuador: the popular epidemiology process Soc Sci Med 60: 799–807. |
| [64] |
Finer M, Jenkins CN, Pimm SL, et al. (2008) Oil and Gas Projects in the Western Amazon: Threats to Wilderness, Biodiversity, and Indigenous Peoples. PLoS One 3: 2932–2008. doi: 10.1371/journal.pone.0002932
|
| [65] |
Da silva EM, Pes-Aguiar MC, Navarro MDFT, et al. (1997) Impact of petroleum pollution on aquatic coastal ecosystems in Brazil. Environ Toxicol Chem 16: 112–118. doi: 10.1002/etc.5620160112
|
| [66] |
Zanardi E, Bı́cego MC, Weber RR (1999) Dissolved/dispersed Petroleum Aromatic Hydrocarbons in the São Sebastião Channel, São Paulo, Brazil. Mar Pollut Bull 38: 410–413. doi: 10.1016/S0025-326X(97)00194-X
|
| [67] |
Netalieva I, Wesseler J, Heijman W (2005) Health costs caused by oil extraction air emissions and the benefits from abatement: The case of Kazakhstan. Energy Policy 33: 1169–1177. doi: 10.1016/j.enpol.2003.11.014
|
| [68] |
Dahl C, Kuralbayeva K (2001) Energy and the environment in kazakhstan. Energy Policy 29: 429–440. doi: 10.1016/S0301-4215(00)00137-3
|
| [69] | Rigina O (2002) Environmental impact assessment of the mining and concentration activities in the Kola Peninsula, Russia by multidate remote sensing. Environ Monit Assess 75: 11–31. |
| [70] |
WalkerTR, Crittenden PD, Dauvalter VA, et al. (2009) Multiple indicators of human impacts on the environment in the Pechora Basin, north-eastern European Russia. Ecol Indic 9: 765–779. doi: 10.1016/j.ecolind.2008.09.008
|
| [71] |
Walker TR, Crittenden PD, Young SD, et al. (2006) An assessment of pollution impacts due to the oil and gas industries in the Pechora basin, north-eastern European Russia. Ecol Indic 6: 369–387. doi: 10.1016/j.ecolind.2005.03.015
|
| [72] |
Reimann C, Halleraker JH, Kashulina G, et al. (1999) Comparison of plant and precipitation chemistry in catchments with different levels of pollution on the Kola Peninsula, Russia. Sci Total Environ 243/244: 169–191. doi: 10.1016/S0048-9697(99)00390-3
|
| [73] |
Perreault T, Valdivia G (2010) Hydrocarbons, popular protest and national imaginaries: Ecuador and Bolivia in comparative context. Geoforum 41: 689–699. doi: 10.1016/j.geoforum.2010.04.004
|
| [74] |
Bebbington AJ (2014) Socio-environmental conflict: an opportunity for mining companies. J Clean Prod 84: 34. doi: 10.1016/j.jclepro.2014.08.108
|
| [75] |
Farrán A, Grimalt J, Albaigés J, et al. (1987) Assessment of petroleum pollution in a Mexican river by molecular markers and carbon isotope ratios. Mar Pollut Bull 18: 284–289. doi: 10.1016/0025-326X(87)90506-6
|
| [76] | Vázquez-Luna D (2012) Environmental bases on the exploitation of crude oil in Mexico. In: Crude Oil Exploration in the World. M.Younes, Ed. InTech, Available from: http://www.intechopen.com/books/crude-oil-exploration-in-the-world/environmental-bases-on-the-exploitationof-crude-oil-in-mexico. |
| [77] |
Hettler J, Irion G, Lehmann B (1997) Environmental impact of mining waste disposal on a tropical lowland river system: A case study on the Ok Tedi Mine, Papua New Guinea. Miner Depos 32: 280–291. doi: 10.1007/s001260050093
|
| [78] |
Swales S, Storey AW, Bakowa KA (2000) Temporal and spatial variations in fish catches in the Fly River system in Papua New Guinea and the possible effects of the Ok Tedi copper mine. Environ Biol Fishes 57: 75–95. doi: 10.1023/A:1007513906281
|
| [79] |
Hilson G (2002) An overview of land use conflicts in mining communities. Land use policy 19: 65–73. doi: 10.1016/S0264-8377(01)00043-6
|
| [80] |
Walton G, Barnett J (2007) The Ambiguities of 'Environmental' Conflict: Insights from the Tolukuma Gold Mine, Papua New Guinea. Soc Nat Resour 21: 1–16. doi: 10.1080/08941920701655635
|
| [81] | O'Faircheallaigh C (2008) Negotiating cultural heritage? Aboriginal-mining company agreements in Australia. Dev Change 39: 25–51. |
| [82] |
Wright IA,Wright S, Graham K, et al. (2011) Environmental protection and management: A water pollution case study within the Greater Blue Mountains World Heritage Area, Australia. Land use policy 28: 353–360. doi: 10.1016/j.landusepol.2010.07.002
|
| [83] | Weng Z, Mudd GM, Martin T, et al. (2012) Pollutant loads from coal mining in Australia: Discerning trends from the National Pollutant Inventory (NPI). Environ Sci Policy 19: 78–89. |
| [84] |
Ostrom E (2009) A general framework for analyzing sustainability of social-ecological systems. Science 325: 419–422. doi: 10.1126/science.1172133
|
| [85] | Barma NH, Kaiser K, MinhLe T, et al. (2012) Rents to Riches? The political economy of Natural Resource–Led Development. International Bank for Reconstruction and Development/International Development Association, The World Bank: Washington, DC, 298p. |
| [86] |
Acheson JM (2006) Institutional Failure in Resource Management. Annu Rev Anthropol 35: 117–134. doi: 10.1146/annurev.anthro.35.081705.123238
|
| [87] | Tarui N (2015) The Role of Institutions in Natural Resource Use. In: Balisacan et al., editors, Sustainable Economic Development. Oxford: Academic Press, pp. 123–136. |
| [88] |
Kumah A (2006) Sustainability and gold mining in the developing world. J Clean Prod 14: 315–323. doi: 10.1016/j.jclepro.2004.08.007
|
| [89] | Spitz K, Trudinger J (2008) Minerals, Wealth and Progress. In: Mining and the Environment From Ore to Metal, CRC Press, Boca Raton, FL.: CRC Press Taylor & Francis Group. 1–25. |
| [90] |
Kaup BZ (2008) Negotiating through nature: The resistant materiality and materiality of resistance in Bolivia's natural gas sector. Geoforum 39: 1734–1742. doi: 10.1016/j.geoforum.2008.04.007
|
| [91] |
Ako RT (2009) Nigeria's Land Use Act: An Anti-Thesis to Environmental Justice. J Afr Law 53: 289. doi: 10.1017/S0021855309990076
|
| [92] |
Aaron KK (2012) New corporate social responsibility models for oil companies in Nigeria's delta region: What challenges for sustainability? Prog Dev Stud 12: 259–27.. doi: 10.1177/146499341201200401
|
| [93] | Nzimande Z, Chauke H (2012) Sustainability through responsible environmental mining. J South. African Inst Min Metall 112: 135–139. |
| [94] | Punam CP, Dabalen D, Land BC, et al. (2015) Socioeconomic Impact of Mining on Local Communities in Africa. Washington, DC, 2015 Working Paper. Report No: ACS14621. |
| [95] | Puppim de Oliveira JA (2000) Implementing Environmental Policies in Developing Countries: Responding to the Environmental Impacts of Tourism Development by Creating Environmntally Protected Areas in Bahia, Brazil. Ph.D. dissertation, Massachusetts Institute of Technology. |
| [96] |
Ascher W (2001) Coping with complexity and organizational interests in natural resource management. Ecosystems 4: 742–757. doi: 10.1007/s10021-001-0043-y
|
| [97] |
Rodela R (2012) Advancing the deliberative turn in natural resource management: An analysis of discourses on the use of local resources. J Environ Manage 96: 26–34. doi: 10.1016/j.jenvman.2011.10.013
|
| [98] |
Bossel H (2001) Assessing Viability and Sustainability: a Systems-based Approach for Deriving Comprehensive Indicator Sets. Conserv Ecol 5: 12. doi: 10.5751/ES-00228-050112
|
| [99] |
Coria J, Sterner T (2011) Natural Resource Management: Challenges and Policy Options. Annu Rev Resour Econ 3: 203–230. doi: 10.1146/annurev-resource-083110-120131
|
| [100] | Franks D, Brereton D, Moran CJ (2009) Surrounded by Change–Collective Strategies for Managing the Cumulative Impacts of Multiple Mines. SDIMI Conf, 6–8. |
| [101] | Bellamy JA, Mcdonald GT, Syme GJ (1999) Policy Review Evaluating Integrated Resource Management Soc Nat Resour 12: 337–353. |
| [102] | Hoyle D (2009) ISO 9000 Quality Systems Handbook - updated for the ISO 9001:2008 standard, Sixth Edit. Elsevier Ltd, pp. 111–135. |
| [103] |
Moser T (2001) MNCs and sustainable business practice: The case of the Colombian and Peruvian petroleum industries. World Dev 29: 291–309. doi: 10.1016/S0305-750X(00)00094-2
|
| [104] |
Rammel C, Stagl S, Wilfing H (2007) Managing complex adaptive systems - A co-evolutionary perspective on natural resource management. Ecol Econ 63: 9–21. doi: 10.1016/j.ecolecon.2006.12.014
|
| [105] | Axelrod R, Cohen MD (1999) Harnessing Complexity: Organizational Implications of a Scientifc Frontier. New York: Free Press,184p. |
| [106] | Collier P, Hoeffler A (2005) Democracy and Natural Resource Rents. Working Paper Series. ESRC Global Poverty Research Group. |
| [107] |
Petak WJ (1981) Environmental management: A system approach. Environ Manage 5: 213–224. doi: 10.1007/BF01873280
|
| [108] |
Stave K (2010) Participatory system dynamics modeling for sustainable environmental management: Observations from four cases. Sustainability 2: 2762–2784. doi: 10.3390/su2092762
|
| [109] |
Mingers J, Rosenhead J (2004) Problem structuring methods in action. Eur J Oper Res 152: 530–554. doi: 10.1016/S0377-2217(03)00056-0
|
| [110] | Ludwig D, Hilborn R, Walters C (1993) Uncertainty, Resource Exploitation, and Conservation: Lessons from History. Ecol Appl 3: 547–549. |
| [111] |
Ballard C, Banks G (2003) The Anthropology of Mining. Annu Rev Anthropol 32: 287–313. doi: 10.1146/annurev.anthro.32.061002.093116
|
| [112] |
Hodge A (2014) Mining company performance and community conflict: moving beyond a seeming paradox. J Clean Prod 84: 27–33. doi: 10.1016/j.jclepro.2014.09.007
|
| [113] |
Jike VT (2004) Environmental Degradation, Social Disequilibrium, and the Dilemma of Sustainable Development in the Niger-Delta of Nigeria. J Black Stud 34: 686–701. doi: 10.1177/0021934703261934
|
| [114] |
Obinaju BE, Alaoma A, Martin FL (2014) Novel sensor technologies towards environmental health monitoring in urban environments: A case study in the Niger Delta (Nigeria). Environ Pollut 192: 222–231. doi: 10.1016/j.envpol.2014.02.004
|
| [115] |
Salako A, Sholeye O, Ayankoya S (2012) Oil spills and community health: Implications for resource limited settings. J Toxicol Environ Heal Sci 4: 145–150. doi: 10.5897/JTEHS12.056
|
| [116] |
Adesopo AA, Asaju AS (2004) Natural Resource Distribution , Agitation for Resource Control Right and the Practice of Federalism in Nigeria. J Human Ecol 15: 277–289. doi: 10.1080/09709274.2004.11905705
|
| [117] |
Ipingbemi O (2009) Socio-economic implications and environmental effects of oil spillage in some communities in the Niger delta. J Integr Environ Sci 6: 7–23. doi: 10.1080/15693430802650449
|
| [118] |
Lockwood M, Davidson J, Curtis A, et al. (2010) Governance Principles for Natural Resource Management. Soc Nat Resour 23: 986–1001. doi: 10.1080/08941920802178214
|
| [119] |
Mccubbin M, Cohen D (1999) A systemic and value-based approach to strategic reform of the mental health system. Heal Care Anal 7: 57–77. doi: 10.1023/A:1009443902415
|
| [120] |
McSherry R (2004) Practice development and health care governance: A recipe for modernization. J Nurs Manag 12: 137–146. doi: 10.1111/j.1365-2834.2004.00461.x
|
| [121] |
Bunch MJ (2003) Soft systems methodology and the ecosystem approach: A system study of the Cooum River and environs in Chennai, India. Environ Manage. 31: 182–197. doi: 10.1007/s00267-002-2721-8
|
| [122] |
Marshall R, Brown D (2003) The Strategy of Sustainability: A Systems Perspective of Environmental Initiatives. California Manage Rev 46: 101–126. doi: 10.2307/41166234
|
| [123] |
Kropff MJ, Bouma J, Jones JW (2001) Systems approaches for the design of sustainable agro-ecosystems. Agric Syst 70: 369–393. doi: 10.1016/S0308-521X(01)00052-X
|
| [124] |
Tejeda J, Ferreira S (2014) Applying systems thinking to analyze wind energy sustainability. Procedia Comput Sci 28: 213–220. doi: 10.1016/j.procs.2014.03.027
|
| [125] |
Checkland PB (1988) Information systems and systems thinking: Time to unite? Int J Inf Manage. 8: 239–248. doi: 10.1016/0268-4012(88)90031-X
|
| [126] |
Checkland PB, Haynes MG (1994) Varieties of Systems Thinking: the Case of Soft Systems Methodology. Syst Dyn Rev 10: 189–197. doi: 10.1002/sdr.4260100207
|
| [127] | Reed MS, Graves A, Dandy N, et al. (2009) Who's in and why? A typology of stakeholder analysis methods for natural resource management. J Environ Manage 90: 1933–1949. |
| [128] |
Magner J (2011) Tailored Watershed Assessment and Integrated Management (TWAIM): A Systems Thinking Approach. Water 3: 590–603. doi: 10.3390/w3020590
|
| 1. | Lu Xu, Chunlai Mu, Qiao Xin, Global boundedness and asymptotic behavior of solutions for a quasilinear chemotaxis model of multiple sclerosis with nonlinear signal secretion, 2023, 28, 1531-3492, 1215, 10.3934/dcdsb.2022118 |
Figures(2) / Tables(2)
Alozie Alaoma, Nikolaos Voulvoulis. Mineral resource active regions: The need for systems thinking in management[J]. AIMS Environmental Science, 2018, 5(2): 78-95. doi: 10.3934/environsci.2018.2.78
DownLoad:
