METHOD FOR DETERMINING AT LEAST ONE METEOROLOGICAL VARIABLE FOR DESCRIBING A STATE OF ATMOSPHERIC WATER

Abstract

The invention relates to a method of determining at least one meteorological quantity (iwvc, ilwc, rr) for describing a state of atmospheric water, in particular of water vapor (wv), condensed water (lw) and/or precipitation (r), comprising the following steps: (a) providing meteorological input data (d1), (b) step in whichfor at least one state (wv, lw) of the atmospheric watersecond data (d2) are calculated from the meteorological input data (d1), said second data (d2) representing a measure of the signal attenuation of signals transmitted though the atmosphere (A), said signal attenuation being caused by said at least one state (wv, lw), (c) providing third data (d3), which represent a measure of the signal attenuation of signals caused by atmospheric water, wherein the third data (d3) are obtained from the measurement of signals (3) which have been transmitted through the atmosphere (A) between signal transmission devices (1,2), in particular between signal transmission devices on satellites (1) and terrestrial signal transmission devices (2), (d) comparing the second data (d2) with the third data (d3), and (e) step in whichfor at least one further state (lw, r) of the atmospheric water, which preferably has not been taken into account in the calculation of step (b)a meteorological quantity (ilwc, rr) for describing said at least one further state (lw, r) is calculated in dependence of the deviation between the second data (d2) and the third data (d3).

Claims

1. A method for determining at least one meteorological quantity for describing a state of atmospheric water, in particular at least one of water vapor, condensed water or precipitation, comprising the following steps: (a) providing meteorological input data, (b) step in whichfor at least one state of the atmospheric watersecond data is calculated from the meteorological input data, said second data representing a measure of the signal attenuation of signals transmitted though the atmosphere (A), said signal attenuation being caused by said at least one state, (c) providing third data, which represent a measure of the signal attenuation of signals caused by atmospheric water, wherein the third data is obtained from the measurement of signals which have been transmitted through the atmosphere between signal transmission devices, in particular between signal transmission devices on satellites and terrestrial signal transmission devices, or between terrestrial signal transmission devices, or both, (d) comparing the second data with the third data, and (e) step in whichfor at least one further state of the atmospheric watera meteorological quantity for describing said at least one further state is calculated in dependence of the deviation between the second data and the third data.

2. The method of claim 1, wherein the calculation of the second data in step (b) is carried out for the state of water vapor, and the second data represents a measure of the signal attenuation caused by water vapor.

3. The method of claim 1, wherein the state for which a meteorological quantity is calculated in step (e) is condensed water in clouds and/or fog.

4. The method according to claim 1, wherein the calculation of the second data in step (b) is carried out for the states of water vapor and condensed water in clouds and/or fog, and the second data represents a measure of the signal attenuation caused by water vapor and condensed water in clouds and/or fog.

5. The method according to claim 1, wherein the state for which a meteorological quantity is calculated in step (e) is precipitation, wherein preferably the meteorological quantity calculated in step (e) is the precipitation rate, in particular the rain rate.

6. The method according to claim 1, wherein the meteorological input data at least comprises values for the air temperature, the air pressure and/or the humidity at a particular level, preferably at the earth's surface.

7. The method according to claim 1, wherein comparing the second data with the third data according to step (d) is carried out by calculating a difference or a ratio between the signal attenuation according to the second data and the signal attenuation according to the third data.

8. The method according to claim 1, wherein the value 0 is assigned to the meteorological quantity according to step (e) if the signal attenuation according to the third data is less than the signal attenuation according to the second data.

9. The method according to claim 1, wherein step (b) further comprises: (b1) for the state of water vapor, water vapor related data is calculated based on the meteorological input data, said water vapor related data representing a measure of the signal attenuation of signals transmitted through the atmosphere, said signal attenuation being caused by water vapor, and (b2) said water vapor related data is compared with the third data, and (b3) if the signal attenuation according to the third data is greater than the signal attenuation according to the water vapor related data, condensed water related data is calculated, which represents a measure of the signal attenuation caused by condensed water in clouds and/or fog, wherein the second data represents a measure of the signal attenuation of signals caused by water vapor and condensed water in clouds and/or fog.

10. The method according to claim 1, wherein the calculation according to step (b) is carried out by means of an atmosphere model which assignsin dependence of the meteorological input data at least one value of at least one meteorological quantity for describing a state of atmospheric water, preferably the water vapor density and/or the condensed water content, to each height segment and/or each volume unit and/or each point or grid point of a model atmosphere.

11. The method according to claim 10, wherein a signal attenuation is calculated from each of the values contained in the model atmosphere, wherein the second data is calculated by integration or summation along, preferably vertical, signal paths (P) through the model atmosphere.

12. The method according to claim 1, wherein a cloud upper boundary is taken into account in the calculation of condensed water in clouds, wherein preferably the cloud upper boundary is derived from satellite images.

13. The method according to claim 1, wherein the calculation of the rain rate is carried out according to the relation $\mathrm{rr}=\sqrt[]{\frac{}{k}},$ with rr representing the rain rate, representing the difference between the signal attenuation according to the second data and the signal attenuation according to the third data, and and k representing parameters, wherein preferably the parameters and k have been selected in accordance with Recommendation ITU-R P.838-3.

14. The method according to claim 1, wherein the calculations of at least one of steps (b), (d), or (e) are performed according to one or more of the following: in dependence of the signal frequency, in dependence of the signal polarization, in dependence of the elevation of the signal propagation direction, for one or several signal frequency/ies, wherein preferably the signal frequency/ies lie between 5 GHz and 100 GHz, in particular between 10 GHz and 50 GHz, or only below a specified height boundary, preferably below 20,000 m, particularly preferably below 10,000 m, even more preferably below 7,500 m.

15-16. (canceled)

17. The method according to claim 1, wherein the at least one further state of the atmospheric water, for which in step (e) a meteorological quantity for describing said at least one further state is calculated, has not been taken into account in the calculation of step (b).

18. The method of claim 1, wherein the meteorological quantity calculated in step (e) is the condensed water content in clouds and/or fog.

19. The method of claim 1, wherein the meteorological quantity calculated in step (e) is the integrated condensed water content in clouds and/or fog.

20. The method of claim 1, wherein calculating the meteorological quantity according to step (e) is carried out from the difference and/or the ratio between the signal attenuation according to the second data and the signal attenuation according to the third data.

21. A data processing system operable to perform the method of claim 1.

22. A computer program product stored on a non-transitory computer-readable medium, including instructions for performing the method of claim 1.

Description

[0055] To provide better understanding of the invention, the invention is explained in more detail with the aid of the following figures.

[0056] The following is shown in a strongly simplified, schematic representation:

[0057] FIG. 1 an atmosphere which only contains water vapor;

[0058] FIG. 2 an atmosphere which only contains water vapor and condensed water (cloud water);

[0059] FIG. 3 an atmosphere which contains water vapor, condensed water (cloud water) and precipitation (rain);

[0060] FIG. 4 a schematic representation of the signal transmission through the atmosphere between signal transmission devices on satellites and terrestrial signal transmission devices;

[0061] FIG. 5 a flow diagram for describing one embodiment of the invention;

[0062] FIG. 6 a flow diagram for describing one embodiment of the invention;

[0063] FIG. 7 a schematic representation of the signal attenuations caused by the individual states in the total attenuation;

[0064] FIG. 8 the individual contributions to the signal attenuation.

[0065] As an introduction, it should be noted that in the embodiments described in different ways, identical parts or method steps are indicated with identical reference numbers or identical component designations; at the same time, the disclosures contained in the entire description may be analogously applied to identical parts with identical reference numbers or identical component designations. Moreover, the position indications chosen in the description, such as at the top, at the bottom, laterally, etc. refer to the figure which is directly represented and described; and if a position changes, said position indications are to be applied analogously to the new position.

[0066] The embodiments show possible variants; however, it should be noted at this point that the invention is not limited to the variants specifically shown; rather, various combinations of the individual variants are possible as well, and, due to the technical information provided by the present invention, this variation possibility is subject to the skills of the person skilled in the art who works in this technical field.

[0067] The scope of protection is determined by the claims. However, the description and the drawings are to be used for construing the claims. Individual features or feature combinations from the different embodiments that are shown and described may per se constitute independent solutions according to the invention. The object underlying the independent solutions according to the invention may be gathered from the description.

[0068] Any indications of value ranges in the present description are to be understood as including any and all subranges thereof; for example, the indication 1 to 10 is to be understood as including all subranges, starting out from the lower boundary 1 to the upper boundary 10; i.e. all subranges start with a lower boundary of 1 or greater and stop at an upper boundary of 10 or less, e.g. 1 to 1.7 or 3.2 to 8.1 or 5.5 to 10.

[0069] For the sake of good order, it should finally be noted that, for better understanding, some of the facts shown in the figures have been represented unscaled and/or enlarged and/or in reduced size.

[0070] FIGS. 1-3 show possible weather constellations that may be described with the method according to the invention. In FIG. 1, a part of an atmosphere A, the troposphere, is shown, where water is only present in the (invisible) state of water vapor wv. E refers to the earth's surface. Preferably, the calculations underlying the invention are only performed up to a specified upper boundary (uppermost dashed line). The boundary may be set in the range between 6,000 hm (from sea level) and the tropopause (which depends on the latitude), for example. This way, on the one hand, the range relevant for weather patterns is taken into consideration, and, on the other hand, the required arithmetic operations are kept at a manageable level. FIG. 2 shows a weather situation where clouds C and thus condensed water lw are present in addition to water vapor wv. The clouds C are located between a lower boundary and an upper boundary CT (cloud top). In FIG. 3, water is also present in the state of precipitation r (rain).

[0071] When signals 3 are transmitted along signal paths P between data transmission devices 1, 2 through the atmosphere A, the present states of atmospheric water cause a signal attenuation, which may be captured and provided in decibel, for example. The different elevations of the signal paths P (see FIG. 4) may be converted, so that the edited data, which represent a measure of the signal attenuation, only refer to one elevation, preferably to the vertical direction.

[0072] By means of FIG. 5, an embodiment of the method according to the invention is now described in detail. The object is to determine a meteorological quantity for describing a state of atmospheric water. As already mentioned, water is particularly present in the states of water vapor wv, condensed water lw (in clouds or fog) and/or precipitation r. Examples of meteorological quantities for such states include: water vapor density (water vapor content), integrated water vapor density iwvc, condensed water content (liquid water content) lwc, integrated condensed water content ilwc, precipitation rate (rain rate) rr.

[0073] The method now comprises the following steps:

[0074] (a) providing meteorological input data (d1). Said data preferably comprise values for the air temperature, the air pressure and/or the humidity at a particular level, preferably at the earth's surface E.

[0075] (b) a step in which, for the state(s) of water vapor wv and/or condensed water lw, second data d2 are calculated from the meteorological input data d1, said second data representing a measure of the signal attenuation of signals transmitted through the atmosphere A, said signal attenuation being caused by water vapor wv and/or condensed water lw. In this process, for example, a quantity describing the statee.g. the water vapor density wvc (e.g. in the unit kg/m.sup.3) and/or the condensed water content lwc (e.g. in the unit kg/m.sup.3)are calculated. From this quantity, the signal attenuation may then be derived. Further quantities, such as the integrated water vapor density iwvc (in the unit kg/m.sup.2) and/or the integrated condensed water content ilwc (in the unit kg/m.sup.2) may then be calculated as a by-product.

[0076] (c) providing third data d3, which represent a measure of the signal attenuation of signals caused by atmospheric water, wherein the third data d3 have been obtained from the measurement of signals 3 transmitted through the atmosphere A between signal transmission devices 1, 2, in particular between signal transmission devices on satellites 1 and terrestrial signal transmission devices 2 (see FIG. 4); additionally or alternatively, measurements of signals that were or are transmitted through the atmosphere A from terrestrial signal transmission devices to other terrestrial signal transmission devices may also be used for the third data d3.

[0077] (d) comparing the second data d2 with the third data d3. The data d2 and d3 may, for example, include the signal attenuation in the unit dB. Comparing the second data d2 with the third data d3 may be carried out by calculating a difference or a ratio between the signal attenuation according to the second data d2 and the signal attenuation according to the third data d3.

[0078] (e) step in which the precipitation rate rr is calculated for the state of precipitation r (said state has not been taken into account in the calculation according to step (b)) in dependence of the deviation between the second data d2 and the third data d3. The precipitation rate rr may be calculated from the difference and/or the ratio between the signal attenuation according to the second data d2 and the signal attenuation according to the third data d3. This means that the precipitation rate rr (in particular the rain rate) may, for example, be represented as a function of this difference: rr=f (d3d2). The value 0 is assigned to the precipitation rate if the signal attenuation according to the third data d3 is less than the signal attenuation according to the second data d2.

[0079] By means of this method, it is ensured that the signal attenuation is not only assigned to the state calculated in step (e) (here: precipitation r), but that contributions of other states (here: water vapor wv and/or condensed water lw) as such are already taken into account beforehand.

[0080] If two states are already taken into account in step (b), step (b) may be carried out according to the preferred embodiment of FIG. 6. The remaining steps may then be carried out in the same way as FIG. 5, as described above.

[0081] According to FIG. 6, the following steps are carried out:

[0082] (b1) for the state of water vapor wv water vapor related data dwv are calculated based on meteorological input data d1, said water vapor related data representing a measure of the signal attenuation of signals transmitted through the atmosphere A, said signal attenuation being caused by water vapor (vw).

[0083] (b2) said water vapor related data dwv are compared with the third data d3. The third data d3 may be the same data that were already used in the method of FIG. 5.

[0084] (b3) if the signal attenuation according to the third data d3 is greater than the signal attenuation according to the water vapor related data dwv, condensed water related data dlw are calculated, which represent a measure of the signal attenuation caused by condensed water lw in clouds C and/or fog.

[0085] The second data d2 now represent a measure of the signal attenuation caused by water vapor wv and condensed water lw in clouds C and/or fog. In other words, the second data are the sum of the water vapor related data dwv and the condensed water related data dlw.

[0086] The second data d2 may then, as shown in FIG. 5, be used further. FIG. 7 shows another embodiment in which the calculation of the second data d2 in step (b) is (only) carried out for the state of water vapor wv. The second data d2 then represent a measure of the signal attenuation caused by water vapor wv. The state for which a meteorological quantity is calculated in step (e) is condensed water lw in clouds C and/or fog. In this regard, the meteorological quantity may be the condensed water content lwc in clouds C and/or fog, preferably the integrated condensed water content ilwc. Such a method is suitable, for example, if it is clear from other observations that no precipitation is present in the evaluated area.

[0087] FIG. 8 illustrates the contributions of the individual states of atmospheric water to the signal attenuation. The principle described above is based on the assumption that the contribution of water vapor dwv, the contribution of condensed water dlw and the contribution of precipitation dr=d3d2 together result in the signal attenuation d3 obtained from measurements. It becomes apparent from this presentation that the method according to the invention has great advantages as not only the presence of a state may be determined, but, due to the gradual assignment with respect to the signal attenuation, precise quantitative statements on a state may be made as well.

[0088] Below, preferred embodiments are described with reference to the individual steps.

[0089] The calculation according to step (b) may be carried out by means of an atmosphere model which, in dependence of the meteorological input data d1, assigns at least one value of at least one meteorological quantity for describing a state wv, lw of atmospheric water, preferably the water vapor density and/or the condensed water content lwc, to each height segment and/or each volume unit and/or each point of a model atmosphere.

[0090] A variety of different atmosphere models are known. They are of fundamentally different types and may, for example, represent adiabatic (dry-adiabatic or moist-adiabatic) layering or inversion (a partial increase in temperature with the height in a particular layer). By means of mathematical relations, based on data obtained from measurements in the area of the earth's surface, such as air temperature, air pressure and humidity, a three-dimensional model may be created. A (linear) temperature gradient (e.g. in the unit K/m) may be assumed for the heightdependent temperature.

[0091] Calculating the water vapor density may be carried out by calculating the saturated vapor pressure at the given temperature and deriving the current vapor pressure from the saturated vapor pressure and the relative humidity. For converting the vapor pressure into water vapor density (e.g. in the unit kg/m.sup.3), the ideal gas equation of thermodynamics may be used.

[0092] Calculating the water vapor density is carried out, for example, by applying the August Roche Magnus formula for calculating saturated vapor pressure from the temperature, calculating the current vapor pressure from the relative humidity and the saturated vapor pressure, and applying the ideal gas equation for converting the vapor pressure into water vapor density.

[0093] In this case, the calculation is made up of the following formulas:

[00002]$e_{\mathrm{sat}}=6.Math.1.Math.0.9.Math.4.Math.\exp \left(\frac{1.Math.7.Math.\mathrm{.25}.Math.t}{t+243.04}\right).Math..Math.e=r_H.Math.e_{s.Math.a.Math.t}$$=\frac{e}{R_V.Math.T},$

[0094] with e designating the current vapor pressure in Pa, e.sub.sat designating the saturated vapor pressure in Pa, t designating the temperature in C., T designating the temperature in K, r.sub.H designating the relative humidity in %, designating the water vapor density in kg/m.sup.3, and R.sub.V=461.5 J/(kg K) designating the individual gas constant of water. In the method according to the invention, the water vapor density (which is also referred to as wvc in the present application) designated with in the last formula may also be used for calculating the second data d2.

[0095] If the condensed water content lwc is to be calculated according to FIG. 6 as well, additional assumptions may be made, such as the cloud lower boundary and the cloud upper boundary, or they may be derived from additional measurements or (satellite) images. For example, the calculation according to step (b3) may be carried out in such a way that wherever (within the cloud lower boundary and the cloud upper boundary) the relative humidity is 1 or greater than a specified threshold value (e.g. >0.99), a (constant) value of the condensed water content is assumed. Preferably, in the cloud layer concerned, a (constant) value for the condensed water content between 0.05 g/m.sup.3 and 0.5 g/m.sup.3, preferably approximately 0.1 g/m.sup.3 is assumed. The signal attenuation may now in turn be calculated from the values that have been (analytically) calculated or assumed (based on empirical values) for the condensed water content.

[0096] This is preferably carried out by calculating the second data d2 by integration or summation along preferably vertical signal paths P (slant paths) through the model atmosphere.

[0097] The calculation of the rain rate rr is preferably carried out according to the relation

[00003]$\mathrm{rr}=\sqrt[]{\frac{}{k}},$

[0098] with representing the difference between the signal attenuation according to the second data d2 and the signal attenuation according to the third data d3, and and k representing parameters. The latter have preferably been selected in accordance with Recommendation ITU-R P.838-3.

[0099] The calculations described above may be carried out in dependence of the signal frequency and/or signal polarization and/or in dependence of the elevation of the signal propagation direction. They may also be carried out for one or several signal frequency/ies, with the signal frequency/ies preferably lying between 5 GHz and 100 GHz, particularly between 10 GHz and 50 GHz.

[0100] As already mentioned, the calculations may be limited to areas below a specified height boundary, preferably below 20.000 m, preferably below 10.000 m, particularly preferably below 7.500 m (above sea level). Usually, the tropopause (which depends on the latitude) forms an uppermost boundary for the above calculations.

LIST OF REFERENCE NUMBERS

[0101] 1 signal transmission device on a satellite [0102] 2 terrestrial signal transmission device [0103] 3 signal [0104] A atmosphere [0105] C clouds [0106] CT cloud upper boundary [0107] E earth's surface [0108] P signal path [0109] wv water vapor [0110] wvc water vapor density [0111] iwvc integrated water vapor density [0112] lw condensed water in clouds and/or fog [0113] lwc condensed water content [0114] ilwc integrated condensed water content [0115] r precipitation [0116] rr rain rate [0117] d1 meteorological input data [0118] d2 second data [0119] d3 third data [0120] dwv water vapor related data [0121] dlw condensed water related data