Large-scale travelling ionospheric disturbances (LSTIDs) are detected using the critical frequency of the $F_{2}$ layer $( f_{o}F_{2})$ of the ionosphere. The HF interferometry (HF-Int) technique is applied to the network of ionosondes in the Australian region to detect and estimate propagation parameters. Here we present traveling ionospheric disturbances (TIDs) characteristics (period, velocity, and intensity) associated with selected geomagnetic storms. The TID response during geomagnetic storms varies significantly owing to the differing storm evolution patterns. Results show the morphology of parameters detected continually changes with time corresponding to highly dynamic response effects on the ionosphere. TIDs detected had a wide range of velocity magnitudes (mean ± standard deviation) such as $687 \pm 16\mathrm{m} /\mathrm{s}$ during onset phase, $516 \pm 125\mathrm{m} /\mathrm{s}$ during main phase and $569 \pm 91\mathrm{m} /\mathrm{s}$ during the recovery phase of the storm presented here. Propagation directions are predominantly equatorward with east and west deviations owing to higher latitude source generations associated with storm induced intensifications. The HF-Int system uses spectral techniques to estimate the dominant period and outputs the spectral energy contribution (SEC) which is a measure of the contribution of the TID to the total variability of the time series. TID activity levels are based on SEC criteria mostly indicated moderate and weak events with a range of periods 60-140 min.
High frequency (HF) radio wave propagation is sensitive to space weather-induced ionospheric disturbances that result from enhanced photoionization and energetic particle precipitation. Recognizing the potential risk to HF radio communication systems used by the aviation industry, as well as potential impacts on GNSS navigation and the risk of elevated radiation levels, the International Civil Aviation Organization (ICAO) initiated the development of a space weather advisory service. For HF systems, this service specifically identifies shortwave fadeout, auroral absorption, polar cap absorption, and post-storm maximum useable frequency depression (PSD) as phenomena impacting HF radio communication and specifies moderate and severe event thresholds to describe event severity. This paper examines the occurrence rate and duration of events crossing the moderate and severe thresholds. Shortwave fadeout was evaluated based on thresholds in the solar X-ray flux. Analysis of 40-years of solar X-ray flux data showed that moderate and severe level solar X-ray flares were observed, on average, 123 and 5 times per 11-year solar cycle, respectively. The mean event duration was 68 min for moderate level events and 132 min for severe level events. Auroral absorption events crossed the moderate threshold for 40 events per solar cycle, with a mean event duration of 5.1 h. The severe threshold was crossed for 3 events per solar cycle with a mean event duration of 12 h. Polar cap absorption had the longest mean duration at ~8 h for moderate events and 1.6 days for severe events; on average, 24 moderate and 13 severe events were observed per solar cycle. Moderate and severe thresholds for shortwave fadeout, auroral absorption, and polar cap absorption were used to determine the expected impacts on HF radio communication. Results for polar cap absorption and shortwave fadeout were consistent with each other, but the expected impact for auroral absorption was shown to be 2–3 times higher. Analysis of 22 years of ionosonde data showed moderate, and severe PSD events occurred, on average, 200 and 56 times per 11-year solar cycle, respectively. The mean event duration was 5.5 h for moderate-level events and 8.5 h for severe-level events. During solar cycles 22 and 23, HF radio communication was expected to experience moderate or severe impacts due to the ionospheric disturbances caused by space weather, a maximum of 163 and 78 days per year, respectively, due to the combined effect of absorption and PSD. The distribution of events is highly non-uniform with respect to the solar cycle: 70% of moderate or severe events were observed during solar maximum compared to solar minimum.
The changes in critical frequency of the F2 layer (foF2) and foF2 deviation (ΔfoF2) have been determined for three geomagnetic storms in March of the years 2012, 2013, and 2015 at low‐latitude stations, Darwin (geomag. lat. 21.96°S) and Townsville (28.95°S), and midlatitude stations, Brisbane (36.73°S), Canberra (45.65°S), and Hobart (54.17°S). The moderate storm during 15–16 March 2012 (Dst = −87 nT) showed a decrease in foF2 at midlatitude and no effect at low‐latitude stations. For the intense storm of 17–18 March 2013 (Dst = −132 nT) and the super storm of 17–18 March 2015 (Dst = −222 nT), some middle‐ to low‐latitude stations showed a short‐duration increase in foF2, but all stations showed a long‐duration decrease in foF2 during the recovery phases with ΔfoF2% varying from 26% (Darwin) to 36.6% at Hobart for the March 2013 storm and above 40% for the March 2015 storm at all of the stations. Short‐duration (~2–4 hr) increase in foF2 seems to be associated with the prompt penetrating electric fields. Long‐duration (>6 hr) decrease in foF2 is mainly accounted to the decrease in thermospheric O/N2 density ratio because of storm‐induced high‐latitude circulation of gas with depleted O/N2 density ratio to lower latitudes and partly due to disturbance dynamo electric fields. A comparison of ionosonde given foF2 for equinoctial storms (March 2013 and 2015) with similar strength Southern Hemisphere winter storms (July 2012 and June 2015) has been made with the IRI‐2016 model foF2 for Darwin, Brisbane, and Canberra stations.
Changes in ionospheric plasma densities can affect society more than ever because of our increasing reliance on communication, surveillance, navigation, and timing technology. Models struggle to predict changes in ionospheric densities at nearly all temporal and spatial scales, especially during geomagnetic storms. Here we combine a 50 year (1965–2015) geomagnetic disturbance storm time (Dst) index with plasma density measurements from a worldwide network of ~132 vertical incidence ionosondes to develop a picture of global scale changes in peak plasma density due to geomagnetic storms. Vertical incidence ionosondes provide measurements of the critical frequency of the ionospheric F2 layer (foF2), a direct measure of the peak electron density (NmF2) of the ionosphere. By dissecting the NmF2 perturbations with respect to the local time at storm onset, season, and storm intensity, it is found that (i) the storm‐associated depletions (negative storm effects) and enhancements (positive storm effects) are driven by different but related physical mechanisms, and (ii) the depletion mechanism tends to dominate over the enhancement mechanism. The negative storm effects, which are detrimental to HF radio links, are found to start immediately after geomagnetic storm onset in the nightside high‐latitude ionosphere. The depletions in the dayside high‐latitude ionosphere are delayed by a few hours. The equatorward expansion of negative storm effects is found to be regulated by storm intensity (farthest equatorward and deepest during intense storms), season (largest in summer), and time of day (generally deeper on the nightside). In contrast, positive storm effects typically occur on the dayside midlatitude and low‐latitude ionospheric regions when the storms are in the main phase, regardless of the season. Closer to the magnetic equator, moderate density enhancements last up to 40 h during the recovery phase of equinox storms, regardless of the local time. Strikingly, high‐latitude plasma densities are moderately enhanced for up to 60 h prior to the actual onset of storms during the equinoxes and summer; a potential precursor of a geomagnetic storm.
Cumulus parameterizations in general circulation models (GCMs) frequently apply mass-flux schemes in their description of tropical convection. Mass flux constitutes the product of the fractional area covered by cumulus clouds in a model grid box and the vertical velocity within the cumulus clouds. The cumulus area fraction profiles can be derived from precipitating radar reflectivity volumes. However, the vertical velocities are difficult to observe, making the evaluation of mass-flux schemes difficult. In this paper, the authors develop and evaluate a parameterization of vertical velocity in convective (cumulus) clouds using only radar reflectivities collected by a C-band polarimetric research radar (CPOL), operating at Darwin, Australia. The parameterization is trained using vertical velocity retrievals from a dual-frequency wind profiler pair located within the field of view of CPOL. The parametric model uses two inputs derived from CPOL reflectivities: the 0-dB Z echo-top height (0-dB Z ETH) and a height-weighted column reflectivity index ( Z HWT ). The 0-dB Z ETH determines the shape of the vertical velocity profile, while Z HWT determines its strength. The evaluation of these parameterized vertical velocities using (i) the training dataset, (ii) an independent wind-profiler-based dataset, and (iii) 1 month of dual-Doppler vertical velocity retrievals indicates that the statistical representation of vertical velocity is reasonably accurate up to the 75th percentile. However, the parametric model underestimates the extreme velocities. The method allows for the derivation of cumulus mass flux and its variability on current GCM scales based only on reflectivities from precipitating radar, which could be valuable to modelers.
This study evaluates the ability of the regional high‐resolution version of the Australian Community Climate and Earth‐System Simulator (ACCESS‐A) model (horizontal resolution of 12 km, 50 height levels with a 37.5 km top) to simulate statistical rainfall properties during three monsoon seasons using the Darwin C‐band polarimetric (CPOL) radar observation as a reference. The comparison consists of characteristics of precipitation, including the total‐domain accumulated daily rainfall, the daily rain intensity and frequency and the probability distribution function of hourly rain rate. These characteristics are then binned using five large‐scale atmospheric regimes that represent the main weather states in that region to assess the ability of the model to reproduce the variability of the statistical rainfall properties as a function of these defined regimes. The main results are that although the model reproduces the accumulated daily rainfall amount reasonably well, albeit tending to slightly overestimate it, the mean rain rate is too weak, compensated by too large a rain frequency of occurrence. In terms of intensity the model tends to perform worse in the dry regimes but improves in the wet regimes, especially during the ‘Deep Westerly’ regime, when most of the monsoonal precipitation is produced. Probability distribution functions of rain rates reveal that although there are small differences amongst regimes, marked discrepancies are found between model and observations. The ACCESS‐A overestimates light rain rates, underestimates moderate to heavy rain rates and produces spurious very heavy rain rates that are not observed at all by the radar. The evaluation of the diurnal cycle of rainfall properties also reveals a general underestimate of rain rates, compensated by an overestimated rain frequency over the whole cycle, resulting in an overestimated amplitude and a wrong timing of the diurnal cycle. In all regimes, the overestimation of total‐domain daily rainfall and rain frequency is found to be due to an overestimation of low rain rates produced by the ACCESS‐A convective scheme, while the underestimation of mean rain rate is due to an underestimate of high rain rates and high rain rates frequency by the ACCESS‐A convective scheme.
Cumulus parameterizations in weather and climate models frequently apply mass-flux schemes in their description of tropical convection. Mass flux constitutes the product of the fractional area covered by convection in a model grid box and the vertical velocity in cumulus clouds. However, vertical velocities are difficult to observe on GCM scales, making the evaluation of mass-flux schemes difficult. Here, the authors combine high-temporal-resolution observations of in-cloud vertical velocities derived from a pair of wind profilers over two wet seasons at Darwin with physical properties of precipitating clouds [cloud-top heights (CTH), convective-stratiform classification] derived from the Darwin C-band polarimetric radar to provide estimates of cumulus mass flux and its constituents. The length of this dataset allows for investigations of the contributions from different cumulus cloud types-namely, congestus, deep, and overshooting convection-to the overall mass flux and of the influence of large-scale conditions on mass flux. The authors found that mass flux was dominated by updrafts and, in particular, the updraft area fraction, with updraft vertical velocity playing a secondary role. The updraft vertical velocities peaked above 10 km where both the updraft area fractions and air densities were small, resulting in a marginal effect on mass-flux values. Downdraft area fractions are much smaller and velocities aremuch weaker than those in updrafts. The area fraction responded strongly to changes in midlevel large-scale vertical motion and convective inhibition (CIN). In contrast, changes in the lower-tropospheric relative humidity and convective available potential energy (CAPE) strongly modulate in-cloud vertical velocities but have moderate impacts on area fractions. Although average mass flux is found to increase with increasing CTH, it is the environmental conditions that seem to dictate the magnitude of mass flux produced by convection through a combination of effects on area fraction and velocity.
In this paper, statistical properties of rainfall are derived from 14 years of Tropical Rainfall Measuring Mission data to optimize the use of flight hours for the upcoming High Altitude Ice Crystals (HAIC)/High Ice Water Content (HIWC) program. This program aims to investigate the convective processes responsible for the generation of the high ice water content that has been recognized as a threat to civil aviation. The probability that convective cells are conducive to HIWC is also further investigated using three years of C-band polarimetric radar data. Further insights into the variability of convective rainfall and favorable conditions for HIWC are also gained using two different methods to characterize the large-scale atmospheric conditions around Darwin, Australia (the Madden-Julian oscillation and the Darwin atmospheric regimes), and the underlying surface type (oceanic vs continental). The main results from the climatology relevant to flight-plan decision making are (i) convective cells conducive to HIWC should be found close to Darwin, (ii) at least 90% of convective cells are conducive to HIWC at 10- and 12-km flight levels, (iii) multiple flights per day in favorable large-scale conditions will be needed so as to utilize the 150 project flight hours, (iv) the largest numbers of HIWC radar pixels are found around 0300 and 1500 local time, and (v) to fulfill the requirement to fly 90 h in oceanic convection and 60 h in or around continental convection, a minimum acceptable size of the convective area has been derived and should serve as a guideline for flight-decision purposes.
Systematic analysis of the surface fair-weather potential gradient (PG) measured for 13 months (July 2005– July 2006) at 10 s resolution over a small tropical island station Suva (18.08◦S, 178.45◦E), Fiji is presented. Based on the solar radiation (>800 W m−2) and surface wind speed (<4.5 m s−1) conditions, 63 fair-weather days were selected. After sorting the data into a range of 0–1000 V m−1, the average PG was ∼139 V m−1. The measured fair-weather PG had a semi-diurnal structure, with a more pronounced peak at 0730 LT (1930 UT) and a reduced peak at 2200 LT (1000 UT). The time of occurrence of the morning peak and the noontime minima agreed well with the Carnegie curve. However, the variation about the global mean in the measured PG was 150% whereas for the Carnegie curve it was only 35%. The effects of the local meteorological parameters on the PG measurements were found to be small. On average, the PG during the dry season (May–October) was always greater than in the wet season (November–April). In contrast, analysis of regional (0–60◦S and 100◦E–160◦W) lightning activity on fair-weather days showed a peak at ∼2000 LT and higher lightning activity during the wet season. These results indicate that the regional thunderstorm activity has no direct connection with the local fair-weather PG at the site.
In this study, methods of convective/stratiform precipitation classification and surface rain-rate estimation based on the Atmospheric Radiation Measurement Program (ARM) cloud radar measurements were developed and evaluated. Simultaneous and collocated observations of the Ka-band ARM zenith radar (KAZR), two scanning precipitation radars [NCAR S-band/Ka-band Dual Polarization, Dual Wavelength Doppler Radar (S-PolKa) and Texas A&M University Shared Mobile Atmospheric Research and Teaching Radar (SMART-R)], and surface precipitation during the Dynamics of the Madden-Julian Oscillation/ARM MJO Investigation Experiment (DYNAMO/AMIE) field campaign were used. The motivation of this study is to apply the unique long-term ARM cloud radar observations without accompanying precipitation radars to the study of cloud life cycle and precipitation features under different weather and climate regimes. The resulting convective/stratiform classification from KAZR was evaluated against precipitation radars. Precipitation occurrence and classified convective/stratiform rain fractions from KAZR compared favorably to the collocated SMART-R and S-PolKa observations. Both KAZR and S-PolKa radars observed about 5% precipitation occurrence. The convective (stratiform) precipitation fraction is about 18% (82%). Collocated disdrometer observations of two days showed an increased number concentration of small and large raindrops in convective rain relative to dominant small raindrops in stratiform rain. The composite distributions of KAZR reflectivity and Doppler velocity also showed distinct structures for convective and stratiform rain. These evidences indicate that the method produces physically consistent results for the two types of rain. A new KAZR-based, two-parameter [the gradient of accumulative radar reflectivity Z(e) (GAZ) below 1 km and near-surface Z(e)] rain-rate estimation procedure was developed for both convective and stratiform rain. This estimate was compared with the exponential Z-R (reflectivity-rain rate) relation. The relative difference between the estimated and surface-measured rainfall rates showed that the two-parameter relation can improve rainfall estimation relative to the Z-R relation.
The technique of radio occultation (RO) is demonstrated to be a powerful tool for studying equatorial F-region irregularities (EFIs) associated with equatorial plasma bubbles. The extensive 4.9 year RO dataset of the Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC) satellites was employed in this study and contains EFI observations under a wide variety of solar and geomagnetic conditions. From an analysis of the EFI occurrence dependence on season/longitude, it is found that the EFI occurrence statistics largely match those reported previously, with the exception of an equinoctial EFI occurrence maximum in the American sector that is absent from previous studies. It is revealed that this maximum is due to enhanced EFI occurrence near the South Atlantic anomaly, where EFIs are expected to be suppressed by particle precipitation. An investigation into the solar activity dependence of the EFI occurrence characteristics revealed significant increases in the range of local times and latitudes with solar activity for most longitude sectors and seasons. Finally, the EFI suppression and enhancement effects of storm-time electric fields are also investigated using the COSMIC data. Citation: Carter, B. A., K. Zhang, R. Norman, V. V. Kumar, and S. Kumar (2013), On the occurrence of equatorial F-region irregularities during solar minimum using radio occultation measurements, J. Geophys. Res. Space Physics, 118, 892-904, doi: 10.1002/jgra.50089.
There is no objective definition to separate cumulus congestus clouds from the shallow cumulus and deep clouds. This has generated misinterpretation about the role of congestus clouds to promote deep convection through the potential of moistening the middle troposphere. In this study, an objective identification for the different tropical cumulus modes is found by examining the occurrence frequency of the cloud cell top heights (CTHs) and near‐ground (at 2.5 km height) rainfall properties of these cells using a three‐season database of the Darwin C‐band polarimetric radar. Four cumulus modes were identified, namely a shallow cumulus mode with CTH in the trade inversion layer (1–3 km), a congestus mode with tops in the highly stable middle troposphere (3–6.5 km), a deep convective mode with tops in the region of free convection (6.5–15 km), and an overshooting convection mode with tops in the tropical tropopause layer (CTH >15 km). The study also investigates the connections between these cumulus modes during heavy rainfall events. The congestus mode occurs predominantly from ~10 h prior to the peak rainfall event to ~2 h past the event. The deep cloud populations (Modes 3 and 4) have their maxima at and shortly after the time of the rainfall peak, with maximum occurrence just below the tropical tropopause layer. A comparison of the heavy rainfall events occurring in morning (oceanic) conditions against the afternoon (continental) conditions revealed a higher ratio of the shallow to the deep cloud population and a shorter transition time from the shallow to the onset of deep population in the morning‐oceanic conditions than the afternoon‐land conditions. It is also found through the analysis of the large‐scale moisture budget data set that for both the morning and afternoon events, the moistening peaked before the peak in the congestus populations.
C-band polarimetric radar measurements spanning two wet seasons are used to perform a critical evaluation of two algorithms for the classification of stratiform and convective precipitation. The first approach is based on the horizontal texture of the radar reflectivity field (two classes: stratiform, convective), and the second approach is based on the properties of the drop size distribution (DSD) parameters as derived from a set of polarimetric variables (three classes: stratiform, mixed, convective). To investigate how well those two methods compare quantitatively, probability density functions of reflectivity, rain rate, 5-dBZ echo top height, and DSD parameters (namely, the median volume diameter and the generalized intercept parameter) are built. The study found that while the two methods agree well on the identification of stratiform precipitation, large differences are obtained for convective rainfall. The texture-based approach seems to classify too many points as being of convective nature compared to the DSD-based method. Among the points that are classified as convective by the texture-based approach, 25% correspond to low concentration of relatively small particles associated with rain rates below 10 mm h(-1). This large proportion of unrealistically low convective rain rates is not produced by the DSD-based approach, which only classifies 4% of the convective points with rain rates below 10 mm h(-1). These points were found to be mainly isolated points embedded within stratiform precipitation and associated with low cloud-top height, suggesting a misclassification of the texture-based approach. Thus, to improve the statistics of the convective class, three modified equations of the peakedness criterion used in the radar-based algorithm are proposed to decrease the number of misclassified points.
C-band polarimetric radar measurements spanning two wet seasons are used to study the effects of the large-scale environment on the statistical properties of stratiform and convective rainfall around Darwin, Australia. The rainfall physical properties presented herein are the reflectivity fields, daily rainfall accumulations and raining area, rain rates, and drop size distribution (DSD) parameters (median volume diameter and normalized intercept parameter). Each of these properties is then analyzed according to five different atmospheric regimes and further separated into stratiform or convective rain categories following a DSD-based approach. The regimes, objectively identified by radiosonde thermodynamic and wind measurements, represent typical wet-season atmospheric conditions: the active monsoon regime, the break periods, the buildup regime, the trade wind regime, and a mixture of inactive/break periods. The large-scale context is found to strongly modulate rainfall and cloud microphysical properties. For example, during the active monsoon regime, the daily rain accumulation is higher than in the other regimes, while this regime is associated with the lowest rain rates. Precipitation in this active monsoon regime is found to be widespread and mainly composed of small particles in high concentration compared to the other regimes. Vertical profiles of reflectivity and DSD parameters suggest that warm rain processes are dominant during this regime. In contrast, rainfall properties in the drier regimes (trade wind/buildup regimes) are mostly of continental origin, with rain rates higher than in the moister regimes. In these drier regimes, precipitation is mainly formed of large raindrops in relatively low concentration due to a larger contribution of the ice microphysical processes on the rainfall formation.
Some cumulus clouds with tops between 3 and 7 km (Cu3km-7km) remain in this height region throughout their lifetime (congestus) while others develop into deeper clouds (cumulonimbus). This study describes two techniques to identify the congestus and cumulonimbus cloud types using data from scanning weather radar and identifies the atmospheric conditions that regulate these two modes. A two-wet-season cumulus cloud database of the Darwin C-band polarimetric radar is analyzed and the two modes are identified by examining the 0-dBZ cloud-top height (CTH) of the Cu3km-7km cells over a sequence of radar scans. It is found that similar to 26% of the classified Cu3km-7km population grow into cumulonimbus clouds. The cumulonimbus cells exhibit reflectivities, rain rates, and drop sizes larger than the congestus cells. The occurrence frequency of cumulonimbus cells peak in the afternoon at similar to 1500 local timea few hours after the peak in congestus cells. The analysis of Darwin International Airport radiosonde profiles associated with the two types of cells shows no noticeable difference in the thermal stability rates, but a significant difference in midtropospheric (5-10 km) relative humidity. Moister conditions are found in the hours preceding the cumulonimbus cells when compared with the congestus cells. Using a moisture budget dataset derived for the Darwin region, it is shown that the existence of cumulonimbus cells, and hence deep convection, is mainly determined by the presence of the midtroposphere large-scale upward motion and not merely by the presence of congestus clouds prior to deep convection. This contradicts the thermodynamic viewpoint that the midtroposphere moistening prior to deep convection is solely due to the preceding cumulus congestus cells.
A persistent problem for numerical weather and climate models is the representation of tropical convective precipitation which for the most part occurs on spatial and temporal scales too small and too short to be explicitly resolved. Given that model parameterizations represent this subgrid convection as a function of the large‐scale atmospheric state, an understanding of the strongest relationships between the two scales is needed. This study introduces a method to create two concurrent long‐term data sets that describe both the large‐scale atmosphere and the characteristics of the small‐scale convection. Important relationships between these two scales are then investigated. It is found that convective precipitation, through convective precipitation area, has the strongest relationship with dynamical variables such as moisture convergence and vertical velocity at midlevels. The magnitude of the fluctuations of convective strength about the mean is found to be anticorrelated with the strength of the large‐scale variables, indicating a more stochastic behavior of tropical convection in weakly than strongly forced regimes, respectively. Atmospheric stability related variables are not found to be positively related to either convective precipitation area or convective precipitation intensity, which is often an assumption made in convective parameterization. On the contrary, in a more unstable atmosphere, there is lower convective precipitation.