Search

A record-breaking east Asian dust storm over recent years occurred in March 2021. The Ångström exponent (AE), which measures the wavelength dependence of aerosol optical thickness (AOT), is significantly sensitive to large aerosols such as dust. Due to the lack of observations during dust storms and the accuracy of the satellite-retrieved AE depending on the instrument and retrieval algorithm, it is possible to estimate the dust storm emission using the time-lagged ground-based AE observations. In this study, the hourly AEs observed by the Aerosol Robotic Network (AERONET) are assimilated with the ensemble Kalman smoother (EnKS) and Weather Research and Forecasting model coupled with Chemistry (WRF-Chem) to optimize simulated dust emissions from 14 to 23 March 2021. The results demonstrate that the additional inclusion of AE can optimize the size distribution of dust emissions and the associated total flux depending on the covariance between time-lagged AE observations and simulated dust emissions in each size bin. Compared to the experiment only assimilating AOT, validation by independent observations from the Skynet Observation NETwork (SONET) shows that assimilating additional AE information reduces the root mean square error (RMSE) of simulated AOT and AE by approximately 17 % and 61 %, respectively. The temporal variation in both simulated AOT and AE is improved through assimilating additional AE information. The assimilation of AOT and AE also makes the magnitude and variations in aerosol vertical extinctions more comparable to the independent Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) observations in both westward and eastward pathways of dust transport. The optimized dust emissions in the Gobi Desert during this period is estimated to be 52.63 Tg and reached a peak value of 3837 kt h−1 at 07:00 UTC on 14 March.

To retrieve columnar intensive aerosol properties from sun–sky photometers, both irradiance and radiance calibration factors are needed. For the irradiance the solar calibration constant, V0, which denotes the instrument counts for a direct normal solar flux extrapolated to the top of the atmosphere, must be determined. The solid view angle, SVA, is a measure of the field of view of the instrument, and it is important for obtaining the radiance from sky diffuse irradiance measurements. Each of the three sun-photometer networks considered in the present study (SKYNET, AERONET, WMO GAW) adopts different protocols of calibration, and we evaluate the performance of the on-site calibration procedures, applicable to every kind of sun–sky photometer but tested in this analysis only on SKYNET Prede POM01 instruments, during intercomparison campaigns and laboratory calibrations held in the framework of the Metrology for Aerosol Optical Properties (MAPP) European Metrology Programme for Innovation and Research (EMPIR) project. The on-site calibration, performed as frequently as possible (ideally monthly) to monitor changes in the device conditions, allows operators to track and evaluate the calibration status on a continuous basis, considerably reducing the data gaps incurred by the periodic shipments for performing centralized calibrations. The performance of the on-site calibration procedures for V0 was very good at sites with low turbidity, showing agreement with a reference calibration between 0.5 % and 1.5 % depending on wavelengths. In the urban area, the agreement decreases between 1.7 % and 2.5 %. For the SVA the difference varied from a minimum of 0.03 % to a maximum of 3.46 %.

In this study, we explored the impacts of radiation damage and chemical composition on the etching time of fission tracks in monazite. Despite the potential of monazite fission-track (MFT) dating as an ultralow-temperature thermochronology, the comprehensive effects of radiation damage and non-formula elements, especially on the etching rate of MFTs, remain unexplored, and established analytical procedures are lacking. We quantified the degree of radiation damage (ΔFHWM) in Cretaceous to Quaternary monazites distributed in the Japan arc through Raman spectroscopy and chemical composition analyses. Subsequently, MFT etching was performed to examine the correlation between these parameters and the etching time. Estimation of the degree of radiation damage showed an increase in radiation damage corresponding to the cooling age of each geological unit. For example, the Toya ignimbrite (ca. 0.1 Ma) and the Kurobegawa granodiorite (< 0.8 Ma), both of which are types of monazite from Quaternary geological units, have ΔFHWM values of 0.27 and 0.55 cm−1, respectively. In contrast, the Muro ignimbrite (ca. 15 Ma) has a ΔFHWM value of 4.01 cm−1, while the Kibe granite and the Sagawa granite, both of which are Cretaceous granitoids, yielded 7.35 and 6.31 cm−1, respectively. MFT etching of these samples according to the existing recipe (6 M HCl at 90 °C for 60–90 min) was completed at 1200, 860, 210, 120, and 90 min for the Toya ignimbrite, Kurobegawa granodiorite, Muro ignimbrite, Sagawa granite, and Kibe granite, respectively. These outcomes highlight an inverse relationship between MFT etching time and the degree of radiation damage in monazite, while the correlation between MFT etching time and chemical composition was unclear. The results affirm earlier considerations that the etching rate of MFTs is strongly influenced by radiation damage. Conversely, young samples with lower levels of radiation damage exhibit higher chemical resistance, suggesting that existing etching recipes may not adequately etch MFTs.

Abstract Earthquakes with magnitudes of 6–7 have been reported even in various active tectonic settings where fault deformation topography have not been detected. Therefore, delineating concealed active faults generating such earthquakes is necessary to reduce earthquake damage; however, few studies exist to provide its clues regarding such faults. The 1984 Western Nagano Earthquake in Japan was a main shock with a magnitude of Mj 6.8 and a depth of 2 km at the source. Solid bedrocks are well‐exposed in the earthquake source region; however, no surface rupture have been identified, and the active fault is known to be concealed. In this study, we collected data on striations observed in fractures by geological survey around the source area of the 1984 Western Nagano Earthquake. Using the collected data, the multiple inverse method was used to estimate the stresses that affected the striation formation. Consequently, stresses similar to acting faults in this area were detected in minor faults around the known concealed active fault. This suggests that the minor faults might be part of the damage zone that has been developed around the concealed active fault. Some minor faults were recognized in Quaternary volcanic rocks, confirming that they experienced displacements recently. This study indicates the possibility of detecting concealed active faults in the bedrock by geological survey.

Abstract Feedback control of qubits is a highly demanded technique for advanced quantum information protocols such as fault-tolerant quantum error correction. Here we demonstrate active reset of a silicon spin qubit using feedback control. The active reset is based on quantum non-demolition (QND) readout of the qubit and feedback according to the readout results, which is enabled by hardware data processing and sequencing. We incorporate a cumulative readout technique to the active reset protocol, enhancing initialization fidelity above a limitation imposed by the single-shot QND readout fidelity. An analysis of the reset protocol implies a pathway to achieve the initialization fidelity sufficient for fault-tolerant quantum computation. These results provide a practical approach to high-fidelity qubit operations in realistic devices.

Abstract Electron spins in semiconductor quantum dots have been intensively studied for implementing quantum computation and high-fidelity single- and two-qubit operations have recently been achieved. Quantum teleportation is a three-qubit protocol exploiting quantum entanglement and it serves as an essential primitive for more sophisticated quantum algorithms. Here we demonstrate a scheme for quantum teleportation based on direct Bell measurement for a single-electron spin qubit in a triple quantum dot utilizing the Pauli exclusion principle to create and detect maximally entangled states. The single spin polarization is teleported from the input qubit to the output qubit. We find this fidelity is primarily limited by singlet–triplet mixing, which can be improved by optimizing the device parameters. Our results may be extended to quantum algorithms with a larger number of semiconductor spin qubits.

The anthropogenic emission of sulfur dioxide (SO2) over China has significantly declined as a consequence of the clean air actions. In this study, we have developed a new emission inversion system based on a four-dimensional local ensemble transform Kalman filter (4D-LETKF) and the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem) to dynamically update the SO2 emission grid by grid over China by assimilating the ground-based hourly SO2 observations. Sensitivity tests for the assimilation system have been conducted firstly to tune four system parameters: ensemble size, horizontal and temporal localization lengths, and perturbation size. Our results reveal that the same random perturbation factors used throughout the whole model grids with assimilating observations within about 180 km can efficiently optimize the SO2 emission, whereas the ensemble size has only little effect. The temporal localization by assimilating only the subsequent hourly observations can reveal the diurnal variation of the SO2 emission, which is better than updating the magnitude of SO2 emission every 12 h by assimilating all the observations within the 12 h window. The inverted SO2 emission over China in November 2016 has declined by an average of 49.4 % since 2010, which is well in agreement with the bottom-up estimation of 48.0 %. Larger reductions of SO2 emission are found over the a priori higher source regions such as the Yangtze River Delta (YRD). The simulated SO2 surface mass concentrations using two distinguished chemical reaction mechanisms are both much more comparable to the observations with the newly inverted SO2 emission than those with the a priori emission. These indicate that the newly developed emission inversion system can efficiently update the SO2 emissions based on the routine surface SO2 observations. The reduced SO2 emission induces the sulfate and PM2.5 surface concentrations to decrease by up to 10 µg m−3 over central China.

This paper is an overview of the progress in sky radiometer technology and the development of the network called SKYNET. It is found that the technology has produced useful on-site calibration methods, retrieval algorithms, and data analyses from sky radiometer observations of aerosol, cloud, water vapor, and ozone. A formula was proposed for estimating the accuracy of the sky radiometer calibration constant F0 using the improved Langley (IL) method, which was found to be a good approximation to observed monthly mean uncertainty in F0, around 0.5 % to 2.4 % at the Tokyo and Rome sites and smaller values of around 0.3 % to 0.5 % at the mountain sites at Mt. Saraswati and Davos. A new cross IL (XIL) method was also developed to correct an underestimation by the IL method in cases with large aerosol retrieval errors. The root-mean-square difference (RMSD) in aerosol optical thickness (AOT) comparisons with other networks took values of less than 0.02 for λ≥500 nm and a larger value of about 0.03 for shorter wavelengths in city areas and smaller values of less than 0.01 in mountain comparisons. Accuracies of single-scattering albedo (SSA) and size distribution retrievals are affected by the propagation of errors in measurement, calibrations for direct solar and diffuse sky radiation, ground albedo, cloud screening, and the version of the analysis software called the Skyrad pack. SSA values from SKYNET were up to 0.07 larger than those from AERONET, and the major error sources were identified as an underestimation of solid viewing angle (SVA) and cloud contamination. Correction of these known error factors reduced the SSA difference to less than 0.03. Retrievals of other atmospheric constituents by the sky radiometer were also reviewed. Retrieval accuracies were found to be about 0.2 cm for precipitable water vapor amount and 13 DU (Dobson Unit) for column ozone amount. Retrieved cloud optical properties still showed large deviations from validation data, suggesting a need to study the causes of the differences. It is important that these recent studies on improvements presented in the present paper are introduced into the existing operational systems and future systems of the International SKYNET Data Center.

High-performance computing resources allow us to conduct numerical simulations with a horizontal grid spacing that is sufficiently high to resolve cloud systems on a global scale, and high-resolution models (HRMs) generally provide better simulation performance than low-resolution models (LRMs). In this study, we execute a next-generation model that is capable of simulating global aerosols using version 16 of the Nonhydrostatic Icosahedral Atmospheric Model (NICAM.16). The simulated aerosol distributions are obtained for 3 years with an HRM using a global 14 km grid spacing, an unprecedentedly high horizontal resolution and long integration period. For comparison, a NICAM with a 56 km grid spacing is also run as an LRM, although this horizontal resolution is still high among current global aerosol climate models. The comparison elucidated that the differences in the various variables of meteorological fields, including the wind speed, precipitation, clouds, radiation fluxes and total aerosols, are generally within 10 % of their annual averages, but most of the variables related to aerosols simulated by the HRM are slightly closer to the observations than are those simulated by the LRM. Upon investigating the aerosol components, the differences in the water-insoluble black carbon and sulfate concentrations between the HRM and LRM are large (up to 32 %), even in the annual averages. This finding is attributed to the differences in the aerosol wet deposition flux, which is determined by the conversion rate of cloud to precipitation, and the difference between the HRM and LRM is approximately 20 %. Additionally, the differences in the simulated aerosol concentrations at polluted sites during polluted months between the HRM and LRM are estimated with normalized mean biases of −19 % for black carbon (BC), −5 % for sulfate and −3 % for the aerosol optical thickness (AOT). These findings indicate that the impacts of higher horizontal grid spacings on model performance for secondary products such as sulfate, and complex products such as the AOT, are weaker than those for primary products, such as BC. On a global scale, the subgrid variabilities in the simulated AOT and cloud optical thickness (COT) in the 1∘×1∘ domain using 6-hourly data are estimated to be 28.5 % and 80.0 %, respectively, in the HRM, whereas the corresponding differences are 16.6 % and 22.9 % in the LRM. Over the Arctic, both the HRM and the LRM generally reproduce the observed aerosols, but the largest difference in the surface BC mass concentrations between the HRM and LRM reaches 30 % in spring (the HRM-simulated results are closer to the observations). The vertical distributions of the HRM- and LRM-simulated aerosols are generally close to the measurements, but the differences between the HRM and LRM results are large above a height of approximately 3 km, mainly due to differences in the wet deposition of aerosols. The global annual averages of the effective radiative forcings due to aerosol–radiation and aerosol–cloud interactions (ERFari and ERFaci) attributed to anthropogenic aerosols in the HRM are estimated to be -0.293±0.001 and -0.919±0.004 W m−2, respectively, whereas those in the LRM are -0.239±0.002 and -1.101±0.013 W m−2. The differences in the ERFari between the HRM and LRM are primarily caused by those in the aerosol burden, whereas the differences in the ERFaci are primarily caused by those in the cloud expression and performance, which are attributed to the grid spacing. The analysis of interannual variability revealed that the difference in reproducibility of both sulfate and carbonaceous aerosols at different horizontal resolution is greater than their interannual variability over 3 years, but those of dust and sea salt AOT and possibly clouds were the opposite. Because at least 10 times the computer resources are required for the HRM (14 km grid) compared to the LRM (56 km grid), these findings in this study help modelers decide whether the objectives can be achieved using such higher resolution or not under the limitation of available computational resources.

The Prede sky radiometer measures direct solar irradiance and the angular distribution of diffuse radiances at the ultraviolet, visible, and near-infrared wavelengths. These data are utilized for the remote sensing of aerosols, water vapor, ozone, and clouds, but the calibration constant, which is the sensor output current of the extraterrestrial solar irradiance at the mean distance between Earth and the Sun, is needed. The aerosol channels, which are the weak gas absorption wavelengths of 340, 380, 400, 500, 675, 870, and 1020 nm, can be calibrated by an on-site self-calibration method, the Improved Langley method. This on-site self-calibration method is useful for the continuous long-term observation of aerosol properties. However, the continuous long-term observation of precipitable water vapor (PWV) by the sky radiometer remains challenging because calibrating the water vapor absorption channel of 940 nm generally relies on the standard Langley (SL) method at limited observation sites (e.g., the Mauna Loa Observatory) and the transfer of the calibration constant by a side-by-side comparison with the reference sky radiometer calibrated by the SL method. In this study, we developed the SKYMAP algorithm, a new on-site method of self-calibrating the water vapor channel of the sky radiometer using diffuse radiances normalized by direct solar irradiance (normalized radiances). Because the sky radiometer measures direct solar irradiance and diffuse radiance using the same sensor, the normalization cancels the calibration constant included in the measurements. The SKYMAP algorithm consists of three steps. First, aerosol optical and microphysical properties are retrieved using direct solar irradiances and normalized radiances at aerosol channels. The aerosol optical properties at the water vapor channel are interpolated from those at aerosol channels. Second, PWV is retrieved using the angular distribution of the normalized radiances at the water vapor channel. Third, the calibration constant at the water vapor channel is estimated from the transmittance of PWV and aerosol optical properties. Intensive sensitivity tests of the SKYMAP algorithm using simulated data of the sky radiometer showed that the calibration constant is retrieved reasonably well for PWV cm, which indicates that the SKYMAP algorithm can calibrate the water vapor channel on-site in dry conditions. Next, the SKYMAP algorithm was applied to actual measurements under the clear-sky and low-PWV ( cm) conditions at two sites, Tsukuba and Chiba, Japan, and the annual mean calibration constants at the two sites were determined. The SKYMAP-derived calibration constants were 10.1 % and 3.2 % lower, respectively, than those determined by a side-by-side comparison with the reference sky radiometer. After determining the calibration constant, we obtained PWV from the direct solar irradiances in both the dry and wet seasons. The retrieved PWV values corresponded well to those derived from a global-navigation-satellite-system–global-positioning-system receiver, a microwave radiometer, and an AERONET (Aerosol Robotic Network) sun–sky radiometer at both sites. The correlation coefficients were greater than 0.96. We calculated the bias errors and the root mean square errors by comparing PWV between the DSRAD (direct solar irradiance) algorithm and other instruments. The magnitude of the bias error and the root mean square error were and  cm for PWV cm, respectively. However, our method tended to underestimate PWV in the wet conditions, and the magnitude of the bias error and the root mean square error became large, and  cm for PWV>3 cm, respectively. This problem was mainly due to the overestimation of the aerosol optical thickness before the retrieval of PWV. These results show that the SKYMAP algorithm enables us to observe PWV over the long term, based on its unique on-site self-calibration method.

Conventional qubit readout methods in silicon spin qubits destroy the quantum state, precluding any further computations based on the outcome. Here, the authors demonstrate quantum non-demolition readout using a second qubit of the same kind, making for a scalable approach.

This study analyzed the aerosol optical properties derived by SKYRAD.pack versions 5.0 and 4.2 (referred to as V5.0 and V4.2) using the radiometer measurements over Qionghai and Yucheng in China, two new sites of the sky radiometer network (SKYNET). As V5.0 uses an a priori size distribution function (SDF) of a bimodal log-normal function, the volume size distribution retrieved by V5.0 presented bimodal patterns with a 0.1–0.2 µm fine particle mode and a 3.0–6.0 µm coarse particle mode both over Qionghai and Yucheng. The differences in the volume size distributions between the two versions were very large for the coarse mode with a radius of over 5 µm. The single scattering albedos (SSAs) by V5.0 correlated with SSAs by V4.2 with R=0.88, 0.87, 0.90, 0.88, and 0.92 at wavelengths of 400, 500, 670, 870, and 1020 nm over Qionghai, respectively. The correlation coefficients were around 0.95, 0.95, 0.96, 0.94, and 0.91 at the five channels in Yucheng. An error of ±5 % for the solid view angle (SVA) introduced about ±2 % differences in retrieved SSA values both by V4.2 and V5.0. An error of ±50 % for ground surface albedo (Ag) caused about 1 % averaged differences in retrieved SSA values by the two versions. With the atmospheric pressure (PRS) increased by 1 %, 2 %, 3 %, and 4 %, the averaged changes in SSAs did not exceed 0.8 % both by V4.2 and V5.0. The SSA differences at 500 nm between the two versions decreased, while aerosol optical depths (AODs) increased over both sites. The seasonal variability of the aerosol properties over Qionghai and Yucheng was investigated based on SKYRAD.pack V5.0. The seasonal averaged AOD over Qionghai had higher values in spring, winter, and autumn and lower values in summer. The AOD averages were commonly higher in summer and spring than in winter and autumn in Yucheng. The lowest seasonal averaged SSAs were both observed in winter at the two sites. The fraction of the fine aerosol particles was much smaller in summer than in other seasons over Qionghai; the volume fraction of the coarse-mode particle in Yucheng had much larger values compared to the fine-mode particle in all seasons. The validation results provide valuable references for continued improvement of the retrieval algorithms of SKYNET and other aerosol observational networks.

Great efforts have been made to simulate atmospheric pollutants, but their spatial and temporal distributions are still highly uncertain. Observations can measure their concentrations with high accuracy but cannot estimate their spatial distributions due to the sporadic locations of sites. Here, we propose an ensemble method by applying a linear minimum variance estimation (LMVE) between multi-model ensemble (MME) simulations and measurements to derive a more realistic distribution of atmospheric pollutants. The LMVE is a classical and basic version of data assimilation, although the estimation itself is still useful for obtaining the best estimates by combining simulations and observations without a large amount of computer resources, even for high-resolution models. In this study, we adopt the proposed methodology for atmospheric radioactive caesium (Cs-137) in atmospheric particles emitted from the Fukushima Daiichi Nuclear Power Station (FDNPS) accident in March 2011. The uniqueness of this approach includes (1) the availability of observed Cs-137 concentrations near the surface at approximately 100 sites, thus providing dense coverage over eastern Japan; (2) the simplicity of identifying the emission source of Cs-137 due to the point source of FDNPS; (3) the novelty of MME with the high-resolution model (3 km horizontal grid) over complex terrain in eastern Japan; and (4) the strong need to better estimate the Cs-137 distribution due to its inhalation exposure among residents in Japan. The ensemble size is six, including two atmospheric transport models: the Weather Research and Forecasting – Community Multi-scale Air Quality (WRF-CMAQ) model and non-hydrostatic icosahedral atmospheric model (NICAM). The results showed that the MME that estimated Cs-137 concentrations using all available sites had the lowest geometric mean bias (GMB) against the observations (GMB =1.53), the lowest uncertainties based on the root mean square error (RMSE) against the observations (RMSE =9.12 Bq m−3), the highest Pearson correlation coefficient (PCC) with the observations (PCC =0.59) and the highest fraction of data within a factor of 2 (FAC2) with the observations (FAC2 =54 %) compared to the single-model members, which provided higher biases (GMB =1.83–4.29, except for 1.20 obtained from one member), higher uncertainties (RMSE =19.2–51.2 Bq m−3), lower correlation coefficients (PCC =0.29–0.45) and lower precision (FAC2 =10 %–29 %). At the model grid, excluding the measurements, the MME-estimated Cs-137 concentration was estimated by a spatial interpolation of the variance used in the LMVE equation using the inverse distance weights between the nearest two sites. To test this assumption, the available measurements were divided into two categories, i.e. learning and validation data; thus, the assumption for the spatial interpolation was found to guarantee a moderate PCC value (> 0.4) within an approximate distance of at least 70 km. Extra sensitivity tests for several parameters, i.e. the site number and the weighting coefficients in the spatial interpolation, the time window in the LMVE and the ensemble size, were performed. In conclusion, the important assumptions were the time window and the ensemble size; i.e. a shorter time window (the minimum in this study was 1 h, which is the observation interval) and a larger ensemble size (the maximum in this study was six, but five is also acceptable if the members are effectively selected) generated better results.

Aerosol vertical information is critical to quantify the influences of aerosol on the climate and environment; however, large uncertainties still persist in model simulations. In this study, the vertical aerosol extinction coefficients from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) onboard the Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) are assimilated to optimize the hourly aerosol fields of the Non-hydrostatic ICosahedral Atmospheric Model (NICAM) online coupled with the Spectral Radiation Transport Model for Aerosol Species (SPRINTARS) using a four-dimensional local ensemble transform Kalman filter (4-D LETKF). A parallel assimilation experiment using bias-corrected aerosol optical thicknesses (AOTs) from the Moderate Resolution Imaging Spectroradiometer (MODIS) is conducted to investigate the effects of assimilating the observations (and whether to include vertical information) on the model performances. Additionally, an experiment simultaneously assimilating both CALIOP and MODIS observations is conducted. The assimilation experiments are successfully performed for 1 month, making it possible to evaluate the results in a statistical sense. The hourly analyses are validated via both the CALIOP-observed aerosol vertical extinction coefficients and the AOT observations from MODIS and the AErosol RObotic NETwork (AERONET). Our results reveal that both the CALIOP and MODIS assimilations can improve the model simulations. The CALIOP assimilation is superior to the MODIS assimilation in modifying the incorrect aerosol vertical distributions and reproducing the real magnitudes and variations, and the joint CALIOP and MODIS assimilation can further improve the simulated aerosol vertical distribution. However, the MODIS assimilation can better reproduce the AOT distributions than the CALIOP assimilation, and the inclusion of the CALIOP observations has an insignificant impact on the AOT analysis. This is probably due to the nadir-viewing CALIOP having much sparser coverage than MODIS. The assimilation efficiencies of CALIOP decrease with increasing distances of the overpass time, indicating that more aerosol vertical observation platforms are required to fill the sensor-specific observation gaps and hence improve the aerosol vertical data assimilation.

Manipulating topological spin textures are demanded for future spintronic devices, but knowledge about phase transitions among different spin textures remain limited. Here, Fujishiro and Kanazawa et al. report chemical-pressure-controlled phase transitions between different topological spin textures in chiral magnets MnSi1−x Ge x .

Mutual control of the electric polarization and magnetization promises for low power consumption spintronic devices but remains challenging. Here the authors show reversal of non-volatile magnetization by electric field as well as the polarization switching by magnetic field in a single-component material, close to room temperature.

In this study, we investigate the feasibility of a multi-pixel scheme in the inversion of aerosol optical properties for multispectral satellite instruments over the ocean. Different from the traditional satellite aerosol retrievals conducted pixel by pixel, we derive the aerosol optical thickness (AOT) of multiple pixels simultaneously by adding a smoothness constraint on the spatial variation of aerosols and oceanic substances, which helps the satellite retrieval, with higher consistency from pixel to pixel. Simulations are performed for two representative oceanic circumstances, open and coastal waters, as well as the land–ocean interface region. We retrieve the AOT for fine, sea spray, and dust aerosols simultaneously using synthetic spectral measurements, which are from the Greenhouse Gases Observing Satellite and Thermal and Near Infrared Sensor for Carbon Observation – Cloud and Aerosol Imager (GOSAT∕TANSO-CAI), with four wavelengths ranging from the ultraviolet to shortwave infrared bands. The forward radiation calculation is performed by a coupled atmosphere–ocean radiative transfer model combined with a three-component bio-optical oceanic module, where the chlorophyll a concentration, sediment, and colored dissolved organic matter are considered. Results show that accuracies of the derived AOT and spectral remote-sensing reflectance are both improved by applying smoothness constraints on the spatial variation of aerosol and oceanic substances in homogeneous or inhomogeneous surface conditions. The multi-pixel scheme can be effective in compensating for the retrieval biases induced by measurement errors and improving the retrieval sensitivity, particularly for the fine aerosols over the coastal water. We then apply the algorithm to derive AOTs using real satellite measurements. Results indicate that the multi-pixel method helps to polish the irregular retrieved results of the satellite imagery and is potentially promising for the aerosol retrieval over highly turbid waters by benefiting from the coincident retrieval of neighboring pixels. A comparison of retrieved AOTs from satellite measurements with those from the Aerosol Robotic Network (AERONET) also indicates that retrievals conducted by the multi-pixel scheme are more consistent with the AERONET observations.

The race to produce a quantum computer has driven the development of many different qubit designs with different benefits and drawbacks. Noiri et al. demonstrate a hybrid device with two coupled semiconductor spin qubits of different designs, which should allow each qubit’s advantages to be exploited.

Retrieval of aerosol optical properties and water-leaving radiance over ocean is challenging since the latter mostly accounts for ∼ 10 % of the satellite-observed signal and can be easily influenced by the atmospheric scattering. Such an effort would be more difficult in turbid coastal waters due to the existence of optically complex oceanic substances or high aerosol loading. In an effort to solve such problems, we present an optimization approach for the simultaneous determination of aerosol optical thickness (AOT) and normalized water-leaving radiance (nLw) from multispectral satellite measurements. In this algorithm, a coupled atmosphere–ocean radiative transfer model combined with a comprehensive bio-optical oceanic module is used to jointly simulate the satellite-observed reflectance at the top of atmosphere and water-leaving radiance just above the ocean surface. Then, an optimal estimation method is adopted to retrieve AOT and nLw iteratively. The algorithm is validated using Aerosol Robotic Network – Ocean Color (AERONET-OC) products selected from eight OC sites distributed over different waters, consisting of observations that covered glint and non-glint conditions from the Moderate Resolution Imaging Spectroradiometer (MODIS) instrument. Results show a good consistency between retrieved and in situ measurements at each site. It is demonstrated that more accurate AOTs are determined based on the simultaneous retrieval method, particularly in shorter wavelengths and sunglint conditions, where the averaged percentage difference (APD) of retrieved AOT is generally reduced by approximate 10 % in visible bands compared with those derived from the standard atmospheric correction (AC) scheme, since all the spectral measurements can be used jointly to increase the information content in the inversion of AOT, and the wind speed is also simultaneously retrieved to compensate the specular reflectance error estimated from the rough ocean surface model. For the retrieval of nLw, atmospheric overcorrection can be avoided in order to have a significant improvement of the inversion of nLw at 412 nm. Furthermore, generally better estimates of band ratios of nLw(443) / nLw(554) and nLw(488) / nLw(554) are obtained using the simultaneous retrieval approach with lower root mean square errors and relative differences than those derived from the standard AC approach in comparison to the AERONET-OC products, as well as the APD values of retrieved Chl which decreased by about 5 %. On the other hand, the standard AC scheme yields a more accurate retrieval of nLw at 488 nm, prompting a further optimization of the oceanic bio-optical module of the current model.

The Tibetan Plateau (TP) is located at the juncture of several important natural and anthropogenic aerosol sources. Satellites have observed substantial dust and anthropogenic aerosols in the atmosphere during summer over the TP. These aerosols have distinct effects on the earth's energy balance, microphysical cloud properties, and precipitation rates. To investigate the transport of summer dust and anthropogenic aerosols over the TP, we combined the Spectral Radiation-Transport Model for Aerosol Species (SPRINTARS) with a non-hydrostatic regional model (NHM). The model simulation shows heavily loaded dust aerosols over the northern slope and anthropogenic aerosols over the southern slope and the east of the TP. The dust aerosols are primarily mobilized around the Taklimakan Desert, where a portion of the aerosols are transported eastward due to the northwesterly current; simultaneously, a portion of the particles are transported southward when a second northwesterly current becomes northeasterly because of the topographic blocking of the northern slope of the TP. Because of the strong upward current, dust plumes can extend upward to approximately 7–8 km a.s.l. over the northern slope of the TP. When a dust event occurs, anthropogenic aerosols that entrained into the southwesterly current via the Indian summer monsoon are transported from India to the southern slope of the TP. Simultaneously, a large amount of anthropogenic aerosol is also transported from eastern China to the east of the TP by easterly winds. An investigation on the transport of dust and anthropogenic aerosols over the plateau may provide the basis for determining aerosol impacts on summer monsoons and climate systems.

An aerosol-coupled global nonhydrostatic model with a stretched-grid system has been developed. Circulations over the global and target domains are simulated with a single model, which includes fine meshes covering the target region to calculate meso-scale circulations. The stretched global model involves lower computational costs to simulate atmospheric aerosols with fine horizontal resolutions compared with a global uniform nonhydrostatic model, whereas it may require higher computational costs compared with the general regional models, because the stretched-grid system calculates inside and outside the target domain. As opposed to general regional models, the stretched-grid system requires neither a nesting technique nor lateral boundary conditions. In this study, we developed a new-type regional model for the simulation of aerosols over Japan, especially in the Kanto areas surrounding Tokyo, with a maximum horizontal resolution of approximately 10 km. This model usually reproduces temporal variations and their averages of the observed weather around Japan. This model generally reproduces monthly mean distributions of the observed sulfate and SO2 over East Asia, with high correlations (R > 0.6), but the underestimation of the simulated concentrations by 40% (sulfate) and 50% (SO2). Their underestimation of the simulated sulfate and SO2 concentrations over East Asia are strongly affected by their underestimation in China and possibly by the uncertainty of the simulated precipitation around Japan. In the Kanto area, this model succeeds in simulating the wind patterns and the diurnal transitions around the center of the Kanto area, although it is inadequate to simulate the wind patterns and the diurnal transitions at some sites located at the edge of the Kanto area and surrounded on three sides by mountains, e.g., Maebashi, mainly due to the insufficient horizontal resolution. This model also generally reproduces both diurnal and synoptic variations of the observed and/or a regional aerosol-transport model, WRF-CMAQ, simulated EC, sulfate, and SO2 concentrations in the Kanto area, especially with their high correlation (R > 0.5) at Komae/Tokyo. Although the aerosol module used in this study is relatively simplified compared to the general regional aerosol models, this study reveals that our proposed model with the stretched-grid system can be applicable for the regional aerosol simulation.

Sun–sky radiometers are instruments created for aerosol study, but they can measure in the water vapour absorption band allowing the estimation of columnar water vapour in clear sky simultaneously with aerosol characteristics, with high temporal resolution. A new methodology is presented for estimating calibration parameters (i.e. characteristic parameters of the atmospheric transmittance and solar calibration constant) directly from the sun–sky radiometer measurements. The methodology is based on the hypothesis that characteristic parameters of the atmospheric transmittance are dependent on vertical profiles of pressure, temperature and moisture occurring at each site of measurement. To obtain the parameters from the proposed methodology some seasonal independent measurements of columnar water vapour taken over a large range of solar zenith angle simultaneously with the sun–sky radiometer measurements, are needed. In this work high time resolution columnar water vapour measurements by GPS were used as independent data set, but also the case when such measurements are not available was considered by developing the surface humidity method (SHM). This methodology makes it possible to retrieve the needed independent data set of columnar water vapour using the standard surface meteorological observations (temperature, pressure and relative humidity) more readily available. The time pattern of columnar water vapour from sun–sky radiometer retrieved using both the methodologies was compared with simultaneous measurements from microwave radiometer, radiosondings and GPS. Water vapour from sun–sky radiometer, obtained using GPS independent measurements, was characterized by an error varying from 1% up to 5%, whereas water vapour from SHM showed an error from 1% up to 11%, depending on the local columnar water occurring at the site during the year. These errors were estimated by comparing water vapour series from sun–sky radiometer against measurements taken by GPS at a nearby station. The accordance between retrievals from sun–sky radiometer and simultaneous measurements from the other instruments was found always within the error both in the case of SHM and of the GPS independent data set. Water vapour obtained using characteristic parameters of the atmospheric transmittance dependent on water vapour was also compared against GPS retrievals, showing a clear improvement with respect to the case when these parameters are kept fixed.

There is great uncertainty regarding the role of mineral dust aerosols in Earth's climate system. One reason for this uncertainty is that the optical properties of mineral dust, such as its single scattering albedo (the ratio of scattering to total extinction), are poorly constrained because ground observations are limited to a few locations and satellite standard products are not available due to the excessively bright surface of the desert in the visible wavelength, which makes robust retrievals difficult. Here, we develop a method to estimate the spatial distributions of the aerosol single scattering albedo (ω0) and optical depth (τa), with daily 1°×1° spatial resolution using data from the Moderate Resolution Imaging Spectroradiometer (MODIS) as well as model simulations of radiative transfer. This approach is based on the "critical surface reflectance" method developed in the literature, which estimates ω0 from the top of the atmospheric radiance. We estimate the uncertainties in ω0 over the Sahara (Asia) to be approximately 0.020 and 0.010 (0.023 and 0.017) for bands 9 and 1, respectively, while the uncertainty in τa is approximately 0.235 and 0.228 (0.464 and 0.370) for bands 9 and 1, respectively. The 5–95% range of the spatial distribution of ω0 over the Sahara (Asia) is approximately 0.90–0.94 and 0.96–0.99 (0.87–0.94 and 0.89–0.97) for bands 9 and 1, respectively, and that of τa over the Sahara (Asia) is approximately 0.8–1.4 and 0.8–1.7 (0.7–2.0 and 0.7–1.9) for bands 9 and 1, respectively. The results for the Sahara indicate a good correlation between ω0 and the surface reflectance, and between ω0 and τa. However, the relationships between ω0, τa, and surface reflectance are less clear in Asia than in the Sahara, and the ω0 values are smaller than those in the Sahara. The regions with small ω0 values are consistent with the regions where coal-burning smoke and carbonaceous aerosols are reported to be transported in previous studies. Because the coal-burning and carbonaceous aerosols are known to be more absorptive and have smaller ω0 values than dust aerosols, our results indicate that the dust aerosols in Asia are contaminated by these anthropogenic aerosols. The spatial distribution of dust optical properties obtained in our work could be useful in understanding the role of dust aerosols in Earth's climate system, most likely through future collaboration with regional and global modelling studies.

We present a validation study of Collection 5 MODIS level 2 Aqua and Terra AOT (aerosol optical thickness) and AE (Ångström exponent) over ocean by comparison to coastal and island AERONET (AErosol RObotic NETwork) sites for the years 2003–2009. We show that MODIS (MODerate-resolution Imaging Spectroradiometer) AOT exhibits significant biases due to wind speed and cloudiness of the observed scene, while MODIS AE, although overall unbiased, exhibits less spatial contrast on global scales than the AERONET observations. The same behaviour can be seen when MODIS AOT is compared against Maritime Aerosol Network (MAN) data, suggesting that the spatial coverage of our datasets does not preclude global conclusions. Thus, we develop empirical correction formulae for MODIS AOT and AE that significantly improve agreement of MODIS and AERONET observations. We show these correction formulae to be robust. Finally, we study random errors in the corrected MODIS AOT and AE and show that they mainly depend on AOT itself, although small contributions are present due to wind speed and cloud fraction in AOT random errors and due to AE and cloud fraction in AE random errors. Our analysis yields significantly higher random AOT errors than the official MODIS error estimate (0.03 + 0.05 τ), while random AE errors are smaller than might be expected. This new dataset of bias-corrected MODIS AOT and AE over ocean is intended for aerosol model validation and assimilation studies, but also has consequences as a stand-alone observational product. For instance, the corrected dataset suggests that much less fine mode aerosol is transported across the Pacific and Atlantic oceans.

In order to reduce uncertainty in the estimation of Direct Aerosol Radiative Forcing (DARF), it is important to improve the estimation of the single scattering albedo (SSA). In this study, we propose a new data processing method to improve SSA retrievals for the SKYNET sky radiometer network, which is one of the growing number of networks of sun-sky photometers, such as NASA AERONET and others. There are several reports that SSA values from SKYNET have a bias compared to those from AERONET, which is regarded to be the most accurate due to its rigorous calibration routines and data quality and cloud screening algorithms. We investigated possible causes of errors in SSA that might explain the known biases through sensitivity experiments using a numerical model, and also using real data at the SKYNET sites at Pune (18.616° N/73.800° E) in India and Beijing (39.586° N/116.229° E) in China. Sensitivity experiments showed that an uncertainty of the order of ±0.03 in the SSA value can be caused by a possible error in the ground surface albedo or solid view angle assumed for each observation site. Another candidate for possible error in the SSA was found in cirrus contamination generated by imperfect cloud screening in the SKYNET data processing. Therefore, we developed a new data quality control method to get rid of low quality or cloud contamination data, and we applied this method to the real observation data at the Pune site in SKYNET. After applying this method to the observation data, we were able to screen out a large amount of cirrus-contaminated data and to reduce the deviation in the SSA value from that of AERONET. We then estimated DARF using data screened by our new method. The result showed that the method significantly reduced the difference of 5 W m−2 that existed between the SKYNET and AERONET values of DARF before screening. The present study also suggests the necessity of preparing suitable a priori information on the distribution of coarse particles ranging in radius between 10 μm and 30 μm for the analysis of heavily dust-laden atmospheric cases.

The direct radiative forcing by sulfate aerosols is still uncertain, mainly because the uncertainties are largely derived from differences in sulfate column burdens and its vertical distributions among global aerosol models. One possible reason for the large difference in the computed values is that the radiative forcing delicately depends on various simplifications of the sulfur processes made in the models. In this study, therefore, we investigated impacts of different parts of the sulfur chemistry module in a global aerosol model, SPRINTARS, on the sulfate distribution and its radiative forcing. Important studies were effects of simplified and more physical-based sulfur processes in terms of treatment of sulfur chemistry, oxidant chemistry, and dry deposition process of sulfur components. The results showed that the difference in the aqueous-phase sulfur chemistry among these treatments has the largest impact on the sulfate distribution. Introduction of all the improvements mentioned above brought the model values noticeably closer to in-situ measurements than those in the simplified methods used in the original SPRINTARS model. At the same time, these improvements also brought the computed sulfate column burdens and its vertical distributions into good agreement with other AEROCOM model values. The global annual mean radiative forcing due to the direct effect of anthropogenic sulfate aerosol was thus estimated to be −0.26 W m−2 (−0.30 W m−2 with a different SO2 inventory), whereas the original SPRINTARS model showed −0.18 W m−2 (−0.21 W m−2 with a different SO2 inventory). The magnitude of the difference between original and improved methods was approximately 50% of the uncertainty among estimates by the world's global aerosol models reported by the IPCC-AR4 assessment report. Findings in the present study, therefore, may suggest that the model differences in the simplifications of the sulfur processes are still a part of the large uncertainty in their simulated radiative forcings.

Aerosols have great impacts on atmospheric environment, human health, and earth's climate. Therefore, information on their spatial and temporal distribution is of paramount importance. Despite numerous studies have examined the variation and trends of BC and AOD over India, only very few have focused on their spatial distribution or even correlating the observations with model simulations. In the present study, a three-dimensional aerosol transport-radiation model coupled with a general circulation model. SPRINTARS, simulated atmospheric aerosol distributions including BC and aerosol optical properties, i.e., aerosol optical thickness (AOT), Ångström Exponent (AE), and single scattering albedo (SSA). The simulated results are compared with both BC measurements by aethalometer and aerosol optical properties measured by ground-based skyradiometer and by satellite sensor, MODIS/Terra over Hyderabad, which is a tropical urban area of India, for the year 2008. The simulated AOT and AE in Hyderabad are found to be comparable to ground-based measured ones. The simulated SSA tends to be higher than the ground-based measurements. Both these comparisons of aerosol optical properties between the simulations with different emission inventories and the measurements indicate that, firstly the model uncertainties derived from aerosol emission inventory cannot explain the gaps between the simulations and the measurements and secondly the vertical transport of BC and the treatment of BC-containing particles can be the main issue in the global model to solve the gap.

Using one month of the cloud-resolving Nonhydrostatic Icosahedral Atmospheric Model (NICAM) simulations, we examined the impact of different definitions of clear-sky flux on the determination of longwave cloud radiative forcing (CRF). Because the satellite-like cloud-free composite preferentially samples drier conditions relative to the all-sky mean state, the conventional clear-sky flux calculation using the all-sky mean state in the model may represent a more humid atmospheric state in comparison to the cloud-free state. The drier bias is evident for the cloud-free composite in the NICAM simulations, causing an overestimation of the longwave CRF by about 10% compared to the NICAM simulated longwave CRF. Overall, water vapor contributions of up to 10% of the total longwave CRF should be taken account for making model-generated cloud forcing comparable to the satellite measurements.

We present sensitivity tests for a global aerosol assimilation system utilizing AERONET observations of AOT (aerosol optical thickness) and AAE (aerosol Ångström exponent). The assimilation system employs an ensemble Kalman filter which requires tuning of three numerical parameters: ensemble size nens, local patch size npatch and inflation factor ρ. In addition, experiments are performed to test the impact of various implementations of the system. For instance, we use a different prescription of the emission ensemble or a different combination of observations. The various experiments are compared against one-another and against independent AERONET and MODIS/Aqua observations. The assimilation leads to significant improvements in modelled AOT and AAE fields. Moreover remaining errors are mostly random while they are mostly systematic for an experiment without assimilation. In addition, these results do not depend much on our parameter or design choices. It appears that the value of the local patch size has by far the biggest impact on the assimilation, which has sufficiently converged for an ensemble size of nens=20. Assimilating AOT and AAE is clearly preferential to assimilating AOT at two different wavelengths. In contrast, initial conditions or a description of aerosol beyond two modes (coarse and fine) have only little effect. We also discuss the use of the ensemble spread as an error estimate of the analysed AOT and AAE fields. We show that a very common prescription of the emission ensemble (independent random modification in each grid cell) can have trouble generating sufficient spread in the forecast ensemble.

We present a global aerosol assimilation system based on an Ensemble Kalman filter, which we believe leads to a significant improvement in aerosol fields. The ensemble allows realistic, spatially and temporally variable model covariances (unlike other assimilation schemes). As the analyzed variables are mixing ratios (prognostic variables of the aerosol transport model), there is no need for the extra assumptions required by previous assimilation schemes analyzing aerosol optical thickness (AOT). We describe the implementation of this assimilation system and in particular the construction of the ensemble. This ensemble should represent our estimate of current model uncertainties. Consequently, we construct the ensemble around randomly modified emission scenarios. The system is tested with AERONET observations of AOT and Angström exponent (AE). Particular care is taken in prescribing the observational errors. The assimilated fields (AOT and AE) are validated through independent AERONET, SKYNET and MODIS Aqua observations. We show that, in general, assimilation of AOT observations leads to improved modelling of global AOT, while assimilation of AE only improves modelling when the AOT is high.

As a result of increasing attention paid to aerosols in climate studies, numerous global satellite aerosol products have been generated. Aerosol parameters and underlining physical processes are now incorporated in many general circulation models (GCMs) in order to account for their direct and indirect effects on the earth's climate, through their interactions with the energy and water cycles. There exists, however, an outstanding problem that these satellite products have substantial discrepancies, that must be lowered substantially for narrowing the range of the estimates of aerosol's climate effects. In this paper, numerous key uncertain factors in the retrieval of aerosol optical depth (AOD) are articulated for some widely used and relatively long satellite aerosol products including the AVHRR, TOMS, MODIS, MISR, and SeaWiFS. We systematically review the algorithms developed for these sensors in terms of four key elements that influence the quality of passive satellite aerosol retrieval: calibration, cloud screening, classification of aerosol types, and surface effects. To gain further insights into these uncertain factors, the NOAA AVHRR data are employed to conduct various tests, which help estimate the ranges of uncertainties incurred by each of the factors. At the end, recommendations are made to cope with these issues and to produce a consistent and unified aerosol database of high quality for both environment monitoring and climate studies.

A parametric study of the instantaneous radiative impact of contrails is presented using three different radiative transfer models for a series of model atmospheres and cloud parameters. Contrails are treated as geometrically and optically thin plane parallel homogeneous cirrus layers in a static atmosphere. The ice water content is varied as a function of ambient temperature. The model atmospheres include tropical, mid-latitude, and subarctic summer and winter atmospheres. Optically thin contrails cause a positive net forcing at top of the atmosphere. At the surface the radiative forcing is negative during daytime. The forcing increases with the optical depth and the amount of contrail cover. At the top of the atmosphere, a mean contrail cover of 0.1% with average optical depth of 0.2 to 0.5 causes about 0.01 to 0.03 Wm-2 daily mean instantaneous radiative forcing. Contrails cool the surface during the day and heat the surface during the night, and hence reduce the daily temperature amplitude. The net effect depends strongly on the daily variation of contrail cloud cover. The indirect radiative forcing due to particle changes in natural cirrus clouds may be of the same magnitude as the direct one due to additional cover.Key words. Atmospheric composition and structure (aerosols and particles) · Meteorology and atmospheric dynamics (climatology · radiative processes)

Radiative properties of cirrus clouds are one of the major unsolved problems in climate studies and global radiation budget. These clouds are generally composed of various ice-crystal shapes, so we tried to evaluate effects of the ice-crystal shape on radiative fluxes. We calculated radiative fluxes of cirrus clouds with a constant geometrical depth, composed of ice crystals with different shapes (hexagonal columns, bullets, bullet-rosettes), sizes and various concentrations. We considered ice particles randomly oriented in space (3D case) and their scattering phase functions were calculated by a ray-tracing method. We calculated radiative fluxes for cirrus layers for different microphysical characteristics by using a discrete-ordinate radiative code. Results showed that the foremost effect of the ice-crystal shape on radiative properties of cirrus clouds was that on the optical thickness, while the variation of the scattering phase function with the ice shape remained less than 3% for our computations. The ice-water content may be a better choice to parameterize the optical properties of cirrus, but the shape effect must be included.

Tumour necrosis factor-α (TNF-α) was measured by enzyme-linked immunosorbent assay and eosinophil cationic protein (ECP) by radio-immunoassay to evaluate TNF-α in nasal allergy. There was no significant difference either between the mean concentrations of TNF-α in nasal secretions from the patients with perennial nasal allergy and those of normal subjects, or between the TNF-α and ECP concentrations. However, reverse transcription polymerase chain reaction showed a specific increase of TNF-α mRNA and IFN-γ mRNA in allergic nasal mucosa after allergen challenge in vitro. These findings suggest a possibility that T cell-derived IFN-γ up-regulates macrophages to elaborate TNF-α, which may play a role in amplifying allergic inflammation in the nose through the cytokine network.

The effect of lyso-PAF on ciliated cells was investigated in vitro. Normal mucosa was surgically obtained from human paranasal sinuses and incubated in the form of tissue culture. Ciliated epithelium was magnified under an inverted microscope, and ciliary movement was photo-electrically measured. Ciliary activity was significantly inhibited by 10−8 M lyso-PAF and could be restored. The effect of lyso-PAF was completely blocked by CV-6209, a specific PAF antagonist. The PAF concentration in the incubation medium of lyso-PAF was determined by radioimmunoassay, because PAF is a well known inhibitor of ciliary activity. PAF gradually increased and after 20 min reached its maximal level. These findings indicated the existence of an enzyme in the paranasal sinus mucosa, by which lyso-PAF is converted to PAF, and that lyso-PAF can inhibit ciliary activity by means of PAF.

Radiative properties of cirrus clouds are one of the major unsolved problems in climate studies and global radiation budget. These clouds are generally composed of various ice-crystal shapes, so we tried to evaluate effects of the ice-crystal shape on radiative fluxes. We calculated radiative fluxes of cirrus clouds with a constant geometrical depth, composed of ice crystals with different shapes (hexagonal columns, bullets, bullet-rosettes), sizes and various concentrations. We considered ice particles randomly oriented in space (3D case) and their scattering phase functions were calculated by a ray-tracing method. We calculated radiative fluxes for cirrus layers for different microphysical characteristics by using a discrete-ordinate radiative code. Results showed that the foremost effect of the ice-crystal shape on radiative properties of cirrus clouds was that on the optical thickness, while the variation of the scattering phase function with the ice shape remained less than 3% for our computations. The ice-water content may be a better choice to parameterize the optical properties of cirrus, but the shape effect must be included.