Parsivel2

The OTT Parsivel2 is a laser-based disdrometer for comprehensive measurement of all types of precipitation. There are currently 3 such disdrometers in operation at the BCO. The longest running disdrometer (451995) was in operation since 2023-05-14, was decommissioned 2025-10-14 due to hardware issues, and was reinstalled 2026-01-22. The other two, 451996 and 451997, were added 2025-09-22. Below we provide the locations of the three parsivels accompanied by their serial numbers.

The device works on the extinction principle and measures precipitation particles using the shadowing effects they cause when they pass through a laser band.

Parsivels 451995 and 451997 (1 and 3) are oriented facing east, while 451996 (2) is facing north.

The Parsivel2 captures both the size and the rate of fall in detail of the individual hydrometeors and classifies them into a range of 32 classes each. The temporal resolution covers between 10 seconds and one hour. The raw data is used to calculate for example precipitation type, amount and intensity.

The measurement principle

OTT Parsivel2 uses a laser-optical sensor to measure precipitation. The transmitter unit of the sensor generates a flat, horizontal beam of light that the receiver unit converts to an electrical signal. This signal changes whenever a hydrometeor falls through the beam anywhere within the measuring area (54 cm2). The degree of dimming is a measure of the size of the hydrometeor and the fall velocity is derived from the extinction signal duration.

Uh oh!

Comparison with other precipitation measurements corroborates findings from the literature that the parsivel does not accurately capture rain events surpassing 20 mm h-1 in intensity.

Data Availability

The data from the first parsivel is available as BCO.disdrometer_PARSIVEL451995_c1 in the catalog.

The next two are published under the keys BCO.disdrometer_PARSIVEL451996_c1_v2 and BCO.disdrometer_PARSIVEL451997_c1_v2.

Versions

The dataset from the first parsivel is available in multiple versions.

• Version 1: first version of the dataset

• Version 2: firmware revision is time-dependent, firmware was upgraded during the measurement period on 2025-08-31 (before installation of second and third parsivels.)

Sample Plot

The plot shows rain rate and particle size distribution during September 2024 (the period where the ORCESTRA Campaign visited the BCO).

import intake
import matplotlib.pylab as plt
cat = intake.open_catalog("https://tcodata.mpimet.mpg.de/catalog.yaml")
ds = cat.BCO.disdrometer_PARSIVEL451995_c1.to_dask()
subset = ds.sel(time="2024-09")
fix, (ax1, ax2) = plt.subplots(1, 2, figsize=(12, 4))
subset["rain_rate"].plot(ax=ax1)
subset["raw"].mean("time").plot(cmap="Blues", ax=ax2);

/builds/tco/bco/docs/.venv/lib/python3.12/site-packages/intake_xarray/base.py:21: FutureWarning: The return type of `Dataset.dims` will be changed to return a set of dimension names in future, in order to be more consistent with `DataArray.dims`. To access a mapping from dimension names to lengths, please use `Dataset.sizes`.
  'dims': dict(self._ds.dims),

The full dataset:

ds

<xarray.Dataset> Size: 3GB
Dimensions:              (d: 32, bnd: 2, time: 1288323, v: 32)
Coordinates:
    alt                  float64 8B ...
  * d                    (d) float32 128B 0.0625 0.1875 0.3125 ... 21.5 24.5
    lat                  float64 8B ...
    lon                  float64 8B ...
  * time                 (time) datetime64[ns] 10MB 2023-05-14T00:46:28 ... 2...
  * v                    (v) float32 128B 0.05 0.15 0.25 0.35 ... 15.2 17.6 20.8
Dimensions without coordinates: bnd
Data variables:
    d_bounds             (d, bnd) float32 256B dask.array<chunksize=(32, 2), meta=np.ndarray>
    error_code           (time) int16 3MB dask.array<chunksize=(262144,), meta=np.ndarray>
    field_n_ved          (time, d) float32 165MB dask.array<chunksize=(8192, 32), meta=np.ndarray>
    field_v              (time, v) float32 165MB dask.array<chunksize=(8192, 32), meta=np.ndarray>
    firmware_iop         (time) uint16 3MB dask.array<chunksize=(262144,), meta=np.ndarray>
    number_of_particles  (time) int32 5MB dask.array<chunksize=(262144,), meta=np.ndarray>
    rain_rate            (time) float32 5MB dask.array<chunksize=(262144,), meta=np.ndarray>
    raw                  (time, v, d) int16 3GB dask.array<chunksize=(256, 32, 32), meta=np.ndarray>
    sample_interval      (time) int16 3MB dask.array<chunksize=(262144,), meta=np.ndarray>
    sensor_status        (time) int16 3MB dask.array<chunksize=(262144,), meta=np.ndarray>
    v_bounds             (v, bnd) float32 256B dask.array<chunksize=(32, 2), meta=np.ndarray>
Attributes: (12/15)
    Conventions:           CF-1.12
    _logical_cutoff_date:  2026-03-07T00:00:00Z
    bcoproc_version:       0.0.0.post1185.dev0+c830e01
    featureType:           timeSeries
    history:               processed with bcoproc l1
    institution:           Max Planck Institute for Meteorology, Hamburg
    ...                    ...
    sensor_serial_number:  451995
    source:                OTT Parsivel2 451995
    station_name:          BCODIS
    summary:               The OTT Parsivel² is a modern laser disdrometer fo...
    title:                 Parsivel 2 Disdrometer data from BCO (Level 1)
    tool_versions:         {"Python": "3.11.2 (main, Apr 28 2025, 14:11:48) [...

xarray.Dataset

• Dimensions:
  • d: 32
  • bnd: 2
  • time: 1288323
  • v: 32
• Coordinates: (6)
  • alt
    ()
    float64
    ...

    long_name :

    height above mean sea level

    standard_name :

    altitude

    units :

    m

    [1 values with dtype=float64]

  • d
    (d)
    float32
    0.0625 0.1875 0.3125 ... 21.5 24.5

    bounds :

    d_bounds

    long_name :

    nominal centers of particle diameter bins

    units :

    mm

    array([ 0.0625,  0.1875,  0.3125,  0.4375,  0.5625,  0.6875,  0.8125,  0.9375,
            1.0625,  1.1875,  1.375 ,  1.625 ,  1.875 ,  2.125 ,  2.375 ,  2.75  ,
            3.25  ,  3.75  ,  4.25  ,  4.75  ,  5.5   ,  6.5   ,  7.5   ,  8.5   ,
            9.5   , 11.    , 13.    , 15.    , 17.    , 19.    , 21.5   , 24.5   ],
          dtype=float32)

  • lat
    ()
    float64
    ...

    long_name :

    north latitude

    units :

    degree_north

    valid_max :

    90.0

    valid_min :

    -90.0

    [1 values with dtype=float64]

  • lon
    ()
    float64
    ...

    long_name :

    east longitude

    units :

    degree_east

    valid_max :

    180.0

    valid_min :

    -180.0

    [1 values with dtype=float64]

  • time
    (time)
    datetime64[ns]
    2023-05-14T00:46:28 ... 2026-03-...

    axis :

    T

    standard_name :

    time

    units_metadata :

    leap_seconds: utc

    array(['2023-05-14T00:46:28.000000000', '2023-05-14T00:47:28.000000000',
           '2023-05-14T00:48:28.000000000', ..., '2026-03-06T23:57:59.000000000',
           '2026-03-06T23:58:59.000000000', '2026-03-06T23:59:59.000000000'],
          shape=(1288323,), dtype='datetime64[ns]')

  • v
    (v)
    float32
    0.05 0.15 0.25 ... 15.2 17.6 20.8

    bounds :

    v_bounds

    long_name :

    nominal centers of particle velocity bins

    units :

    m s-1

    array([ 0.05,  0.15,  0.25,  0.35,  0.45,  0.55,  0.65,  0.75,  0.85,  0.95,
            1.1 ,  1.3 ,  1.5 ,  1.7 ,  1.9 ,  2.2 ,  2.6 ,  3.  ,  3.4 ,  3.8 ,
            4.4 ,  5.2 ,  6.  ,  6.8 ,  7.6 ,  8.8 , 10.4 , 12.  , 13.6 , 15.2 ,
           17.6 , 20.8 ], dtype=float32)

• Data variables: (11)
  • d_bounds
    (d, bnd)
    float32
    dask.array<chunksize=(32, 2), meta=np.ndarray>

    long_name :

    bounds of particle diameter bins

    Array Chunk Bytes 256 B 256 B Shape (32, 2) (32, 2) Dask graph 1 chunks in 2 graph layers Data type float32 numpy.ndarray

  • error_code
    (time)
    int16
    dask.array<chunksize=(262144,), meta=np.ndarray>

    long_name :

    Error Code

    Array Chunk Bytes 2.46 MiB 512.00 kiB Shape (1288323,) (262144,) Dask graph 5 chunks in 2 graph layers Data type int16 numpy.ndarray

  • field_n_ved
    (time, d)
    float32
    dask.array<chunksize=(8192, 32), meta=np.ndarray>

    long_name :

    mean volume equivalent diameter per class

    units :

    0.1 lg(re 1 m3 mm)

    Array Chunk Bytes 157.27 MiB 1.00 MiB Shape (1288323, 32) (8192, 32) Dask graph 158 chunks in 2 graph layers Data type float32 numpy.ndarray

  • field_v
    (time, v)
    float32
    dask.array<chunksize=(8192, 32), meta=np.ndarray>

    long_name :

    mean velocity per class

    units :

    mm s-1

    Array Chunk Bytes 157.27 MiB 1.00 MiB Shape (1288323, 32) (8192, 32) Dask graph 158 chunks in 2 graph layers Data type float32 numpy.ndarray

  • firmware_iop
    (time)
    uint16
    dask.array<chunksize=(262144,), meta=np.ndarray>

    flag_meanings :

    2.11.6 2.11.11

    flag_values :

    [1, 2]

    Array Chunk Bytes 2.46 MiB 512.00 kiB Shape (1288323,) (262144,) Dask graph 5 chunks in 2 graph layers Data type uint16 numpy.ndarray

  • number_of_particles
    (time)
    int32
    dask.array<chunksize=(262144,), meta=np.ndarray>

    long_name :

    Number of particles

    units :

    1

    Array Chunk Bytes 4.91 MiB 1.00 MiB Shape (1288323,) (262144,) Dask graph 5 chunks in 2 graph layers Data type int32 numpy.ndarray

  • rain_rate
    (time)
    float32
    dask.array<chunksize=(262144,), meta=np.ndarray>

    long_name :

    Rain Intensity

    standard_name :

    rainfall_rate

    units :

    mm h-1

    Array Chunk Bytes 4.91 MiB 1.00 MiB Shape (1288323,) (262144,) Dask graph 5 chunks in 2 graph layers Data type float32 numpy.ndarray

  • raw
    (time, v, d)
    int16
    dask.array<chunksize=(256, 32, 32), meta=np.ndarray>

    long_name :

    particle count per class

    Array Chunk Bytes 2.46 GiB 512.00 kiB Shape (1288323, 32, 32) (256, 32, 32) Dask graph 5033 chunks in 2 graph layers Data type int16 numpy.ndarray

  • sample_interval
    (time)
    int16
    dask.array<chunksize=(262144,), meta=np.ndarray>

    long_name :

    Sample Interval

    units :

    s

    Array Chunk Bytes 2.46 MiB 512.00 kiB Shape (1288323,) (262144,) Dask graph 5 chunks in 2 graph layers Data type int16 numpy.ndarray

  • sensor_status
    (time)
    int16
    dask.array<chunksize=(262144,), meta=np.ndarray>

    long_name :

    Sensor status

    Array Chunk Bytes 2.46 MiB 512.00 kiB Shape (1288323,) (262144,) Dask graph 5 chunks in 2 graph layers Data type int16 numpy.ndarray

  • v_bounds
    (v, bnd)
    float32
    dask.array<chunksize=(32, 2), meta=np.ndarray>

    long_name :

    bounds of particle velocity bins

    Array Chunk Bytes 256 B 256 B Shape (32, 2) (32, 2) Dask graph 1 chunks in 2 graph layers Data type float32 numpy.ndarray

• Indexes: (3)
  • d
    PandasIndex

    PandasIndex(Index([0.0625, 0.1875, 0.3125, 0.4375, 0.5625, 0.6875, 0.8125, 0.9375, 1.0625,
           1.1875,  1.375,  1.625,  1.875,  2.125,  2.375,   2.75,   3.25,   3.75,
             4.25,   4.75,    5.5,    6.5,    7.5,    8.5,    9.5,   11.0,   13.0,
             15.0,   17.0,   19.0,   21.5,   24.5],
          dtype='float32', name='d'))

  • time
    PandasIndex

    PandasIndex(DatetimeIndex(['2023-05-14 00:46:28', '2023-05-14 00:47:28',
                   '2023-05-14 00:48:28', '2023-05-14 00:49:28',
                   '2023-05-14 00:50:28', '2023-05-14 00:51:28',
                   '2023-05-14 00:52:28', '2023-05-14 00:53:28',
                   '2023-05-14 00:54:28', '2023-05-14 00:55:28',
                   ...
                   '2026-03-06 23:50:59', '2026-03-06 23:51:59',
                   '2026-03-06 23:52:59', '2026-03-06 23:53:59',
                   '2026-03-06 23:54:59', '2026-03-06 23:55:59',
                   '2026-03-06 23:56:59', '2026-03-06 23:57:59',
                   '2026-03-06 23:58:59', '2026-03-06 23:59:59'],
                  dtype='datetime64[ns]', name='time', length=1288323, freq=None))

  • v
    PandasIndex

    PandasIndex(Index([0.05000000074505806, 0.15000000596046448,                0.25,
            0.3500000238418579, 0.44999998807907104,   0.550000011920929,
            0.6499999761581421,                0.75,  0.8500000238418579,
             0.949999988079071,   1.100000023841858,  1.2999999523162842,
                           1.5,  1.7000000476837158,   1.899999976158142,
             2.200000047683716,  2.5999999046325684,                 3.0,
            3.4000000953674316,   3.799999952316284,   4.400000095367432,
             5.199999809265137,                 6.0,   6.800000190734863,
             7.599999904632568,   8.800000190734863,  10.399999618530273,
                          12.0,  13.600000381469727,  15.199999809265137,
            17.600000381469727,  20.799999237060547],
          dtype='float32', name='v'))

• Attributes: (15)

  Conventions :

  CF-1.12

  _logical_cutoff_date :

  2026-03-07T00:00:00Z

  bcoproc_version :

  0.0.0.post1185.dev0+c830e01

  featureType :

  timeSeries

  history :

  processed with bcoproc l1

  institution :

  Max Planck Institute for Meteorology, Hamburg

  license :

  CC0-1.0

  location :

  The Barbados Cloud Observatory (BCO), Deebles Point, Barbados, West Indies

  platform :

  BCO

  sensor_serial_number :

  451995

  source :

  OTT Parsivel2 451995

  station_name :

  BCODIS

  summary :

  The OTT Parsivel² is a modern laser disdrometer for comprehensive measurement of all precipitation types. The Parsivel² captures both the size and speed of falling particles, classifying them into one of 32 separate size and velocity classes.

  title :

  Parsivel 2 Disdrometer data from BCO (Level 1)

  tool_versions :

  {"Python": "3.11.2 (main, Apr 28 2025, 14:11:48) [GCC 12.2.0]", "frozenlist": "1.6.0", "bcologs": "0.0.0.post1185.dev0+c830e01", "pluggy": "1.5.0", "pyyaml": "6.0.2", "toolz": "1.0.0", "idna": "3.10", "six": "1.17.0", "dask": "2025.4.1", "partd": "1.4.2", "aiosignal": "1.3.2", "pytz": "2025.2", "tzdata": "2025.2", "iniconfig": "2.1.0", "python-dateutil": "2.9.0.post0", "ruff": "0.13.0", "cftime": "1.6.4.post1", "fsspec": "2025.3.2", "packaging": "25.0", "xarray": "2025.4.0", "pytest": "8.3.5", "yarl": "1.20.0", "aiohappyeyeballs": "2.6.1", "pandas": "2.2.3", "importlib_metadata": "8.7.0", "locket": "1.0.0", "numpy": "2.2.5", "attrs": "25.3.0", "netcdf4": "1.7.2", "cloudpickle": "3.1.1", "click": "8.1.8", "certifi": "2025.4.26", "datarecords": "0.0.0.post1185.dev0+c830e01", "zipp": "3.21.0", "pendulum": "3.1.0", "propcache": "0.3.1", "multidict": "6.4.3", "bcoproc": "0.0.0.post1185.dev0+c830e01", "aiohttp": "3.11.18"}