Weather Ελληνικά Sign in
Home › Article

Soil moisture: what the station shows, when the field needs irrigating, and how it binds to evapotranspiration and VPD

Irrigation Soil Climate
Diagram of the ladder of soil water, from saturation to the permanent wilting point, with volumetric water content and tension
The same quantity of water is described by two entirely different numbers. Volumetric water content says how much water is there, tension says how easily the root can take it. The station usually displays the first. The plant feels the second.

Savvas Pastopoulos, Agronomist MSc, Pastopoulos Agronomics, Neos Mylotopos, Pella, Greece

The agrometeorological station shows a percentage for soil moisture, and the grower is asked to decide whether to irrigate. The percentage on its own is not a decision. The same "22%" means a saturated sand and a dry clay. And the point at which a field starts to run short of water is not fixed, because it moves with how much water the atmosphere is demanding that day. What follows is what the sensor measures, how the percentage becomes millimetres of water, where the irrigation limit sits, what evapotranspiration is and how it is computed step by step, what vapour pressure deficit is, and two calculators for the field.

A note on units

The Greek edition of this article works in stremmata; this English edition uses hectares throughout. 1 ha = 10 stremmata. Depths of water in millimetres are the same in both: 1 mm = 1 litre per square metre = 10 m³ per hectare.

The number that is missing from the screen

A soil moisture reading means something only with three further pieces of information, and those are rarely published. Which quantity is measured, at what depth, and whether the sensor was calibrated on that particular soil. Without them, the percentage rises and falls correctly but does not correspond to a quantity of water. How large the error becomes is shown in section 3.

-33 kPathe tension at field capacity, and -1,500 kPa at the permanent wilting point [6]
+115%error in volumetric water content on a clay soil with factory calibration [4]
2.1 kPathe VPD at which kiwifruit stomatal conductance falls to half [11]
4.2 kPathe air VPD of an afternoon at 35 °C and 25% relative humidity

1. Four different quantities, one word

The words "soil moisture" are used for four different quantities. Two say how much water is present. One says how tightly the soil holds it. The fourth is not measured at all, it is computed.

QuantityUnitWhat it actually says
Volumetric water content, θv m³ of water per m³ of soil, or the same times one hundred as % The volume of water over the volume of soil. It measures amount. Multiplied by depth it gives millimetres of water [2].
Gravimetric water content, w kg of water per kg of dry soil, or % The mass of water over the mass of dry soil. It is what the laboratory oven produces. It converts to volumetric with the bulk density, θv = w × ρb [2].
Tension or water potential kPa, centibar, or pF How much force the root must spend to extract the water. It measures availability, not amount [3]. One centibar equals one kPa.
Percent of field capacity or percent available water % dimensionless Where today's moisture sits between the wilting point and field capacity. It is not measured. It needs both ends of the range for that particular soil.

The distinction has a practical consequence. Two soils may both read 25%, one of them soaking and the other one wilting the plant. At the same percentage, the tension in a sand and in a clay is entirely different. The plant does not perceive amount, it perceives tension.

The pF, for anyone who meets it in older texts

The pF is the tension written logarithmically, by analogy with pH [30]. Field capacity corresponds to about pF 2 and the permanent wilting point to pF 4.2.

Sources: 2 3 30

2. What the percentage on the screen means

The stations installed in Greek fields carry sensors from different manufacturers. The word on the screen is the same. The quantity behind it is not.

The catalogue of a single manufacturer contains five soil moisture instruments returning four different units. Some return volumetric water content in m³/m³, one returns a percentage from 0 to 100, and two return tension in kPa [10]. All five are sold as "soil moisture sensors".

What is not published, and that is a finding

A search of the public documentation of the agrometeorological station networks operating in Greece found no statement at all on three things, without which the percentage does not become a number with meaning.

  • Which physical quantity the percentage is. Volumetric water content, percent of field capacity, or an index between the two ends of the available water. The three give entirely different numbers for the same soil.
  • At what depths the sensors sit. The wording that is found is "at various depths", without centimetres.
  • Whether a soil-specific calibration was carried out. This is the most important of the three, for the reason the next section explains.

This finding is not an accusation against any provider, it is the state of the public documentation as it was found. The consequence is that these three questions are put to the provider in writing before the reading is used for an irrigation decision.

Until that answer exists, the reading remains useful as relative monitoring. It shows when the soil wets, how fast it dries, how deep the irrigation water reached, and when the curve stopped falling, which means the root stopped drawing. It remains useless as an absolute quantity, that is for computing millimetres of water or for comparison with a limit from a textbook table.

What the curve says, even without calibration

After an irrigation the curve rises sharply, stays high for a few hours and then falls at two rates. The sharp fall of the first hours is drainage. The slow fall that follows is evapotranspiration. When the slow fall stops and the line goes flat while the weather stays hot, the plant has closed its stomata and is no longer drawing water from that depth. All of this is visible without knowing the unit.

Sources: 10

3. How the sensor works and why it needs calibration

The capacitance sensor, the most widespread type in stations, does not measure water. It measures an electrical property of the mixture of soil, water and air, the apparent dielectric permittivity. Water is a strongly polar molecule, minerals and air are not, so the more water is present the higher the permittivity [3]. Between permittivity and water content sits a calibration equation, and that is where the problem is.

What shifts the curve

The factory equation is an average for mineral soils. Four properties of the particular soil displace it, bulk density, mineralogy, texture and salinity [3]. The manufacturer quotes an accuracy of 3 to 4 percentage points with the factory curve and 1 to 2 points after a calibration in the field [3]. The figures measured in real soil are worse.

How large the error is in practice

On a clay soil, with the factory settings, relative errors reached +115% for sensors that rely on the real part of the permittivity and +245% for sensors that rely on the modulus. The average overestimation was 0.10 m³/m³ in the drier part of the range. After a calibration on the same soil, errors fell below 0.05 m³/m³ [4].

An overestimation of 0.10 m³/m³ means the station reads 25% while the soil holds 15%. For a sandy loam with a permanent wilting point at 11%, that difference is the whole decision. At 25% the field is comfortable. At 15% it is four points above wilting.

The same was found in a Greek study. An evaluation of two commercial capacitance sensors found significant deviations from the behaviour the manufacturer states, and derived calibration equations for the particular soil [5].

What this means for the grower

It does not mean the sensor is useless. It means the relative change is reliable and the absolute value is not. The curve rises and falls correctly. The number on the vertical axis may be displaced by ten points.

The cheap correction is two reference points taken from the field itself. The upper point is taken two days after a deep irrigation or a heavy rain, once drainage has stopped. The lower one is taken at the end of a long dry period without irrigation, once the curve has gone flat. The range between them is the available water of that field, expressed in the language of its own sensor. After that, the percentage becomes interpretable even if what it measures is still unknown.

Tensiometers do not have this problem

The tensiometer and the granular matrix sensor measure tension directly, that is the quantity the plant feels, and they need no per-soil calibration. The granular matrix sensor is not affected by salinity either [10]. The price is that a tensiometer works only up to about 80 kPa, needs maintenance, and does not give millimetres of water. For a field where the question is when to irrigate, tension is the better quantity. For how much, volumetric water content is needed.

Sources: 3 4 5 10

4. Field capacity and permanent wilting point

The available water of a soil has two ends. The upper one is field capacity, the moisture that remains once the gravitational water has drained, one to three days after saturation. The lower one is the permanent wilting point, the moisture at which the plant does not recover even in a saturated atmosphere.

The tension values

The classic limits are -33 kPa for field capacity and -1,500 kPa for the permanent wilting point, that is one third of an atmosphere and fifteen atmospheres [2]. In the applied literature they are given as ranges, 10 to 33 kPa and 500 to 3,000 kPa [6].

Texture moves the upper end. For coarse-textured soils field capacity sits near -10 kPa, because the large pores empty at a much lower tension [7]. If -33 kPa is used on a sandy soil, the available water comes out smaller than the real one.

Two values that are repeated without having been measured

The engineer who built the instruments with which those values were established writes that in his own work the permanent wilting point never once came out at fifteen atmospheres, and that field capacity as a steady state a few days after rain is something he has never seen [7]. He adds something that bears directly on the alluvial soils of the plain. When a fine-textured horizon sits on a coarse-textured one, drainage stops at the interface and the profile holds twice as much water as the -33 kPa value predicts.

The values remain useful as reference points. The consequence is that a measurement in the field overrides the table wherever one exists, and that in a layered soil the table underestimates.

The volumetric values by texture

The table below gives the typical ranges by texture class, as FAO-56 lists them [1].

Soil textureField capacity θFCWilting point θWPDifference, available water
Sand0.07 - 0.170.02 - 0.070.05 - 0.11
Loamy sand0.11 - 0.190.03 - 0.100.06 - 0.12
Sandy loam0.18 - 0.280.06 - 0.160.11 - 0.15
Loam0.20 - 0.300.07 - 0.170.13 - 0.18
Silt loam0.22 - 0.360.09 - 0.210.13 - 0.19
Silt0.28 - 0.360.12 - 0.220.16 - 0.20
Silty clay loam0.30 - 0.370.17 - 0.240.13 - 0.18
Silty clay0.30 - 0.420.17 - 0.290.13 - 0.19
Clay0.32 - 0.400.20 - 0.240.12 - 0.20

Values in m³ of water per m³ of soil. For a percentage, multiply by one hundred. The FAO-56 table has no rows for clay loam and sandy clay loam, which are common classes in the Greek plain.

For the missing classes, and as a second independent source, the following are mean values from American data, in percentages [6].

Soil textureField capacityWilting pointAvailable water
Sand10%4%6%
Loamy sand16%7%9%
Sandy loam21%9%12%
Loam27%12%15%
Silt loam30%15%15%
Sandy clay loam36%16%20%
Sandy clay32%18%14%
Clay loam29%18%11%
Silty clay loam28%15%13%
Silty clay40%20%20%
Clay40%22%18%

Averages from American soils, which have not been measured in Greek ones. The same source also gives measurements from real Oklahoma fields, where loam came out at 25% and 13% instead of 27% and 12%. The departure of the measured values from the averages is itself the argument for measuring in the field.

Sources: 1 2 6 7

5. From the percentage to millimetres, the available water

The percentage becomes useful the moment it is converted into millimetres. Millimetres are the common language of rainfall, irrigation and evapotranspiration. One millimetre of water equals one litre per square metre, that is ten cubic metres per hectare.

The two quantities

TAW = 1,000 × (θFC - θWP) × Zr

where TAW is the total available water of the root zone in millimetres, θFC and θWP are field capacity and wilting point in m³/m³, and Zr is the rooting depth in metres [1].
RAW = p × TAW

where RAW is the readily available water, that is the part of the available water that may be depleted without reducing evapotranspiration, and p is the depletion fraction [1].

The p is the number that does the work. It is the fraction of the available water that can be depleted before water stress appears, and FAO-56 tabulates it by crop. It ranges from 0.30 to 0.70, with 0.50 as the common value [1].

CropRooting depth ZrDepletion fraction p
Peach, apricot, plum1.0 - 2.0 m0.50
Cherry, apple, pear1.0 - 2.0 m0.50
Kiwifruit0.7 - 1.3 m0.35
Olive1.2 - 1.7 m0.65
Almond1.0 - 2.0 m0.40
Walnut1.7 - 2.4 m0.50
Grape, table1.0 - 2.0 m0.35
Grape, wine1.0 - 2.0 m0.45
Citrus, 70% canopy1.2 - 1.5 m0.50
Cotton1.0 - 1.7 m0.65
Maize, grain1.0 - 1.7 m0.55
Tomato0.7 - 1.5 m0.40
Potato0.4 - 0.6 m0.35
Alfalfa1.0 - 2.0 m0.55
Sugar beet0.7 - 1.2 m0.55

From FAO-56 Table 22 [1]. Rooting depth is given as a range because it depends on the rootstock, the soil and the irrigation system. Under drip with closely spaced emitters the active root is shallower than the range of the table, and the available water smaller.

A worked example

A peach orchard on a sandy loam. From the table, θFC = 0.23 and θWP = 0.11, a difference of 0.12. The effective rooting depth is 1.2 metres.

TAW = 1,000 × 0.12 × 1.2 = 144 millimetres, that is 1,440 cubic metres of water per hectare.

With p = 0.50 for stone fruit, RAW = 0.50 × 144 = 72 millimetres. That is what may be consumed before the tree begins to suffer.

If the crop is using 6 millimetres a day, the 72 millimetres run out in twelve days. If it is using 3, in twenty-four.

Why one sensor at one depth is not enough

The root zone does not empty evenly. The root draws first from the upper layers, where the water is given up most easily, and moves down as those dry. A sensor at 20 centimetres shows a fall long before the profile is empty. One at 60 centimetres shows stability while the upper part has already dried out.

The recommendation for tree crops is two to three measurements at different depths per position, typically around 30, 60 and 90 centimetres [8]. The deepest sensor is the one that judges the dose. It says whether the water reached that far, that is whether the dose was sufficient, or whether it passed through, that is whether it is draining below the root.

Sources: 1 8

6. When the soil is considered to need irrigating

There are three ways to answer the question, and they are not equivalent.

The first way, with tension

It is the simplest, because the soil need not be known. For orchards, the general recommendation is to begin irrigating at 30 kPa, that is 30 centibar [8]. In field and vegetable crops the limits are more differentiated.

CropTension limit for starting irrigation
Celery20 - 30 kPa
Potato30 - 50 kPa
Lettuce40 - 60 kPa
Onion45 - 65 kPa
Citrus50 - 70 kPa
Carrot55 - 65 kPa
Cabbage, cauliflower60 - 70 kPa
Tomato60 - 150 kPa
Sweet corn50 - 80 kPa
Alfalfa80 - 150 kPa
Cotton100 - 120 kPa

Values from American applied literature [6]. No corresponding Greek table of tension limits was found, neither in Ministry guidance nor in a published Greek study. They are given as an order of magnitude and not as a recommendation for a Greek field, because the limit depends on the irrigation system and on how frequently it runs.

The same tension limit corresponds to an entirely different depletion depending on the soil. The tension at which half the available water is gone is 20 kPa in sand, 30 in loamy sand, 40 in sandy loam, 65 in loam and 90 in silt loam [8]. A limit of 30 kPa empties half the available water of a sand, and barely one fifth of a silt loam.

The second way, with the depletion fraction

This is the FAO-56 way. Irrigation starts when the depletion of the root zone reaches RAW, that is p times the available water. It requires the volumetric water content and both ends of the range to be known. It has one large advantage, though, the limit is not fixed, it adapts to the demand of the atmosphere. Section 11 develops this.

The third way, with the water balance

This uses no soil sensor. It starts from a full root zone, subtracts the crop evapotranspiration every day, adds rainfall and irrigation, and irrigates when the balance reaches the limit [9]. It is bookkeeping, and it works. Its weakness is that errors accumulate. Every rainfall the station mismeasured and every irrigation that went unrecorded shifts the balance for the rest of the season.

The best practice is the combination

The water balance says how much to apply, because it counts in millimetres. The soil sensor says whether it happened, because it shows whether the water reached the depth it should have and whether the root is taking it. A sensor without a balance never says how much water to apply. A balance without a sensor drifts silently.

Sources: 1 6 8 9

7. Evapotranspiration, what it is exactly

Diagram showing evapotranspiration as the sum of evaporation and transpiration, and the three steps from ET0 to ETc adj
Evapotranspiration is the sum of two flows with entirely different control. Evaporation is a physical process governed by the weather. Transpiration is a biological one governed by the plant. The calculation always takes three steps, and the station supplies only the first.

The two flows

Water leaves a field by two routes that happen at the same time and are not easily separated in a measurement.

Evaporation is the conversion of liquid water into vapour from the soil surface, from the wet surfaces of leaves after rain or spraying, and from free water. It is purely physical. It is fuelled by energy, mainly solar radiation, and governed by temperature, by the wind that removes the saturated vapour from the surface, and by how wet the surface itself is. As the upper soil layer dries, evaporation falls sharply, because the water must now rise from deeper down.

Transpiration is the loss of water from inside the plant to the atmosphere, through the stomata. Water reaching the leaf evaporates from the cell walls inside the mesophyll, and the vapour leaves through the stoma. It is biological in the sense that the plant controls it: it opens the stomata to take in carbon dioxide and closes them when the loss of water becomes threatening. The important point is that transpiration is not waste. It is the unavoidable price of photosynthesis, because the same opening that lets carbon dioxide in lets water out.

In an early crop, with small plants and bare soil, most of the evapotranspiration is evaporation. In a fully developed canopy that shades the ground, most of it is transpiration. The same value of six millimetres means entirely different things in the two cases, and calls for entirely different measures: in the first, evaporation is reduced by covering the soil and by irrigating without wetting the surface; in the second, transpiration is not reduced without reducing yield as well.

Why it is measured in millimetres and not in litres

Evapotranspiration is expressed in millimetres of water depth per day, exactly like rainfall, so that the two can be subtracted directly. One millimetre corresponds to one litre per square metre or ten cubic metres per hectare. A day with an evapotranspiration of 7 millimetres empties seventy cubic metres per hectare. Over ten hectares of peach, seven hundred cubic metres a day.

The energy behind the number

Evaporating one kilogram of water takes about 2.45 megajoules of energy, which is the latent heat of vaporisation at ordinary temperatures [1]. That explains two things that look unrelated.

The first is that evapotranspiration has an upper limit set by the sun. As much energy as reaches the surface, that much water can at most leave. This is why a cloudy day has a low evapotranspiration however warm it is, and why December evapotranspiration in the Peloponnese is around one millimetre a day while August approaches five and a half [22].

The second is that transpiration cools the leaf. The energy spent evaporating the water no longer heats the tissue. A leaf transpiring freely is cooler than a leaf that has closed its stomata, and the difference can be several degrees. This is the basis of thermal detection of water stress, and it returns in section 10.

Sources: 1 22

8. Computing evapotranspiration, step by step

Nobody measures the evapotranspiration of a field directly. Direct measurement needs a lysimeter or an eddy covariance tower, and is done at research installations. In practice evapotranspiration is computed from meteorological data, in three steps.

8.1 First step, the reference surface

The problem to be solved is that evapotranspiration depends at once on the weather and on the crop itself. The solution was to define a hypothetical reference surface, the same everywhere, so that all the weight falls on the weather. FAO-56 defines it as a hypothetical reference crop with an assumed crop height of 0.12 metres, a fixed surface resistance of 70 seconds per metre and an albedo of 0.23 [1]. It is, in practice, a well-watered grass that is never short of water.

The evapotranspiration of that surface is called reference evapotranspiration and is written ET0. It is a purely meteorological quantity. It describes no real crop, and that is the point: two fields with different crops in the same area have the same ET0.

8.2 The Penman-Monteith equation

ET0 = [ 0.408 · Δ · (Rn - G) + γ · 900/(T+273) · u2 · (es - ea) ] / [ Δ + γ · (1 + 0.34 · u2) ]
SymbolWhat it isUnit and comment
ET0Reference evapotranspirationmillimetres per day
RnNet radiation at the surfaceMJ per m² per day. It is the energy left once what is reflected and what is re-emitted have been subtracted.
GSoil heat flux densityMJ per m² per day. On a daily basis it is very small and is often taken as zero.
TMean daily air temperature at 2 metres°C
u2Wind speed at 2 metresmetres per second. If the anemometer sits at another height, the value is adjusted to 2 metres.
es - eaSaturation vapour pressure deficitkPa. This is the VPD. It enters the evapotranspiration equation directly.
ΔSlope of the saturation vapour pressure curvekPa per °C. Computed from temperature.
γPsychrometric constantkPa per °C. Depends on atmospheric pressure, therefore on elevation.
0.408Conversion factorConverts MJ per m² into millimetres of water. It is the reciprocal of the latent heat of vaporisation.

The structure of the equation says something worth reading without the mathematics. The first term of the numerator, the one with the radiation, is the energy side: how much energy there is to evaporate water. The second term, the one with the wind and the vapour pressure deficit, is the aerodynamic side: how fast the atmosphere removes the vapour. A day with strong sun and no wind is driven by the first. A day with a dry wind is driven by the second, and this is why the etesian winds dry a field faster than the temperature alone would justify.

8.3 Second step, the crop coefficient

ETc = Kc × ET0

The crop evapotranspiration under standard conditions, that is, in disease-free, well-fertilised crops, in large fields and with optimum soil water [1].

The crop coefficient Kc absorbs all the differences between the real crop and the reference grass: height, albedo, canopy resistance, and how much soil stays exposed. It changes through the season, which is why it is given as three values, for the initial, the mid and the late stage.

CropKc initialKc midKc endHeight
Stone fruit, no ground cover, killing frost0.450.900.653.0 m
Stone fruit, active ground cover, killing frost0.501.150.903.0 m
Cherry, apple, pear, no ground cover, killing frost0.450.950.704.0 m
Cherry, apple, pear, active ground cover, killing frost0.501.200.954.0 m
Kiwifruit0.401.051.053.0 m
Olive, 40 to 60% ground cover by canopy0.650.700.703 - 5 m
Almond, no ground cover0.400.900.655.0 m
Walnut orchard0.501.100.654 - 5 m
Grape, table0.300.850.452.0 m
Grape, wine0.300.700.451.5 - 2.0 m
Cotton0.351.15 - 1.200.70 - 0.501.2 - 1.5 m
Maize, grain0.301.200.60 - 0.352.0 m
Tomato, processing0.601.150.70 - 0.900.6 m
Potato0.501.150.750.6 m
Sugar beet0.351.200.700.5 m

From FAO-56 Table 12 [1]. FAO-56 does not separate peach from apricot and plum, it gives them as one group, and likewise cherry together with apple and pear. The pair of values in the late Kc refers to the condition of the crop at harvest.

The difference the ground cover makes is very large

In stone fruit, the mid-season Kc is 0.90 with no ground cover and 1.15 with an active one. At an ET0 of 7 millimetres, that is the difference between 6.3 and 8.1 millimetres a day, nearly eighteen cubic metres per hectare every day. A ground cover is not a neutral choice for the water balance, and an irrigation programme that ignores it falls short consistently.

8.4 Adjusting Kc for the Greek climate

The table values are calibrated for a sub-humid climate, with a minimum relative humidity around 45% and wind around 2 metres per second. The Greek summer is not that. FAO-56 gives an explicit adjustment.

Kc = Kc(table) + [ 0.04 · (u2 - 2) - 0.004 · (RHmin - 45) ] × (h/3)0.3

with u2 from 1 to 6 m/s, RHmin from 20 to 80%, and h the mean crop height from 0.1 to 10 metres [1].

The adjustment raises Kc whenever the air is drier than 45% or the wind stronger than 2 m/s, which is almost always the case in the Greek plain in July. For a 3-metre peach orchard with an active ground cover, a minimum relative humidity of 25% and a wind of 2 m/s, the term (h/3)0.3 equals one and the adjustment gives:

Kc = 1.15 + [ 0 - 0.004 × (25 - 45) ] × 1 = 1.15 + 0.08 = 1.23

At an ET0 of 7 millimetres, the difference between 1.15 and 1.23 is half a millimetre a day, that is 150 cubic metres per hectare over a month. For the late Kc the same adjustment is applied only when the table value exceeds 0.45 [1].

8.5 Third step, the stress coefficient

ETc adj = Ks × Kc × ET0

The crop evapotranspiration under non-standard conditions, that is, the field as it really is [1].
Ks = (TAW - Dr) / (TAW - RAW)

where Dr is the current root zone depletion in millimetres. As long as depletion has not reached RAW, Ks equals 1 and there is no stress. Beyond RAW, Ks falls linearly to zero at the wilting point [1].

This step is often skipped, and the error runs in one direction. Once the soil has dried, actual consumption is lower than ETc, because the plant has closed its stomata. A balance that subtracts ETc every day without the Ks coefficient overestimates depletion and leads to more irrigation than the real requirement.

Sources: 1

9. Greek evapotranspiration figures

The general FAO-56 classification places Northern Greece in July in the category "temperate region, arid and semi-arid, warm", with a range of 6 to 9 millimetres per day [1]. That is a generic bracket and not a measurement. The Greek measurements that were found are the following.

LocationData periodAnnual ET0Summer value
Kozani, elevation 625 m48 years, 1962 to 2010952 millimetresJuly 166.8 millimetres, about 5.4 a day [19]
Krania, Grevena, elevation 952 m34 years, 1961 to 1994822.8 millimetresJuly 143.9 millimetres, about 4.6 a day [19]
Northern Greece, national estimatenot stated1,058 millimetres [20]not stated
South-eastern Greece, arid and semi-aridnot stated2,015 millimetres [20]not stated
Thessaly plain, simulation1961 to 1990, six-model ensemble900 to 1,100 millimetres [21]June to August, 420 to 540 millimetres
Peloponnese, 62 stationsAugusts 2016 to 2019not statedAugust 5.42 to 5.62 millimetres a day [22]
There is no published series for Pella or Imathia

The search found no published reference evapotranspiration series for the Giannitsa plain, for Imathia, or for Thessaloniki with a daily July figure. The nearest documented one is Kozani, at an elevation of 625 metres. A lowland field at twenty metres would be expected to have a higher evapotranspiration than Kozani, but no measurement supports that and no figure is given here. The gap is a research question, not something to be filled with an estimate.

Projected increase

For the Thessaly plain, the six-model ensemble projects an increase in annual reference evapotranspiration of 6 to 7% by the 2021 to 2050 period and of 10.5% by 2071 to 2100, with the summer increase at 10.5 to 12% [21]. These are simulations and not observations.

How much water Greek crops consume

The published Greek figures are fewer than one would expect, and most of them concern field crops.

CropQuantityWhere the figure comes from
Maize500 to 800 millimetres in total, of which about 300 from irrigation, in 5 to 6 applicationsIntegrated pest management guidance of the Greek Ministry of Rural Development and Food [27]
CottonAbout 700 millimetres in total, peaking around 10 millimetres a day at flowering and boll setGreek publication on rational irrigation [24]
Olive2,000 to 4,000 cubic metres per hectare per year, with a Mediterranean Kc of 0.65 to 0.75Teaching material, Agricultural University of Athens [26]
Kiwifruit, PieriaApplied irrigation of 599 to 1,350 millimetres per year, depending on plot and seasonCase study on two plots, 2020 and 2021 [23]
The kiwifruit figure is consumption, not requirement

The values of 599 to 1,350 millimetres are what was applied and not what the crop needed. The authors themselves conclude that in most cases the water applied was higher than the crop irrigation water needs, in certain months higher than double the needs [23]. It is one of the few published Greek datasets that directly show over-irrigation in a tree crop.

The table Kc does not fit a sparse planting

In eddy covariance measurements in a peach orchard in central Portugal, measured evapotranspiration was 1.2 to 3.6 millimetres a day with a mean of 2.5 from June to September, and the measured Kc came out at about 0.5 against the 0.9 of the table, that is, FAO-56 overestimated by about 35% [25]. The reason is not the species but the sparse, incompletely developed canopy.

The consequence for Greek fields is immediate: in a young planting, in a planting with low ground cover, or in a season with heavy summer pruning, the table Kc overestimates, and the irrigation computed with it is larger than the requirement.

The simpler methods, what they give in Greece

When a station does not measure solar radiation or wind, the Penman-Monteith equation cannot be fully applied and simpler methods are used. Two Greek comparative studies give a picture.

At two stations in Western Macedonia, nine methods were compared against FAO-56. The Penman-type methods gave the closest agreement, with a coefficient of determination above 0.97, while Hargreaves-Samani gave 0.923 at Krania and 0.854 at Kozani. All methods perform better in winter and diverge most in the warm season, that is, exactly when the accuracy is needed [19].

At three Greek cultivation sites, ERA5-Land reanalysis data underestimated in situ evapotranspiration by about 35%, while among the empirical methods Valiantzas overestimated by 32 to 62% and Copais by 24 to 56%. Hargreaves-Samani was found to be the most efficient of the simple methods [29]. In Peloponnese Augusts, the same method underestimated systematically with a root mean square difference of 0.72 to 0.76 millimetres a day [22].

What this means in practice

An error of 0.75 millimetres a day, over a thirty-day August, is 22.5 millimetres, that is 225 cubic metres per hectare. If the method underestimates, those cubic metres are missing from the programme and the field falls short without it showing. Before a station's evapotranspiration reading is used, it is worth establishing which method computes it, and whether the station actually measures solar radiation and wind or estimates them.

Sources: 1 19 20 21 22 23 24 25 26 27 29

10. Vapour pressure deficit

Air can hold a maximum amount of water vapour, and that maximum rises exponentially with temperature. The vapour pressure deficit, abbreviated VPD, is the difference between what the air can hold and what it already holds. It is, literally, how empty the air is.

10.1 The calculation

e°(T) = 0.6108 × exp[ 17.27 × T / (T + 237.3) ]

the saturation vapour pressure in kPa, with the temperature T in degrees Celsius [1].
VPD = e°(T) × (1 - RH/100)

the vapour pressure deficit in kPa, from the temperature and relative humidity of a single moment.

What matters in the first equation is the exp. Saturation vapour pressure does not rise in proportion to temperature, it rises exponentially. From 25 to 35 degrees, that is ten degrees, it nearly doubles. This is why a day at 38 degrees is not "a little worse" than one at 33, it is far worse.

Air temperatureRelative humiditySaturation vapour pressure e°(T)VPD
25 °C50%3.168 kPa1.584 kPa
30 °C40%4.243 kPa2.546 kPa
33 °C35%5.030 kPa3.270 kPa
35 °C30%5.623 kPa3.936 kPa
35 °C25%5.623 kPa4.217 kPa
38 °C20%6.625 kPa5.300 kPa
40 °C15%7.376 kPa6.269 kPa

Computed from the FAO-56 equations [1]. No published VPD measurements from a Greek field were found, in Greek or in English. The values in the table are computed and not measured, and are given so that the order of magnitude of a Greek summer afternoon is visible.

10.2 The error of using mean humidity

The usual way to compute a daily VPD is to take the mean temperature and the mean relative humidity of the day. It is wrong, and FAO-56 states so explicitly: using mean air temperature instead of daily minimum and maximum results in a lower mean saturation vapour pressure, a smaller deficit, and therefore an underestimation of reference evapotranspiration [1].

The correct way uses the extremes of the day.

es = [ e°(Tmax) + e°(Tmin) ] / 2
ea = [ e°(Tmin) · RHmax/100 + e°(Tmax) · RHmin/100 ] / 2

and the deficit is their difference [1].
DayCorrect VPD, FAO-56With mean humidityWith mean temperature and mean humidity
35 / 20 °C, humidity 80 / 25%2.342 kPa1.891 kPa, that is -19.3%1.744 kPa, that is -25.6%
38 / 22 °C, humidity 70 / 20%3.046 kPa2.549 kPa, that is -16.3%2.334 kPa, that is -23.4%
32 / 18 °C, humidity 85 / 35%1.700 kPa1.364 kPa, that is -19.8%1.267 kPa, that is -25.5%

Computed from the FAO-56 equations for three typical summer days of Northern Greece. The days are chosen as examples and do not come from a particular station.

The error is large and always in the same direction. The rough calculation underestimates the daily deficit by 23 to 26%. In the third row of the table, 1.70 kPa becomes 1.27 kPa, so the day moves from "above the point at which kiwifruit stomata are closing" to "comfortably below it". The decision changes because of the way the number was computed.

10.3 The leaf does not have the temperature of the air

The driving force of transpiration is not the deficit of the air but the deficit between the inside of the leaf, which is saturated at leaf temperature, and the air. A leaf in the sun is warmer than the air whenever the radiation it absorbs exceeds what transpiration and convection remove.

In measurements of leaf temperature in tree canopies, sunlit leaves were found to be 8 to 10 degrees warmer than the air at midday, with a maximum recorded difference of 12.06 degrees. The crossover sits near 25 degrees of air temperature: above it leaves run warmer, below it cooler [15].

Those measurements were made in forest canopies in Australia, in a temperate eucalypt woodland and in a tropical rainforest. No equivalent measurement in a Mediterranean orchard was found. They are carried over as the order of magnitude of the phenomenon and not as a value for a Greek tree.

The consequence is easy to compute. With air at 35 degrees and 25% relative humidity, the actual vapour pressure of the air is fixed at 1.406 kPa.

Leaf above airLeaf temperatureSaturation vapour pressure at the leafLeaf to air deficitAgainst the station VPD
035 °C5.623 kPa4.217 kPa1.00
+2 °C37 °C6.275 kPa4.869 kPa1.15
+5 °C40 °C7.376 kPa5.970 kPa1.42
+8 °C43 °C8.640 kPa7.234 kPa1.72
+10 °C45 °C9.582 kPa8.177 kPa1.94
The station figure is the lower bound

At the leaf to air temperature differences that have been measured in sunlit leaves, the deficit the leaf actually faces is 1.4 to 1.9 times the air deficit the station displays. There is also the reverse side, and it is useful: a leaf transpiring freely cools itself and its difference from the air is small, while a leaf that has closed its stomata heats up. The leaf to air temperature difference is, on its own, a measure of water stress.

10.4 Where the stomata close

There is no single number for all species. There are measurements by species, and they are instructive.

Kiwifruit

In gas exchange measurements on young vines, stomatal conductance declined steadily from 0.8 to 2.5 kPa, in both irrigated and droughted plants. There is therefore no safe plateau: the decline starts very low. Conductance reached half of its maximum at 2.1 to 2.2 kPa, in both treatments [11].

The work was done in Italy. The species is the same as in Greek plantings.

Olive

In an olive orchard at Córdoba, leaf conductance declined linearly with the vapour pressure deficit, with a slope that depends on water status. In well-watered trees with a potential above -1.65 MPa, the relationship was gs = 206.4 - 20.7 × VPD, that is a loss of about 21 units of conductance per kPa. In trees below -4 MPa there was no relationship left, because the stomata were already shut [12].

The work was done in Spain, in a Mediterranean olive orchard.

For field crops, a review collected the deficit above which transpiration stops rising linearly, that is where the stomata begin to restrict: durum wheat 1.05 to 1.10 kPa, soybean 1.10 to 2.19, sorghum 1.17 to 4.42, maize 1.30 to 4.14, pearl millet 1.45 to 4.35, chickpea 1.69 to 3.07, wheat 1.86 to 3.89 [13]. The ranges are wide because they depend on genotype and conditions, but the lower end sits consistently around 1 kPa.

For trees, a census of five species in Sicily gave the deficit each species experiences during the growth of its fruit: peach, mean 2.1 and maximum 3.84 kPa, mango 1.11 and 1.97, olive 1.1 and 2.15, orange 1.1 and 1.8, loquat 0.6 and 1.4 [14]. Peach therefore develops its fruit at a deficit that is permanently above the limits of the field crops and above the point at which kiwifruit has already closed half its stomata.

A very low deficit is a problem too

The discussion centres on high deficit, but there is an opposite end. In a review of greenhouse microclimate, a deficit below 0.3 kPa is associated with mineral deficiencies, soft growth, guttation and fungal disease, and specifically below 0.43 kPa with calcium deficiency, because transpiration is restricted and the calcium does not reach the fruit. Above 1.5 kPa wilting, leaf roll and stunted growth appear [28].

The link to calcium is worth keeping, because it explains why nights with a saturated atmosphere are not harmless for the fruit. The subject is developed in the article on calcium and postharvest physiology in sweet cherry.

The values concern a greenhouse and not an open field crop. The range of 0.45 to 1.25 kPa that circulates widely as the "greenhouse optimum" is not used here, because the review itself attributes it to a commercial weblog rather than to a peer-reviewed study [28].

10.5 The midday depression

The photosynthesis of a tree does not follow the sun. It peaks in the morning and falls towards midday, even when the soil holds plenty of water. In the Spanish olive orchard, assimilation fell to 6 to 11 units on days of high evaporative demand, against 15 units on a day of low demand, a depression of 27 to 60% [12]. In kiwifruit, the authors attribute the midday depression explicitly to the combination of the deficit arising inside the leaf from transpiration and the direct effect of the rising vapour pressure deficit [11].

What it means for irrigation

The midday depression is not corrected by water. It is a response to the atmosphere and not to the soil. What adequate irrigation achieves is that no second limitation is added from the soil side, so that the depression stays as shallow as possible and the leaf recovers quickly in the afternoon. A field that is adequately irrigated in a heatwave does not avoid the depression. It avoids the permanent damage.

Sources: 1 11 12 13 14 15 28

11. How the three quantities bind together

Comparison of the irrigation limit on a mild day and in a heatwave, for the same field
The same soil, the same crop, the same available water. What changes is how much of it may be depleted, and the atmosphere is what changes it. The adjustment is not an opinion, it is part of FAO-56 itself.

11.1 The same moisture, two different answers

The central finding of the recent literature is that the availability of soil water is not a property of the soil alone. It is defined jointly with the demand of the atmosphere. When the vapour pressure deficit is high, the transpiration rate required to satisfy the atmosphere creates a hydraulic gradient that the soil and root pathway cannot sustain. Stomata close at a soil water content that would have been perfectly adequate on a humid day [13].

The mechanism is hydraulic. As the rhizosphere dries, the hydraulic conductance of the soil and of the soil to root interface collapses, and that collapse interacts with the flow rate the deficit demands. It is not that the plant "senses" less water. It is that the water does not arrive in time.

Two measurements that show the two sides

The first side, the deficit limits even with enough water. In semi-arid ecosystems, within periods immediately after rain, that is with adequate soil moisture, high vapour pressure deficit alone reduced gross photosynthesis by 7 to 27% and surface conductance by 22 to 34% [18]. The measurements were made in Arizona, in grassland, savanna and shrubland.

The second side, the same deficit costs more when water is short. In kiwifruit, conductance fell to half at the same deficit, 2.1 to 2.2 kPa, in both irrigated and droughted vines. The droughted vines, however, started from a lower level and ended the day with 2.7 against 4.7 grams of net carbon gain per square metre [11].

11.2 The controlled experiment

The cleanest demonstration that the two act independently comes from a greenhouse, where one can be held constant while the other changes. Tomato and cucumber were tested at two levels of soil moisture, 100% and 60% of field capacity, crossed with two levels of vapour pressure deficit, high at 3 to 6 kPa and low at 1.5 to 2.5 kPa [16].

Three findings matter for the field.

  • At the same water shortage, the low deficit rescued the plant. At 60% of field capacity, the relative water content of the leaves was 11% higher at the low deficit than at the high one [16]. Same soil, different atmosphere, different plant water status.
  • The low deficit restored photosynthesis. The reduction in stomatal conductance caused by drought was largely reversed when the deficit was kept low [16].
  • Soil moisture remained the dominant factor for biomass. The fresh and dry weights of tomato at 100% of field capacity were 40.25% and 36.79% higher than at 60% [16]. The low deficit helps, but it does not replace water.

The experiment was carried out in a greenhouse, in pots, in China, with fogging systems controlling humidity. It does not transfer as a practice to open field crops, where the vapour pressure deficit is not controlled. It transfers as a mechanism, that is, as proof that the two factors are separate and additive.

11.3 Why the two coincide more often than chance would suggest

The two adverse conditions do not arrive independently. High atmospheric demand accelerates evapotranspiration, which reduces soil moisture, and the continuing loss of soil water in turn limits evapotranspiration, increasing sensible heat and making the near-surface air hotter and drier still [17].

In an analysis of 37 years over Chinese territory, the probability of high vapour pressure deficit and low soil moisture occurring together was found to be significantly greater than the product of the two separate probabilities, and the difference grew as the stress threshold became more severe. The coupling of the two was strongest in the semi-arid and sub-humid regions, that is, in the transitional zones, and weakest in the extremely arid and in the humid ones [17].

The study concerns China and rests on remote sensing and reanalysis data, with soil moisture representing the surface layer of 0 to 7 centimetres. The percentages and the regions do not transfer. What transfers is the structure of the phenomenon, and the fact that the Greek agricultural zone belongs climatically to exactly the transitional categories where the coupling was strongest.

11.4 The quantitative expression, inside FAO-56 itself

FAO-56 does not leave the matter in qualitative terms. The depletion fraction p is tabulated for an evapotranspiration of about 5 millimetres a day, and is explicitly adjusted for demand.

p = p(table) + 0.04 × (5 - ETc)

with the adjusted p always limited between 0.1 and 0.8, and ETc in millimetres per day [1].

That single line is the translation of this whole section into practice. The irrigation limit moves with the demand of the atmosphere, without anything changing in the soil.

Cropp from tableMild day, ETc 3 mmAverage day, ETc 5 mmHeatwave, ETc 8 mm
Stone fruit, cherry0.500.580.500.38
Kiwifruit0.350.430.350.23
Olive0.650.730.650.53
Cotton0.650.730.650.53
Maize0.550.630.550.43

The peach orchard example of section 5, with an available water of 144 millimetres, gives a readily available water of 83.5 millimetres on the mild day and 54.7 in the heatwave. In days the difference is larger still, because the heatwave empties a smaller reserve faster: twenty-eight days against seven.

The sentence that the whole chapter comes down to

A fixed irrigation limit, of the kind "irrigate when the sensor falls to 18%", is correct only for one particular level of atmospheric demand. In May it is too strict and leads to unnecessary applications. In a July heatwave it is too lax and lets the tree reach stress before any alarm sounds. The limit is a function, not a number.

For stone fruit there is also the reverse case, where controlled water stress is deliberate rather than a mistake. The subject is developed in the article on deficit irrigation in peach and nectarine, where the stress is applied at a specific phenological stage and for a specific purpose.

Sources: 1 11 13 16 17 18

12. The two calculators

The first converts the sensor reading into millimetres of water and into days. The second computes the vapour pressure deficit from temperature and relative humidity.

Calculator 1: from the percentage to millimetres and days

The field capacity and wilting point values are filled in automatically from the soil texture, and are replaced with measured values from the field wherever those exist.

QuantityValueWhat it means
Total available water, TAW-All the water between field capacity and wilting, for the stated depth
Adjusted depletion fraction p-From the table, adjusted for today's demand
Readily available water, RAW-What may be depleted without stress
Water present right now-Above the wilting point
Depletion from field capacity-What has already gone
Percent of available water remaining-The reading translated onto the scale of this field
Remaining until the irrigation limit-In millimetres and in days, at the stated evapotranspiration
Dose to refill to field capacity-Net dose, before the efficiency of the system
Fill in the fields.

The net dose is not the dose that will be applied. It is divided by the efficiency of the irrigation system, which for well-adjusted drip runs around 0.90 and lower for sprinkler. The calculation assumes that the sensor reading is volumetric water content. If it is not, the result does not hold, and section 2 explains why.

Calculator 2: vapour pressure deficit

From the temperature and relative humidity of a single moment. For the daily FAO-56 deficit the maximum and minimum values of the day are needed, as section 10.2 explains.

QuantityValueWhat it means
Saturation vapour pressure e°(T)-How much vapour the air can hold at this temperature
Actual vapour pressure ea-How much it already holds
Vapour pressure deficit of the air-What the station displays
Leaf to air deficit-What a sunlit leaf actually faces
Dew point-The temperature at which dew will form
Fill in the fields.

The reference limits in the comment come from the published breakpoint deficits of field crops and from the measured response of kiwifruit [11] [13]. They are indicative and concern those particular species and the particular conditions of those studies.

13. What the station provider needs to be asked

The questions that follow are not technical detail. Without their answers, the soil moisture reading cannot enter the calculations of this article.

The seven questions

  1. What the percentage actually is. Volumetric water content in m³ per m³, percent of field capacity, or a normalised index between the two ends of the available water.
  2. At which depths, in centimetres, the sensors sit, and whether the value displayed comes from one depth or is an average of several.
  3. What technology the sensor uses. Capacitance or FDR, TDR, granular matrix, or tensiometer. Whether the value is an amount or a tension depends on this.
  4. Whether a soil-specific calibration was carried out or the factory curve is used. If it was, by what method and when.
  5. Exactly where the sensor sits in relation to the emitter and to the trunk. Half a metre of difference changes the curve entirely under drip.
  6. By what method the displayed evapotranspiration is computed, and whether the station measures solar radiation and wind speed or estimates them from temperature.
  7. Whether the irrigation recommendation incorporates a crop coefficient, and which one. A recommendation that gives cubic metres without knowing the crop, the stage and the ground cover cannot be right except by accident.
What to do until the answers arrive

The reading is used as relative and not as absolute. The value two days after a deep irrigation is recorded, which is the upper end for that field, and the value at the end of a long dry spell, which is the lower one. The limit is then set as a percentage of the distance between them, and not from a textbook table. That is an empirical two-point calibration, it is free, and it solves most of the problem.

14. The points that remain open

Writing this article brought up gaps that are not covered by estimation and are stated as they are.

  • The unit, the depth and the calibration of the soil moisture sensors of the agrometeorological station networks operating in Greece are not published. Without them, the reading does not convert into millimetres.
  • No published reference evapotranspiration series was found for Pella, Imathia or Thessaloniki with a daily July value. The nearest one is Kozani, at 625 metres of elevation.
  • No Greek table of tension limits for starting irrigation by crop was found, neither in Ministry guidance nor in a published study. The values in circulation come from American sources.
  • No measured seasonal crop evapotranspiration was found for peach, cherry or kiwifruit in Greece, by lysimeter or eddy covariance. For kiwifruit there are applied irrigation amounts, which is not the same thing.
  • No published vapour pressure deficit measurements from a Greek field were found. The values given here are computed from the equations.
  • No measurement of the leaf to air temperature difference in a Mediterranean orchard was found. The values used come from forest canopies in Australia.
  • No measurement of the midday depression of photosynthesis in a Greek orchard was found. The nearest ones are from Spain and Italy.
The two gaps that close most easily

The first two need no research. The first closes with a written question to the station provider. The second closes with a station recording solar radiation, wind, temperature and humidity, and with the computation of FAO-56 reference evapotranspiration over two or three years. The rest need an experimental field.

15. Sources

The numbers in the text refer to the list below.

Balance and limits of the sources The list contains 30 sources: five Greek, of which two carry primary Greek evapotranspiration data, one is a doctoral thesis and one an official Ministry document; one international FAO technical handbook; five university extension fact sheets; two educational or commercial documents from manufacturers; and the remainder publications in peer-reviewed journals.

Three things are stated explicitly as non-Greek. The tension limits for starting irrigation come from American sources. The leaf temperature measurements come from forest canopies in Australia. The stomatal response measurements to the vapour pressure deficit come from Italy, Spain and Sicily. All three are carried over as mechanism and as order of magnitude, and not as values for a Greek field.

FAO-56 accounts for more than one third of the citations, and that is stated. The reason is that it is the computation standard itself and not one research paper among others. Wherever a measured Greek value exists, it takes precedence.
Disclaimer

The values of field capacity, permanent wilting point, crop coefficients and irrigation limits given here are indicative and come from the international literature. They are not an irrigation recommendation for a particular field. The final dose and frequency are determined after soil analysis, measurement in that particular field, and on-site assessment by a licensed agronomist.

The two calculators in this article are tools for understanding and assume that the data entered are correct. The first in particular assumes that the sensor reading is volumetric water content. If the reading is something else, or if the sensor has not been calibrated for that soil, the result does not hold.

Plant protection products are used solely on the basis of the approvals of the Greek Ministry of Rural Development and Food and of the product label. Those approvals are valid in Greece only and do not apply to other countries. Foliar spraying over trees carrying fruit can cause marking on the fruit and scorching of the foliage, and a prior test on a limited number of trees is required, with the same product and the same rate, and a wait of 5 to 7 days before general application. Spraying is done in the late afternoon or early morning, never in extreme heat and never on foliage under water stress.

Every intervention is made at the user's responsibility, following on-site assessment by a licensed agronomist. Pastopoulos Agronomics accepts no liability for any damage arising from the use of the information in this article.

Related articles

Growing avocado in Greece: cultivars, soil, climate, diseases and pests

Growing avocado in Greece: cultivars, soil, climate, diseases and pests

Three botanical races, and the Mexican one supplies frost tolerance while the West Indian one supplies tolerance of lime and salt: Greece asks for both at once. With Greek data from ELGO Chania: across ten orchards over six years Fuerte met the dry matter threshold in only two seasons and Zutano in none. In detail: synchronous dichogamy and the cultivars one by one, the temperature limits, soils at pH up to 8.3 and calcium carbonate up to 20.7 per cent, eleven rootstocks, planting distances and training.

Attapulgite: what it actually does to soil, at what rate, and what its use costs

Attapulgite: what it actually does to soil, at what rate, and what its use costs

A Greek mineral with a specific surface area of 220 m2/g and an exchange capacity of 57 meq/100 g. What the field trials on sandy soils showed, why the rate threshold is 6,000 kg per hectare, what the IARC fibre classification means, which particle grade for which use, and why the rise in pH is the main argument against it on Mediterranean calcareous soils.