Most of the quantitative information about the magnetic field vector in solar prominences comes from the analysis of the Hanle effect acting on lines formed by scattering. As these lines can be of non-negligible optical thickness, it is of interest to study the line formation process further. We investigate the multidimensional effects on the interpretation of spectropolarimetric observations, particularly on the inference of the magnetic field vector. We do this by analyzing the differences between multidimensional models, which involve fully self-consistent radiative transfer computations in the presence of spatial inhomogeneities and velocity fields, and those which rely on simple one-dimensional geometry. We study the formation of a prototype line in ad hoc inhomogeneous, isothermal 2D prominence models. We solve the NLTE polarized line formation problem in the presence of a large-scale oriented magnetic field. The resulting polarized line profiles are then interpreted (i.e. inverted) assuming a simple 1D slab model. We find that differences between input and the inferred magnetic field vector are non-negligible. Namely, we almost universally find that the inferred field is weaker and more horizontal than the input field. Spatial inhomogeneities and radiative transfer have a strong effect on scattering line polarization in the optically thick lines. In real-life situations, ignoring these effects could lead to a serious misinterpretation of spectropolarimetric observations of chromospheric objects such as prominences.

State-of-the-art methods in multidimensional NLTE radiative transfer are based on the use of local approximate lambda operator within either Jacobi or Gauss-Seidel iterative schemes. Here we propose another approach to the solution of 2D NLTE RT problems, Forth-and-Back Implicit Lambda Iteration (FBILI), developed earlier for 1D geometry. In order to present the method and examine its convergence properties we use the well-known instance of the two-level atom line formation with complete frequency redistribution. In the formal solution of the RT equation we employ short characteristics with two-point algorithm. Using an implicit representation of the source function in the computation of specific intensities, we compute and store the coefficients of the linear relations J = a + bS between the mean intensity J and the corresponding source function S. The use of iteration factors in the 'local' coefficients of these implicit relations in two 'inward' directions, along with the update of the source function in other two, 'outward', directions leads to four times faster solution than the Jacobi's one. Moreover, the update made in all four consecutive sweeps of the grid leads to an acceleration by a factor of 6-7 compared to the Jacobi iterative scheme.

Polarized NLTE radiative transfer in the presence of scattering in spectral lines and/or in continua may be cast in a so-called reduced form for six reduced components of the radiation field. In this formalism the six components of the reduced source function are angle-independent quantities. It thus reduces drastically the storage requirement of numerical codes. This approach encounters a fundamental problem when the medium is illuminated by a polarized incident radiation, because there is a priori no way of relating the known (and measurable) Stokes parameters of the incident radiation to boundary conditions for the reduced equations. The origin of this problem is that there is no unique way of deriving the radiation reduced components from its Stokes parameters (only the inverse operation is clearly defined). The method proposed here aims at enabling to work with arbitrary incident radiation field (polarized or unpolarized). In previous works an ad-hoc treatment of the boundary conditions, applying to cases where the incident radiation is unpolarized, has been used. In this note we show that it is possible to account for the incident radiation in a rigorous way, without any assumption on its properties, by expressing the radiation field as the sum of a directly transmitted radiation and of a so-called diffuse radiation. The diffuse radiation field obeys a transfer equation with no incident radiation that may be solved in the reduced form. The first scattering of the incident radiation introduces primary creation terms in the six components of the reduced source function. Once the reduced polarized transfer problem is solved for the diffuse radiation field, its Stokes parameters can be computed. We perform numerical computations of such cases, showing that the emergent line-polarization may be significantly affected by the polarization of the incident radiation.