Micro-climate

The micro-climatic/meteorological conditions measured close to the object or component are given as the second argument of the simulation functions shown earlier.

PlantBiophysics usually uses a special data structure from the PlantMeteo package to declare those conditions, and to pre-compute other required variables. This data structure is a type called Atmosphere.

The mandatory variables to provide are: T (air temperature in °C), Rh (relative humidity, 0-1), Wind (the wind speed in m s-1) and P (the air pressure in kPa).

We can declare such conditions using Atmosphere such as:

using PlantMeteo
meteo = Atmosphere(T = 20.0, Wind = 1.0, P = 101.3, Rh = 0.65)
Atmosphere(date = Dates.DateTime("2024-07-11T15:00:24.236"), duration = Dates.Second(1), T = 20.0, Wind = 1.0, P = 101.3, Rh = 0.65, Precipitations = 0.0, Cₐ = 400.0, e = 1.5255470730405223, eₛ = 2.3469954969854188, VPD = 0.8214484239448965, ρ = 1.2037851579511918, λ = 2.4537e6, γ = 0.06723680111943287, ε = 0.5848056484857892, Δ = 0.14573378083416522, clearness = Inf, Ri_SW_f = Inf, Ri_PAR_f = Inf, Ri_NIR_f = Inf, Ri_TIR_f = Inf, Ri_custom_f = Inf)

The Atmosphere also computes other variables based on the provided conditions, such as the vapor pressure deficit (VPD) or the air density (ρ). You can also provide those variables as inputs if necessary. For example if you need another way of computing the VPD, you can provide it as follows:

using PlantMeteo
Atmosphere(T = 20.0, Wind = 1.0, P = 101.3, Rh = 0.65, VPD = 0.82)
Atmosphere(date = Dates.DateTime("2024-07-11T15:00:24.569"), duration = Dates.Second(1), T = 20.0, Wind = 1.0, P = 101.3, Rh = 0.65, Precipitations = 0.0, Cₐ = 400.0, e = 1.5255470730405223, eₛ = 2.3469954969854188, VPD = 0.82, ρ = 1.2037851579511918, λ = 2.4537e6, γ = 0.06723680111943287, ε = 0.5848056484857892, Δ = 0.14573378083416522, clearness = Inf, Ri_SW_f = Inf, Ri_PAR_f = Inf, Ri_NIR_f = Inf, Ri_TIR_f = Inf, Ri_custom_f = Inf)

To access the values of the variables after instantiation, we can use the dot syntax. For example if we need the vapor pressure at saturation, we would do as follows:

meteo.eₛ
2.3469954969854188

See the documentation of the function if you need more information about the variables: Atmosphere.

If you want to simulate several time-steps with varying conditions, you can do so by using Weather instead of Atmosphere.

Weather is just an array of Atmosphere along with some optional metadata. For example for three time-steps, we can declare it like so:

using PlantMeteo
w = Weather(
    [
        Atmosphere(T = 20.0, Wind = 1.0, P = 101.3, Rh = 0.65),
        Atmosphere(T = 23.0, Wind = 1.5, P = 101.3, Rh = 0.60),
        Atmosphere(T = 25.0, Wind = 3.0, P = 101.3, Rh = 0.55)
    ],
    (
        site = "Montpellier",
        other_info = "another crucial metadata"
    )
)
TimeStepTable{Atmosphere{(:date, :duration,...}(3 x 22):
╭─────┬─────────────────────────┬──────────────┬─────────┬─────────┬─────────┬──
│ Row │                    date │     duration │       T │    Wind │       P │ ⋯
│     │          Dates.DateTime │ Dates.Second │ Float64 │ Float64 │ Float64 │ ⋯
├─────┼─────────────────────────┼──────────────┼─────────┼─────────┼─────────┼──
│   1 │ 2024-07-11T15:00:24.607 │     1 second │    20.0 │     1.0 │   101.3 │ ⋯
│   2 │ 2024-07-11T15:00:24.607 │     1 second │    23.0 │     1.5 │   101.3 │ ⋯
│   3 │ 2024-07-11T15:00:24.607 │     1 second │    25.0 │     3.0 │   101.3 │ ⋯
╰─────┴─────────────────────────┴──────────────┴─────────┴─────────┴─────────┴──
                                                              17 columns omitted
Metadata: `Dict("site" => "Montpellier", "other_info" => "another crucial metadata")`

As you see the first argument is an array of Atmosphere, and the second is a named tuple of optional metadata such as the site or whatever you think is important.

A Weather can also be declared from a DataFrame, provided each row is an observation from a time-step, and each column is a variable needed for Atmosphere (see the help of Atmosphere for more details on the possible variables and their units).

Here's an example of using a DataFrame as input:

using CSV, DataFrames, PlantMeteo
file = joinpath(dirname(dirname(pathof(PlantMeteo))),"test","data","meteo.csv")
df = CSV.read(file, DataFrame; header=5, skipto = 6, dateformat = "yyyy/mm/dd")
# Select and rename the variables:
select!(df, :date, :temperature => :T, :relativeHumidity => (x -> x ./ 100 ) => :Rh, :wind => :Wind, :atmosphereCO2_ppm => :Cₐ)
df[!,:duration] .= 1800 # Add the time-step duration, 30min

# Make the weather, and add some metadata:
Weather(df, (site = "Aquiares", file = file))
TimeStepTable{Atmosphere{(:date, :duration,...}(3 x 22):
╭─────┬────────────┬──────────┬─────────┬─────────┬─────────┬─────────┬─────────
│ Row │       date │ duration │       T │    Wind │       P │      Rh │ Precip ⋯
│     │ Dates.Date │    Int64 │ Float64 │ Float64 │ Float64 │ Float64 │        ⋯
├─────┼────────────┼──────────┼─────────┼─────────┼─────────┼─────────┼─────────
│   1 │ 2016-06-12 │     1800 │    25.0 │     1.0 │ 101.325 │     0.6 │        ⋯
│   2 │ 2016-06-12 │     1800 │    26.0 │     1.5 │ 101.325 │    0.62 │        ⋯
│   3 │ 2016-06-12 │     1800 │    25.3 │     1.5 │ 101.325 │    0.58 │        ⋯
╰─────┴────────────┴──────────┴─────────┴─────────┴─────────┴─────────┴─────────
                                                              16 columns omitted
Metadata: `Dict("file" => "/home/runner/.julia/packages/PlantMeteo/8tGGz/test/data/meteo.csv", "site" => "Aquiares")`

One can also directly import the Weather from an Archimed-ϕ-formatted meteorology file (a csv file optionally enriched with some metadata). In this case, the user can rename and transform the variables from the file to match the names and units needed in PlantBiophysics using a DataFrame.jl-alike syntax:

using Dates, PlantMeteo

meteo = read_weather(
    joinpath(dirname(dirname(pathof(PlantMeteo))),"test","data","meteo.csv"),
    :temperature => :T,
    :relativeHumidity => (x -> x ./100) => :Rh,
    :wind => :Wind,
    :atmosphereCO2_ppm => :Cₐ,
    date_format = DateFormat("yyyy/mm/dd")
)
TimeStepTable{Atmosphere{(:date, :duration,...}(3 x 29):
╭─────┬─────────────────────┬──────────────────────┬─────────┬─────────┬────────
│ Row │                date │             duration │       T │    Wind │       ⋯
│     │      Dates.DateTime │ Dates.CompoundPeriod │ Float64 │ Float64 │ Float ⋯
├─────┼─────────────────────┼──────────────────────┼─────────┼─────────┼────────
│   1 │ 2016-06-12T12:00:00 │           30 minutes │    25.0 │     1.0 │ 101.3 ⋯
│   2 │ 2016-06-12T12:30:00 │           30 minutes │    26.0 │     1.5 │ 101.3 ⋯
│   3 │ 2016-06-12T13:00:00 │           30 minutes │    25.3 │     1.5 │ 101.3 ⋯
╰─────┴─────────────────────┴──────────────────────┴─────────┴─────────┴────────
                                                              25 columns omitted
Metadata: `Dict{String, Any}("name" => "Aquiares", "latitude" => 15.0, "altitude" => 100.0, "use" => [:clearness], "file" => "/home/runner/.julia/packages/PlantMeteo/8tGGz/test/data/meteo.csv")`

Helper functions

PlantBiophysics provides some helper functions to compute some micro-climate related variables.

Here is a complete list of these functions:

  • vapor_pressure computes e (kPa), the vapor pressure from the air temperature and the relative humidity
  • e_sat computes eₛ (kPa), the saturated vapor pressure from the air temperature
  • air_density computes ρ (kg m-3), the air density from the air temperature, the pressure, and some constants
  • latent_heat_vaporization computes λ (J kg-1), the latent heat of vaporization from the air temperature and a constant
  • psychrometer_constant computes γ (kPa K−1), the psychrometer "constant" from the air pressure, the latent heat of vaporization and some constants
  • atmosphere_emissivity(T,e,constants.K₀) computes ε (0-1), the atmosphere emissivity from the air temperature, the vapor pressure and a constant
  • e_sat_slope computes Δ (0-1), the slope of the saturation vapor pressure at air temperature, from the air temperature
Note

All constants are found in Constants