Research article Special Issues

Summer cooling potential of urban vegetation—a modeling study for Melbourne, Australia

  • Received: 17 March 2015 Accepted: 11 May 2015 Published: 24 June 2015
  • The summer cooling potential of urban vegetation is investigated using an urban climate model for the current and future climates in the Melbourne central business district (CBD) area with various urban forms and vegetation schemes. Simulation results suggest that the average seasonal summer temperatures can be reduced in the range of around 0.5 and 2°C if the Melbourne CBD were replaced by vegetated suburbs and planted parklands, respectively, benefiting a reduction in the number of hot days. It was also found that despite the projected warming in the future and variations in the climate projections among different climate models, the average seasonal cooling potential due to various urban vegetation schemes may not change significantly in comparison with those predicted for the current climate, indicating little dependency on climate change. This finding suggests that the average seasonal cooling potential as a result of urban vegetation in future climates may be empirically quantified in similar amounts to those under the current climate. When urban climate models are used, the cooling potential of urban vegetation in future climates may be quantified by modeling several selected years with one or a few climate models.

    Citation: Dong Chen, Marcus Thatcher, Xiaoming Wang, Guy Barnett, Anthony Kachenko, Robert Prince. Summer cooling potential of urban vegetation—a modeling study for Melbourne, Australia[J]. AIMS Environmental Science, 2015, 2(3): 648-667. doi: 10.3934/environsci.2015.3.648

    Related Papers:

    [1] Nalin Fonseka, Jerome Goddard Ⅱ, Alketa Henderson, Dustin Nichols, Ratnasingham Shivaji . Modeling effects of matrix heterogeneity on population persistence at the patch-level. Mathematical Biosciences and Engineering, 2022, 19(12): 13675-13709. doi: 10.3934/mbe.2022638
    [2] Nazanin Zaker, Christina A. Cobbold, Frithjof Lutscher . The effect of landscape fragmentation on Turing-pattern formation. Mathematical Biosciences and Engineering, 2022, 19(3): 2506-2537. doi: 10.3934/mbe.2022116
    [3] Nancy Azer, P. van den Driessche . Competition and Dispersal Delays in Patchy Environments. Mathematical Biosciences and Engineering, 2006, 3(2): 283-296. doi: 10.3934/mbe.2006.3.283
    [4] James T. Cronin, Jerome Goddard II, Amila Muthunayake, Ratnasingham Shivaji . Modeling the effects of trait-mediated dispersal on coexistence of mutualists. Mathematical Biosciences and Engineering, 2020, 17(6): 7838-7861. doi: 10.3934/mbe.2020399
    [5] Eric Ruggieri, Sebastian J. Schreiber . The Dynamics of the Schoener-Polis-Holt model of Intra-Guild Predation. Mathematical Biosciences and Engineering, 2005, 2(2): 279-288. doi: 10.3934/mbe.2005.2.279
    [6] Robert Stephen Cantrell, Chris Cosner, William F. Fagan . The implications of model formulation when transitioning from spatial to landscape ecology. Mathematical Biosciences and Engineering, 2012, 9(1): 27-60. doi: 10.3934/mbe.2012.9.27
    [7] Chang-Yuan Cheng, Kuang-Hui Lin, Chih-Wen Shih . Coexistence and extinction for two competing species in patchy environments. Mathematical Biosciences and Engineering, 2019, 16(2): 909-946. doi: 10.3934/mbe.2019043
    [8] Kuang-Hui Lin, Yuan Lou, Chih-Wen Shih, Tze-Hung Tsai . Global dynamics for two-species competition in patchy environment. Mathematical Biosciences and Engineering, 2014, 11(4): 947-970. doi: 10.3934/mbe.2014.11.947
    [9] Jesus R. Artalejo, A. Economou, M.J. Lopez-Herrero . Evaluating growth measures in an immigration process subject to binomial and geometric catastrophes. Mathematical Biosciences and Engineering, 2007, 4(4): 573-594. doi: 10.3934/mbe.2007.4.573
    [10] Jacques A. L. Silva, Flávia T. Giordani . Density-dependent dispersal in multiple species metapopulations. Mathematical Biosciences and Engineering, 2008, 5(4): 843-857. doi: 10.3934/mbe.2008.5.843
  • The summer cooling potential of urban vegetation is investigated using an urban climate model for the current and future climates in the Melbourne central business district (CBD) area with various urban forms and vegetation schemes. Simulation results suggest that the average seasonal summer temperatures can be reduced in the range of around 0.5 and 2°C if the Melbourne CBD were replaced by vegetated suburbs and planted parklands, respectively, benefiting a reduction in the number of hot days. It was also found that despite the projected warming in the future and variations in the climate projections among different climate models, the average seasonal cooling potential due to various urban vegetation schemes may not change significantly in comparison with those predicted for the current climate, indicating little dependency on climate change. This finding suggests that the average seasonal cooling potential as a result of urban vegetation in future climates may be empirically quantified in similar amounts to those under the current climate. When urban climate models are used, the cooling potential of urban vegetation in future climates may be quantified by modeling several selected years with one or a few climate models.


    There has been significant progress in conducting research in the photovoltaic (PV) field over the last decade [1,2,3,4]. PV energy is generated by the conversion of sunlight and temperature. The particular interest in PV energy has driven researchers to enhance and optimize the efficiency of energy conversion within PV systems [5,6,7,8]. The autonomous functioning of a PV system in optimal conditions requires the optimization of several physical parameters related to the PV generator used [9,10,11]. Therefore, introducing an electronic converter between the PV generator and the load is essential. A step-up converter is suitable to achieve the goal of the studied work [12,13,14]. Indeed, as defined by the duty cycle, the factor characterizing the step-up device should be controlled by a specific technique of the PV power maximization [15,16].

    In the literature, the numerous documented power maximization techniques are classified into three categories namely: Classical iterative methods, revolutionary methods based on artificial intelligence and hybrid methods [17,18,19,20]. The first class is mainly based on processing instantaneous values of power, currents and voltages, such as Perturb & Observe (P & O) and Incremental Conductance (INC). The second class includes techniques based on artificial intelligence, such as Artificial Neural Networks (ANN) and Fuzzy Logic (FL). On the other hand, the third class includes the techniques defined by the combination of two or more maximum power point tracker (MPPT) methods, namely the following: The hybrid two-stage adaptive method, the novel Lyapunov-based rapid and ripple-free method, and the adjustable variable step based MRAC method. In [21], the P & O method was observed with an oscillatory convergence around the maximum power point (MPP), thus causing instability and a margin of error in the found power value. These oscillations are explained by the difference between two instantaneous power values. Loukriz et al. [22] presented the INC method, known by its damped oscillation due to the small convergence step generated by the ratio of the instantaneous variation of current to the instantaneous variation of voltage. On the other hand, Goel et al. [23] displayed the ability of the ANN tool to rapidly and intelligently predict the value of the duty cycle to control the used converter. Moreover, [24] describes the converter control strategy as a set of linguistic rules that allow the value to reach the desired MPP quickly with a better performance. The hybrid, two-stage, adaptive MPPT method proposed by [25] is formed by two stages: The first stage presents a control block to find the reference voltage for each MPP; and the second stage is an adaptive model reference controller block that determines concise values of the duty cycle to hold on stable the found MPP. This strategy shows high performances under varying irradiances and temperatures. Manna et al. [26] suggested a Lyapunov-based robust model reference adaptive controller to quickly find the MPP for rapid variations in irradiance, temperature and output load. Singh et al. [27] experimented with an efficient technique for different weather scenarios defined by an adjustable P & O variable step.

    In a detailed study, we separately applied P & O and ANN as an MPPT method to control a boost converter and closely observed their efficiencies and performances. We found that P & O is accurate, unstable, and sensitive to initial values while ANN is less accurate, stable, and did not require initial values. We combined these two techniques with the aim to start the convergence by ANN trained by the Levenberg-Marquardt (LM) method to approach the optimum of the duty cycle and then hand it over to P & O to continue the convergence. Therefore, the obtained method is sensitive to the irradiance changes considerably affecting the precision. These results can be explained by the high sensitivity of P & O toward the changing climatic conditions [28].

    While trying to take advantage of the intelligent prediction quality of ANN and its low sensitivity towards initial values, we thought of improving its training process to enhance its low precision. In principle, the LM is the most used method to ensure the learning of ANNs. In the literature, the convergence of gradient descent methods, such as LM and Gauss-Newton (GN) tends towards a local minimum knowing that the shape of the error function to be optimized can have either a single global minimum or several local minimums [29]. This can explain the low precision observed when searching for MPP by an ANN trained by LM. On the other hand, in [30], we proposed an LM combined with simulated annealing (SA) for the parameter identification of the single diode model of the solar cell and we obtained the most accurate results compared to those reported in the literature. The interest behind the use of SA is interpreted by the existence of several local minima in the shape of the error function according to the damping factor, thus characterizing the equation of the LM method. In the present paper, for the first time, we proposed to combine the conventional GN method with the heuristic technique SA (GNSA) to guarantee global convergence during the optimization process treated in the training of ANN. Indeed, GNSA is proposed to ensure the training of the ANN network by adjusting the characteristic weights and bias wij, wjm, bij and bjm to predict the correct value of the duty cycle corresponding to a defined value of temperature and irradiance.

    The three techniques (i.e., P & O, ANN and ANN-P & O) were studied to establish a comparative study and highlight the effectiveness of the proposed ANN trained by the GNSA (ANN-GNSA) approach. After an in-depth study under variable irradiance and temperature, it is observed that the proposed ANN-GNSA method demonstrates a significant efficiency, a good accuracy during the search for the MPP and a high robustness towards the meteorological condition's variation.

    The main contributions of this work are presented as follows:

    √ For the first time, this paper proposed the novel ANN-GNSA approach to precisely track the right MPP into a PV system based on a step-up converter with a high speed of convergence under variable irradiance and temperature values.

    √ The idea behind the proposition of the improved GNSA to train the ANN network is explained by its descent directions quality guided by the globally convergent SA method to adjust correctly the weights and the bias of the hidden and the output layer of ANN tool at each iteration.

    √ A fair comparison is carried out between the outcomes of the proposed method and those obtained by the P & O, ANN and ANN-P & O approaches to demonstrate the high performances of ANN-GNSA in terms of its outstanding accuracy and the speed of convergence toward the desired MPP.

    • First, the ANN-GNSA is studied at fixed meteorological conditions (G = 1000 W/m2 and T = 25℃), where it reaches the right MPP = 18.59 W during a convergence time of 0.04832 s.

    • The proposed method controls the step-up converter under variable irradiances G = [200 W/m2, 300 W/m2, 700 W/m2, 1000 W/m2, 800 W/m2, 400 W/m2], with an efficiency of 99.54% and under variable temperatures T = [15℃, 35℃, 45℃, 5℃] with an efficiency of 99.98%. The tracking of MPP in the both cases is ensured with a high speed compared to P & O, ANN and ANN-P & O.

    The functioning within a direct connection between a PV generator (GPV) and a load is carried out according to a non-optimal power. Certainly, the power transmitted by a GPV to the output might not align with the maximum power [31].

    Figure 1.  Illustration of a photovoltaic system.

    For this purpose, it is essential to utilize an impedance adaptation stage in direct-direct mode (DC-DC) in order to optimize the power at the output of the GPV. In this context, a step-up (Boost) chopper-type static electronic converter presents the most appropriate device for this work. A suitable MPPT settling the duty cycle of the boost device to its optimal values presents the appropriate solution to maximize the power and enable the optimal operation of the PV system.

    A PV module is a collection of solar cells generating continuous electrical energy. The mode of connection of the cells, either in series or in parallel, is carried out according to the energy requirement of the PV voltage or the PV current [32].

    The equivalent circuit of a PV module based on a single diode is simply the same circuit of the PV cell but scaled by Ns. The latter presents the count of solar cells linked in series. The multiplication by Ns is applied for series resistance (Rs), parallel resistance (Rsh) and the diode's ideality factor (n) [33].

    Figure 2.  The electrical model of PV module with one diode.

    Iph is the photocurrent, Is denotes the diode's saturation current, IPV indicates the PV current and VPV presents the PV voltage.

    The mathematical equation deduced from the electrical characterization of the PV module, formulating Rs, Rsh, Iph, Is and n parameters as a function of IPV and VPV is presented by Eq 1 [34].

    $ {\mathrm{I}}_{\mathrm{P}\mathrm{V}} = {\mathrm{I}}_{\mathrm{p}\mathrm{h}}-{\mathrm{I}}_{\mathrm{s}}\left[\mathrm{exp}\left(\frac{\frac{{\mathrm{V}}_{\mathrm{P}\mathrm{V}}}{{\mathrm{N}}_{\mathrm{s}}}+{\mathrm{R}}_{\mathrm{s}}{\mathrm{I}}_{\mathrm{P}\mathrm{V}}}{\mathrm{n}{\mathrm{V}}_{\mathrm{t}\mathrm{h}}}\right)-1\right]-\frac{\frac{{\mathrm{V}}_{\mathrm{P}\mathrm{V}}}{{\mathrm{N}}_{\mathrm{s}}}+{\mathrm{R}}_{\mathrm{s}}{\mathrm{I}}_{\mathrm{P}\mathrm{V}}}{{\mathrm{R}}_{\mathrm{s}\mathrm{h}}} $ (1)

    where Vth is the thermal voltage calculated using (A. T)/q, A presents the Boltzmann constant equal to 1.3806503 × 10-23 J/°K, q denotes the charge of the electron equal to 1.60217646 × 10-19 C and T is the temperature of the cell measured in Kelvin.

    The expression of the five intrinsic parameters of the PV module in terms of climatic factors is crucial to establish a faithful representative model. Eqs (2) and (3) allow the modelling of a GPV, and show the relation with the irradiance (G) and the temperature (T) [35,36]:

    $ {I}_{PV} = \left({I}_{PV, n}+{K}_{i}\mathrm{\Delta }T\right)\frac{G}{{G}_{n}} $ (2)
    $ {I}_{s} = \frac{{I}_{sc, n}+{K}_{i}\mathrm{\Delta }T}{exp\left(\frac{q\left({V}_{oc, n}+{K}_{i}\mathrm{\Delta }T\right)}{n{V}_{th}}\right)-1} $ (3)

    where ΔT denotes the variation between the given T and Tn which is the nominal value of T, Voc, n is the nominal value of the open circuit voltage, Gn presents the nominal value of the irradiation, IPV, n indicates the nominal value of the photovoltaic current, Ki is the coefficient of the current and Isc, n presents the short circuit current of the solar cell.

    In this study, a boost type static electronic converter was used to increase the energy coming from the module intended to supply the output load. Figure 3 illustrates an electrical circuit of the boost converter based on an inductance L, an input capacitance Cint, an output capacitance Cout, a diode D, and a Mosfet K controlled by the duty cycle at its trigger [37].

    Figure 3.  The electrical circuit of the boost converter.

    VPV is the generated voltage by the PV module, Vs is the voltage at the output load and ipv is the input current.

    The duty cycle D is calculated as follows:

    $ D = \frac{{{V_D} - {V_{PV}} + {V_s}}}{{{V_D} + {V_s} - {V_{sw}}}} $ (4)

    where VD is the diode voltage, Vs is the output voltage at the load, Vsw is the Mosfet voltage and VPV is the voltage at the output of the PV module.

    The inductance L is found by Eq 5 using ripple current inductors ΔiL [38]:

    $ L = \frac{{{V_{PV}}\left( {{V_s} - {V_{PV}}} \right)}}{{\Delta {i_L} \times {f_s} \times {V_s}}} $ (5)

    where fs is the switching frequency of the step-up device.

    The value of the input capacitance can be found by the following equation [38]:

    $ C = \frac{{{I_{s(\max )}} \times D}}{{{f_s} \times \Delta {V_s}}} $ (6)

    To guide the module towards its operating MPP, the power optimization process requires a specific mechanism called the "maximum power point tracker" (MPPT) settling an adaptation stage [39,40,41].

    Figure 4.  Position of MPP in the IPV(VPV) curve.

    To achieve this, the implementation of a specific MPPT method is crucial. The most prevalent power maximization methods reported in the literature include the following:

    • The conventional perturbation & observation (P & O) method,

    • The classical artificial neural network (ANN) method,

    • The ANN in conjunction with the P & O method.

    These methods have been studied to underscore the performance of the novel ANN-GNSA methods.

    As its name indicates, the P & O method applies perturbations at the VPV voltage and observes the sign of the instantaneous power variation in order to decide the convergence direction [21].

    Figure 5.  Presentation of P & O algorithm.

    The perturbation is ensured by applying either positive or negative adjustments, depending on the convergence orientation towards the MPP point. This process first begins by incrementing the voltage by adding a positive step. This increment addition operation is adopted if the adjusted value of the generated power is always positive. Once this variation turns negative, the voltage is either decreased or increased by either a negative step or by a positive step, respectively, depending on the position of the functioning point, placed either on the left or on the right of the MPP, respectively.

    Due to its black box characteristic, this tool allows us to model problems without any prior knowledge of the physical or theoretical functioning of a system. The ANN helps to establish a significant liaison between its inputs and outputs. This connection is ensured by updating the characteristic parameters, named "weights" (wij and wjm) and "bias" (bij, bjm), thus marking its internal structure formed by layers of neurons [23].

    The neural network operates as a MPPT command to automatically generate the exact duty cycle according to the temperature (T) and the irradiance (G) values given at the input. The effectiveness of an ANN highly depends on the success of the learning process as well as on the optimal number of hidden neurons, which was fixed at six neurons in this study. The learning process is generally ensured by the LM method whose training principle is based on the backpropagation technique. Typically, the inputs of the ANN network are T and G, and the output is simply represented by the duty cycle D, which is intended to control the activation of the chopper Mosfet. The proposed ANN network aimed to find the PPM point, which is illustrated as follows:

    Figure 6.  The suggested ANN architecture.

    where i refers to the index of the output layer, j pertains to the index of the hidden layer, m denotes the index of the input layer, wij signifies the weight linking the hidden neurons to the output, wjm represents the weight linking the input neurons to the hidden ones, bij indicates the bias of the output neurons and bjm denotes the bias of the hidden neurons.

    The generation of the exact value of the duty cycle strongly depends on the learning accuracy of the network and the optimal number of the hidden neurons. The training process mainly requires an experimental database presenting the examples (T, G), as well as the corresponding values of the duty cycle D. The database D = f (G, T) used in this paper includes 600 examples sweeping a large margin of T and G, thus allowing the network to be able to predict the exact value of D regardless of the value of T and G.

    The training process is divided into two stages: the learning stage and the testing stage. During the learning stage, almost 500 examples are provided to the network during the convergence of the classical LM algorithm. During this phase, the network intelligently learns how to predict the value of the D factor according to T and G, which are given as inputs. The prediction capacity and efficiency of the network are demonstrated during the test phase using 100 examples of data D = f (G, T), which are completely different from the examples already processed during the learning stage.

    Figure 7.  The learning process of the ANN by the optimization of MSSE.

    The studied objective function "mean of sum squared errors (MSSE)" is defined as follows [42]:

    $ MSSE = \frac{1}{N}{\sum }_{r = 1}^{N}\left[{\sum }_{i = 1}^{Ns}\left[{\left({D}_{predicted}\left(i, r\right)-{D}_{target}\left(i, r\right)\right)}^{2}\right]\right] $ (7)

    where N is the number of examples, Ns is the number of output neurons, r denotes the index of the used target data, I indicates the index of the used output, Dpredicted presents the output value computed by the network and Dtarget is the target value of the output.

    The principle of this approach consists of starting the execution of the algorithm with the ANN model. The predicted value of D, which is close to the optimal one, is taken by the P & O method as an initial value to continue its convergence and judiciously seek the exact value of D corresponding to the MPP [43].

    After correctly training the ANN network with approximately 600 examples of D = f (T, G) during the two phases (i.e., learning and testing in the optimization process and identification of the correct values of the weights and the bias), the ANN tool provides the well-determined wij, wjm, bij and bjm, which are able to predict the correct value of the duty D for any value of T and G.

    The ANN executed first in the ANN-P & O algorithm is merely presented as a feedforward equation Eq 8 based on the well-identified wij, wjm, bij and bjm, thus resulting from the training process [44].

    $ D = {F}^{s}\left[{\sum }_{i = 1}^{Ns}\left[{w}_{ij}{F}^{c}\left({\sum }_{j = 1}^{Nc}\left({w}_{jm}{x}_{m}+{b}_{jm}\right)\right)+{b}_{ij}\right]\right] $ (8)

    I denotes the index of the output neurons, j presents the index of the hidden neurons, m is the index of the input neurons, yi indicates the output of the network, xm is the inputs of the network [T, G], Fs is the activation function of the output neurons, Fc is the activation function of the hidden neurons and Nc is the number of hidden neurons.

    The GN method is a conventional iterative technique aimed to ensure a parameter identification using descent directions dn. In this paper, GN is used to train the ANN network through the adjustment of the weights and the biases wij, wjm, bij and bjm at each iteration n during the optimization process of the MSSE error function, according to the following expressions cited in [45]:

    $ {\theta }_{n+1} = {\theta }_{n}+\alpha {d}_{n} $ (9)
    $ {d}_{n} = -\frac{{J}^{\text{'}}ϵ}{{J}^{\text{'}}J} $ (10)

    where dn is the descent direction at an iteration n, J is a matrix formed by four Jacobeans or the first four derivatives functions of Eq 8, J' presents the transpose of the Jacobean, $ \alpha $ indicates either the step size factor or the length of descent direction, $ ϵ $ denotes the generated error between the computed Eq 8 and the experimental Dtarget, θ is the vector of the fours weights and biases to be adjusted [wij, wjm, bij, bjm] and n is the number of iterations.

    Figure 8.  Presentation of ANN-P & O.

    The convergence of the GNSA algorithm requires declares the initial value of θ and the expression of either the four first derivatives or the Jacobeans "J" of Eq 8 at the start of its code. The latter presents D(θ). Each Jacobean among the four ones is obtained by deriving D(θ) according to a parameter from the four (wij, wjm, bij, bjm) parameters, such as: J = [$

    \begin{array}{ccc}\delta D/\delta {w}_{ij} & \delta D/\delta {w}_{jm} & \begin{array}{cc}\delta D/\delta {b}_{ij} & \delta D/\delta {b}_{jm}\end{array}
    \end{array} $]. During its execution, at an iteration n, GNSA determines the optimal value of αn by the SA technique and multiplies this value to the calculated value of the descent direction dn (Eq 10). The resulting value is added to the previously found value of θn to generate the adjusted value θn+1 to be processed during the next iteration. Just after, the iteration of n increments and this principle of GNSA functioning loops until finding the good value of θ, corresponding to the minimum error of MSSE.

    The SA approach was chosen after a further study of the behavior of MSSE according to large marge of α values at each iteration n of GNSA convergence during the training of the ANN. From this, the MSSE is observed with many local minima. To reach the global minimum, the SA is executed at each run of GN to find the most precise α value that leads to a good identification of the weights and the biases wij, wjm, bij and bjm. This combination between the GN and SA techniques defined by GNSA is used in this research as a new neural network training approach presented by the enhanced ANN-GNSA method.

    The training process of the ANN by the novel optimization GNSA technique is illustrated in the following flowchart and then explained in detail step by step.

    Step 1: - Give the experimental samples of the inputs (T, G) and the output D,

    - iteration n = 1,

    - Simultaneously start the two stages of learning and the test,

    Step 2: - Learn the ANN model by the improved GNSA method using the characteristic equations Eq. 9 and Eq. 10 to adjust the four parameters θ = [wij, wjm, bij, bjm],

    -The learning phase is carried out using its own experimental data of D and (T, G), fixed at 500 examples,

    Step 3: - Compute the MSSElearning criterion at an iteration n based on 100 examples of test.

    Step 4: - Give the adjusted weights θn at an iteration n,

    Step 5: - Compute the feedforward equation Eq 8 based on the learning examples and by using the adjusted weights θn at a fixed iteration n,

    - Compute the feedforward equation Eq 8 based on the 100 examples of the test data and by using the adjusted weights θn obtained at the same iteration n,

    Step 6: - Compute the H(n) = MSSEtest criterion,

    Step 7: - If H(n) < H(n - 1)

    • Reset a counter variable "c = 0",

    • Increment the iteration n = n + 1,

    • Loop to step 2.

    - Otherwise, increment the counter variable c = c + 1

    Step 8: - If c different to 6

    • Loop to step 2,

    - Otherwise, if c equal to 6

    • Stop the training process,

    Step 9: Give the optimal and the well-adjusted value of the θ parameter.

    Figure 9.  Flowchart of the proposed ANN-GNSA.

    To emphasize the efficiency of the new method, three techniques (P & O, ANN and ANN-P & O) have been examined in a detailed study to compare their results with those obtained by the new ANN-GNSA method.

    The used PV module is Solarex MSX-20L. It is comprised of 36 cells attached in series and characterized by the parameter's values shown in Table 1. On the other hand, the overall parameters values characterizing the structure of the step-up converter is shown in Table 2.

    Table 1.  The settings of Solarex MSX-20L module.
    Parameters Values
    Maximum power Pmax 20 W
    Current at maximum power IMPP 1.17 A
    Voltage at maximum power VMPP 17.1 V
    Short-circuit current Icc 1.27 A
    Open-circuit voltage Vcc 20.8 V

     | Show Table
    DownLoad: CSV
    Table 2.  Calculated parameters of the boost converter.
    Parameters Values
    Input capacitor Cin 771.6 µF
    Inductance L 189.8 µH
    Switching frequency f 30 kHz
    Output capacitor Cout 57.16 µF

     | Show Table
    DownLoad: CSV

    The control of a PV system is given by the "MPPT" subsystem which contains the code of the algorithm being processed (P & O, ANN, ANN-P & O or ANN-GNSA). The load connected to the output is a 30 Ω resistor.

    Figure 10.  Illustration under Simulink of a photovoltaic system based on MPPT.

    Equations (1), (2) and (3) have been implemented as blocks in the "Solarex module" subsystem to emulate the PV module, as shown in Figure 11.

    Figure 11.  Schema under Simulink of PV module Solarex MSX-20L.

    The graphical presentation of the boost converter is shown in the model illustrated in Figure 4. The corresponding model is presented by the hereby illustration.

    Figure 12.  Schema under Simulink of boost converter.

    A fixed value of the pair (G, T) corresponds to a given value of the response IPV(VPV). In this study, several responses measured experimentally for 600 different values of the pair (G, T) of the Solarex module were taken into consideration. Figures 13 and 14 show some experimental curves of IPV(VPV) for the variables G and T.

    Figure 13.  Three IPV(VPV) characteristic under variable irradiance.
    Figure 14.  Three IPV(VPV) characteristic under variable temperature.

    In order to extract all 600 values of D to be called out during the training of ANN at its output, we proposed to execute a P & O algorithm to generate the corresponding value of D for a precise curve of PPV (VPV). PPV is obtained by multiplying the VPV and the IPV. The P & O algorithm used is performed 600 times for 600 different PPV (VPV) curves obtained for 600 different values of T and G. The examples considered during the training process are formed by (G, T) and D.

    Five hundred examples of D = f (G, T) are used during the learning of ANN-GNSA and 100 examples D = f (G, T) are employed during the test of its effectiveness.

    The number of hidden neurons have been judiciously studied by observing MSSEtest for different number of hidden neurons, where 15 is the optimal number at which the obtained MSSEtest is the lowest one, as shown in Figure 15.

    Figure 15.  Study of the number of hidden neurons.

    The learning process of ANN-GNSA is based on the minimization of the objective function MSSElearning. Figure 16 shows a decreasing curve until the iteration 270. The learning stop is decided at iteration 223 instead of 270 to avoid over-learning caused by the memorization of the network of the set of examples D = f (G, T), which leads to an erroneous prediction of the values of D at the output.

    Figure 16.  Convergence of MSSElearning.

    A good decision of the learning stop is made based on the convergence of ANN-GNSA during the test phase. Indeed, when the MSSEtest curve successively increases six times during a decrease of MSSElearning, the learning of the ANN should stop (Figure 17).

    Figure 17.  Convergence of MSSEtest.

    The values of D found by ANN-GNSA are shown in Figure 18. An illustration of the target values and the predicted values of D in the same figure shows the compatibility amongst the two, which proves the effectiveness of ANN-GNSA.

    Figure 18.  Different samples of duty cycle D according to number of examples.

    To shed light on the exact maximum power value reached by each MPPT method (P & O, ANN, ANN-P & O or ANN-GNSA), the PPV have been presented according to time in Figure 19. The real target maximum power PPV corresponding to G = 1000 W/m2 and T = 25℃ is set at 18.59 W. As it is clearly seen, the ANN-GNSA method runs into an MPP equal to 18.59 W during 0.04382 s as compared to 18.55 W at 0.06667 s for the ANN-P & O, 18.28 W at 0.08008 s for the ANN and a range of power [18.31 W – 18.57 W] during [0.05808 s – 0.09156 s] for the P & O. Therefore, the proposed MPPT technique is the fastest and the most precise to track the right MPP value.

    Figure 19.  Evolution of the PV power according to time of P & O, ANN, ANN-P & O and of the proposed ANN-GNSA technique.

    The generated value of the control factor D of the step-up converter varies depending on the used MPPT method. Indeed, Figure 20 shows that the values of D given by the P & O are changeable, thus leading to oscillations at the PPV curve (Figure 19). D is fixed at 0.3103 for the ANN and at 0.28 for the ANN-GNSA, where the second one is best explained by the achieved power PPV = 18.59 W, which is considered the desired value. On the other hand, the ANN-P & O starts from the value of D = 0.3103 generated by the ANN and converges to a value of D = 0.2863 close to that provided by the ANN-GNSA. From these observations, we deduce that the value of D = 0.28 is the most optimal value for G = 1000 W/m2 and T = 25℃.

    Figure 20.  Evolution of the duty cycle D according to time of P & O, ANN, ANN-P & O and of the proposed ANN-GNSA technique.

    In Figure 21, the evolutions of PPV(VPV) obtained by the four studied MPPT methods were also taken into consideration. In the latter a vertical offset of 0.5 W is established between the four curves to shed light on the tracking accuracy provided by the improved ANN-GNSA method.

    Figure 21.  Evolution of the PV power according to PV voltage of P & O, ANN, ANN-P & O and of the proposed ANN-GNSA technique.

    The P & O, ANN and ANN-P & O methods can't correctly achieve the right MPPs, as is the case of G = 300 W/m2 and G = 400 W/m2. On other hand, the ANN method is less efficient at tracking the desired MPPs; it remains stuck in its neighborhood, especially in the case of G = 200 W/m2, G = 300 W/m2, G = 1000 W/m2 and G = 400 W/m2. However, the ANN-P & O and the proposed ANN-GNSA are more efficient in reaching the correct MPP values.

    To evaluate the performance of the P & O, ANN, ANN-P & O and ANN-GNSA techniques to reach to optimal value of MPP, we are based on the convergence time to be concluded from Figure 22 and on the tracking efficiency η to be calculated from the obtained maximum PPV (MPP) of each MPPT method through the following expression [46]:

    $ \eta = \frac{{P}_{PVmean}}{{P}_{PVmax}} $ (11)
    Figure 22.  Evolution of the PV power according to time under variable irradiance of P & O, ANN, ANN-P & O and of the proposed ANN-GNSA technique.

    where PPVmean is the mean PV power delivered by the photovoltaic module, and PPVmax is either the maximum PV power or the MPP.

    Table 3 shows the optimal power MPP to achieve at each irradiance, as well as the maximum power generated by each method among the four methods. The quality of the tracking of MPP by each technique is illustrated by the calculated efficiency values and by the reached convergence time (tc).

    Table 3.  The efficiencies and the tracking times obtained by the P & O, ANN, ANN-P & O and by the proposed ANN-GNSA method for variable irradiance.


    G(W/m2)
    Target MPP (W)
    P & O ANN ANN-P & O ANN-GNSA
    Obtained
    MPP
    (W)
    tc(s) ObtainedMPP
    (W)
    tc(s) Obtained MPP
    (W)
    tc(s) Obtained MPP
    (W)
    tc(s)
    200 1.858 1.586 0.1095 - 0.1152 0.4517 0.0239 1.553 0.1766 1.613 0.1
    300 4.026 3.437 0.5772 - 0.5831 2.243 0.5659 3.371 0.6028 4.026 0.5887
    700 12.72 12.68 1.045 - 1.047 12.7 1.038 12.72 1.105 12.72 1.04
    1000 18.59 18.55 1.541 - 1.549 18.28 1.547 18.55 1.66 18.59 1.517
    800 14.69 14.61 2.025 - 2.039 14.63 2.032 14.67 2.082 14.67 2.005
    400 6.676 5.9 2.598 - 2.604 5.766 2.606 6.076 2.624 6.676 2.594
    η (%) 96.93 92.34 97.28 99.54

     | Show Table
    DownLoad: CSV

    The P & O, ANN and ANN-P & O approaches take the longest time to reach the MPPs with a significant gap compared to their target values. The proposed ANN-GNSA is faster and tracks the MPP with high efficiency. The convergence time margins for the P & O are [0.1095 s-0.1152 s], [0.5772 s-0.5831 s], [1.045 s-1.047 s], [1.541 s-1.549 s], [2.025 s-2.039 s] and [2.598 s-2.604 s]. For the ANN method, the convergence times are 0.0239 s, 0.5659 s, 1.038 s, 1.547 s, 2.032 s and 2.606 s and for the ANN-P & O method, the convergence times are 0.1766 s, 0.6028 s, 1.105 s, 1.66 s, 2.082 s and 2.624 s. The speed of tracking the MPP by the proposed ANN-GNSA method is considered faster than the other three techniques for G set at (1000 W/m2, 800 W/m2 and 400 W/m2), as it is shown by the convergence times: 1.517 s, 2.005 s, 2.594 s. For a G defined at (200 W/m2, 300 W/m2 and 700 W/m2), the ANN is observed faster than the ANN-GNSA, though less precise in the reached MPP values. The MPP search efficiencies are as follows: 96.93% for the P & O, 92.34% for the ANN, 97.28% for the ANN-P & O and 99.54% for the proposed approach. From these all outcomes, the ANN-GNSA demonstrates its outstanding and its superiority to track the optimal MPP with a high efficiency during less convergence time.

    From Figure 24 and Table 4, the ANN-GNSA outperforms the other comparative methods in terms of the convergence times and the efficiency. Indeed, the proposed method has an efficiency of 99.98% as compared to 99.87% for the P & O, 95.98% for the ANN and 99.16% for the ANN-P & O. Moreover, the ANN-GNSA demonstrates its high speed in tracking the MPP for different temperature values: For T = 15℃ the tc = 0.2096 s, for T = 35℃ the tc = 0.8198 s, for T = 45℃ the tc 1.549 s and for T = 5℃ the tc = 2.301 s.

    Figure 23.  Evolution of the PV power according to voltage under variable irradiance of P & O, ANN, ANN-P & O and of the proposed technique ANN-GNSA.
    Figure 24.  Evolution of the PV power according to voltage under variable temperature of P & O, ANN, ANN-P & O and of the proposed technique ANN-GNSA.
    Table 4.  The efficiencies and the tracking times obtained by the P & O, ANN, ANN-P & O and by the proposed ANN-GNSA method for variable temperature.


    T(℃)

    Target MPP
    (W)
    P & O ANN ANN-P & O ANN-GNSA
    Obtained
    MPP
    (W)
    tc(s) ObtainedMPP
    (W)
    tc(s) Obtained MPP
    (W)
    tc(s) Obtained MPP
    (W)
    tc(s)
    15 19.39 19.32 0.2589 17.890 0.2657 19.16 0.2266 19.37 0.2096
    35 17.77 17.77 0.8724 17.61 0.8741 17.61 0.8849 17.77 0.8198
    45 17.35 17.32 1.559 16.84 1.608 17.31 1.558 17.35 1.549
    5 20.19 20.2 2.329 19.36 2.404 20.02 2.418 20.2 2.301
    η (%) 99.87 95.98 99.19 99.98

     | Show Table
    DownLoad: CSV

    A novel and enhanced method, named the ANN trained by GNSA, is proposed in this paper to control a step-up converter within a PV system under varying irradiances and temperatures. The ANN-GNSA is evolved to ensure the fast and accurate tracking of the MPP, under variable irradiance and temperature values, without oscillations. Some MPPT methods such as the P & O, ANN and ANN-P & O are used to shed light on the performances of the proposed controller. For a fixed G = 1000 W/m2 and T = 25℃, the ANN-GNSA has the smallest tracking time (0.04382 s) compared to other MPPT methods [0.05808 s – 0.09156 s] for the P & O, 0.08008 s for the ANN and 0.06667 s for the ANN-P & O. Under a variable irradiance, the obtained efficiency is 99.54% compared to 96.93% for the P & O, 92.34% for the ANN and 97.28% for the ANN-P & O; on other side, the variable temperature is equal to 99.98% as compared to 99.87% for the P & O, 95.98% for the ANN and 99.19% for the ANN-P & O. By varying the irradiance in the margin [200 W/m2, 300 W/m2, 700 W/m2, 1000 W/m2, 800 W/m2, 400 W/m2], the proposed MPPT controller converges faster toward the MPP with convergence times of 0.1 s, 0.5887 s, 1.04 s, 1.517 s, 2.005 s and 2.594 s, respectively. A similar behavior is observed by varying the temperature in the margin [15℃, 35℃, 45℃, 5℃]; the convergence times are 0.2096 s, 0.8198 s, 1.549 s and 2.301 s, respectively. From these outcomes, the enhanced ANN-GNSA is considered a promising MPPT technique due to its outstanding performance under variable irradiances and temperatures.

    In a future study, we will lead a comparative study of the behavior and the tracking time of the ANN-GNSA implemented on FPGA, Arduino and on a microcontroller.

    The author declares that Artificial Intelligence (AI) tools are not used in the creation of this article.

    The author declares that there is no conflict of interest in this paper.

    [1] United Nations, World Urbanization Prospects, the 2014 Revision. 2014. Available from: http://esa.un.org/unpd/ accessed March 2015.
    [2] DSEWPaC (Department of Sustainability, Environment, Water, Population and Communities), State of the Environment 2011, 2011. Available from: http://www.environment.gov.au/science/soe/2011-report/contents
    [3] Morris CJG, Simmonds I, Plummer N (2001) Quantification of the Influences of Wind and Cloud on the Nocturnal Urban Heat Island of a Large City. J Appl Meteorol 40: 169-182. doi: 10.1175/1520-0450(2001)040<0169:QOTIOW>2.0.CO;2
    [4] Coutts A, Beringer J, Tapper N (2010) Changing urban climate and CO2 emissions: implications for the development of policies for sustainable cities. Urban Policy Res 28: 27-47. doi: 10.1080/08111140903437716
    [5] Wu K, Yang XQ (2013) Urbanization and heterogeneous surface warming in eastern China. Chin Sci Bull 58: 1363-1373. doi: 10.1007/s11434-012-5627-8
    [6] DSE (Department of Sustainability and Environment), Melbourne 2030: Planning for Sustainable Growth. State Government of Victoria, 2002. Available from: http://www.dtpli.vic.gov.au/planning/plans-and-policies/planning-for-melbourne/melbournes-strategic-planning-history/melbourne-2030-planning-for-sustainable-growth
    [7] Konopacki S, Akbari H (2002) Energy savings for heat island reduction strategies in Chicago and Houston (including updates for Baton Rouge, Sacramento, and Salt Lake City). Draft Final Report, LBNL-49638, University of California, Berkeley.
    [8] Kolokotroni M, Giannitsaris I, Watkins R (2006) The effect of the London Urban Heat Island on building summer cooling demand and night ventilation strategies. Sol Energ 80: 383-392. doi: 10.1016/j.solener.2005.03.010
    [9] Priyadarsini R (2011) Urban Heat Island and its Impact on Building Energy Consumption. Adv Build Energ Res 3: 261-270.
    [10] Changnon SA, Kunkel KE, Reinke BC (1996) Impacts and responses to the 1995 heat wave: A call to action. Bull Am Meteorol Soc 77: 1497-1505. doi: 10.1175/1520-0477(1996)077<1497:IARTTH>2.0.CO;2
    [11] Haines A, Kovats RS, Campbell-Lendrum D, et al. (2006) Climate change and human health: impacts, vulnerability, and mitigation. Lancet 367: 2101-2109. doi: 10.1016/S0140-6736(06)68933-2
    [12] Cadot E, Rodwin VG, Spira A (2007) In the heat of the summer: lessons from the heat waves in Paris. J Urban Health 84: 466-468. doi: 10.1007/s11524-007-9161-y
    [13] Luber G (2008) Climate Change and Extreme Heat Events, Am J Prev Med 35: 429-435.
    [14] Nicholls N, Skinner C, Loughnan M, et al. (2008) A simple heat alert system for Melbourne, Australia, Int J Biometeorol 52: 375-384.
    [15] Kovats RS, Hajat S (2008) Heat stress and public health: a critical review. Annu Rev Public Health 29: 41-55. doi: 10.1146/annurev.publhealth.29.020907.090843
    [16] Loughnan ME, Nicholls N, Tapper NJ (2010) The effects of summer temperature, age and socioeconomic circumstance on acute myocardial infarction admissions in Melbourne, Australia. Int J Health Geogr 9: 41. doi: 10.1186/1476-072X-9-41
    [17] DHS, January 2009 Heatwave in Victoria: an assessment of health impacts. Victorian Government Department of Human Services, Victoria, 2009. Available from: http://docs.health.vic.gov.au/docs/doc/F7EEA4050981101ACA257AD80074AE8B/$FILE/heat_health_impact_rpt_Vic2009.pdf
    [18] Alexander LV, Arblaster J (2008) Assessing trends in observed and modelled climate extremes over Australia in relation to future projections. Int J Climatol 29: 417-435.
    [19] Solecki WD, Rosenzweig C, Parshall L, et al. (2005) Mitigation of the heat island effect in urban New Jersey. Global Environ Chang Part B: Environ Hazards 6: 39-49.
    [20] Gill SE, Handley JF, Ennos AR, et al. (2008) Characterising the urban environment of UK cities and towns: A template for landscape planning. Landscape Urban Plann 87: 210-222. doi: 10.1016/j.landurbplan.2008.06.008
    [21] Gaitani N, Spanou A, Saliari M, et al. (2011) Improving the microclimate in urban areas: a case study in the centre of Athens. Build Serv Eng Res Technol 32: 53-71. doi: 10.1177/0143624410394518
    [22] Kleerekoper L, van Esch M, Salcedo TB (2012) How to make a city climate proof, addressing the urban heat island effect. Resour Conserv Recycl 64: 30-38. doi: 10.1016/j.resconrec.2011.06.004
    [23] DCLG, Building a greener future: towards zero carbon development. UK: Department of Communities and Local Government, 2007. Available from: http://webarchive.nationalarchives.gov.uk/20120919132719/http://www.communities.gov.uk/documents/planningandbuilding/pdf/153125.pdf
    [24] European Commission, Low energy buildings in Europe: current state of play, definition and best practice. 2009. Available from:http: //regions202020.eu/cms/inspiration/resources/PublicResource/330/view
    [25] Wang XM, Chen D, Ren ZG (2010) Assessment of climate change impact on residential building heating and cooling energy requirement in Australia. Build Environ 45: 1663-1682. doi: 10.1016/j.buildenv.2010.01.022
    [26] Rosenfeld AH, Akbari H, Romm JJ, et al. (1998) Cool communities: strategies for heat island mitigation and smog reduction. Energ Build 28: 51-62. doi: 10.1016/S0378-7788(97)00063-7
    [27] Akbari H, Pomerantz M, Taha H (2001) Cool surfaces and shade trees to reduce energy use and improve air quality in urban areas. Sol Energ 70: 295-310. doi: 10.1016/S0038-092X(00)00089-X
    [28] U.S. Environmental Protection Agency, Reducing Urban Heat Islands: Compendium of Strategies. 2008. Available from: http: //www.epa.gov/hiri/resources/compendium.htm
    [29] Memon RA, Leung DYC, Liu CH (2008) A review on the generation, determination and mitigation of Urban Heat Island. J Environ Sci 20: 120-128. doi: 10.1016/S1001-0742(08)60019-4
    [30] Santamouris M (2012) Cooling the cities - a review of reflective and green roof mitigation technologies to fight heat island and improve comfort in urban environments. Sol Energ 103: 680-703.
    [31] Liu K, Bass B (2005) Performance of Green Roof Systems. National Research Council Canada, Report No. NRCC-47705, Toronto, Canada.
    [32] Rosenzweig C, Solecki W, Parshall L, et al. (2006) Mitigating New York City's Heat Island with Urban Forestry, Living Roofs, and Light Surfaces. Sixth Symposium on the Urban Environment and Forum on Managing our Physical and Natural Resources, American Meteorological Society, Atlanta, GA.
    [33] Alexandri E, Jones P (2008) Temperature decreases in an urban canyon due to green walls and green roofs in diverse climates. Build Environ 43: 480-493. doi: 10.1016/j.buildenv.2006.10.055
    [34] Wong JKW, Lau LSK (2013) From the ‘urban heat island’ to the ‘green island’? A preliminary investigation into the potential of retrofitting green roofs in Mongkok district of Hong Kong. Habit Int 39: 25-35.
    [35] Susca T, Gaffin SR, Dell'Osso GR (2011) Positive effects of vegetation: Urban heat island and green roofs, Environ Pollut 159: 2119-2126.
    [36] Xu T, Sathaye J, Akbari H, et al. (2012) Quantifying the direct benefits of cool roofs in an urban setting: Reduced cooling energy use and lowered greenhouse gas emissions. Build Environ 48: 1-6.
    [37] Smith GB, Aguilar JLC, Gentle AR, et al. (2012) Multi-parameter sensitivity analysis: A design methodology applied to energy efficiency in temperate climate houses. Energy Build 55: 668-673. doi: 10.1016/j.enbuild.2012.09.007
    [38] Bozonnet E, Musy M, Calmet I, et al. (2013) Modeling methods to assess urban fluxes and heat island mitigation measures from street to city scale. Int J Low Carbon Technol 0: 1-16.
    [39] Sailor DJ (1995) Simulated urban response to modifications in surface albedo and vegetative cover. J Appl Meteor 34: 1694- 1704. doi: 10.1175/1520-0450-34.7.1694
    [40] Taha H, Douglas S, Haney J (1997) Mesoscale meteorological and air quality impacts of increased urban albedo and vegetation. Energ Build 25: 169-177. doi: 10.1016/S0378-7788(96)01006-7
    [41] Shashua-Bar L, Hoffman ME (2000) Vegetation as a climatic component in the design of an urban street: An empirical model for predicting the cooling effect of urban green areas with trees. Energy Build 31: 221-35. doi: 10.1016/S0378-7788(99)00018-3
    [42] Kong FH, Yin HW, Wang CZ, et al. (2014) A satellite image-based analysis of factors contributing to the green-space cool island intensity on a city scale. Urban For Urban Gree 13: 846-853. doi: 10.1016/j.ufug.2014.09.009
    [43] Perini K, Magliocco A (2014) Effects of vegetation, urban density, building height, and atmospheric conditions on local temperatures and thermal comfort. Urban For Urban Gree 13:495-506. doi: 10.1016/j.ufug.2014.03.003
    [44] Chen C, Yue Y, Jiang W (2014) Numerical Simulation on Cooling Effects of Greening for Alleviating Urban Heat Island Effect in North China. Appl Mech Mater 675-677: 1227-1233. doi: 10.4028/www.scientific.net/AMM.675-677.1227
    [45] Yang JB, Liu HN, Sun JN, et al. (2015) Further Development of the Regional Boundary Layer Model to Study the Impacts of Greenery on the Urban Thermal Environment. J Appl Meteor Climatol 54: 137-152. doi: 10.1175/JAMC-D-14-0057.1
    [46] Gromke C, Blocken B, Janssen W, et al. (2015) CFD analysis of transpirational cooling by vegetation: Case study for specific meteorological conditions during a heat wave in Arnhem, Netherlands, Build Environ 83: 11-26.
    [47] Middel A, Chhetri N, Quay R (2015) Urban forestry and cool roofs: Assessment of heat mitigation strategies in Phoenix residential neighborhoods. Urban For Urban Gree 14: 178-186. doi: 10.1016/j.ufug.2014.09.010
    [48] Thatcher M, Hurley P (2012) Simulating Australian urban climate in a mesoscale atmospheric numerical model. Boundary Layer Meteorol 142: 149-175. doi: 10.1007/s10546-011-9663-8
    [49] Hurley P, Physick W, Luhar A (2005) TAPM - A practical approach to prognostic meteorological and air pollution modeling. Environ Model Sotfw 20: 737-752. doi: 10.1016/j.envsoft.2004.04.006
    [50] Masson V (2000) A physically-based scheme for the urban energy budget in atmospheric models. Boundary Layer Meteorol 94: 357-397. doi: 10.1023/A:1002463829265
    [51] Harman I, Barlow J, Belcher S (2004) Scalar fluxes from urban street canyons. Part II: Model. Boundary Layer Meteorol 113: 387-409. doi: 10.1007/s10546-004-6205-7
    [52] Coutts A, Beringer J, Tapper N (2007) Impact of increasing urban density on local climate: spatial and temporal variations in the surface energy balance in Melbourne, Australia. J Appl Meteor Climatol 46: 477-493. doi: 10.1175/JAM2462.1
    [53] McGregor J, C-CAM: Geometric aspects and dynamical formulation. CSIRO Atmospheric Research, 2005. Available from: www.cmar.csiro.au/e-print/open/mcgregor_2005a.pdf
    [54] McGregor J, Dix M (2008) An updated description of the conformal-cubic atmospheric model. High Resolution Simulation of the Atmosphere and Ocean, Hamilton, K. and Ohfuchi, W., Eds., Springer 51-76.
    [55] Katzfey J, McGregor J, Nguyen K, et al. Dynamical downscaling techniques: impacts on regional climate change signals. In: Anderssen RS, Braddock RD and Newham LTH (ed) 18th World IMACS congress and MODSIM09 international congress on modelling and simulation. 2009, July; Cairns. Modelling and Simulation Society of Australia and New Zealand. pp 2377-2383.
    [56] Nguyen K, Katzfey J, McGregor J (2011) Global 60 km simulations with the Conformal Cubic Atmospheric Model: Evaluation over the tropics. Clim Dyn 39: 637-654.
    [57] City of Melbourne, Urban Forest Strategy Make a Great City Greener 2013-2032. 2015. Available from: http://www.melbourne.vic.gov.au/Sustainability/UrbanForest/Documents/Urban_Forest_Strategy.pdf
  • This article has been cited by:

    1. Daniel Coombs, Michael A. Gilchrist, Colleen L. Ball, Evaluating the importance of within- and between-host selection pressures on the evolution of chronic pathogens, 2007, 72, 00405809, 576, 10.1016/j.tpb.2007.08.005
    2. Kazuhiko Nara, Five questions about Streptococcus mutans: Theoretical study of its transmission and colonisation, 2012, 57, 00039969, 1498, 10.1016/j.archoralbio.2012.03.008
  • 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(7024) PDF downloads(1291) Cited by(6)

Figures and Tables

Figures(8)  /  Tables(2)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog