calculate_temperature_response: Calculate temperature-dependent parameter values

In PhotoGEA: Photosynthetic Gas Exchange Analysis

Calculate temperature-dependent parameter values

Description

Calculate leaf-temperature-dependent values of various parameters using various temperature response functions.

Usage

  calculate_temperature_response(
    exdf_obj,
    temperature_response_parameters,
    tleaf_column_name = 'TleafCnd'
  )

Arguments

exdf_obj	An exdf object representing data from a Licor gas exchange measurement system.
temperature_response_parameters	A list, where each element describes the temperature response of a parameter value. The name of each element must be the name of the parameter. Each element must be a list itself, whose named elements must include the type of temperature response function to use (type), thee units of the parameter (units), and the values of necessary temperature response parameters. See below for more details.
tleaf_column_name	The name of the column in exdf_obj that contains the leaf temperature in units of degrees C.

Details

Some key photosynthetic parameters are known to vary with temperature according to well-established temperature response functions such as the Arrhenius equation. The calculate_temperature_response function can be used to calculate such temperature-dependent parameter values at leaf temperature.

Depending on the type value supplied in each element of temperature_response_parameters, one of several possible functions will be used to calculate the temperature response:

• When type is 'Arrhenius', the calculate_temperature_response_arrhenius function will be used.

• When type is 'Gaussian', the calculate_temperature_response_gaussian function will be used.

• When type is 'Johnson', the calculate_temperature_response_johnson function will be used.

• When type is 'Polynomial', the calculate_temperature_response_polynomial function will be used.

Values of type are not case-sensitive.

It is rare to directly specify these parameters; instead, it is more typical to use one of the pre-set values such as those included in c3_temperature_param_sharkey.

Value

An exdf object based on exdf_obj that includes one new column for each element of temperature_response_parameters, where the temperature-dependent values of these new columns are determined using the temperature values specified by the tleaf_column_name column. The category of each of these new columns is calculate_temperature_response to indicate that they were created using this function.

Examples

# Read an example Licor file included in the PhotoGEA package
licor_file <- read_gasex_file(
  PhotoGEA_example_file_path('ball_berry_1.xlsx')
)

# In this example we will calculate temperature-dependent values of two
# parameters:
#
# - The `Kc` parameter (in units of `micromol mol^(-1)`) will be calculated
#   using an Arrhenius function with scaling constant `c` = 38.05 and activation
#   energy `Ea` = 79.43 kJ / mol.
#
# - The `Jmax` parameter (in units of `micromol m^(-2) s^(-1)) will be
#   using a Gaussian function with optimal temperature `t_opt` = 43 degrees C
#   and width `sigma` = 16 degrees C.
#
# So the `temperature_response_parameters` list will contain two elements,
# defined as follows:

trp <- list(
  Kc = list(
    type = 'Arrhenius',
    c = 38.05,
    Ea = 79.43,
    units = 'micromol mol^(-1)'
  ),
  Jmax = list(
    type = 'Gaussian',
    optimum_rate = 4,
    t_opt = 43,
    sigma = 16,
    units = 'micromol m^(-2) s^(-1)'
  )
)

# Now we can calculate the values of Kc and Jmax at the measured leaf
# temperatures recorded in the log file
licor_file <- calculate_temperature_response(licor_file, trp)

licor_file$units$Kc      # View the units of the new `Kc` column
licor_file$categories$Kc # View the category of the new `Kc` column
licor_file[,'Kc']        # View the values of the new `Kc` column

licor_file$units$Jmax      # View the units of the new `Jmax` column
licor_file$categories$Jmax # View the category of the new `Jmax` column
licor_file[,'Jmax']        # View the values of the new `Jmax` column

PhotoGEA documentation built on April 11, 2025, 5:48 p.m.