Research article Special Issues

Greening and “un”greening Adelaide, South Australia

  • The original design for Adelaide, the capital city of the state of South Australia, incorporated a green belt (known as the Park Lands) around the city centre, itself laid out on a one square mile (2.59 km2) grid and including five large public squares. The Park Lands provided a barrier to urban sprawl and covered approximately 9.31 km2, of which 1.53 km2 has been used subsequently for cultural institutions, railways, cemeteries, sporting facilities and other constructions. In addressing issues of greening pertaining to Adelaide, the Park Lands and its management represents a core element in the evolving history of the city's growth. This paper will consider some of the contradictions within this growth, examining the changing attitudes of government and the populace to the Park Lands and also to the increasing sprawl of the city. It can be argued that this sprawl has been antithetical to maintenance of biodiversity and principles of “greening”, not only during the main phase of expansion in the 1960s and 1970s but also in recent years when planned development on prime farmland and other “green” areas is contributing to problems for provision of transport infrastructure and generally reducing capacity for sustainability. The potential for conflict between the desire to maintain biodiversity versus protection for the growing number of people moving into bushfire risk areas is just one of several examples of problems arising as a result of a relaxed attitude to low-density expansion. In examining these problems the paper will present maps of the changing footprint of Adelaide and will elaborate new “greening” initiatives that include green roofs, new systems of water harvesting, community-supported agriculture and schemes directly aimed at creating low-carbon living. A consistent theme will be the contradictions within plans for the city between greening and “un”greening.

    Citation: Guy M. Robinson, Zhiling Liu. Greening and “un”greening Adelaide, South Australia[J]. AIMS Environmental Science, 2015, 2(3): 511-532. doi: 10.3934/environsci.2015.3.511

    Related Papers:

    [1] Simphiwe M. Simelane, Justin B. Munyakazi, Phumlani G. Dlamini, Oluwaseun F. Egbelowo . Projections of human papillomavirus vaccination and its impact on cervical cancer using the Caputo fractional derivative. Mathematical Biosciences and Engineering, 2023, 20(7): 11605-11626. doi: 10.3934/mbe.2023515
    [2] Emmanuel Addai, Mercy Ngungu, Musibau Abayomi Omoloye, Edmore Marinda . Modelling the impact of vaccination and environmental transmission on the dynamics of monkeypox virus under Caputo operator. Mathematical Biosciences and Engineering, 2023, 20(6): 10174-10199. doi: 10.3934/mbe.2023446
    [3] Rahat Zarin, Usa Wannasingha Humphries, Amir Khan, Aeshah A. Raezah . Computational modeling of fractional COVID-19 model by Haar wavelet collocation Methods with real data. Mathematical Biosciences and Engineering, 2023, 20(6): 11281-11312. doi: 10.3934/mbe.2023500
    [4] Adnan Sami, Amir Ali, Ramsha Shafqat, Nuttapol Pakkaranang, Mati ur Rahmamn . Analysis of food chain mathematical model under fractal fractional Caputo derivative. Mathematical Biosciences and Engineering, 2023, 20(2): 2094-2109. doi: 10.3934/mbe.2023097
    [5] Ritu Agarwal, Pooja Airan, Mohammad Sajid . Numerical and graphical simulation of the non-linear fractional dynamical system of bone mineralization. Mathematical Biosciences and Engineering, 2024, 21(4): 5138-5163. doi: 10.3934/mbe.2024227
    [6] Somayeh Fouladi, Mohammad Kohandel, Brydon Eastman . A comparison and calibration of integer and fractional-order models of COVID-19 with stratified public response. Mathematical Biosciences and Engineering, 2022, 19(12): 12792-12813. doi: 10.3934/mbe.2022597
    [7] Shuang-Hong Ma, Hai-Feng Huo . Global dynamics for a multi-group alcoholism model with public health education and alcoholism age. Mathematical Biosciences and Engineering, 2019, 16(3): 1683-1708. doi: 10.3934/mbe.2019080
    [8] Alia M. Alzubaidi, Hakeem A. Othman, Saif Ullah, Nisar Ahmad, Mohammad Mahtab Alam . Analysis of Monkeypox viral infection with human to animal transmission via a fractional and Fractal-fractional operators with power law kernel. Mathematical Biosciences and Engineering, 2023, 20(4): 6666-6690. doi: 10.3934/mbe.2023287
    [9] Elamin H. Elbasha . Model for hepatitis C virus transmissions. Mathematical Biosciences and Engineering, 2013, 10(4): 1045-1065. doi: 10.3934/mbe.2013.10.1045
    [10] Hardik Joshi, Brajesh Kumar Jha, Mehmet Yavuz . Modelling and analysis of fractional-order vaccination model for control of COVID-19 outbreak using real data. Mathematical Biosciences and Engineering, 2023, 20(1): 213-240. doi: 10.3934/mbe.2023010
  • The original design for Adelaide, the capital city of the state of South Australia, incorporated a green belt (known as the Park Lands) around the city centre, itself laid out on a one square mile (2.59 km2) grid and including five large public squares. The Park Lands provided a barrier to urban sprawl and covered approximately 9.31 km2, of which 1.53 km2 has been used subsequently for cultural institutions, railways, cemeteries, sporting facilities and other constructions. In addressing issues of greening pertaining to Adelaide, the Park Lands and its management represents a core element in the evolving history of the city's growth. This paper will consider some of the contradictions within this growth, examining the changing attitudes of government and the populace to the Park Lands and also to the increasing sprawl of the city. It can be argued that this sprawl has been antithetical to maintenance of biodiversity and principles of “greening”, not only during the main phase of expansion in the 1960s and 1970s but also in recent years when planned development on prime farmland and other “green” areas is contributing to problems for provision of transport infrastructure and generally reducing capacity for sustainability. The potential for conflict between the desire to maintain biodiversity versus protection for the growing number of people moving into bushfire risk areas is just one of several examples of problems arising as a result of a relaxed attitude to low-density expansion. In examining these problems the paper will present maps of the changing footprint of Adelaide and will elaborate new “greening” initiatives that include green roofs, new systems of water harvesting, community-supported agriculture and schemes directly aimed at creating low-carbon living. A consistent theme will be the contradictions within plans for the city between greening and “un”greening.


    Mathematical modeling and analysis can provide us with information about the characteristics of important phenomena and processes that arise in real-world situations, as well as predict possible outcomes. The study of mathematical models describing infectious diseases, in particular, has become one of the most powerful and effective approaches to analyzing [1], understanding [2], and predicting transmission mechanisms as well as infectious disease characteristics [3,4,5] Tungiasis is a parasitic skin disease caused by female sand flea Tunga penetrans also known as the jigger flea, an insect that lives mostly in dry sand or soil, and is found mostly in tropical and sub-tropical regions [6]. Jigger fleas generally attack feet or hands by making a burrow through the skin of humans or a variety of mammals that serves as reservoir hosts. While in the skin of the host, the impregnated female parasite continues to feed on blood by inserting its proboscis into dermal capillaries, releases her eggs and eventually dies [7]. Tungiasis can be transmitted by the infestation of human or animal reservoirs. Mammalian species that can act as reservoirs include cats, dogs, rats, bovines and pigs [8]. Mostly, the infection can be transmitted in the absence of animal reservoirs. This happens inside houses and classrooms without paved or sealed floors when skin touches floors or soil wherein adult sand fleas have developed [8].

    Unlike other diseases which attract attention and funding, not much is known about the effects of tungiasis due to minimal investigation into its dynamics and lack of awareness [9]. During the $ 66^{th} $ World Health Assembly in May 2013, a resolution to intensify and consolidate measures against neglected tropical diseases like tungiasis was agreed upon by the World Health Organization and its member states. The resolution contained in [10], resolves to plan investments and implement ways of improving the health and social well-being of the communities affected by these diseases. The burden of Tungiasis is mostly felt in rural communities and resource-poor urban neighbourhoods. In the Americas alone, it is estimated that more than 20 million people are at risk. If not attended to, Tungiasis disfigure and mutilate the feet. This then results in loss of morbidity, which has an adverse effect on the quality of life and on family economics if the affected person is the breadwinner.

    There are currently no known effective drugs for the treatment of tungiasis [9]. The current treatment involves identifying the parasite and mechanically removing it with a sharp-pointed sterile object like pins or needles, then applying an antiseptic dressing on the lesion. Tungiasis can also be effectively treated in medical facilities by surgically extracting submerged fleas under sterile conditions [11]. Using repellents, anti-inflammatory creams, applying anti-parasitic agents topically, and using disinfectants like potassium permanganate to wash affected areas are other techniques to deter sand fleas [7,11]. Moreover, there have been several attempts to verify the use of medicated lotion or ointments in the treatment of tungiasis. Heukelbach et. al. [12] showed that topical ivermectin, metrifonate or thiabendazole can each significantly reduce the number of lesions caused by embedded sand fleas. However, they identified that more research is required to ascertain the optimal doses and administration of these medicines. The importance of maintaining good personal cleanliness and wearing shoes in preventing jigger infection cannot be overstated. Children and the elderly need to be well-cared for, and low-income areas need to have access to cheap basic amenities and healthcare facilities [11].

    Exploring existing mathematical models on Tungiasis, we found that few mathematical models have been used to understand the disease transmission dynamics and various ways to control the epidemic. These models are integer-order models. Researchers have presented models incorporating hygiene as a control strategy and treatments [13], optimal control and Jigger flea reservoirs [14], public health education [11], and impact of protection [15] to model the Tungiasis. In the study by Muehlen et. al. [16], illiteracy was found to correlate with high tungiasis infestation. Hence, the impact of an educational campaign will be highly important in reducing the transmission dynamics of the disease. Public health education targeted to people with precarious living conditions in rural areas and in shanty towns may be a useful approach to curbing disease spread. Public health education is pivotal in the management and control of infectious and non-infectious diseases. This manuscript presents and analyses a nonlinear fractional model of tungiasis propagation dynamics with the impact of public health education for the first time.

    Of late, fractional calculus has been employed for numerous problems in engineering and science. Valuable monographs of fractional differential equations (FDE) and their utilisation are discussed in [17,18,19]. Fractional Caputo derivative has been used to model various infectious and non-infectious diseases to predict the different types of disease transmission dynamics and possible control strategies (see for example [20,21,22]) and many others. The appropriate capture of real-life dynamics at any time is one of the main reasons researchers opt to use fractional operators. Also, dealing with FDE systems allows us to model real-life circumstances incorporating history, memory and heredity. Hence, FDE gives us more realistic ways of describing disease dynamics because they capture complicated behavioural patterns of these biological systems. There are other fractional derivative definitions such as Riemann-Liouville fractional derivative, Caputo–Fabrizio derivative, Atangana-Baleanu fractional derivative, and many more. These derivative operators have additional properties compared to the Caputo derivative used in this study.

    In literature, there are not many mathematical models using the fractional derivative to predict or model the spread of tungiasis in the human population. This, to a large extent, is influenced by the fact that tungiasis is also a neglected disease according to WHO [11]. In this study, we use the Caputo fractional derivative to model tungiasis dynamics. We are motivated by the fact that following studies in literature, modelling disease dynamics using fractional derivatives is moreover fitting and intuitive. We describe the model using the Caputo derivative because of its ease of implementation as well as the property that for constants, the Caputo derivative is zero. Another attribute of the Caputo operator is that, in deriving the model, it allows using standard initial conditions described as derivatives of integer order [23]. We first show that the fractional model governing equations have non-negative bounded and distinctive solutions. Steady-state existence and stability are investigated on the basis of the basic reproduction number. The tungiasis FDE model is solved by using the generalised Adams–Bashforth-Moulton mechanism. The obtained results indicated that investigating the disease dynamics in fractional order provides additional insights into tungiasis propagation. In particular, we discuss how public health education, treatment and the contact rate between individuals affect overall disease propagation in the population.

    The manuscript is organised as follows: Section 2 provides basics and preliminaries of the Caputo fractional derivative. Development of the model is presented in Section 3. Model properties such as its positivity, boundedness, steady states and stability are discussed in Section 4. In Section 5, the generalised Adams-Bashforth-Moulton method is implemented in the developed fractional model. Section 6 presents a discussion of the numerical results obtained to corroborate the analytical investigation. Lastly, Section 7 provides a conclusion of the findings of the present investigation.

    Definition 1. The fractional integral of order $ \alpha > 0 $ of a function $ x \in C ([0, b], \mathbb{R}^+) $, where $ b > 0 $, is given by [24,25,26,27]

    $ Iαx(t)=1Γ(α)t0(tτ)α1x(τ)dτ,t>0.
    $
    (2.1)

    Definition 2. Suppose $ x \in C^1 ([0, b], \mathbb{R}^+) $, where $ b > 0 $. The Caputo fractional derivative of order $ \alpha > 0 $ of $ x(t) $ is given by [18,23,25]

    $ CDαtx(t)=InαDnx(t)={1Γ(nα)t0(tτ)nα1dndτnx(τ)dτ,n1<α<n,(ddt)nx(t),α=n,
    $
    (2.2)

    where $ n $ is a positive integer.

    Operator $ ^{C} D _{t} ^{\alpha} $ satisfies the following properties:

    (i). ${} ^{C} D _{t} ^{\alpha} I^{\alpha} x(t) = x(t) $

    (ii). $ {I^\alpha }^CD_t^\alpha x(t) = x(t) - \sum\limits_{\vartheta = 0}^{n - 1} {\frac{{{x^{(\vartheta)}}(a)}}{{\vartheta!}}} {(t - a)^\vartheta }, t > a. $

    For $ 0 < \alpha \le 1 $, the Caputo fractional derivative (2.2) of order $ \alpha > 0 $ reduces to

    $ CDαtx(t)=1Γ(1α)t0(tτ)αx(τ)dτ.
    $
    (2.3)

    Next, we implement definition (2.3) to the non-linear tungiasis model.

    The human population is divided into five sub-populations according to their disease status: susceptible population ($ S(t) $), infected people but are unaware that they are infected ($ I_u (t) $), infected people and know that they are infected ($ I_a (t) $), treated population ($ T(t) $) and recovered population ($ R(t) $) such that

    $ N(t)=S(t)+Iu(t)+Ia(t)+T(t)+R(t).
    $
    (3.1)

    The model considers tungiasis dynamics, which is predominantly transmitted by sand fleas. People enter the susceptible sub-population by recruitment, $ \pi $, either by birth or immigration. $ \beta $ is the contact rate, $ \mu $ is the natural mortality rate, $ \nu $ is the recovery rate due to treatment, $ \gamma $ is treatment waning $ \phi $ is progression rate from educated to infectious class, $ \epsilon $ is the literacy level on tungiasis dynamics, $ \omega $ is the treatment rate, and $ d_u $ and $ d_a $ are the disease induced death rate [28]. The incidence rate is defined as the number of new cases generated per unit time during the disease dynamics [29,30]. In mathematical modeling, the incidence is known for playing an essential role in driving the infection during an epidemic. Numerous other studies have used a bilinear incidence rate to study disease models (see [32,33]). The bilinear incidence implies that the number of new cases per time becomes saturated within the total population. To bypass this complexity, we combine the bilinear and the saturation incidence function, which we define to be $ \frac{\omega I_a}{1+\omega I_a} $ in our model, as it gives a better insight into the dynamics of the disease progression and reduction in disease incidence rate [29,30]. The state variables and parameters are summarized in Table 1. Figure 1 provides a graphical representation of tungiasis dynamics.

    Table 1.  Model state variable and parameter descriptions.
    State variables Description Units
    S(t) susceptible to jigger infestation people per time
    $ I_u (t) $ individuals infected with jigger but unaware people per time
    $ I_a (t) $ infected individuals with jigger and aware people per time
    T(t) treated individuals people per time
    R(t) individuals who recovered from jigger infestation people per time
    Parameters Description Units
    $ \pi $ recruitment rate Individual per time
    $ \beta $ contact rate per time
    $ \mu $ natural mortality rate per time
    $ \nu $ recovery rate from treatment per time
    $ \gamma $ progression rate of recovered to susceptible per time
    $ \phi $ rate of progression from $ I_a $ to $ I_b $ class per time
    $ \epsilon $ literacy level per time
    $ \omega $ treatment rate per time
    $ d_u, d_a $ disease-induced death rate per time

     | Show Table
    DownLoad: CSV
    Figure 1.  Transmission diagram of tungiasis with saturated treatment.

    The model assumptions are:

    1. We assume that infected individuals undergo a treatment course before they recover [13] i.e. no natural recovery.

    2. Define $ \theta = 1 - \epsilon $, where $ \epsilon $ is the literacy level on Tungiasis dynamics. In essence, $ \theta $ can be viewed as the ignorance level to educational campaigns carried out by health officials to educate people about Tungiasis transmission dynamics. Hence, as more people get educated on the disease dynamics, $ \epsilon $ increases and thus $ \theta $ decreases, which in turn reduces the number of people flowing to class $ I_a $.

    The combination of the model diagram Figure 1, Table 1, model description and assumptions yields the following system of non-linear ODEs:

    $ dSdt=πθβSIuβSIa+γRμSdIudt=θβSIuϕIu(μ+du)IudIadt=βSIa+ϕIuωIa1+ωIa(μ+da)IadTdt=ωIa1+ωIa(μ+ν)TdRdt=νT(μ+γ)R},
    $
    (3.2)

    subject to

    $ S(0)=S00,Iu(0)=Iu00,Ia(0)=Ia00,T(0)=T00,R(0)=R00.
    $
    (3.3)

    Tungiasis model (3.2) is then converted into an FDE model of order $ \alpha $ given by

    $ CDαtS(t)=πβSIaβθSIu+γRμSCDαtIu(t)=βθSIu(μ+ϕ+du)IuCDαtIa(t)=βSIa+ϕIuωIa1+ωIa(μ+da)IaCDαtT(t)=ωIa1+ωIa(μ+ν)TCDαtR(t)=νT(μ+γ)R},
    $
    (3.4)

    where $ 0 < \alpha \le 1 $. Since model (3.2) traces the human beings, all variables and parameters are positive for all $ t $. In the following sections, we explore the dynamics of model (3.4).

    Let $ \Phi = \left\lbrace x(t) \in \mathbb{R}^5_{+} : x(t) \ge 0 \right\rbrace $ and $ x(t) = \left[S(t), I_u(t), I_a(t), T(t), R(t) \right] ^T $. To show that the model solutions are positive, recall the following theorem and consider the corollary below.

    Lemma 1 (Generalized Mean Value Theorem [34]). Suppose that $ x(t) \in C[a, b] $ and the Caputo derivative $ ^{C} D _{t} ^{\alpha} x(t) \in C(a, b] $ for $ 0 < \alpha \le 1 $, then we have

    $ x(t)=x(a)+1Γ(α)CDαtx(ξ)(ta)α
    $
    (4.1)

    with $ 0 \le \xi \le t $ for $ t \in (a, b] $.

    Corollary 1. Suppose $ x(t) \in C[0, b] $ and the Caputo derivative $ ^{C} D _{t} ^{\alpha} x(t) \in C(0, b] $ for $ 0 < \alpha \le 1 $. From Lemma 1, if $ ^{C} D _{t} ^{\alpha} x(t) \ge 0 $ for $ t \in (0, b) $, then the function $ x(t) $ is non-decreasing and if $ ^{C} D _{t} ^{\alpha} x(t) \le 0 $ for $ t \in (0, b) $, then the function $ x(t) $ is non-increasing for all $ t \in [0, b] $ [24].

    Proof. The proof to the Corollary follows from Lemma 1.

    Theorem 1. Let $ (S, I_u, I_a, T, R) $ be a solution of system (3.4) subject to initial conditions (3.3) on $ t > 0 $, and the closed set $ \Phi = \left\lbrace (S, I_u, I_a, T, R) \in \mathbb{R}^5_{+} : 0 < N (t) \leq \frac{\pi }{\mu } \right\rbrace $. The set $ \Phi $ is positively invariant plus attracting for the dynamics described by the system (3.4).

    Proof. According to Theorem 3.1 and Remark 3.2 by Lin [35], we can ascertain the solution of system (3.4) subject to initial conditions (3.3) on $ (0, \infty) $. The solution exists and is distinctive. Thereafter, we have shown that $ \Phi $ is positively invariant. From (3.4), we have that

    $ CDαtS(t)|S=0=π+αR>0,CDαtIu(t)|Iu=0=0,CDαtIa(t)|Ia=0=ϕIu0,CDαtT(t)|T=0=ωIa1+ωIa0,CDαtR(t)|R=0=νT0.
    $

    The solution will remain in $ \Phi $ according to Corollary 1. $ \Phi $ is attracting because the vector field gravitates into $ \Phi $ at each hyperplane bounding the non-negative orthant. From the model system (3.4), the entire human population dynamics are governed by

    $ CDαtN(t)=πμNduIudaIa.
    $
    (4.2)

    It then follows from Eq. (4.2) that

    $ CDαtN(t)πμN.
    $
    (4.3)

    Hence,

    $ limtSupNπμ.
    $
    (4.4)

    Thus, by Lemma 1, solutions of (3.4) are all in $ \Phi = \left\lbrace (S, I_u, I_a, T, R) \in \mathbb{R}^5_{+} : 0 < N (t) \leq \frac{\pi }{\mu } \right\rbrace $. This means that $ N(t) $ is bounded by $ \frac{\pi}{\mu} $.

    Particularly, disease models have a disease-free equilibrium $ (DFE) $ plus an endemic equilibrium $ (EE) $ (at least). Their stability is expressed in terms of the reproduction number. To find these equilibrium points, set the derivatives of (3.4) to zero. This gives a disease-free steady state ($ \mathbb{E}_{0}^{*} $) given as

    $ E0=(S0,Iu0,Ia0,T0,R0)=(πμ,0,0,0,0).
    $
    (4.5)

    The endemic equilibrium ($ \mathbb{E}_{1}^{*} $) is

    $ E1=(S,Iu,Ia,T,R),
    $
    (4.6)

    where

    $ S=Q0βθ,Iu=Iaω(α(μQ2Q0βθν)+βIaQ3Q2Q0+μ2Q2Q0)πβθQ3(Iaω+1)Q2+Q3Q2Q0(βIa+μ)βθQ3(Iaω+1)Q2Q0,T=ωIaQ2(1+ωIa),R=νωIaQ3(1+ωIa)Q2,
    $
    (4.7)

    with $ Q_{0} = (\mu+\phi+d_u) $, $ Q_{1} = (\mu+d_a) $, $ Q_{2} = (\mu+\nu) $ and $ Q_{3} = (\mu+\gamma) $.

    The Tungiasis basic reproduction number, which is the average number of humans that can be infected by an infectious individual, is thus calculated using the methodology in [31] as follows. Considering the transfer and the transition matrix, we have that

    $ {F} = \left( βθπμ000βπμ0000
    \right)\; \; \; \text{and}\; \; \; {V} = \left( Q000ϕω+Q100ωQ2
    \right). $

    Therefore, the Tungiasis basic reproduction number represented by $ \mathcal{R}_{t} $ is the spectral radius of the next generation matrix generated by computing $ FV^{-1}. $ Hence,

    $ \mathcal{R}_{t} = \rho(FV^{-1}) = \text{max}\biggl\{\mathcal{R}_{t1}, \mathcal{R}_{t2}\biggr\}. $

    where

    $ \mathscr{R}_{t1} = \dfrac{\beta \pi \theta}{\mu Q_{0}}\; \; \text{and}\; \; \mathscr{R}_{t2} = \dfrac{\beta \pi }{\mu(\omega+Q_{1}) }. $

    Theorem 1. The Tungiasis model (3.4) is locally asymptotically stable (LAS) at $ \mathbb{E}_{0}^{*} $ whenever $ \mathcal{R}_{t} < 1 $ and unstable otherwise.

    Proof. Evaluating the Jacobian matrix of system (3.4) gives the following matrix

    $ JE0=(μβθπμβπμ0α0βθπμQ00000ϕβπμωQ10000ωQ20000νQ3),
    $
    (4.8)

    Matrix (4.8) has five distinct negative eigenvalues given by $ -\mu $, $ -Q_{2} $, $ -Q_{3} $, $ -Q_{0}\Big(1-\mathscr{R}_{t1}\Big) $, and $ -(\omega+Q_{1})(1-\mathscr{R}_{t2}) $. Hence, the local stability of $ \mathbb{E}_{0}^{*} $ is established if $ \mathscr{R}_{t1} < 1 $ and $ \mathscr{R}_{t2} < 1 $, i.e, $ E_{0}^{*} $ is asymptotically stable since $ \mathcal{R}_{t1}, \mathscr{R}_{t2}\subset\mathscr{R}_{t}. $

    Naik et al. [23] gives a brief description of the generalised Adams-Bashford-Moulton method. We use this method to solve model the nonlinear model (3.4) numerically. We utilize Matlab R2019b [36] and Wolfram 9.0 [37] for the codes and simulations.

    We first briefly describe the generalised Adam-Bashforth-Moulton technique utilised in solving the model (3.4). Consider

    $ CDαtx(t)=f(t,x(t)),0tT,
    $
    (5.1)

    subject to

    $ xj(0)=xj0,j=0,1,2,[α]1.
    $
    (5.2)

    From property $ (ii.) $ and the fractional integral definition, as defined earlier, operate Eq. (5.1) with this integral operator setting the lower limit $ a = 0 $. To obtain solution $ x(t) $, solve the equation:

    $ x(t)=[α]1j=0xj0j!tj+1Γ(α)t0(tτ)α1f(τ,x(τ))dτ.
    $
    (5.3)

    Set the step-size $ h = \frac{T}{N} $ so that $ t_n = n h $ and $ n = 0, 1, 2, \ldots , N \in \mathbb{Z}^{+} $ for $ 0 \le t \le T $. The Adams-Bashforth-Moulton scheme is utilised to integrate Eq. (5.3) (see [18,23]). Equation (5.3) is then discretized as

    $ xh(tn+1)=[α]1j=0xj0j!tjn+1+hαΓ(α+2)f(tn+1,xph(tn+1))+hαΓ(α+2)nq=0aq,n+1f(tq,xh(tq)),
    $
    (5.4)

    where

    $ aq,n+1={nα+1(nα)(n+α)α,ifq=0,(nq+2)α+1+(nq)α+12(nq+1)α+1,if0<qn,1,ifq=n+1,
    $
    (5.5)

    and the prediction $ x_h ^p (t_{n+1}) $ is given by

    $ xph(tn+1)=[α]1j=0xj0j!tjn+1+1Γ(α)nq=0bq,n+1f(tq,xh(tq)),
    $
    (5.6)

    where

    $ bq,n+1=hαα((n+1q)α(nq)α).
    $
    (5.7)

    Since $ x_h (t_p) $ approximates $ x (t_p) $, the error estimate is

    $ maxq=0,1,2,,m|x(tp)xh(tp)|=O(hp),
    $
    (5.8)

    with $ p = \mbox{min} (2, 1 + \alpha) $.

    We then use the technique described above to numerically solve (3.4). From Eq. (5.4), the scheme for (3.4) is

    $ S(tn+1)=S0+hαΓ(α+2)fS(tn+1,Sp(tn+1),Ipu(tn+1),Ipa(tn+1),Tp(tn+1),Rp(tn+1))+hαΓ(α+2)nq=0aq,n+1fS(tq,S(tq),Iu(tq),Ia(tq),T(tq),R(tq)),Iu(tn+1)=Iu0+hαΓ(α+2)fIu(tn+1,Sp(tn+1),Ipu(tn+1),Ipa(tn+1),Tp(tn+1),Rp(tn+1))+hαΓ(α+2)nq=0aq,n+1fIu(tq,S(tq),Iu(tq),Ia(tq),T(tq),R(tq)),Ia(tn+1)=Ia0+hαΓ(α+2)fIa(tn+1,Sp(tn+1),Ipu(tn+1),Ipa(tn+1),Tp(tn+1),Rp(tn+1))+hαΓ(α+2)nq=0aq,n+1fIa(tq,S(tq),Iu(tq),Ia(tq),T(tq),R(tq)),T(tn+1)=T0+hαΓ(α+2)fT(tn+1,Sp(tn+1),Ipu(tn+1),Ipa(tn+1),Tp(tn+1),Rp(tn+1))+hαΓ(α+2)nq=0aq,n+1fT(tq,S(tq),Iu(tq),Ia(tq),T(tq),R(tq)),R(tn+1)=R0+hαΓ(α+2)fR(tn+1,Sp(tn+1),Ipu(tn+1),Ipa(tn+1),Tp(tn+1),Rp(tn+1))+hαΓ(α+2)nq=0aq,n+1fR(tq,S(tq),Iu(tq),Ia(tq),T(tq),R(tq)),
    $

    where

    $ Sp(tn+1)=S0+1Γ(α)nq=0bq,n+1fS(tq,S(tq),Iu(tq),Ia(tq),T(tq),R(tq)),Ipu(tn+1)=Iu0+1Γ(α)nq=0bq,n+1fIu(tq,S(tq),Iu(tq),Ia(tq),T(tq),R(tq)),Ipa(tn+1)=Ia0+1Γ(α)nq=0bq,n+1fIa(tq,S(tq),Iu(tq),Ia(tq),T(tq),R(tq)),Tp(tn+1)=T0+1Γ(α)nq=0bq,n+1fT(tq,S(tq),Iu(tq),Ia(tq),T(tq),R(tq)),Rp(tn+1)=R0+1Γ(α)nq=0bq,n+1fR(tq,S(tq),Iu(tq),Ia(tq),T(tq),R(tq)).
    $

    Furthermore, the functions $ f_i \left(t_q, S (t_q), I_u (t_q), I_a (t_q), T (t_q), R (t_q) \right) $, $ i = S, I_u, I_a, T, R $ are computed as follows

    $ fS(tq,S(tq),Iu(tq),Ia(tq),T(tq),R(tq))=πβSIaβθSIu+γRμS,
    $
    (5.9)
    $ fIu(tq,S(tq),Iu(tq),Ia(tq),T(tq),R(tq))=βθSIu(μ+ϕ+du)Iu,
    $
    (5.10)
    $ fIa(tq,S(tq),Iu(tq),Ia(tq),T(tq),R(tq))=βSIa+ϕIuωIa1+ωIa(μ+da)Ia,
    $
    (5.11)
    $ fT(tq,S(tq),Iu(tq),Ia(tq),T(tq),R(tq))=ωIa1+ωIa(μ+ν)T,
    $
    (5.12)
    $ fR(tq,S(tq),Iu(tq),Ia(tq),T(tq),R(tq))=νT(μ+γ)R.
    $
    (5.13)

    Likewise, the functions

    $ fS(tn+1,Sp(tn+1),Ipu(tn+1),Ipa(tn+1),Tp(tn+1),Rp(tn+1)),fIu(tn+1,Sp(tn+1),Ipu(tn+1),Ipa(tn+1),Tp(tn+1),Rp(tn+1)),fIa(tn+1,Sp(tn+1),Ipu(tn+1),Ipa(tn+1),Tp(tn+1),Rp(tn+1)),fT(tn+1,Sp(tn+1),Ipu(tn+1),Ipa(tn+1),Tp(tn+1),Rp(tn+1)),fR(tn+1,Sp(tn+1),Ipu(tn+1),Ipa(tn+1),Tp(tn+1),Rp(tn+1)),
    $

    are computed from Eq. (5.9)-(5.13), respectively, at the points $ t_{n+1} $ for $ n = 1, 2, 3, \ldots m $. We thus present numerical simulation in the next section.

    Here, we parameterise the model using the available published literature. Our parameterisation uses the historical emergence of tungiasis in sub-Saharan Africa. The parameter values, at which the model is in its endemic state, are summarised in Table 3.

    Using the initial data in Table 2 and the parameters in Table 3, we obtain trajectories for the fractional model (3.4) as illustrated in Figure 2. Considering $ 0\le t \le 100 $ and initial values as described in Table 2, we give plots of the sub-populations for various fractional orders $ \alpha $. From the values in Table 1, we have $ \mathcal{R}_{t1} = 0.0007356 < 1 $, $ \mathcal{R}_{t2} = 0.0699499 < 1 $ and subsequently $ \mathcal{R}_0 < 1 $, implying that tungiasis eventually die out even with no intervention in place. Theoretically, if $ \mathcal{R}_0 > 1 $, then tungiasis will persist in the population until the introduction of intervention strategies to reduce $ \mathcal{R}_0 $. The fractional order $ \alpha $ considered varies from $ \alpha = 0.75 $ to $ \alpha = 1 $ as depicted in Figure 2. When $ \alpha = 1 $, we have the classical integer order model.

    Table 2.  The initial data numerical simulation.
    $ S(0) $ $ I_u (0) $ $ I_a (0) $ $ T(0) $ $ R(0) $
    900 300 120 20 5

     | Show Table
    DownLoad: CSV
    Table 3.  Model parameter values used for numerical simulations.
    Symbol Value Source Symbol Value Source
    $ \pi $ $ 4.4 \times 10 ^{-3} $ [11] $ \gamma $ $ 4.27 \times 10 ^{-1} $ [11]
    $ \beta $ $ 1.4989 \times 10 ^{-2} $ [11] $ \phi $ 0.001
    $ \mu $ $ 1.6 \times 10 ^{-2} $ [11] $ \theta $ $ 1.431 \times 10 ^{-2} $ [11]
    $ \nu $ 0.05 $ d_a $ $ 4.2 \times 10 ^{-2} $ [11]
    $ \omega $ 0.01 $ d_u $ $ 6.2 \times 10 ^{-2} $ Estimated

     | Show Table
    DownLoad: CSV
    Figure 2.  Projections of tungiasis model (3.4) for different $ \alpha $ using parameters in Table 1. We observe that population classes converge quicker to their steady states when $ \alpha $ is increased.

    The graphs demonstrate that fractional order significantly affects the dynamical behaviour of all the sub-populations. The memory effect of the model system grows as $ \alpha $ decreases from 1 to 0.75. Over time, the amount of educated and infected people increases when $ \alpha $ decreases from 1 to 0.75. In a general sense, low education levels on tungiasis and its transmission dynamics thereof result in a delay in identifying tungiasis-exposed people, a rise in the tungiasis undiagnosed individuals, accelerated progression of the disease in the population, as well as a rise in individuals eventually diagnosed with the disease. Conversely, being educated on tungiasis transmission dynamics trigger the susceptibles to take precautionary measures such as a change in behaviour and regular application of a repellent based on coconut oil to reduce the probability of infection. This consequently results in minimal infection growth in society. We may then conclude from Figure 2 that the derivative order $ \alpha $ doubles as a precautionary measure in the fight against infection propagation.

    Worth noting is that each sub-population projection maintains the same trend when $ \alpha $ is varied, although their actual values are moderately different. In Figure 2a, we observe that the $ S(t) $ curves are decreasing over time and converge to equilibrium point $ S_0 = 0.271605 $. Figures 2b-2e indicate that all $ I_u (t) $, $ I_a (t) $, $ T(t) $ and $ R(t) $ decrease with time and finally converge to the equilibrium point ($ I_{u_0}, I_{a_0}, T_0, R_0 $) = (0, 0, 0, 0). Furthermore, the existence of attractors for different fractional order $ \alpha $ are given in Figure 3 for the case $ \mathcal{R}_0 < 1 $. Figure 3 indicates that each population class exists and enters a state of remaining unchanged indefinitely.

    Figure 3.  Existence of attractors for different fractional order $ \alpha $ in model (3.4). We observe decreasing the $ \alpha $ does not impact the dynamics of the system.

    In summary, we see from Figure 2 and 3 that population classes converge quicker to their steady states when $ \alpha $ is increased. Figure 2 indicates the numerical solution converging to the DFE $ \mathbb{E}_{0}^{*} = (0.271605, 0, 0, 0, 0) $. This displays a positive outcome that the cases of Tungiasis-infected individuals tend to zero as $ t \to \infty $. Thus, we can conclude that the derivative order $ \alpha $ captures the role of experience or knowledge that individuals have on the disease's history. The obtained approximate results confirm that FDEs provide plentiful dynamics and tend to better express biological systems when compared to their integer model versions. As a result, the presence of the fractional order $ \alpha $ makes our results more logical than perhaps what integer order-only models would present.

    The impact of public health education is quantified by simulating the different sub-populations for different values of $ \theta $. Recall that $ \theta = 1 - \epsilon $. If $ \epsilon $ is quantifies literacy on public health issues, then $ \theta $ is the ignorance or illiteracy level. Figures 4 and 5 show plots of $ I_u $ and $ I_a $ for $ \theta = 0.3 $ and $ \theta = 0.9 $, respectively, for different fractional order values. We see that increasing the fractional order from 0.75 to 1 results in a decrease in equilibrium values of the infected but unaware and infected but aware populations. Figure 6 shows a plot of the $ I_u $ and $ I_a $ populations for different values of $ \theta $ and a fixed $ \alpha $. The figure indicates that the equilibrium number of infected and unaware individuals generally increases when $ \theta $ increases.

    Figure 4.  Results indicating the effect of $ \theta = 0.3 $ on the $ I_u $ and $ I_a $ fractional dynamics for different fractional order $ \alpha $.
    Figure 5.  Numerical results showing the effect of $ \theta = 0.9 $ on the $ I_u $ and $ I_a $ fractional dynamics for different fractional order $ \alpha $.
    Figure 6.  Results showing the effect of $ \theta $ on the $ I_u $ and $ I_a $ fractional dynamics for a fixed fractional order $ \alpha $.

    The relationship between $ \theta $ and $ \epsilon $ implies that $ I_u $ decrease as $ \epsilon $ increases. However, varying $ \theta $ has very little to no effect on the equilibrium number of those infected and aware individuals. We also plotted the recovered and treated populations for varying values of $ \theta $ and observed that these sub-populations do not change as we increase $ \theta $. We conclude that increasing public health education campaigns/literacy on how Tungiasis spreads (reducing $ \theta $) plays a significant role in reducing the number of people who eventually end up in class $ I_u $. Ideally, we want to reduce the number of people who end up in this class. If people are educated on Tungiasis dynamics and preventative measures, then they will adhere to all precautionary measures in preventing possible infection. Also, individuals literate on the disease may easily identify symptoms if infected. Then, they can quickly seek medical help and avoid complications brought about by prolonged disease exposure.

    The effect of the treatment rate $ \omega $ on the population dynamics is presented in Figure 7, 8 and 9. Figures 7 and 8 show plots of $ I_a $, $ T $ and $ R $ for $ \omega = 0.1 $ and $ \omega = 0.3 $, respectively, and different fractional order values. We observe that for a fixed $ \omega $ and different fractional order, the equilibrium values of $ I_a $, $ T $ and $ R $ decrease as the fractional order increases from 0.57 to 1. In Figure 9, we also observe that for a given fractional order, increasing $ \omega $ results in a decrease in equilibrium values of the infected class. Equilibrium values of $ T $ and $ R $ increase on average due to the increase in $ \omega $. This is consistent with what is expected of treatment. As more people get treatment and recover from disease, the quality of life of the people improves such that the effects of the disease burden are decreased. If more infected people recover quickly, these people resume their duties (jobs) sooner that otherwise would be the case in the absence of treatment. Treatment negates any potential adverse economic effects due to prolonged sickness. Hence, research into effective drugs for the treatment of tungiasis should be heightened.

    Figure 7.  Numerical results showing the effect of $ \omega = 0.1 $ on the $ I_a $, $ T $ and $ R $ fractional dynamics for different fractional order $ \alpha $.
    Figure 8.  Numerical results showing the effect of $ \omega = 0.3 $ on the $ I_a $, $ T $ and $ R $ fractional dynamics for different fractional order $ \alpha $.
    Figure 9.  Results showing the effect of $ \omega $ on the $ I_a $, $ T $ and $ R $ fractional dynamics for a fixed fractional order $ \alpha $.

    The effect of the contact or transmission rate $ \beta $ is analysed in Figure 10, 11 and 12. Figures 10 and 11 show dynamics of $ I_u $, $ I_a $, $ T $ and $ R $ for $ \beta = 0.1 $ and $ \beta = 0.3 $ with varying fractional order, respectively. Similarly, the equilibrium values of $ I_u $, $ I_a $, $ T $ and $ R $ decrease as the fractional order increases from 0.75 to 1. Qualitatively, the trends of the curves are the same for different fractional values. Figure 12 indicate plots of $ I_u $, $ I_a $, $ T $ and $ R $ for different $ \beta $ for a fixed fractional order. The results show that decreasing $ \beta $ results in a decrease in equilibrium values of $ I_u $. For $ I_a $, we observe a similar trend even though there are times when a lower $ \beta $ results in a higher $ I_a $. For instance, for approximately $ 32.0565 < t < 95.9409 $, $ \beta = 0.2 $ results in higher equilibrium values of $ I_a $ than the case for $ \beta = 0.3 $. However, in the long run, we observe that as $ \beta $ decreases, $ I_a $ decreases, similarly to the $ I_u $ case. If a decrease in $ \beta $ results in a decrease in both $ I_u $ and $ I_a $, then we expect a decrease in both $ T $ and $ R $. This is because $ I_u $ and $ I_a $ feed into $ T $ and subsequently, $ R $. Figure 12 depicts this expected outcome, that is, as $ \beta $ decreases, the equilibrium values of $ T $ and $ R $ decrease too. As a public health intervention, officials should devise strategies that target reducing the contact or transmission rate. Things like wearing shoes, cleaning/disinfection ground surfaces where people walk barefoot, and reducing contact between human beings and stray animals like rats, bovines, pigs, wild cats and dogs, etc can reduce the transmission or spread of tungiasis.

    Figure 10.  Numerical results showing the effect of $ \beta = 0.1 $ on $ I_u $, $ I_a $, $ T $ and $ R $ fractional dynamics for different fractional order $ \alpha $.
    Figure 11.  Numerical results showing the effect of $ \beta = 0.3 $ on $ I_u $, $ I_a $, $ T $ and $ R $ fractional dynamics for different fractional order $ \alpha $.
    Figure 12.  Results showing the effect of $ \beta $ on $ I_u $, $ I_a $, $ T $ and $ R $ fractional dynamics for a fixed fractional order $ \alpha $.

    In this manuscript, we developed and analysed a fractional model for the dynamics of tungiasis, with the aim of establishing the potential impact of public health education and treatment in reducing disease burden. Using the Caputo fractional derivative, we first established if the developed model is well-posed. The boundedness and positivity of the solution of the model are determined using the generalised mean-value theorem. We also determined the model steady-states and presented their stability based on the basic reproduction number. To solve the resulting system of fractional differential equations, we used the Adam-Bashford-Moulton method. The numerical results are presented graphically for different fractional order $ \alpha $. We also studied the impact of public health education, treatment and the contact rate on the overall disease dynamics. This is an attempt to identify key parameters that public health officials can target in their attempts to reduce disease burden in communities. Our numerical results envisage that an increase in media campaigns will result in an increase in the number of educated, thus reducing the number of infected but unaware individuals. If people are infected and aware of their disease status, then they can quickly seek medical help and avoid complications brought about by prolonged disease exposure. Also, programs such as disinfection of public floors in endemic areas, provision of shoes for the disadvantaged, and reduced contact with stray animals, targeted at reducing disease transmission should be prioritised by public health officials. Research into effective tungiasis treatment drugs should also be heightened to reduce disease burden and suffering.

    We acknowledged that the work presented in this paper has some limitations. The fact that the model was not fitted to any epidemiological data in a particular setting is an issue we still need to address. Our results agree with that found by other mathematical models that have previously analysed tungiasis dynamics in sub-Saharan African settings, see, for instance, Nyang'inja et al. [11].

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this manuscript.

    We thank the reviewers for carefully checking and reading our manuscript. The manuscript has improved immensely because of their insightful comments.

    [1] Henderson K (2006) History and myth: the origin of Colonel Light's “Park Grounds” and their universal significance, In: Garnaut C, and Round K, (eds.), The Adelaide Parklands Symposium. A balancing act: past-present-future. Adelaide: Centre for Settlement Studies and the Bob Hawke Prime Ministerial Centre, University of South Australia, 73-89.
    [2] Johnson DL (2006) Foundations for Adelaide's Park Lands, In: Garnaut C, and Round K, (eds.), The Adelaide Parklands Symposium. A balancing act: past-present-future, Adelaide: Centre for Settlement Studies and the Bob Hawke Prime Ministerial Centre, University of South Australia, 106-114.
    [3] Johnson DL (2012) Anticipating municipal parks: London to Adelaide to Garden City, Kent Town SA: Wakefield Press.
    [4] Goldsworthy K (2011) Adelaide, Sydney: NewSouth Publishing, 14.
    [5] Bunker R (1988) Systematic colonization and town planning in Australia and New Zealand. Plann Perspect 3: 59-80. doi: 10.1080/02665438808725652
    [6] Williams M (1966) The parkland towns of Australia and New Zealand. Geogr Rev 56: 67-89. doi: 10.2307/212735
    [7] Budge T, Butt A, Chersterfield M, et al. (2012) Does Australia need a national policy to preserve agricultural land? Surry Hills, NSW: Australian Farm Institute.
    [8] Freestone R (2010) Urban nation: Australia's planning heritage, Clayton, Vic.: CSIRO Publishing.
    [9] Bunker R, Houston P (2003) Prospects for the rural-urban fringe in Australia: observations from a brief history of the landscapes around Sydney and Adelaide. Aust Geogr Stud 41: 303-323. doi: 10.1046/j.1467-8470.2003.00236.x
    [10] Davison A (2006) Stuck in a cul-de-sac: suburban history and urban sustainability in Australia. Urban Policy Res 24: 201-216. doi: 10.1080/08111140600704137
    [11] McLoughlin B (1991) Urban consolidation and urban sprawl: a question of density. Urban Policy Res 9: 148-156. doi: 10.1080/08111149108551499
    [12] Jones D (2000) A gardenesque vision for Adelaide: planning and planting the Adelaide Park Lands in the 1880s, In: Garnaut C, and Hamnett S, (eds.), Proceedings of the Fifth Australian urban history/planning history conference, Adelaide: University of South Australia, 233-244.
    [13] Jones D (2010) Validating planning heritage: Evidence and reasoning for National Heritage recognition of the “City of Adelaide Plan” and Park Lands. T Roy Soc SAust 134: 177-197.
    [14] Hay I, Hughes A, Tutton M (2004) Monuments, memory and marginalisation in Adelaide's Prince Henry Gardens. Geogr Ann B 86: 201-216. doi: 10.1111/j.0435-3684.2004.00162.x
    [15] Hamilton M, Parkland poachers mar vision for city. The Advertiser, April 13, 2010. Available from: http://www.gdp-research.com.au/parks_1.htm
    [16] Sumerling P (2011) The Adelaide Park Lands: a social history, Adelaide: Wakefield Press.
    [17] Planning SA (2001) The future of open space in Adelaide: Parklands 21 strategy, Adelaide: Government of South Australia, 17.
    [18] Mulgavin D (2004) Adelaide's greenway: River Torrens linear park. Landscape Urban Plan 68: 223-240. doi: 10.1016/j.landurbplan.2003.07.002
    [19] Hassell and Partners (1979) River Torrens study: a co-ordinated development scheme, Adelaide: Government Printer.
    [20] Water Resources Branch, Engineering and Water Supply Department South Australia (1980) River Torrens flood mitigation study: flood frequency and hydrograph determinations, Adelaide: EWSDSA.
    [21] Howard E (1902) Garden Cities of tomorrow, London: S. Sonnenschein & Co., Ltd.
    [22] Jones D, Adelaide Park Lands and squares cultural landscape assessment study, prepared for the Corporation of the City of Adelaide, 2007. Available from: http://www.adelaidecitycouncil.com/city-living/sustainable-adelaide/adelaide-park-lands/cultural-landscape-study
    [23] Bell P, Parklands, 2015. Available from: http://adelaidia.sa.gov.au/places/parklands
    [24] Daly J (1987) Decisions and disasters. Alienation of the Adelaide Parklands, Adelaide: Bland House, 137-164.
    [25] Adelaide Parklands Authority (APA) (2010) Adelaide Parklands Management Strategy: Towards 2020, Adelaide: Adelaide City Council and State Government of South Australia.
    [26] Van Vliet N (1991) Adelaide's parklands: the legacy of enlightened planning. Landscape Architect Rev 12: 8-12.
    [27] Chartingtransport.com, Comparing the residential densities of Australian cities 2011. Available from: http://chartingtransport.com/2012/10/19/comparing-the-residential-densities-of-australian-cities-2011/
    [28] Marsden S (2004) Twentieth century heritage survey: Stage 1: Post-Second World War (1946-1959). Overview history, Adelaide: Marsden Russell Historians.
    [29] Bunker R (1990) Urban design in a metropolitan setting: a case study of Adelaide. Town Plann Rev 61: 21-40. doi: 10.3828/tpr.61.1.6v4n661p02575k17
    [30] Department of Premier and Cabinet South Australia, Report on the metropolitan area of Adelaide 2012. Available from: https://www.sa.gov.au/topics/housing-property-and-land/building-and-development/land-supply-and-planning-system/historical-planning-strategies/report-on-the-metropolitan-area-of-adelaide-1962
    [31] Australian Bureau of Statistics (ABS), Australian historical population statistics - 3105.0.65.001, Canberra: Government of Australia, 2008. Available from: http://www.abs.gov.au/AUSSTATS/abs@.nsf/DetailsPage/3105.0.65.0012008?OpenDocument
    [32] Hutchings A (1988) Visions for Adelaide: the persistence of planning ideas. Urban Policy Res 6: 136-139. doi: 10.1080/08111148808551335
    [33] Hutchings A (1993) Consistent vision: the planning of Metropolitan Adelaide 1962-1992, In: Freestone R, (ed.), Spirited cities: urban planning, traffic and environmental management in the Nineties: essays for Hans Westerman, Sydney: Federation Press, 190-208.
    [34] Hutchings A (1998), Planning history in practice: A commentary from South Australia. Australian Planner 35: 122-126.
    [35] Department of Environment and Planning South Australia (DEPSA) (1992) South Australian Planning Review - 2020 Vision: Planning Strategy for Metropolitan Adelaide, Adelaide: DEPSA.
    [36] Hamnett S (2005) Ten years of metropolitan strategic planning in South Australia. Proceedings of the second state of Australian cities conference, city structure section, Brisbane: Urban Research Program, Griffith University.
    [37] Hamnett S, Parham S (1992) 2020 Vision: a planning strategy for metropolitan Adelaide. Urban Futures 2: 85-88.
    [38] Llewellyn-Smith M (2012) Behind the scenes: The politics of planning Adelaide, Adelaide: University of Adelaide Press, 309-330.
    [39] Premier of South Australia (2003) Planning Strategy for the Development of Metropolitan Adelaide, Adelaide: Department of Premier and Cabinet SA.
    [40] Department of Planning and Transport Infrastructure South Australia (DPTISA) (2010) The 30-Year Plan for Greater Adelaide, Adelaide: Government of South Australia.
    [41] Gleeson B (2006) Australian heartlands: making space for hope in the suburbs. Crows Nest, NSW: Allen & Unwin.
    [42] Stretton H (1970) Ideas for Australian cities. Vermont South, Vic: ARRB Group Limited.
    [43] Stretton H (2005) Density, efficiency and equality in Australian cities. In: E Burton, M Jenks and K Williams (eds.) The compact city: A sustainable urban form, 37-43. London: Taylor & Francis e-Library.
    [44] Low N, Gleeson B, Green R, et al. (2005) The green city: sustainable homes, sustainable suburbs. Sydney: NSW Press.
    [45] Troy PN (1996) The perils of urban consolidation: A discussion of Australian housing and urban development policies. Leichhardt, NSW: Federation Press.
    [46] Troy PN (2003) Saving our cities with suburbs. Griffith Rev 2: 80-90.
    [47] Troy PN (2004) The structure and form of the Australian city: Prospects for improved urban planning. Urban Policy Program Issues Paper 1. Brisbane: Griffith University
    [48] Buxton M, Tieman G (2005) Patterns of urban consolidation in Melbourne: Planning policy and the growth of medium density housing. Urban Policy Res 23: 137-157. doi: 10.1080/08111470500135086
    [49] Gleeson B, Dodson J, Spiller M (2010) Metropolitan governance for the Australian city: the case for reform. Urban Policy Program Issues Paper 12. Brisbane: Griffith University.
    [50] Curtis C (2006) Network city: Retrofitting the Perth Metropolitan Region to facilitate sustainable travel. Urban Policy Res 24: 159-180. doi: 10.1080/08111140600703691
    [51] Mees P (2009) How dense are we? Another look at urban density and transport patterns in Australia, Canada and the USA. Road Transport Res 18: 58-67.
    [52] Baum S, Gleeson B (2010) Space and place: Social exclusion in Australia's suburban heartlands. Urban Policy Res 28: 135-159. doi: 10.1080/08111141003693133
    [53] Bostock L, Gleeson B (2004) Contested housing landscapes? Social inclusion, deinstitutionalison and housing policy in Australia. Aust J Soc Issues 39: 41-62.
    [54] Troy PN (2012) Accommodating Australians: Commonwealth Government involvement in housing. Leichhardt, NSW: Federation Press.
    [55] Beer A, Kearins B, Pieters H (2007) Housing affordability and planning in Australia: The challenge of policy under neo-liberalism. Housing Stud 22: 11-24. doi: 10.1080/02673030601024572
    [56] McGuirk P (2005) Neoliberalist planning? Re-thinking and re-casting Sydney's metropolitan planning. Geogr Res 43: 59-70.
    [57] Public Transport Users Association (PTUA) (2007) Moving Australians sustainably: Transport policy in the national interest. Melbourne: PTUA.
    [58] Government of South Australia, Building a stronger South Australia: 3 Public transport, 2013. Available from: http://stronger.sa.gov.au/wp-content/uploads/2013/11/building_a_stronger_sa-public_transport.pdf
    [59] Gleeson B (2008) Critical commentary. Waking from the dream: an Australian perspective on urban resilience. Urban Stud 45: 2653-2668.
    [60] SGS Economics & Planning, The adequacy of greenfield land supply in metro Adelaide. Final report for DPTISA, 2013. Available from: http://dpti.sa.gov.au/__data/assets/pdf_file/0016/123208/20130361_Final_Report_130930.pdf
    [61] Gardner B, SA Government aims to halt urban sprawl. Australian Property Institute News, 2014. Available from: http://news.api.org.au/article/537936/sa_government_aims_halt_urban_sprawl/
    [62] Adelaide City Council (ACC) (2012) The City of Adelaide Strategic Plan 2012-16, Adelaide: ACC.
    [63] Low Carbon Living Co-operative Research Centre (LCLCRC), RP3017: Adelaide Living Laboratory Hub, 2015. Available from: http://www.lowcarbonlivingcrc.com.au/research/program-3-engaged-communities/rp3017-adelaide-living-laboratory-hub
    [64] LCLCRC, Program 2: Low Carbon Precincts, 2015. Available from: http://www.lowcarbonlivingcrc.com.au/research/program-2-low-carbon-precincts
    [65] Government of South Australia, Tonsley, 2015. Available from: http://www.tonsley.com/project_vision
    [66] Garnaut R, Lochiel Park: case study of a green village - South Australia. Garnaut Climate Change Review, 2011. Available from: http://www.garnautreview.org.au/CA25734E0016A131/WebObj/Casestudy-GreenVillage-LochielPark-SouthAustralia/$File/Case%20study%20-%20Green%20Village%20-%20
    [67] Edwards J, Pocock B (2011) Comfort, convenience and cost: The calculus of sustainable living at Lochiel Park, Adelaide: Centre for Work + Life, University of South Australia.
    [68] Barr S, Gilg A (2007) A conceptual framework for understanding and analysing attitudes toward environmental behaviour. Geogr Ann B 89: 361-279. doi: 10.1111/j.1468-0467.2007.00266.x
    [69] Shove E (2003) Users, technologies and expectations of comfort, cleanliness and convenience. Innovation 16: 193-206.
    [70] DPTISA, Bowden urban village and environs Development Plan Amendment, 2012. Available from: http://www.sa.gov.au/__data/assets/pdf_file/0014/22622/Bowden_A3_comm_info_template_DPTI_singles.pdf
    [71] Renewal SA, Bowden: Life more interesting. Redefining city living, 2015. Available from: http://lifemoreinteresting.com.au/vision/
    [72] DPTISA, Integrated transport and land use plan, 2013. Available from: http://www.transportplan.sa.gov.au/
    [73] Renewal SA, Delivering an inspiring urban future: Annual Report 2013-14. Available from: https://renewalsa.sa.gov.au/renewalsa/cms/wp-content/uploads/2014/05/Delivering-an-inspiring-urban-future-Annual-Report-2013-14-30-September-2014-with-Financial-Statements_1.pdf
    [74] Hopgood D (1986) Metropolitan Adelaide and its future, In: Hutchings A, and Bunker R, (eds.), With conscious purpose. A history of town planning in South Australia, Adelaide: Wakefield Press and the Royal Australian Planning Institute South Australian Division.
    [75] Garnaut C (2008) The Adelaide Parklands and the endurance of the green belt idea in South Australia, In: Amati M. (ed.), Urban Green Belts in the Twenty-first Century, Aldershot UK and Burlington VT: Ashgate Publishing, 107-128.
    [76] DPCSA, Metropolitan Open Space System (MOSS), 2013. https://www.sa.gov.au/__data/assets/pdf_file/0015/11274/MOSS_Factsheet_Aug_2013.pdf
    [77] Planning SA (2001) The future of open space in Adelaide: Parklands 21 strategy, Adelaide: Planning SA, 11.
    [78] Mahmoudi P, Hatton MacDonald D, Crossman ND, et al. (2013) Space matters: the importance of amenity in planning metropolitan growth. Aust J Agr Resour Ec 57: 38-59. doi: 10.1111/j.1467-8489.2012.00608.x
    [79] Planning SA (2004) Coast Park 2036. South Aust Geogr 18: 7-8.
    [80] King CA (2008) Community resilience and contemporary agri‐ecological systems: reconnecting people and food, and people with people. Syst Res and Behav Sci 25: 111-124. doi: 10.1002/sres.854
    [81] Giradet H (2005) Urban agriculture and sustainable urban development, In: Viljoen A, (ed.) CPLUs: continuous productive urban landscapes - Designing urban agriculture for sustainable cities, Oxford and Burlington, MA: Architectural Press, 32-39.
    [82] Lea E, Phillips J, Ward M, et al. (2006) Farmers' and consumers' beliefs about Community-Supported Agriculture in Australia: A qualitative study. Ecol Food Nutr 45: 61-86. doi: 10.1080/03670240500530592
    [83] McInerney A, Community Supported Agriculture (CSA) alive and well in Adelaide - Jupiter Creek Farm, 2015. Available from: http://www.lambsearsandhoney.com/2012/03/community-supported-agriculture-csa-alive-well-in-adelaide-jupiter-creek-farm/
    [84] Beach Organics, Beach Organics Community Supported Agriculture system (CSA), 2015. Available from: http://www.beachorganics.com.au/csa.html
    [85] Guitart D, Pickering C, Byrne J (2012) Past results and future directions in urban community gardens research. Urban For Urban Gree 11: 364-373. doi: 10.1016/j.ufug.2012.06.007
    [86] Kingsley JY, Townsend M (2006) Dig in-to social capital: community gardens as mechanisms for growing urban social connectedness. Urban Policy Res 24: 525-537. doi: 10.1080/08111140601035200
    [87] Kitchen Gardens SA, Community gardens, 2015. Available from: http://www.kitchengardenssa.com.au/index.php/running-your-garden
    [88] Markow K, Coveney J, Booth S (2014) Improving access to community-based food systems in Adelaide, South Australia: Strategies to encourage low-socio-economic status groups to participate. J Hunger Environ Nutr 9: 113-134. doi: 10.1080/19320248.2013.840550
    [89] Winter M (2003) Geographies of food: agro-food geographies - making reconnections. Prog Hum Geog 27: 505-513. doi: 10.1191/0309132503ph446pr
    [90] Gaynor A (2006) Harvest of the suburbs: An environmental history of growing food in Australian cities, Crawley, WA: University of Western Australia Press.
    [91] Robinson GM (2009) Towards sustainable agriculture: current debates. Geogr Compass 3: 1757-1773. doi: 10.1111/j.1749-8198.2009.00268.x
    [92] Goodman D, DuPuis EM (2002) Knowing food and growing food: beyond the production-consumption debate in the sociology of agriculture. Sociol Ruralis 42: 5-22. doi: 10.1111/1467-9523.00199
    [93] Renting H, Marsden TK, Banks J (2003) Understanding alternative food networks: exploring the role of short food supply chains in rural development. Environ Plann A 35: 393-412. doi: 10.1068/a3510
    [94] Jarosz L (2008) The city in the country: Growing alternative food networks in Metropolitan areas. J Rural Stud 24: 231-244. doi: 10.1016/j.jrurstud.2007.10.002
    [95] Fielke S, Bardsley DK (2013) South Australian farmers' markets: tools for enhancing the multifunctionality of Australian agriculture. GeoJournal 78: 759-776. doi: 10.1007/s10708-012-9464-8
    [96] Department of the Environment, Water, Heritage and the Arts (DEWHA) (2010) Rain gardens, green roofs and infiltration systems, In: Department of Planning and Local Government, Water sensitive urban design technical manual for the Greater Adelaide region, Adelaide: Government of South Australia.
    [97] Fifth Creek Studio (2012) Green roof trials monitoring report, Adelaide: South Australian Government, prepared for the SA Government's Building Innovation Fund and Aspen Development Fund No. 1.
    [98] Razzaghmanesh M, Beecham S, Introducing green roofs into the urban and built environments of Adelaide', paper presented to the 13th National StreetTree Symposium, 2012. Available from: http://treenetmedia.com/up/pdf/2012/Introducing%20Green%20Roofs_Razzaghmanesh%20and%20Beechman.pdf
    [99] Hopkins K, Goodwin C (2011) Living architecture: Green roofs and walls, Clayton, Vic.: CSIRO Publishing.
    [100] Mibus R (2006) Green roofs in the changing Australian landscape, paper presented to the International Landcare Conference, Melbourne, 8-11 October.
    [101] Bolan N, Laurenson S, Kunhikrishnan A, et al. Northern Adelaide Plains recycling scheme - champion in the management of recycled water for sustainable production and environmental protection, 2008. Available from: http://irrigation.org.au/documents/publications-resources/conference_papers_2008/200508_Stream3_Bolan1.pdf
    [102] Miller J, Roots J (2012) Changes in Australian agriculture and land use: implications for future food security. Int J Agric Sustain 10: 25-39. doi: 10.1080/14735903.2012.646731
    [103] Thorpe D, Adelaide Parklands to benefit from high quality recycled water. ReWater Summer 2008/9, 2-3. Available from: http://www.recycledwater.com.au/uploads/File/newsletters/reWater%20Summer%2008.pdf
    [104] Hassanli AM, Kazemi F (2012) An investigation into the tolerance and sensitivity of the Adelaide Parklands landscape plants to the Glenelg recycled wastewater. Adelaide: University of South Australia.
    [105] Naughton D, Recycled water damages park lands: report. South Australian Policy Online, 2012. Available from: http://www.sapo.org.au/opin/opin24284.html
    [106] Weatherill J, Rau J, News release: Making the Parklands a place for everyone, 2013. Available from: http://dpti.sa.gov.au/__data/assets/pdf_file/0007/118429/sundayparklands.pdf
    [107] Capital City Committee Adelaide, Major infrastructure projects, 2015. Available from:
    [108] http://capcity.adelaide.sa.gov.au/priorities/major-infrastructure-projects-2/
    [109] 108. Adelaide Hills Council (AHC) (2013) Adelaide Hills Council biodiversity strategy, Adelaide Hills: AHC, 5.
    [110] 109. Department of the Environment, Australian Government, Australia's 15 National Biodiversity Hotspots, Canberra: Australian Government, 2006. Available from: http://www.environment.gov.au/topics/biodiversity/biodiversity-conservation/biodiversity-hotspots/national-biodiversity-hotspots
    [111] 110. Armstrong DM (2003) Mammals, In: Armstrong DM, Croft SJ, and Foulkes JN (eds.), A biological survey of the Southern Mount Lofty Ranges, South Australia 2000-2001, Adelaide: Biodiversity Survey and Monitoring Section Science and Conservation Directorate Department for Environment and Heritage South Australia, 129-156.
    [112] 111. Gill, AM, McKenna DJ, Wouters MA (2014) Landscape fire, biodiversity decline and a rapidly changing milieu: A microcosm of global issues in an Australian biodiversity hotspot. Land 3: 1091-1136. doi: 10.3390/land3031091
    [113] 112. Tait CJ, Daniels CB, Hill RS (2005) Changes in species assemblages within the Adelaide Metropolitan area, Australia, 1836-2002 Ecol Appl 15: 346-359.
    [114] 113. James SW (2014) Protecting Sydney's peri-urban agriculture: Moving beyond a housing/farming dichotomy. Geogr Res 52: 377-386. doi: 10.1111/1745-5871.12048
    [115] 114. Mason D, Knowd I (2010) The emergence of urban agriculture: Sydney, Australia. Int J Agric Sustain 8: 62-71. doi: 10.3763/ijas.2009.0474
    [116] 115. Merson J, Attwater R, Ampt P, et al. (2010) The challenges to urban agriculture in the Sydney basin and lower Blue Mountains region of Australia. Int J Agric Sustain 8: 72-85. doi: 10.3763/ijas.2009.0464
    [117] 116. Butt A (2013a) Functional change and the peri-urban region: Food systems and agricultural vulnerability. Econ Pap 32: 308-316.
    [118] 117. Butt A (2013b) Development, dilution, and functional change in the peri-urban landscape: what does it really mean for agriculture? In: Farmar-Bowers Q, Higgins V, and Millar J, (eds.), Food security in Australia: Challenges and prospects for the future, New York: Springer, 425-441.
    [119] 118. Houston P (2005) Re-valuing the fringe: Some findings on the value of agricultural production in Australia's peri-urban regions. Geogr Res 43: 209-223. doi: 10.1111/j.1745-5871.2005.00314.x
    [120] 119. Carey R, Krumholz F, Duignan K, et al. (2011) Integrating agriculture and food policy to achieve sustainable peri-urban fruit and vegetable production in Victoria, Australia. J Agr Food Syst Commun Dev 1: 181-195.
    [121] 120. Houston P (2003) The national audit of peri-urban agriculture. Aust Planner 40: 43-45. doi: 10.1080/07293682.2003.9995278
    [122] 121. Pires V, Planning for urban agriculture planning in Australian cities, 2011. Available from: http://soac.fbe.unsw.edu.au/2011/papers/SOAC2011_0055_final(1).pdf
    [123] 122. Penman TD, Bradstock RA, Price O (2013) Modelling the determinants of ignition in the Sydney Basin, Australia: implications for future management. Int J Wildland Fire 22: 469-478. doi: 10.1071/WF12027
    [124] 123. Adelaide Hills Council (2014) Development Plan. Adelaide: DPTISA.
    [125] 124. Government of South Australia, Requirements for bushfire risk areas, 2015. Available from:
    [126] https://www.sa.gov.au/topics/housing-property-and-land/building-and-development/residential-building-regulations/fire-safety-and-bushfire-regulations/building-and-development-requirements-for-bushfire-prone-areas/requirements-for-bushfire-risk-areas
    [127] 125. Government of South Australia (2003) South Australian Development Act 2003. Adelaide, Australia: Government of SA.
    [128] 126. Department of Planning and Local Government (DPLGSA) 2011. Protecting Regulated and Significant Trees: Community Information. Adelaide, Australia: Government of SA.
    [129] 127. Mannix L, Adelaide's density barely changed in 30 years, 2013. Available from: http://indaily.com.au/news/2013/11/07/adelaides-density-barely-changed-in-30-years/
  • This article has been cited by:

    1. Gilberto González-Parra, Hana M. Dobrovolny, Editorial: Mathematical foundations in biological modelling and simulation, 2024, 21, 1551-0018, 7084, 10.3934/mbe.2024311
    2. CW Chukwu, S. Y. Tchoumi, Z. Chazuka, M. L. Juga, G. Obaido, Assessing the impact of human behavior towards preventative measures on COVID-19 dynamics for Gauteng, South Africa: a simulation and forecasting approach, 2024, 9, 2473-6988, 10511, 10.3934/math.2024514
    3. Marlene Thielecke, Hannah McNeilly, Francis Mutebi, Mike B. Banalyaki, Rebecca Arono, Susanne Wiese, Felix Reichert, George Mukone, Hermann Feldmeier, High Level of Knowledge about Tungiasis but Little Translation into Control Practices in Karamoja, Northeastern Uganda, 2023, 8, 2414-6366, 425, 10.3390/tropicalmed8090425
  • Reader Comments
  • © 2015 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(11115) PDF downloads(1537) Cited by(10)

Figures and Tables

Figures(4)  /  Tables(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog