Refine
Year of publication
Document Type
- Article (40)
- Monograph/Edited Volume (3)
- Postprint (3)
- Other (2)
- Preprint (2)
Is part of the Bibliography
- yes (50) (remove)
Keywords
- discrepancy principle (5)
- regularization (5)
- nonlinear operator (4)
- logarithmic source condition (3)
- Aerosols (2)
- Arctic haze (2)
- Boundary value methods (2)
- Eigenvalues (2)
- Finite difference method (2)
- Fourth order Sturm-Liouville problem (2)
Institute
- Institut für Mathematik (50) (remove)
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.
The Levenberg–Marquardt regularization for the backward heat equation with fractional derivative
(2022)
The backward heat problem with time-fractional derivative in Caputo's sense is studied. The inverse problem is severely ill-posed in the case when the fractional order is close to unity. A Levenberg-Marquardt method with a new a posteriori stopping rule is investigated. We show that optimal order can be obtained for the proposed method under a Hölder-type source condition. Numerical examples for one and two dimensions are provided.
The Ill-posed Problem of Multiwavelength Lidar Data by a Hybrid Method of Variable Projection
(1999)
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.
The Runge-Kutta type regularization method was recently proposed as a potent tool for the iterative solution of nonlinear ill-posed problems. In this paper we analyze the applicability of this regularization method for solving inverse problems arising in atmospheric remote sensing, particularly for the retrieval of spheroidal particle distribution. Our numerical simulations reveal that the Runge-Kutta type regularization method is able to retrieve two-dimensional particle distributions using optical backscatter and extinction coefficient profiles, as well as depolarization information.
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.
In this work, a closure experiment for tropospheric aerosol is presented. Aerosol size distributions and single scattering albedo from remote sensing data are compared to those measured in-situ. An aerosol pollution event on 4 April 2009 was observed by ground based and airborne lidar and photometer in and around Ny-Alesund, Spitsbergen, as well as by DMPS, nephelometer and particle soot absorption photometer at the nearby Zeppelin Mountain Research Station.
The presented measurements were conducted in an area of 40 x 20 km around Ny-Alesund as part of the 2009 Polar Airborne Measurements and Arctic Regional Climate Model Simulation Project (PAMARCMiP). Aerosol mainly in the accumulation mode was found in the lower troposphere, however, enhanced backscattering was observed up to the tropopause altitude. A comparison of meteorological data available at different locations reveals a stable multi-layer-structure of the lower troposphere. It is followed by the retrieval of optical and microphysical aerosol parameters. Extinction values have been derived using two different methods, and it was found that extinction (especially in the UV) derived from Raman lidar data significantly surpasses the extinction derived from photometer AOD profiles. Airborne lidar data shows volume depolarization values to be less than 2.5% between 500 m and 2.5 km altitude, hence, particles in this range can be assumed to be of spherical shape. In-situ particle number concentrations measured at the Zeppelin Mountain Research Station at 474 m altitude peak at about 0.18 mu m diameter, which was also found for the microphysical inversion calculations performed at 850 m and 1500 m altitude. Number concentrations depend on the assumed extinction values, and slightly decrease with altitude as well as the effective particle diameter. A low imaginary part in the derived refractive index suggests weakly absorbing aerosols, which is confirmed by low black carbon concentrations, measured at the Zeppelin Mountain as well as on board the Polar 5 aircraft.
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.
In this study we present iterative methods using rational approximations, e.g... Pade approximants, which work very well for strongly ill-conditioned systems. In principle all methods of the family are convergent. One type of those methods has the advantage that their convergence behavior is very fast without additional a-priori information on the optimal relaxation parameter. (c) 2005 Elsevier Inc. All rights reserved
In this study we present iterative regularization methods using rational approximations, in particular, Pade approximants, which work well for ill-posed problems. We prove that the (k,j)-Pade method is a convergent and order optimal iterative regularization method in using the discrepancy principle of Morozov. Furthermore, we present a hybrid Pade method, compare it with other well-known methods and found that it is faster than the Landweber method. It is worth mentioning that this study is a completion of the paper [A. Kirsche, C. Bockmann, Rational approximations for ill-conditioned equation systems, Appl. Math. Comput. 171 (2005) 385-397] where this method was treated to solve ill-conditioned equation systems. (c) 2006 Elsevier Inc. All rights reserved.
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 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.
We present a summary on the current status of two inversion algorithms that are used in EARLINET (European Aerosol Research Lidar Network) for the inversion of data collected with EARLINET multiwavelength Raman lidars. These instruments measure backscatter coefficients at 355, 532, and 1064 nm, and extinction coefficients at 355 and 532 nm. Development of these two algorithms started in 2000 when EARLINET was founded. The algorithms are based on a manually controlled inversion of optical data which allows for detailed sensitivity studies. The algorithms allow us to derive particle effective radius as well as volume and surface area concentration with comparably high confidence. The retrieval of the real and imaginary parts of the complex refractive index still is a challenge in view of the accuracy required for these parameters in climate change studies in which light absorption needs to be known with high accuracy. It is an extreme challenge to retrieve the real part with an accuracy better than 0.05 and the imaginary part with accuracy better than 0.005-0.1 or +/- 50 %. Single-scattering albedo can be computed from the retrieved microphysical parameters and allows us to categorize aerosols into high-and low-absorbing aerosols. On the basis of a few exemplary simulations with synthetic optical data we discuss the current status of these manually operated algorithms, the potentially achievable accuracy of data products, and the goals for future work. One algorithm was used with the purpose of testing how well microphysical parameters can be derived if the real part of the complex refractive index is known to at least 0.05 or 0.1. The other algorithm was used to find out how well microphysical parameters can be derived if this constraint for the real part is not applied. The optical data used in our study cover a range of Angstrom exponents and extinction-to-backscatter (lidar) ratios that are found from lidar measurements of various aerosol types. We also tested aerosol scenarios that are considered highly unlikely, e.g. the lidar ratios fall outside the commonly accepted range of values measured with Raman lidar, even though the underlying microphysical particle properties are not uncommon. The goal of this part of the study is to test the robustness of the algorithms towards their ability to identify aerosol types that have not been measured so far, but cannot be ruled out based on our current knowledge of aerosol physics. We computed the optical data from monomodal logarithmic particle size distributions, i.e. we explicitly excluded the more complicated case of bimodal particle size distributions which is a topic of ongoing research work. Another constraint is that we only considered particles of spherical shape in our simulations. We considered particle radii as large as 7-10 mu m in our simulations where the Potsdam algorithm is limited to the lower value. We considered optical-data errors of 15% in the simulation studies. We target 50% uncertainty as a reasonable threshold for our data products, though we attempt to obtain data products with less uncertainty in future work.