docs/depolarization/depolarization.rst

changeset 67
ae93ea881019
child 68
f697817dad5f
equal deleted inserted replaced
66:4c925b7dd7e5 67:ae93ea881019
1 **Single Calculus Chain **
2
3 **version: 4.0**
4
5 **date: Date (fixed)**
6
7 **DRAFT**
8
9 This document describes the main changes implemented in the SCC v4.0
10 with respect to what already provided in the SCC v3.11. It will be also
11 reported the modifications the users need to perform to run the new
12 version of SCC.
13
14 Table of Contents
15
16 1. Particle Linear Depolarization Ratio Implementation 3
17
18 1.1 Background 3
19
20 1.2 Polarization calibration 4
21
22 1.3 SCC procedure to calculate the PLDRP 4
23
24 2.Changes of the SCC input format 8
25
26 3.Real Example 10
27
28 3.1 Modification of polarization channel parameters 10
29
30 3.2 Definition of new calibration configuration and product 12
31
32 3.3 Definition of “Raman/Elastic backscatter and linear depolarization
33 ratio” 16
34
35 Particle Linear Depolarization Ratio Implementation
36 ===================================================
37
38 The most important improvement included in the SCC v4.0 is the
39 implementation of a new optical product which is the particle linear
40 depolarization ratio.
41
42 **Background**
43 --------------
44
45 The calculation of the volume linear depolarization ratio profile
46 (*VLDR*) and particle linear depolarization ratio profile (*PLDR*) needs
47 two different steps:
48
49 1. the calibration of the polarization sensitive lidar channels;
50
51 2. the calculation of the *VLDR* or *PLDR* itself.
52
53 The SCC allows the user to make both the above points. In particular the
54 calibration step is made by a completely new module called
55 **scc\_calibrator** which computes the *apparent calibration factor*
56 h\ :sup:`\*` out of the pre-processed data provided by the standard
57 **ELPP** (Earlinet Lidar Pre-Processor) module and it records it in the
58 SCC database (SCC\_DB). Once logged into the SCC\_DB this factor can be
59 used whenever it is necessary.
60
61 The raw lidar calibration measurements should be put in a NetCDF file
62 which has the same structure as the “standard” raw SCC NetCDF input file
63 (for more details see sections 2 and 3.2).
64
65 New signal types have been introduced to take into account special
66 channel configurations used for calibration purposes.
67
68 Moreover new product types for both calibration and *PLDR* calculation
69 have been defined. As, in principle, it is possible to calculate the
70 *PLDR* only when the aerosol backscatter coefficient profile is
71 available the following new products have been defined:
72
73 1. *Linear polarization calibration (factor* h) *(product\_type\_id=6);*
74
75 2. *Raman backscatter and linear depolarization ratio
76 (product\_type\_id=7);*
77
78 3. *Elastic backscatter and linear depolarization ratio
79 (product\_type\_id=8).*
80
81 The first product in the above list is used only for calibration while
82 the other two are used for the calculation of *PLDR*. Basically, in most
83 of the cases, the products 2 and 3 are equivalent to the corresponding
84 backscatter product types with the exception that also the following new
85 variables are available:
86
87 double VolumeDepol(Length) ;
88
89 double ErrorVolumeDepol(Length) ;
90
91 ErrorVolumeDepol:long\_name = "absolute error of VolumeDepol" ;
92
93 double ParticleDepol(Length) ;
94
95 double ErrorParticleDepol(Length) ;
96
97 ErrorParticleDepol:long\_name = "absolute error of ParticleDepol" ;
98
99 **Polarization calibration**
100 ----------------------------
101
102 An important point is the definition of reliable *PLDR* calibration
103 procedures. Within EARLINET the following calibration procedures are
104 currently used:
105
106 a) Rayleigh calibration;
107
108 b) +45 calibration method, or D90 calibration method (made by +45 and
109 -45 measurements);
110
111 c) 3 signals (total, cross and parallel).
112
113 It is well known that method a) could produce easily large errors on
114 *PLDR* which cannot be controlled. For this reason only the methods b)
115 and c) can be used to provide reliable polarization calibrations and so
116 only those methods will be implemented in the SCC.
117
118 For what it concerns the method c) it, basically, requires to solve the
119 equation:
120
121 in two different of atmospheric layers with considerably different
122 *VLDR*. So to calibrate in this way the implementation of automatic
123 layer identification in the SCC is required. As at moment this feature
124 is not yet available within the SCC *ONLY* the method b) is considered.
125
126 SCC procedure to calculate the PLDRP
127 ------------------------------------
128
129 According to what mentioned before the SCC calculates the *PLDR* through
130 the following steps:
131
132 1. The user needs to create a new system configuration in the SCC\_DB
133 including only lidar channels used for the calibration. One (or more)
134 *Linear polarization calibration (product\_type\_id=6)* product
135 should be associated to this new configuration (see section 3.2 for
136 more details);
137
138 2. This new system configuration should contain only the polarization
139 channels in the configuration used for the calibration (for example
140 rotated in the polarization plane of +45 degrees). A channel in
141 calibration measurement configuration should have a *DIFFERENT*
142 channel ID from the channel ID corresponding to the same channel in
143 standard measurement configuration. For example, if a system has two
144 polarization channels which in standard measurement configuration
145 correspond to the channel ID=1 and 2 respectively, the same physical
146 channels under calibration measurement configuration should
147 correspond to different channel IDs (let's say ID=3 and 4 for the +45
148 degrees polarization rotated channels and ID=5 and 6 for the -45
149 degrees polarization rotated ones in case D90 calibration method is
150 used). Moreover, the polarization channels should be labeled
151 correctly using the new signal types available (*+45elPT, +45elPR,
152 -45elPT, -45elPR, +45elPTnr, +45elPTfr, +45elPRnr, +45elPRfr,
153 -45elPTnr, -45elPTfr, -45elPRnr, -45elPRfr).* For more details see
154 section 3.2;
155
156 3. In SCC v4.0 the polarization channels are *NOT* labeled on the base
157 of their polarization state (as it was done in the SCC v3.11) but
158 *ALWAYS* as transmitted and reflected channels. So the channels that
159 in SCC v3.11 were labeled as *elCP, elCPnr, elCPfr, elPP, elPPnr
160 elPPfr* will be labeled in SCC v4.0 as *elPR, elPRnr elPRfr elPT,
161 elPTnr elPTfr* where the letter *T* stands from transmitted and the
162 letter *R* for reflected.
163
164 **WARNING:** In switching from the SCC v3.11 to SCC v4.0 the following
165 modifications have been made on *ALL* channels of *ALL* registered
166 configurations:
167
168 *elPP→elPR*
169
170 *elCP→elPT*
171
172 *elPPnr→elPRnr*
173
174 *elPPfr→ elPRfr*
175
176 *elCPnr→ elPTnr*
177
178 *elCPfr→ elPTfr*
179
180 Please be sure these modifications reflect to your actual lidar setup
181 (cross channels are transmitted and parallel channels are reflected);
182
183 1. The user needs to submit a file (same format as raw SCC input file)
184 containing the raw data for the lidar channels defined at the point 1
185 (see section 3.2 for more details);
186
187 2. The file at point 2 is pre-processed by **ELPP** module which applies
188 the standard pre-processing procedures applied to “standard” lidar
189 data;
190
191 3. The pre-processed files are then processed by the new modules
192 **scc\_calibrator** which calculates h\ :sup:`\*` *the apparent
193 calibration factor* and logs it into the SCC\_DB;
194
195 4. The user needs to create a new system configuration in the SCC\_DB
196 (which should be different from the one used for the calibration) and
197 associate it the new product *Raman backscatter and linear
198 depolarization ratio (product\_type\_id=7)* or *Elastic backscatter
199 and linear depolarization ratio (product\_type\_id=8).* Alternatively
200 the calculation of those products can be added to an already existing
201 lidar configuration as long as it is different from the calibration
202 one;
203
204 5. The product defined at point 5 should be linked to the product
205 containing the polarization calibration (defined at point 1) in a way
206 that the *apparent calibration factor* can be selected from the
207 SCC\_DB (see section 3.3 and in particular figure 3.4);
208
209 6. The user needs to submit another SCC raw data file containing the
210 “standard” measurements;
211
212 7. Finally **ELPP** and **ELDA** will produce a b-file containing
213 backscatter coefficient profile and *PLDR*. In particular this
214 calculation is made in two different steps: from the pre-processed
215 lidar polarization signals, and taking into account the *apparent
216 calibration factor* and the *calibration factor correction K*
217 (defined as option of *Linear polarization calibration* product\ *)*
218 written into the SCC\_DB, an “apparent” *VLDR* *d\ :sup:`\*`* is
219 calculated. Even if *d\ :sup:`\*`* is a calibrated quantity it can be
220 still affected by possible systematic errors due to not perfect
221 optics or alignment of the system;
222
223 8. To take into account these errors a corrected *VLDR* (*d)* is
224 calculated using the *polarization cross-talk correction parameters*
225 *G* and *H* calculated on the base of Müller matrix formalism. These
226 cross-talk correction parameters (*G* and *H*) are stored in the
227 SCC\_DB for each lidar channels (see section 3.1 in particular figure
228 3.2). Finally the *PLDR* is calculated using the backscatter
229 coefficient profile and the molecular LDRP calculated by ELPP
230 considering the center wavelength and bandwidth of the channels
231 interference filter.
232
233 The *apparent calibration factor* h\ :sup:`\*` is calculated by the
234 **scc\_calibrator** module as the geometrical mean of the ratio of the
235 +/-45 degrees reflected to the +/- 45 degrees transmitted signals within
236 an altitude calibration range defined by the users in the raw data input
237 files.
238
239 In case of +45 calibration method h\ :sup:`\*` is calculated by:
240
241 (1.1)
242
243 While in case of D90 calibration method:
244
245 (1.2)
246
247 **ELDA** module calculates the “apparent” *VLDR*:
248
249 (1.3)
250
251 the *VLDR*
252
253 (1.4)
254
255 and the *PLDR*
256
257 (1.5)
258
259 where:
260
261 - h\ :sup:`\*` is the *apparent calibration factor* calculated by
262 **scc\_calibrator**
263
264 - *K* is the *calibration factor correction* defined as polarization
265 product option
266
267 - *I\ :sub:`T`* and I\ *:sub:`R`* are the transmitted and the reflected
268 signals in the polarization detection set-up
269
270 - *G\ :sub:`T,R`* and *H\ :sub:`T,R`* are *polarization cross-talk
271 correction parameters* for the transmitted and reflected signals used
272 to correct for systematic errors. Both these factors are defined in
273 the SCC\_DB for each lidar channel.
274
275 - *d\ :sub:`m`* is the molecular linear depolarization ratio calculated
276 by ELPP
277
278 - *R* is the backscatter ratio
279
280 Please note once again that the polarization channels are described in
281 terms of transmitted and reflected signals. This means that according to
282 different lidar instrumental configurations, the transmitted or the
283 reflected channel can contain total, perpendicular or parallel polarized
284 signals.
285
286 In order to retrieve the backscatter profile the total signal must be
287 obtained combining the transmitted and reflected polarized signals. The
288 following formula is used:
289
290 (1.6)
291
292 The formulas above are general and can be adapted to all possible
293 polarization lidar configurations selecting the right polarization
294 cross-talk correction parameters (see Table 1.1).
295
296 Let's suppose, for example, we have the perpendicular polarized lidar
297 signal on the transmitted channel and the parallel polarized on
298 reflected channel. For an ideal system (no diattenuation and cross-talk)
299 we have:
300
301 If, on the other hands, we have the perpendicular polarized lidar signal
302 on reflected channel and the total polarized on the transmitted for and
303 ideal system we have:
304
305 **Table 1.1:** Polarization cross-talk correction parameters for ideal
306 systems
307
308 +----------------------+-----------------------------+-----------------+-----------------+-----------------+
309 | Laser polarization | Detected in lidar channel |
310 +----------------------+-----------------------------+-----------------+-----------------+-----------------+
311 | | Transmitted | Reflected |
312 +----------------------+-----------------------------+-----------------+-----------------+-----------------+
313 | | *G\ :sub:`T`* | *H\ :sub:`T`* | *G\ :sub:`R`* | *H\ :sub:`R`* |
314 +----------------------+-----------------------------+-----------------+-----------------+-----------------+
315 | total | 1 | 0 | 1 | 0 |
316 +----------------------+-----------------------------+-----------------+-----------------+-----------------+
317 | parallel | 1 | 1 | 1 | 1 |
318 +----------------------+-----------------------------+-----------------+-----------------+-----------------+
319 | cross | 1 | -1 | 1 | -1 |
320 +----------------------+-----------------------------+-----------------+-----------------+-----------------+
321
322 The *apparent calibration factor* (h:sup:`\*`), *the calibration factor
323 correction* (*K*) and the *polarization cross-talk correction
324 parameters* are stored by **ELPP** module in the intermediate NetCDF
325 files using the following variables:
326
327 - *Polarization\_Channel\_Gain\_Factor (apparent calibration factor* –
328 h\ :sup:`\*`)
329
330 - *Polarization\_Channel\_Gain\_Factor\_Correction (calib. factor
331 corr.* – *K*)
332
333 - *G\_T*
334
335 - *H\_T*
336
337 - *G\_R*
338
339 - *H\_R*
340
341 Finally new usecases have been defined to take into account all the
342 possible lidar configurations. The details on that are provided as a
343 separate file.
344
345 Changes of the SCC input format
346 ===============================
347
348 The following minor changes have been applied to raw SCC data format:
349
350 1. The optional variable *ID\_Range* has been *REMOVED*;
351
352 2. The *OPTIONAL* variable *int Signal\_Type(channels)* has been added.
353 The possible values are the same available in the SCC\_DB:
354
355 0 *→* elT
356
357 1 *→* elTnr
358
359 2 *→* elTfr
360
361 3 *→* vrRN2
362
363 4 *→* vrRN2nr
364
365 5 *→* vrRN2fr
366
367 6 *→* elPR
368
369 7 *→* elPT
370
371 8 *→* pRRlow
372
373 9 *→* pRRhigh
374
375 10 *→* elPRnr
376
377 11 *→* elPRfr
378
379 12 *→* elPTnr
380
381 13 *→* elPTfr
382
383 14 *→* vrRH2O
384
385 15 *→* pRRhighnr
386
387 16 *→* pRRhighfr
388
389 17 *→* pRRlownr
390
391 18 *→* pRRlowfr
392
393 19 *→* vrRH2Onr
394
395 20 *→* vrRH2Ofr
396
397 21 *→* elTunr
398
399 *22 → +45elPT*
400
401 *23 → +45elPR*
402
403 *24 → -45elPT*
404
405 *25 → -45elPR*
406
407 *26 → +45elPTnr*
408
409 *27 → +45elPTfr*
410
411 *28 → +45elPRnr*
412
413 *29 → +45elPRfr*
414
415 *30 → -45elPTnr*
416
417 *31 → -45elPTfr*
418
419 *32 → -45elPRnr*
420
421 *33 → -45elPRfr*
422
423 **WARNING:** It this variable is found in the SCC input file the
424 corresponding settings in the SCC database will be *overwritten*. Unless
425 you don't have any valid reason to overwrite the database value this
426 variable should not be used.
427
428 1. The variables:
429
430 *double Pol\_Calib\_Range\_Min(channels)*
431
432 *double Pol\_Calib\_Range\_Max(channels) *
433
434 have been added. Both these variable are *mandatory* for any calibration
435 raw dataset.
436
437 These variable should be included only the polarization calibration
438 measurements and should specify the altitude range (meters) in which the
439 polarization calibration should be made. For more details see section
440 3.3;
441
442 1. The variable *Depolarization\_Factor* has been *removed*.
443
444 The SCC v3.11 used this variable to get polarization calibration factor
445 for the calculation of the total signal out of cross and parallels ones.
446 As the SCC v4.0 is able to calculate the same parameter by itself, the
447 use of this variable is *NOT* possible anymore. The recommended way to
448 get a valid and quality assured depolarization calibration factor is to
449 submit to the SCC v4.0 a polarization calibration dataset and let the
450 SCC to calculate such factor.
451
452 To make this change more smooth and to provide the users with the
453 possibility to continue to analyze their data with the SCC v4.0 even if
454 a calibration dataset has not been submitted yet, it will be possible
455 for a *LIMITED* period of time to submit the calibration constant via
456 the SCC web interface. The SCC will keep track of the used calibration
457 method (automatic or manual).
458
459 **WARNING:** After this transition period *only* automatic calibration
460 will be allowed!
461
462 1. The new *optional* variable:
463
464 *string channel\_string\_ID(channels)*
465
466 has been introduced.
467
468 Starting from SCC v4.0 the lidar channel can be identified not only by
469 using integers (as it happened until SCC v3.11) but also by using
470 strings.
471
472 The procedure implemented in the SCC v4.0 to recognize the lidar channel
473 within the raw lidar data is fully backward compatible (old format files
474 are accepted as they are by SCC v4.0).
475
476 **WARNING:** Please note that the definition of the new string variable
477 requires netCDF-4 format! The type *string* is not supported in netCDF-3
478 format!
479
480 Real Example
481 ============
482
483 This section describes all the practical steps the users need to follow
484 to switch from SCC v3.11 to new SCC v4.0.
485
486 **IMPORTANT:**
487
488 If your lidar system is not equipped with any polarization channels *NO*
489 changes are required. In this case, the SCC v4.0 should work using the
490 same input files and the same database configurations you have used with
491 the SCC v3.11. Anyway as in the SCC v4.0 several bugs have been fixed,
492 it is recommended to re-run all the measurement IDs you have submitted.
493 For doing that you just need to reprocess all your data without the need
494 to submit raw data files already uploaded on the server.
495
496 The practical example reported below describes the modifications
497 required to use the SCC v4.0 for lidar systems equipped with
498 polarization channels.
499
500 Modification of polarization channel parameters
501 -----------------------------------------------
502
503 In what it follows it is assumed you already have registered one or more
504 lidar configurations in the SCC database and that such configurations
505 have been already used to produce optical products (aerosol extinction
506 and/or backscatter coefficients) by means of the SCC v3.11.
507
508 Let's assume your 3+2 system is registered in the SCC database and the
509 settings used by the SCC v3.11 are the ones summarized in table 3.1.
510
511 **Table 3.1:** Example of configuration in SCC v3.11
512
513 +----------------+--------------+----------------+-------------+-----------+
514 | Channel Name | Channel ID | Channel Type | nighttime | daytime |
515 +----------------+--------------+----------------+-------------+-----------+
516 | 355 | 1 | elT |  |  |
517 +----------------+--------------+----------------+-------------+-----------+
518 | 387 | 2 | vrRN2 |  | |
519 +----------------+--------------+----------------+-------------+-----------+
520 | 532 cross | 3 | elCP |  |  |
521 +----------------+--------------+----------------+-------------+-----------+
522 | 532 parallel | 4 | elPP |  |  |
523 +----------------+--------------+----------------+-------------+-----------+
524 | 607 | 5 | vrRN2 |  | |
525 +----------------+--------------+----------------+-------------+-----------+
526 | 1064 | 6 | elT |  |  |
527 +----------------+--------------+----------------+-------------+-----------+
528
529 We assume there are 2 system configurations called “nighttime” and
530 “daytime”. The nighttime configuration contains all the available lidar
531 channels (in order to calculate, for example, the aerosol extinction at
532 355 and 532nm and the aerosol backscatter at 355, 532 and 1064nm) while
533 in daytime conditions only elastic channels are used (only elastic
534 backscatter coefficients are generated).
535
536 To make these settings working with SCC v4.0 it is needed to modify
537 *ONLY* the products properties involving the polarization channels (532
538 cross and parallel). All the products not involving the polarization
539 channels *DO NOT* need any modification and should work in the SCC v4.0
540 exactly as they did in SCC v3.11. In the example above the aerosol
541 extinction and backscatter coefficient at 355nm, the extinction at 532nm
542 as well as the backscatter coefficient at 1064nm do not required any
543 modification. Let's focus on the modifications needed for the
544 calculation of backscatter at 532nm.
545
546 |image0| How to select signal types
547
548 The first modification concerns the settings of the channel type for the
549 532 cross and 532 parallel polarization channels. Starting from SCC v4.0
550 polarization channels are identified as transmitted and reflected
551 polarization channels and not on the base of their polarization state.
552 So suppose if we suppose the cross polarized channel is transmitted by a
553 polarizer beam splitter cube and the parallel is reflected the value
554 reported in table 3.1 should be modified as they appear in table 3.2. So
555 using the SCC web interface, the signal type of the 532 cross channel
556 should be changed from elCP to elPT and and in the same way the 532
557 parallel channel should be changed from elPP to elPR (see figure 3.1).
558
559 **Table 3.2:** The same of table 3.1 but with new channel types
560 introduced in SCC v4.0
561
562 +----------------+--------------+----------------+-------------+-----------+
563 | Channel Name | Channel ID | Channel Type | nighttime | daytime |
564 +----------------+--------------+----------------+-------------+-----------+
565 | 355 | 1 | elT |  |  |
566 +----------------+--------------+----------------+-------------+-----------+
567 | 387 | 2 | vrRN2 |  | |
568 +----------------+--------------+----------------+-------------+-----------+
569 | 532 cross | 3 | **elPT** |  |  |
570 +----------------+--------------+----------------+-------------+-----------+
571 | 532 parallel | 4 | **elPR** |  |  |
572 +----------------+--------------+----------------+-------------+-----------+
573 | 607 | 5 | vrRN2 |  | |
574 +----------------+--------------+----------------+-------------+-----------+
575 | 1064 | 6 | elT |  |  |
576 +----------------+--------------+----------------+-------------+-----------+
577
578 The other change about the polarization channels required to run the SCC
579 v4.0 is the definition of the polarization crosstalk parameters for all
580 the polarization channels available. Such parameters can be defined for
581 each polarization channel using the SCC web interface (see figure 3.2).
582 In particular among the channel parameters there is a new tab called
583 *Polarization crosstalk parameters* where it is possible to insert the
584 values from for the parameters *G* and *H* and the corresponding
585 statistical and systematic errors if available. In case you have
586 measured *G* and *H* for your polarization channels please insert the
587 corresponding values there. Otherwise you can insert the ideal values as
588 reported in table 1.1.
589
590 |image1| *Polarization crosstalk parameters* tab in channel properties
591 (SCC v4.0).
592
593 Definition of new calibration configuration and product
594 -------------------------------------------------------
595
596 In this section we will see how to set the polarization calibration
597 parameters: the calibration constant (called h\ :sup:`\*` in section
598 1.3) and the correction to calibration constant (called K in section
599 1.3).
600
601 In order to provide such parameters you need to define a new system
602 configuration to be used *only* for calibration purposes. Such new
603 configuration should include the polarization channels in the
604 measurement configuration used for the calibration. Let's suppose we
605 want to use the D90 calibration method.
606
607 In this case we need to define a new configuration (called for example
608 “depol\_calibration”) as reported in the table 3.3. As you can see the
609 configuration “depol\_calibration” includes 4 “new” channels. Actually
610 the channels “532 cross +45 degrees” (channel ID=10) and “532 cross -45
611 degrees” (channel ID=12) refer to the same physical channel “532 cross”
612 reported with channel ID=3 in table 3.2. Anyway we need to define two
613 new channel IDs to identify the “532 cross” channel in the two
614 polarization rotated configurations (+45 and -45 degrees) needed to
615 apply the D90 calibration method. The same is true for the “532
616 parallel” channel. The polarization rotated channels should be labeled
617 with the corresponding signal type as reported in table 3.3 (see figure
618 3.1).
619
620 **Table 3.3:** Polarization calibration configurations assuming D90
621 calibration method
622
623 +----------------------------+--------------+----------------+----------------------+
624 | Channel Name | Channel ID | Channel Type | depol\_calibration |
625 +----------------------------+--------------+----------------+----------------------+
626 | 532 cross +45 degrees | 10 | +45elPT |  |
627 +----------------------------+--------------+----------------+----------------------+
628 | 532 parallel +45 degrees | 11 | +45elPR |  |
629 +----------------------------+--------------+----------------+----------------------+
630 | 532 cross -45 degrees | 12 | -45elPT |  |
631 +----------------------------+--------------+----------------+----------------------+
632 | 532 parallel -45 degrees | 13 | -45elPR |  |
633 +----------------------------+--------------+----------------+----------------------+
634
635 Finally we should add to the configuration “depol\_calibration” a
636 product “\ *Linear polarization calibration”* to be used for the
637 calibration. According to the example given above and to the usecase
638 document attached we should use an usecase=4 for this example.
639
640 Other “\ *Linear polarization calibration”* options to be specified are
641 reported in figure 3.3. The most important factor you should insert here
642 is the *Pol calibration correction factor* (K). The ideal value for this
643 parameter is 1. Anyway if you have measured the parameter K please fill
644 in the measured value and the corresponding measurement errors.
645
646 |image2| Options for *Linear polarization calibration product*.
647
648 As you can see it is possible to fill in only the K correction factor
649 and not the calibration constant h\ :sup:`\*`.
650
651 Actually for a *LIMITED* period of time it will be possible to fill in
652 also the constant h\ :sup:`\*` using a temporary tab called
653 *Polarization calibration constant*. This has been done to provide the
654 users with the possibility to continue to use the SCC even if an
655 automatic calibration made by the SCC was not submitted yet. Anyway
656 after a transition period it will be *not* possible to provide
657 calibration constant using this procedure and the parameter h\ :sup:`\*`
658 can be calculated *ONLY* by the SCC as result of the submission of a
659 proper calibration raw input dataset. The format of this input file is
660 the same as the standard SCC input file. The only difference is that is
661 should contain calibration measurements instead of standard
662 measurements. Following our example, such file should contain the
663 measurement performed at +45 and -45 degrees at 532nm. Also the channel
664 IDs in the file should reflect the ones reported in table 3.3.
665
666 Moreover this raw input file has to contain the variables:
667
668 *double Pol\_Calib\_Range\_Min(channels)*
669
670 *double Pol\_Calib\_Range\_Max(channels) *
671
672 where to specify the altitude ranges in meters in which the polarization
673 calibration should be done.
674
675 According to the table 3.3 this file should be something similar to:
676
677 dimensions:
678
679 channels = 4 ;
680
681 nb\_of\_time\_scales = 1 ;
682
683 points = 16380 ;
684
685 scan\_angles = 1 ;
686
687 time = UNLIMITED ; // (3 currently)
688
689 variables:
690
691 int channel\_ID(channels) ;
692
693 double Background\_Low(channels) ;
694
695 double Background\_High(channels) ;
696
697 int id\_timescale(channels) ;
698
699 double Laser\_Pointing\_Angle(scan\_angles) ;
700
701 int Molecular\_Calc ;
702
703 int Laser\_Pointing\_Angle\_of\_Profiles(time, nb\_of\_time\_scales) ;
704
705 int Raw\_Data\_Start\_Time(time, nb\_of\_time\_scales) ;
706
707 int Raw\_Data\_Stop\_Time(time, nb\_of\_time\_scales) ;
708
709 int Laser\_Shots(time, channels) ;
710
711 double Raw\_Lidar\_Data(time, channels, points) ;
712
713 double Pressure\_at\_Lidar\_Station ;
714
715 double Temperature\_at\_Lidar\_Station ;
716
717 double Pol\_Calib\_Range\_Min(channels) ;
718
719 double Pol\_Calib\_Range\_Max(channels) ;
720
721 // global attributes:
722
723 :System = "mysystem" ;
724
725 :Longitude\_degrees\_east = 15.723771 ;
726
727 :RawData\_Start\_Time\_UT = "220000" ;
728
729 :RawData\_Start\_Date = "20130620" ;
730
731 :Measurement\_ID = "20130620po00" ;
732
733 :Altitude\_meter\_asl = 760. ;
734
735 :RawData\_Stop\_Time\_UT = "230333" ;
736
737 :Latitude\_degrees\_north = 40.601039 ;
738
739 data:
740
741 channel\_ID = 10, 11, 12, 13 ;
742
743 Background\_Low = 30000, 30000, 30000, 30000 ;
744
745 Background\_High = 50000, 50000, 50000, 50000 ;
746
747 id\_timescale = 0, 0, 0, 0 ;
748
749 Laser\_Pointing\_Angle = 0 ;
750
751 Molecular\_Calc = 0 ;
752
753 Laser\_Pointing\_Angle\_of\_Profiles =
754
755 0,
756
757 0,
758
759 0 ;
760
761 Raw\_Data\_Start\_Time =
762
763 0,
764
765 300,
766
767 600 ;
768
769 Raw\_Data\_Stop\_Time =
770
771 210,
772
773 510,
774
775 810 ;
776
777 Laser\_Shots =
778
779 1200, 1200, 1200, 1200,
780
781 1200, 1200, 1200, 1200,
782
783 1200, 1200, 1200, 1200 ;
784
785 Pressure\_at\_Lidar\_Station = 1010 ;
786
787 Temperature\_at\_Lidar\_Station = 14 ;
788
789 Pol\_Calib\_Range\_Min = 1000, 1000, 1000, 1000 ;
790
791 Pol\_Calib\_Range\_Min = 2000, 2000, 2000, 2000 ;
792
793 Raw\_Lidar\_Data = …...;
794
795 The file above assume the following calibration measurements have been
796 done:
797
798 1. First +45 degrees acquisition followed by a corresponding -45 degrees
799 acquisition
800
801 a. Measurement at +45 degrees
802
803 Start Time: 20130620 22:00:00
804
805 Stop Time: 20130620 22:01:00
806
807 Shots: 1200
808
809 a. Measurement at -45 degrees
810
811 Start Time: 20130620 22:02:30
812
813 Stop Time: 20130620 22:03:30
814
815 Shots: 1200
816
817 1. Second +45 degrees acquisition followed by a corresponding -45
818 degrees acquisition
819
820 a. Measurement at +45 degrees
821
822 Start Time: 20130620 22:05:00
823
824 Stop Time: 20130620 22:06:00
825
826 Shots: 1200
827
828 a. Measurement at -45 degrees
829
830 Start Time: 20130620 22:07:30
831
832 Stop Time: 20130620 22:08:30
833
834 Shots: 1200
835
836 1. Third +45 degrees acquisition followed by a corresponding -45 degrees
837 acquisition
838
839 a. Measurement at +45 degrees
840
841 Start Time: 20130620 22:10:00
842
843 Stop Time: 20130620 22:11:00
844
845 Shots: 1200
846
847 a. Measurement at -45 degrees
848
849 Start Time: 20130620 22:12:30
850
851 Stop Time: 20130620 22:13:30
852
853 Shots: 1200
854
855 As you can see there are 3 cycles of consecutive measurements at +45 and
856 -45 degrees. That's way the dimension time is set to 3.
857
858 The first +/-45 degrees measurement starts at “20130620 22:00:00” (start
859 time of the first +45 measurement) and stops at “20130620 22:03:30”
860 (stop time of the fist -45 measurement). As a consequence, according to
861 the values of the global attributes RawData\_Start\_Date and
862 RawData\_Start\_Time\_UT we have to set:
863
864 Raw\_Data\_Start\_Time[0]=0 (start of the first +45 measurement in
865 seconds since RawData\_Start\_Time\_UT)
866
867 Raw\_Data\_Stop\_Time[0]=210 (stop of the first -45 measurement in
868 seconds since RawData\_Start\_Time\_UT)
869
870 Following a similar procedure for the other 2 cycles we have:
871
872 Raw\_Data\_Start\_Time[1]=300 (start of the second +45 measurement in
873 seconds since RawData\_Start\_Time\_UT)
874
875 Raw\_Data\_Stop\_Time[1]=510 (stop of the second -45 measurement in
876 seconds since RawData\_Start\_Time\_UT)
877
878 Raw\_Data\_Start\_Time[2]=600 (start of the third +45 measurement in
879 seconds since RawData\_Start\_Time\_UT)
880
881 Raw\_Data\_Stop\_Time[2]=810 (stop of the third -45 measurement in
882 seconds since RawData\_Start\_Time\_UT)
883
884 Moreover, according to the order of the channels in the channel\_ID
885 variable, the Raw\_Lidar\_Data array should be filled as it follows:
886
887 Raw\_Lidar\_Data[0][0][points] → 1\ :sup:`st` measured transmitted
888 signal at +45 degrees
889
890 Raw\_Lidar\_Data[0][1][points] → 1\ :sup:`st` measured reflected signal
891 at +45 degrees
892
893 Raw\_Lidar\_Data[0][2][points] → 1\ :sup:`st` measured transmitted
894 signal at -45 degrees
895
896 Raw\_Lidar\_Data[0][3][points] → 1\ :sup:`st` measured reflected signal
897 at -45 degrees
898
899 Raw\_Lidar\_Data[1][0][points] → 2\ :sup:`nd` measured transmitted
900 signal at +45 degrees
901
902 Raw\_Lidar\_Data[1][1][points] → 2\ :sup:`nd` measured reflected signal
903 at +45 degrees
904
905 Raw\_Lidar\_Data[1][2][points] → 2\ :sup:`nd` measured transmitted
906 signal at -45 degrees
907
908 Raw\_Lidar\_Data[1][3][points] → 2\ :sup:`nd` measured reflected signal
909 at -45 degrees
910
911 Raw\_Lidar\_Data[2][0][points] → 3\ :sup:`rd` measured transmitted
912 signal at +45 degrees
913
914 Raw\_Lidar\_Data[2][1][points] → 3\ :sup:`rd` measured reflected signal
915 at +45 degrees
916
917 Raw\_Lidar\_Data[2][2][points] → 3\ :sup:`rd` measured transmitted
918 signal at -45 degrees
919
920 Raw\_Lidar\_Data[2][3][points] → 3\ :sup:`rd` measured reflected signal
921 at -45 degrees
922
923 Once this file has been created it needs to be submitted to the SCC and
924 linked to the configuration “depol\_calibration”. The result of the SCC
925 analysis on this file will be the calculation of the calibration
926 constant h\ :sup:`\*` that will be logged into the SCC database and can
927 be used to calibrate Raman/Elastic backscat ter products (see section
928 3.3).
929
930 **Definition of “Raman/Elastic backscatter and linear depolarization ratio”**
931 -----------------------------------------------------------------------------
932
933 In order to calculate the *PLDR* we need to modify the polarization
934 related products linked to the “standard” measurement configurations
935 (the configuration called “nighttime” and/or “daytime” in table 3.2).
936
937 Let's suppose we have defined the following products (defined already in
938 SCC v3.11):
939
940 **Table 3.4:** Example of products configuration in SCC v3.11
941
942 +-----------------------+--------------+-----------------------+-------------+-----------+
943 | Product Name | Product ID | Product Type | nighttime | daytime |
944 +-----------------------+--------------+-----------------------+-------------+-----------+
945 | Raman backscatter | 1 | Raman backscatter |  | |
946 | | | | | |
947 | 355nm | | | | |
948 +-----------------------+--------------+-----------------------+-------------+-----------+
949 | Extinction | 2 | Extinction |  | |
950 | | | | | |
951 | 387nm | | | | |
952 +-----------------------+--------------+-----------------------+-------------+-----------+
953 | Raman backscatter | 3 | Raman backscatter |  | |
954 | | | | | |
955 | 532nm | | | | |
956 +-----------------------+--------------+-----------------------+-------------+-----------+
957 | Extinction | 4 | Extinction |  | |
958 | | | | | |
959 | 532nm | | | | |
960 +-----------------------+--------------+-----------------------+-------------+-----------+
961 | Elastic backscatter | 5 | Elastic backscatter | |  |
962 | | | | | |
963 | 355nm | | | | |
964 +-----------------------+--------------+-----------------------+-------------+-----------+
965 | Elastic backscatter | 6 | Elastic backscatter | |  |
966 | | | | | |
967 | 532nm | | | | |
968 +-----------------------+--------------+-----------------------+-------------+-----------+
969 | Elastic backscatter | 7 | Elastic backscatter |  |  |
970 | | | | | |
971 | 1064nm | | | | |
972 +-----------------------+--------------+-----------------------+-------------+-----------+
973
974 Product ID=1, 2, 4, 5, 7 do not need any modification as they do not
975 involve polarization channels. The only product that need to be modified
976 are the Product ID=3 and 6. To produce b532 files containing also *PLDR*
977 we need to modify the “nighttime” and “daytime” configurations to
978 include a product of type “Raman bakscatter and linear depolarization
979 ratio” or “Elastic bakscatter and linear depolarization ratio”
980 respectively. So the configuration reported in table 3.4 should be
981 changed to match what is included in table 3.5.
982
983 **Table 3.5:** The same of table 3.4 but with new product types
984 introduced in SCC v4.0
985
986 +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+
987 | Product Name | Product ID | Product Type | nighttime | daytime |
988 +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+
989 | Raman backscatter | 1 | Raman backscatter |  | |
990 | | | | | |
991 | 355nm | | | | |
992 +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+
993 | Extinction | 2 | Extinction |  | |
994 | | | | | |
995 | 387nm | | | | |
996 +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+
997 | Raman backscatter | 10 | **Raman backscatter and linear depolarization ratio** |  | |
998 | | | | | |
999 | 532nm | | | | |
1000 +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+
1001 | Extinction | 4 | Extinction |  | |
1002 | | | | | |
1003 | 532nm | | | | |
1004 +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+
1005 | Elastic backscatter | 5 | Elastic backscatter | |  |
1006 | | | | | |
1007 | 355nm | | | | |
1008 +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+
1009 | Elastic backscatter | 11 | **Elastic backscatter and linear depolarization ratio** | |  |
1010 | | | | | |
1011 | 532nm | | | | |
1012 +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+
1013 | Elastic backscatter | 7 | Elastic backscatter |  |  |
1014 | | | | | |
1015 | 1064nm | | | | |
1016 +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+
1017
1018 As you can see in table 3.5, the old product IDs=3 and 6 (present in
1019 table 3.4) have been replaced with the new product ID=10 and 11 to
1020 guarantee the calculation of *PLDR*.
1021
1022 It is important to set among the product options of the product ID=10
1023 and 11 which calibration product we want to use for calibration (see
1024 section 3.2). This can be done using the SCC web interface setting the
1025 appropriate setting in the tab *Polarization calibration products* (see
1026 figure 3.4). According to the current example you should set here the
1027 calibration product defined in section 3.2.
1028
1029 |image3| How to link a product to calibrate with a calibration product.
1030
1031 **WARNING:** Please not that also *Raman/Elastic backscatter products*
1032 need to be linked to a calibration product because the calibration
1033 constant and the corresponding correction factor is needed to calculate
1034 the total signal out of the two polarization components even if the
1035 *PLDR* is not involved in the product calculation.
1036
1037 .. |image0| image:: ./media/image1.png
1038 :width: 6.69514in
1039 :height: 2.40764in
1040 .. |image1| image:: ./media/image2.png
1041 :width: 6.69306in
1042 :height: 1.71458in
1043 .. |image2| image:: ./media/image3.png
1044 :width: 6.69306in
1045 :height: 1.77431in
1046 .. |image3| image:: ./media/image4.png
1047 :width: 6.69306in
1048 :height: 0.36389in

mercurial