Abstract

Bullet rosette particles are common in ice clouds, and the bullets may often be hollow. Here the single-scattering properties of randomly oriented hollow bullet rosette ice particles are investigated. A bullet, which is an individual branch of a rosette, is defined as a hexagonal column attached to a hexagonal pyramidal tip. For this study, a hollow structure is included at the end of the columnar part of each bullet branch and the shape of the hollow structure is defined as a hexagonal pyramid. A hollow bullet rosette may have between 2 and 12 branches. An improved geometric optics method is used to solve for the scattering of light in the particle. The primary optical effect of incorporating a hollow end in each of the bullets is to decrease the magnitude of backscattering. In terms of the angular distribution of scattered energy, the hollow bullets increase the scattering phase function values within the forward scattering angle region from 1° to 20° but decrease the phase function values at side- and backscattering angles of 60°–180°. As a result, the presence of hollow bullets tends to increase the asymmetry factor. In addition to the scattering phase function, the other elements of the phase matrix are also discussed. The backscattering depolarization ratios for hollow and solid bullet rosettes are found to be very different. This may have an implication for active remote sensing of ice clouds, such as from polarimetric lidar measurements. In a comparison of solid and hollow bullet rosettes, the effect of the differences on the retrieval of both the ice cloud effective particle size and optical thickness is also discussed. It is found that the presence of hollow bullet rosettes acts to decrease the inferred effective particle size and to increase the optical thickness in comparison with the use of solid bullet rosettes.

Abstract

The present study investigates interdecadal modulations of the El Niño–Southern Oscillation (ENSO) influence on the climate of the northwest Pacific (NWP) and East Asia (EA) in early boreal summer following a winter ENSO event, based on 19 simulations from phase 5 of the Coupled Model Intercomparison Project (CMIP5). In the historical run, 8 out of 19 models capture a realistic relationship between ENSO and NWP early summer climate—an anomalous anticyclone develops over the NWP following a winter El Niño event—and the interdecadal modulations of this correlation. During periods when the association between ENSO and NWP early summer climate is strong, ENSO variance and ENSO-induced anomalies of summer sea surface temperature (SST) and tropospheric temperature over the tropical Indian Ocean (TIO) all strengthen relative to periods when the association is weak.

In future projections with representative concentration pathways 4.5 and 8.5, the response of TIO SST, tropospheric temperature, and NWP anomalous anticyclone to ENSO all strengthen regardless of ENSO amplitude change. In a warmer climate, low-level specific humidity response to interannual SST variability strengthens following the Clausius–Clapeyron equation. The resultant intensification of tropospheric temperature response to interannual TIO warming is suggested as the mechanism for the strengthened ENSO effect on NWP–EA summer climate.

Abstract

The Moderate Resolution Imaging Spectroradiometer (MODIS) and the Atmospheric Infrared Sounder (AIRS) measurements from the NASA Earth Observing System Aqua satellite enable global monitoring of the distribution of clouds during day and night. The MODIS is able to provide a high-spatial-resolution (1–5 km) cloud mask, cloud classification mask, cloud-phase mask, cloud-top pressure (CTP), and effective cloud amount during both the daytime and the nighttime, as well as cloud particle size (CPS) and cloud optical thickness (COT) at 0.55 μm during the daytime. The AIRS high-spectral-resolution measurements reveal cloud properties with coarser spatial resolution (13.5 km at nadir). Combined, MODIS and AIRS provide cloud microphysical properties during both the daytime and nighttime. A fast cloudy radiative transfer model for AIRS that accounts for cloud scattering and absorption is described in this paper. One-dimensional variational (1DVAR) and minimum-residual (MR) methods are used to retrieve the CPS and COT from AIRS longwave window region (790–970 cm⁻¹ or 10.31–12.66 μm, and 1050–1130 cm⁻¹ or 8.85–9.52 μm) cloudy radiance measurements. In both 1DVAR and MR procedures, the CTP is derived from the AIRS radiances of carbon dioxide channels while the cloud-phase information is derived from the collocated MODIS 1-km phase mask for AIRS CPS and COT retrievals. In addition, the collocated 1-km MODIS cloud mask refines the AIRS cloud detection in both 1DVAR and MR procedures. The atmospheric temperature profile, moisture profile, and surface skin temperature used in the AIRS cloud retrieval processing are from the European Centre for Medium-Range Weather Forecasts forecast analysis. The results from 1DVAR are compared with the operational MODIS products and MR cloud microphysical property retrieval. A Hurricane Isabel case study shows that 1DVAR retrievals have a high correlation with either the operational MODIS cloud products or MR cloud property retrievals. 1DVAR provides an efficient way for cloud microphysical property retrieval during the daytime, and MR provides the cloud microphysical property retrievals during both the daytime and nighttime.

Abstract

El Niño’s influence on the subtropical northwest (NW) Pacific climate increased after the climate regime shift of the 1970s. This is manifested in well-organized atmospheric anomalies of suppressed convection and a surface anticyclone during the summer (June–August) of the El Niño decay year [JJA(1)], a season when equatorial Pacific sea surface temperature (SST) anomalies have dissipated. In situ observations and ocean–atmospheric reanalyses are used to investigate mechanisms for the interdecadal change. During JJA(1), the influence of the El Niño–Southern Oscillation (ENSO) on the NW Pacific is indirect, being mediated by SST conditions over the tropical Indian Ocean (TIO). The results here show that interdecadal change in this influence is due to changes in the TIO response to ENSO.

During the postregime shift epoch, the El Niño teleconnection excites downwelling Rossby waves in the south TIO by anticyclonic wind curls. These Rossby waves propagate slowly westward, causing persistent SST warming over the thermocline ridge in the southwest TIO. The ocean warming induces an antisymmetric wind pattern across the equator, and the anomalous northeasterlies cause the north Indian Ocean to warm through JJA(1) by reducing the southwesterly monsoon winds. The TIO warming excites a warm Kelvin wave in tropospheric temperature, resulting in robust atmospheric anomalies over the NW Pacific that include the surface anticyclone. During the preregime shift epoch, ENSO is significantly weaker in variance and decays earlier than during the recent epoch. Compared to the epoch after the mid-1970s, SST and wind anomalies over the TIO are similar during the developing and mature phases of ENSO but are very weak during the decay phase. Specifically, the southern TIO Rossby waves are weaker, so are the antisymmetric wind pattern and the North Indian Ocean warming during JJA(1). Without the anchor in the TIO warming, atmospheric anomalies over the NW Pacific fail to develop during JJA(1) prior to the mid-1970s. The relationship of the interdecadal change to global warming and implications for the East Asian summer monsoon are discussed.

Abstract

The 2015 Paris Agreement proposed targets to limit global-mean surface temperature (GMST) rise well below 2°C relative to preindustrial level by 2100, requiring a cease in the radiative forcing (RF) increase in the near future. In response to changing RF, the deep ocean responds slowly (ocean slow response), in contrast to the fast ocean mixed layer adjustment. The role of the ocean slow response under low warming targets is investigated using representative concentration pathway (RCP) 2.6 simulations from phase 5 of the Coupled Model Intercomparison Project. In RCP2.6, the deep ocean continues to warm while RF decreases after reaching a peak. The deep ocean warming helps to shape the trajectories of GMST and fuels persistent thermosteric sea level rise. A diagnostic method is used to decompose further changes after the RF peak into a slow warming component under constant peak RF and a cooling component due to the decreasing RF. Specifically, the slow warming component amounts to 0.2°C (0.6°C) by 2100 (2300), raising the hurdle for achieving the low warming targets. When RF declines, the deep ocean warming takes place in all basins but is the most pronounced in the Southern Ocean and Atlantic Ocean where surface heat uptake is the largest. The climatology and change of meridional overturning circulation are both important for the deep ocean warming. To keep the GMST rise at a low level, substantial decrease in RF is required to offset the warming effect from the ocean slow response.

Abstract

Cloud longwave scattering is generally neglected in general circulation models (GCMs), but it plays a significant and highly uncertain role in the atmospheric energy budget as demonstrated in recent studies. To reduce the errors caused by neglecting cloud longwave scattering, two new radiance adjustment methods are developed that retain the computational efficiency of broadband radiative transfer simulations. In particular, two existing scaling methods and the two new adjustment methods are implemented in the Rapid Radiative Transfer Model (RRTM). The results are then compared with those based on the Discrete Ordinate Radiative Transfer model (DISORT) that explicitly accounts for multiple scattering by clouds. The two scaling methods are shown to improve the accuracy of radiative transfer simulations for optically thin clouds but not effectively for optically thick clouds. However, the adjustment methods reduce computational errors over a wide range, from optically thin to thick clouds. With the adjustment methods, the errors resulting from neglecting cloud longwave scattering are reduced to less than 2 W m⁻² for the upward irradiance at the top of the atmosphere and less than 0.5 W m⁻² for the surface downward irradiance. The adjustment schemes prove to be more accurate and efficient than a four-stream approximation that explicitly accounts for multiple scattering. The neglect of cloud longwave scattering results in an underestimate of the surface downward irradiance (cooling effect), but the errors are almost eliminated by the adjustment methods (warming effect).

Abstract

The FORMOSAT-7/COSMIC-2 Global Navigation Satellite System (GNSS) Radio Occultation (RO) satellite constellation was launched in June 2019 as a successor of the FORMOSAT-3/COSMIC mission. The Central Weather Bureau (CWB) of Taiwan has received FORMOSAT-7/COSMIC-2 GNSS RO data in real time from the Taiwan Analysis Center for COSMIC. With the global numerical prediction system at CWB, a parallel semioperational experiment assimilating the FORMOSAT-7/COSMIC-2 bending angle data with all other operational observation data has been conducted to evaluate the impact of the FORMOSAT-7/COSMIC-2 data. The first seven-month results show that the quality of the early FORMOSAT-7/COSMIC-2 data has been satisfactory for assimilation. Consistent and significant positive impacts on global forecast skills have been observed since the start of the parallel experiment, with the most significant impact found in the tropical region, reflecting the low-inclination orbital design of the satellites. The impact of the FORMOSAT-7/COSMIC-2 RO data is also estimated using the ensemble forecast sensitivity to observation impact (EFSOI) method, showing an average positive impact per observation similar to other existing GNSS RO datasets, while the total impact is impressive by virtue of its large amount. Sensitivity experiments suggest that the quality control processes built in the Gridpoint Statistical Interpolation (GSI) system for RO data work well to achieve a positive impact by the low-level FORMOSAT-7/COSMIC-2 RO data, while more effort on observation error tuning should be focused to obtain an optimal assimilation performance. This study demonstrates the usefulness of the FORMOSAT-7/COSMIC-2 RO data in global numerical weather prediction during the calibration/validation period and leads to the operational use of the data at CWB.

Abstract

The unique characteristics of land–air coupling and troposphere–stratosphere interaction over the Tibetan Plateau (TP), the highest landform in the world, play a vital role in weather and climate on regional and global scales. Although a great deal of research has been carried out, large gaps remain in our understanding of TP land–air coupling and its climate effects, due to a lack of observations and the issue of model biases. To address these obstacles, a 10-yr national research program entitled “Changes in the Land–Air Coupled System over the Tibetan Plateau and its Impacts on Global Climate (LASTPIC)” was launched by the National Natural Science Foundation of China in January 2014. What LASTPIC does revolves around three aspects: TP land–air coupled processes; TP’s influence on global climate; and reanalysis and model. This paper mainly focuses on the data collection, scientific understanding, and model development of LASTPIC in terms of TP land–atmosphere–ocean coupling and its global climate impacts since program’s inception.

Abstract

This study evaluates the impact of dust–radiation–cloud interactions on the development of a mesoscale convective system (MCS) by comparing numerical experiments run with and without dust–radiation and/or dust–cloud interactions. An MCS that developed over North Africa on 4–6 July 2010 is used as a case study. The CloudSat and Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) satellites passed over the center of the MCS after it reached maturity, providing valuable profiles of aerosol backscatter and cloud information for model verification. The model best reproduces the MCS’s observed cloud structure and morphology when both dust–radiation and dust–cloud interactions are included. Our results indicate that the dust–radiation effect has a far greater influence on the MCS’s development than the dust-cloud effect. Results show that the dust-radiative effect, both with and without the dust–cloud interaction, briefly delays the MCS’s formation but ultimately produces a stronger storm with a more extensive anvil cloud. This is caused by dust–radiation-induced changes to the MCS’s environment. The impact of the dust–cloud effect on the MCS, on the other hand, is greatly affected by the presence of the dust–radiation interaction. The dust–cloud effect alone slows initial cloud development but enhances heterogeneous ice nucleation and extends cloud lifetime. When the dust–radiation interaction is added, increased transport of dust into the upper portions of the storm—due to a dust–radiation-driven increase in convective intensity—allows dust–cloud processes to more significantly enhance heterogeneous freezing activity earlier in the storm’s development, increasing updraft strength, hydrometeor growth (particularly for ice particles), and rainfall.

Abstract

The current cloud thermodynamic phase discrimination by Cloud-Aerosol Lidar Pathfinder Satellite Observations (CALIPSO) is based on the depolarization of backscattered light measured by its lidar [Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP)]. It assumes that backscattered light from ice crystals is depolarizing, whereas water clouds, being spherical, result in minimal depolarization. However, because of the relationship between the CALIOP field of view (FOV) and the large distance between the satellite and clouds and because of the frequent presence of oriented ice crystals, there is often a weak correlation between measured depolarization and phase, which thereby creates significant uncertainties in the current CALIOP phase retrieval. For water clouds, the CALIOP-measured depolarization can be large because of multiple scattering, whereas horizontally oriented ice particles depolarize only weakly and behave similarly to water clouds. Because of the nonunique depolarization–cloud phase relationship, more constraints are necessary to uniquely determine cloud phase. Based on theoretical and modeling studies, an improved cloud phase determination algorithm has been developed. Instead of depending primarily on layer-integrated depolarization ratios, this algorithm differentiates cloud phases by using the spatial correlation of layer-integrated attenuated backscatter and layer-integrated particulate depolarization ratio. This approach includes a two-step process: 1) use of a simple two-dimensional threshold method to provide a preliminary identification of ice clouds containing randomly oriented particles, ice clouds with horizontally oriented particles, and possible water clouds and 2) application of a spatial coherence analysis technique to separate water clouds from ice clouds containing horizontally oriented ice particles. Other information, such as temperature, color ratio, and vertical variation of depolarization ratio, is also considered. The algorithm works well for both the 0.3° and 3° off-nadir lidar pointing geometry. When the lidar is pointed at 0.3° off nadir, half of the opaque ice clouds and about one-third of all ice clouds have a significant lidar backscatter contribution from specular reflections from horizontally oriented particles. At 3° off nadir, the lidar backscatter signals for roughly 30% of opaque ice clouds and 20% of all observed ice clouds are contaminated by horizontally oriented crystals.