Refine
Year of publication
Document Type
- Article (44)
- Monograph/Edited Volume (6)
- Postprint (6)
- Preprint (6)
- Part of a Book (2)
- Other (2)
- Doctoral Thesis (1)
- Part of Periodical (1)
Keywords
- discrepancy principle (5)
- regularization (5)
- aerosol size distribution (4)
- nonlinear operator (4)
- inverse ill-posed problem (3)
- inversion (3)
- iterative regularization (3)
- logarithmic source condition (3)
- multiwavelength lidar (3)
- 47A52 (2)
Institute
- Institut für Mathematik (52)
- Interdisziplinäres Zentrum für Dynamik komplexer Systeme (6)
- Extern (5)
- Institut für Physik und Astronomie (4)
- Zentrum für Umweltwissenschaften (4)
- Mathematisch-Naturwissenschaftliche Fakultät (3)
- Institut für Umweltwissenschaften und Geographie (1)
- Referat für Presse- und Öffentlichkeitsarbeit (1)
In this paper we extract the aerosol microphysical properties for a collection of mineral dust cases measured by multi-wavelength depolarization Raman lidar systems located at the National Technical University of Athens (NTUA, Athens, Greece) and the Andalusian Institute for Earth System Research (IISTA-CEAMA, Granada, Spain). The lidar-based retrievals were carried out with the Spheroidal Inversion eXperiments software tool (SphInX) developed at the University of Potsdam (Germany). The software uses regularized inversion of a two-dimensional enhancement of the Mie model based on the spheroid-particle approximation with the aspect ratio determining the particle shape. The selection of the cases was based on the transport time from the source regions to the measuring sites. The aerosol optical depth as measured by AERONET ranged from 0.27 to 0.54 (at 500 nm) depending on the intensity of each event. Our analysis showed the hourly mean particle linear depolarization ratio and particle lidar ratio values at 532 nm ranging from 11 to 34% and from 42 to 79 sr respectively, depending on the mixing status, the corresponding air mass pathways and their transport time. Cases with shorter transport time showed good agreement in terms of the optical and SphInX-retrieved microphysical properties between Athens and Granada providing a complex refractive index value equal to 1.4 + 0.004i. On the other hand, the results for cases with higher transport time deviated from the aforementioned ones as well as from each other, providing, in particular, an imaginary part of the refractive index ranging from 0.002 to 0.005. Reconstructions of two-dimensional shape-size distributions for each selected layer showed that the dominant effective particle shape was prolate with diverse spherical contributions. The retrieved volume concentrations reflect overall the intensity of the episodes.
In this work we extract the microphysical properties of aerosols for a collection of measurement cases with low volume depolarization ratio originating from fire sources captured by the Raman lidar located at the National Institute of Optoelectronics (INOE) in Bucharest. Our algorithm was tested not only for pure smoke but also for mixed smoke and urban aerosols of variable age and growth. Applying a sensitivity analysis on initial parameter settings of our retrieval code was proved vital for producing semi-automatized retrievals with a hybrid regularization method developed at the Institute of Mathematics of Potsdam University. A direct quantitative comparison of the retrieved microphysical properties with measurements from a Compact Time of Flight Aerosol Mass Spectrometer (CToF-AMS) is used to validate our algorithm. Microphysical retrievals performed with sun photometer data are also used to explore our results. Focusing on the fine mode we observed remarkable similarities between the retrieved size distribution and the one measured by the AMS. More complicated atmospheric structures and the factor of absorption appear to depend more on particle radius being subject to variation. A good correlation was found between the aerosol effective radius and particle age, using the ratio of lidar ratios (LR: aerosol extinction to backscatter ratios) as an indicator for the latter. Finally, the dependence on relative humidity of aerosol effective radii measured on the ground and within the layers aloft show similar patterns. (C) 2015 Elsevier Inc. All rights reserved.
In this work, an evaluation of an intense biomass burning event observed over Ny-angstrom lesund (Spitsbergen, European Arctic) in July 2015 is presented. Data from the multi-wavelengths Raman-lidar KARL, a sun photometer and radiosonde measurements are used to derive some microphysical properties of the biomass burning aerosol as size distribution, refractive index and single scattering albedo at different relative humidities. Predominantly particles in the accumulation mode have been found with a bi-modal distribution and dominance of the smaller mode. Above 80% relative humidity, hygroscopic growth in terms of an increase of particle diameter and a slight decrease of the index of refraction (real and imaginary part) has been found. Values of the single scattering albedo around 0.9 both at 355nm and 532nm indicate some absorption by the aerosol. Values of the lidar ratio are around 26sr for 355nm and around 50sr for 532nm, almost independent of the relative humidity. Further, data from the photometer and surface radiation values from the local baseline surface radiation network (BSRN) have been applied to derive the radiative impact of the biomass burning event purely from observational data by comparison with a clear background day. We found a strong cooling for the visible radiation and a slight warming in the infra-red. The net aerosol forcing, derived by comparison with a clear background day purely from observational data, obtained a value of -95 W/m(2) per unit AOD500.
This paper examines and develops matrix methods to approximate the eigenvalues of a fourth order Sturm-Liouville problem subjected to a kind of fixed boundary conditions. Furthermore, it extends the matrix methods for a kind of general boundary conditions. The idea of the methods comes from finite difference and Numerov's methods as well as boundary value methods for second order regular Sturm-Liouville problems. Moreover, the determination of the correction term formulas of the matrix methods is investigated in order to obtain better approximations of the problem with fixed boundary conditions since the exact eigenvalues for q = 0 are known in this case. Finally, some numerical examples are illustrated.
This paper examines and develops matrix methods to approximate the eigenvalues of a fourth order Sturm-Liouville problem subjected to a kind of fixed boundary conditions, furthermore, it extends the matrix methods for a kind of general boundary conditions. The idea of the methods comes from finite difference and Numerov's method as well as boundary value methods for second order regular Sturm-Liouville problems. Moreover, the determination of the correction term formulas of the matrix methods are investigated in order to obtain better approximations of the problem with fixed boundary conditions since the exact eigenvalues for q = 0 are known in this case. Finally, some numerical examples are illustrated.
In this paper, we present the convergence rate analysis of the modified Landweber method under logarithmic source condition for nonlinear ill-posed problems. The regularization parameter is chosen according to the discrepancy principle. The reconstructions of the shape of an unknown domain for an inverse potential problem by using the modified Landweber method are exhibited.
In this paper, we present the convergence rate analysis of the modified Landweber method under logarithmic source condition for nonlinear ill-posed problems. The regularization parameter is chosen according to the discrepancy principle. The reconstructions of the shape of an unknown domain for an inverse potential problem by using the modified Landweber method are exhibited.
In this paper, we investigate the continuous version of modified iterative Runge–Kutta-type methods for nonlinear inverse ill-posed problems proposed in a previous work. The convergence analysis is proved under the tangential cone condition, a modified discrepancy principle, i.e., the stopping time T is a solution of ∥𝐹(𝑥𝛿(𝑇))−𝑦𝛿∥=𝜏𝛿+ for some 𝛿+>𝛿, and an appropriate source condition. We yield the optimal rate of convergence.
In this paper, we investigate the continuous version of modified iterative Runge–Kutta-type methods for nonlinear inverse ill-posed problems proposed in a previous work. The convergence analysis is proved under the tangential cone condition, a modified discrepancy principle, i.e., the stopping time T is a solution of ∥𝐹(𝑥𝛿(𝑇))−𝑦𝛿∥=𝜏𝛿+ for some 𝛿+>𝛿, and an appropriate source condition. We yield the optimal rate of convergence.
Retrieval of aerosol extinction coefficient profiles from Raman lidar data by inversion method
(2012)
We regard the problem of differentiation occurring in the retrieval of aerosol extinction coefficient profiles from inelastic Raman lidar signals by searching for a stable solution of the resulting Volterra integral equation. An algorithm based on a projection method and iterative regularization together with the L-curve method has been performed on synthetic and measured lidar signals. A strategy to choose a suitable range for the integration within the framework of the retrieval of optical properties is proposed here for the first time to our knowledge. The Monte Carlo procedure has been adapted to treat the uncertainty in the retrieval of extinction coefficients.
This work is devoted to the convergence analysis of a modified Runge-Kutta-type iterative regularization method for solving nonlinear ill-posed problems under a priori and a posteriori stopping rules. The convergence rate results of the proposed method can be obtained under a Holder-type sourcewise condition if the Frechet derivative is properly scaled and locally Lipschitz continuous. Numerical results are achieved by using the Levenberg-Marquardt, Lobatto, and Radau methods.
This work is devoted to the convergence analysis of a modified Runge-Kutta-type iterative regularization method for solving nonlinear ill-posed problems under a priori and a posteriori stopping rules. The convergence rate results of the proposed method can be obtained under Hölder-type source-wise condition if the Fréchet derivative is properly scaled and locally Lipschitz continuous. Numerical results are achieved by using the Levenberg-Marquardt and Radau methods.
This paper further improves the Lie group method with Magnus expansion proposed in a previous paper by the authors, to solve some types of direct singular Sturm-Liouville problems. Next, a concrete implementation to the inverse Sturm-Liouville problem algorithm proposed by Barcilon (1974) is provided. Furthermore, computational feasibility and applicability of this algorithm to solve inverse Sturm-Liouville problems of higher order (for n=2,4) are verified successfully. It is observed that the method is successful even in the presence of significant noise, provided that the assumptions of the algorithm are satisfied. In conclusion, this work provides a method that can be adapted successfully for solving a direct (regular/singular) or inverse Sturm-Liouville problem (SLP) of an arbitrary order with arbitrary boundary conditions.
In this paper Lie group method in combination with Magnus expansion is utilized to develop a universal method applicable to solving a Sturm–Liouville problem (SLP) of any order with arbitrary boundary conditions. It is shown that the method has ability to solve direct regular (and some singular) SLPs of even orders (tested for up to eight), with a mix of (including non-separable and finite singular endpoints) boundary conditions, accurately and efficiently. The present technique is successfully applied to overcome the difficulties in finding suitable sets of eigenvalues so that the inverse SLP problem can be effectively solved. The inverse SLP algorithm proposed by Barcilon (1974) is utilized in combination with the Magnus method so that a direct SLP of any (even) order and an inverse SLP of order two can be solved effectively.
In this paper Lie group method in combination with Magnus expansion is utilized to develop a universal method applicable to solving a Sturm–Liouville problem (SLP) of any order with arbitrary boundary conditions. It is shown that the method has ability to solve direct regular (and some singular) SLPs of even orders (tested for up to eight), with a mix of (including non-separable and finite singular endpoints) boundary conditions, accurately and efficiently. The present technique is successfully applied to overcome the difficulties in finding suitable sets of eigenvalues so that the inverse SLP problem can be effectively solved. The inverse SLP algorithm proposed by Barcilon (1974) is utilized in combination with the Magnus method so that a direct SLP of any (even) order and an inverse SLP of order two can be solved effectively.
Regularized inversion of microphysical atmospheric particle parameters - theory and application
(2013)
Retrieving the distribution of aerosols in the atmosphere via remote sensing techniques is a highly complex task that requires dealing with a wide range of different problems stemming both from Physics and Mathematics. We focus on retrieving this distribution from multi-wavelength lidar data for aerosol ensembles consisting of spherical particles via an iterative regularization technique. The optical efficiencies for spherical scatterers are examined to account for the behavior of the underlying integral equation. The ill-posedness of the problem and the conditioning of the discretized problem are analyzed. Some critical points in the model, like the assumed wavelength-independence of the refractive index and the fixed grid of investigated refractive indices, are studied with regard to their expected impact on the regularized solution. A new Monte-Carlo type method is proposed for retrieval of the refractive index. To validate the results, the developed algorithm is applied to two measurement cases of burning biomass gained from multi-wavelength Raman lidar.
Strong events of long-range transported biomass burning aerosol were detected during July 2013 at three EARLINET (European Aerosol Research Lidar Network) stations, namely Granada (Spain), Leipzig (Germany) and Warsaw (Poland). Satellite observations from MODIS (Moderate Resolution Imaging Spectroradiometer) and CALIOP (Cloud-Aerosol Lidar with Orthogonal Polarization) instruments, as well as modeling tools such as HYSPLIT (Hybrid Single-Particle Lagrangian Integrated Trajectory) and NAAPS (Navy Aerosol Analysis and Prediction System), have been used to estimate the sources and transport paths of those North American forest fire smoke particles. A multiwavelength Raman lidar technique was applied to obtain vertically resolved particle optical properties, and further inversion of those properties with a regularization algorithm allowed for retrieving microphysical information on the studied particles. The results highlight the presence of smoke layers of 1-2 km thickness, located at about 5 km a.s.l. altitude over Granada and Leipzig and around 2.5 km a.s.l. at Warsaw. These layers were intense, as they accounted for more than 30% of the total AOD (aerosol optical depth) in all cases, and presented optical and microphysical features typical for different aging degrees: color ratio of lidar ratios (LR532/LR355) around 2, alpha-related angstrom exponents of less than 1, effective radii of 0.3 mu m and large values of single scattering albedos (SSA), nearly spectrally independent. The intensive microphysical properties were compared with columnar retrievals form co-located AERONET (Aerosol Robotic Network) stations. The intensity of the layers was also characterized in terms of particle volume concentration, and then an experimental relationship between this magnitude and the particle extinction coefficient was established.
Strong events of long-range transported biomass burning aerosol were detected during July 2013 at three EARLINET (European Aerosol Research Lidar Network) stations, namely Granada (Spain), Leipzig (Germany) and Warsaw (Poland). Satellite observations from MODIS (Moderate Resolution Imaging Spectroradiometer) and CALIOP (Cloud-Aerosol Lidar with Orthogonal Polarization) instruments, as well as modeling tools such as HYSPLIT (Hybrid Single-Particle Lagrangian Integrated Trajectory) and NAAPS (Navy Aerosol Analysis and Prediction System), have been used to estimate the sources and transport paths of those North American forest fire smoke particles. A multiwavelength Raman lidar technique was applied to obtain vertically resolved particle optical properties, and further inversion of those properties with a regularization algorithm allowed for retrieving microphysical information on the studied particles. The results highlight the presence of smoke layers of 1-2 km thickness, located at about 5 km a.s.l. altitude over Granada and Leipzig and around 2.5 km a.s.l. at Warsaw. These layers were intense, as they accounted for more than 30% of the total AOD (aerosol optical depth) in all cases, and presented optical and microphysical features typical for different aging degrees: color ratio of lidar ratios (LR532/LR355) around 2, alpha-related angstrom exponents of less than 1, effective radii of 0.3 mu m and large values of single scattering albedos (SSA), nearly spectrally independent. The intensive microphysical properties were compared with columnar retrievals form co-located AERONET (Aerosol Robotic Network) stations. The intensity of the layers was also characterized in terms of particle volume concentration, and then an experimental relationship between this magnitude and the particle extinction coefficient was established.