Imaging of structures in the high-latitude ionosphere: model comparisons

The tomographic reconstruction technique generates a two-dimensional latitude versus height electron density distribution from sets of slant total electron content measurements (TEC) along ray paths between beacon satellites and ground-based radio receivers. In this note, the technique is applied to TEC values obtained from data simulated by the Sheffield/UCL/SEL Coupled Thermosphere/Ionosphere/Model (CTIM). A comparison of the resulting reconstructed image with the input modelled data allows for verification of the reconstruction technique. All the features of the high-latitude ionosphere in the model data are reproduced well in the tomographic image. Reconstructed vertical TEC values follow closely the modelled values, with the F-layer maximum density (NmF2) agreeing generally within about 10%. The method has also been able successfully to reproduce underlying auroral-E ionisation over a restricted latitudinal range in part of the image. The height of the F2 peak is generally in agreement to within about the vertical image resolution (25 km).     

Trends in the ionospheric E and F regions over Europe

Continuous observations in the ionospheric E and F regions have been regularly carried out since the fifties of this century at many ionosonde stations. Using these data from 31 European stations long-term trends have been derived for different parameters of the ionospheric E layer (h E, foE), F1 layer (foF1) and F2 layer (hmF2,foF2). The detected trends in the E and F1 layers (lowering of the E region height hE; increase of the peak electron densities of the E and F1 layers, foE and foF1) are in qualitative agreement with model predictions of an increasing atmospheric greenhouse effect. In the F2 region, however, the results are more complex. Whereas in the European region west of 30° E negative trends in hmF2 (peak height of the F2 layer) and in the peak electron density (foF2) have been found, in the eastern part of Europe (east of 30° E) positive trends dominate in both parameters. These marked longitudinal differences cannot be explained by an increasing greenhouse effect only, here probably dynamical effects in the F2 layer seem to play an essential role.     

Annual and semiannual variations in the ionospheric F2-layer. I. Modelling

Annual, seasonal and semiannual variations of F2-layer electron density (NmF2) and height (hmF2) have been compared with the coupled thermosphere-ionosphere-plasmasphere computational model (CTIP), for geomagnetically quiet conditions. Compared with results from ionosonde data from midlatitudes, CTIP reproduces quite well many observed features of NmF2, such as the dominant winter maxima at high midlatitudes in longitude sectors near the magnetic poles, the equinox maxima in sectors remote from the magnetic poles and at lower latitudes generally, and the form of the month-to-month variations at latitudes between about 60°N and 50°S. CTIP also reproduces the seasonal behaviour of NmF2 at midnight and the summer-winter changes of hmF2. Some features of the F2-layer, not reproduced by the present version of CTIP, are attributed to processes not included in the modelling. Examples are the increased prevalence of the winter maxima of noon NmF2 at higher solar activity, which may be a consequence of the increase of F2-layer loss rate in summer by vibrationally excited molecular nitrogen, and the semiannual variation in hmF2, which may be due to tidal effects. An unexpected feature of the computed distributions of NmF2 is an east-west hemisphere difference, which seems to be linked to the geomagnetic field configuration. Physical discussion is reserved to the companion paper by Rishbeth et al.     

Geomagnetic control of the foF2 long-term trends

Further development of the method proposed by Danilov and Mikhailov is presented. The method is applied to reveal the foF2 long-term trends on 30 Northern Hemisphere ionosonde stations. Most of them show significant foF2 trends. A pronounced dependence of trend magnitude on geomagnetic (invariant) latitude is confirmed. Periods of negative/positive foF2 trends corresponding to the periods of long-term increasing/decreasing geomagnetic activity are revealed for the first time. Pronounced diurnal variations of the foF2 trend magnitude are found. Strong positive foF2 trends in the post-midnight-early-morning LT sector and strong negative trends during daytime hours are found on the sub-auroral stations for the period with increasing geomagnetic activity. On the contrary middle and lower latitude stations demonstrate negative trends in the early-morning LT sector and small negative or positive trends during daytime hours for the same period. All the morphological features revealed of the foF2 trends may be explained in the framework of contemporary F2-region storm mechanisms. This newly proposed F2-layer geomagnetic storm concept casts serious doubts on the hypothesis relating the F2-layer parameter long-term trends to the thermosphere cooling due to the greenhouse effect.     

Observations of Pc 3–4 and Pi 2 geomagnetic pulsations in the low-latitude ionosphere

Day-time Pc 3–4 (≃5–60 mHz) and night-time Pi 2 (≃5–20 mHz) ULF waves propagating down through the ionosphere can cause oscillations in the Doppler shift of HF radio transmissions that are correlated with the magnetic pulsations recorded on the ground. In order to examine properties of these correlated signals, we conducted a joint HF Doppler/magnetometer experiment for two six-month intervals at a location near L = 1.8. The magnetic pulsations were best correlated with ionospheric oscillations from near the F region peak. The Doppler oscillations were in phase at two different altitudes, and their amplitude increased in proportion to the radio sounding frequency. The same results were obtained for the O- and X-mode radio signals. A surprising finding was a constant phase difference between the pulsations in the ionosphere and on the ground for all frequencies below the local field line resonance frequency, independent of season or local time. These observations have been compared with theoretical predictions of the amplitude and phase of ionospheric Doppler oscillations driven by downgoing Alfvén mode waves. Our results agree with these predictions at or very near the field line resonance frequency but not at other frequencies. We conclude that the majority of the observations, which are for pulsations below the resonant frequency, are associated with downgoing fast mode waves, and models of the wave-ionosphere interaction need to be modified accordingly.     

Some new features of electron density irregularities over SHAR during strong spread F

An RH-560 rocket flight was conducted from Sriharikota rocket range (SHAR) (14°N, 80°E, dip latitude 5.5°N) to study electron density and electric field irregularities during spread F. The rocket was launched at 2130 local time (LT) and it attained an apogee of 348 km. Results of electron density fluctuations are presented here. Two extremely sharp layers of very high electron density were observed at 105 and 130 km. The electron density increase in these layers was by a factor of 50 in a vertical extent of 10 km. Large depletions in electron density were observed around 175 and 238 km. Both sharp layers as well as depletions were observed also during the descent. The presence of sharp layers and depletions during the ascent and the descent of the rocket as well as an order of magnitude less electron density, in 150/300 km region during the descent, indicate the presence of strong large-scale horizontal gradients in the electron density. Some of the valley region irregularities (165/178 km), in the intermediate scale size range, observed during this flight, show spectral peaks at 2 km and can be interpreted in terms of the image striation theory suggested by Vickrey et al. The irregularities at 176 km do not exhibit any peak at kilometer scales and appear to be of new type. The growth rate of intermediate scale size irregularities, produced through generalized Rayleigh Taylor instability, was calculated for the 200/330 km altitude, using observed values of electron density gradients and an assumed vertically downward wind of 20 ms^–1. These growth rate calculations suggest that the observed irregularities could be produced by the gradient drift instability.     

Annual and semiannual variations in the ionospheric F2-layer: II. Physical discussion

The companion paper by Zou et al. shows that the annual and semiannual variations in the peak F2-layer electron density (NmF2) at midlatitudes can be reproduced by a coupled thermosphere-ionosphere computational model (CTIP), without recourse to external influences such as the solar wind, or waves and tides originating in the lower atmosphere. The present work discusses the physics in greater detail. It shows that noon NmF2 is closely related to the ambient atomic/molecular concentration ratio, and suggests that the variations of NmF2 with geographic and magnetic longitude are largely due to the geometry of the auroral ovals. It also concludes that electric fields play no important part in the dynamics of the midlatitude thermosphere. Our modelling leads to the following picture of the global three-dimensional thermospheric circulation which, as envisaged by Duncan, is the key to explaining the F2-layer variations. At solstice, the almost continuous solar input at high summer latitudes drives a prevailing summer-to-winter wind, with upwelling at low latitudes and throughout most of the summer hemisphere, and a zone of downwelling in the winter hemisphere, just equatorward of the auroral oval. These motions affect thermospheric composition more than do the alternating day/night (up-and-down) motions at equinox. As a result, the thermosphere as a whole is more molecular at solstice than at equinox. Taken in conjunction with the well-known relation of F2-layer electron density to the atomic/molecular ratio in the neutral air, this explains the F2-layer semiannual effect in NmF2 that prevails at low and middle latitudes. At higher midlatitudes, the seasonal behaviour depends on the geographic latitude of the winter downwelling zone, though the effect of the composition changes is modified by the large solar zenith angle at midwinter. The zenith angle effect is especially important in longitudes far from the magnetic poles. Here, the downwelling occurs at high geographic latitudes, where the zenith angle effect becomes overwhelming and causes a midwinter depression of electron density, despite the enhanced atomic/molecular ratio. This leads to a semiannual variation of NmF2. A different situation exists in winter at longitudes near the magnetic poles, where the downwelling occurs at relatively low geographic latitudes so that solar radiation is strong enough to produce large values of NmF2. This circulation-driven mechanism provides a reasonably complete explanation of the observed pattern of F2 layer annual and semiannual quiet-day variations.     

First in situ measurement of electric field fluctuations during strong spread F in the Indian zone

An RH-560 rocket flight was conducted from Sriharikota rocket range (SHAR) (14°N, 80°E, dip 14°N) along with other experiments, as a part of equatorial spread F (ESF) campaign, to study the nature of irregularities in electric field and electron density. The rocket was launched at 2130 local time (LT) and it attained an apogee of 348 km. Results of vertical and horizontal electric field fluctuations are presented here. Scale sizes of electric field fluctuations were measured in the vertical direction only. Strong ESF irregularities were observed in three regions, viz., 160/190 km, 210/257 km and 290/330 km. Some of the valley region vertical electric field irregularities (at 165 km and 168 km), in the intermediate-scale size range, observed during this flight, show spectral peak at kilometer scales and can be interpreted in terms of the image striation theory suggested by Vickrey et al. The irregularities at 176 km do not exhibit any peak at kilometer scales and appear to be of a new type. Scale sizes of vertical electric field fluctuations showed a decrease with increasing altitude. The most prominent scales were of the order of a few kilometers around 170 km and a few hundred meters around 310 km. Spectra of intermediate-scale vertical electric field fluctuations below the base of the F region (210/257 km) showed a tendency to become slightly flatter (spectral index n = –2.1 ± 0.7) as compared to the valley region (n = –3.6 ± 0.8) and the region below the F peak (n = –2.8 ± 0.5). Correlation analysis of the electron density and vertical electric field fluctuations suggests the presence of a sheared flow of current in 160/330 km region.     

A supplementary approach for estimating reaeration rate coefficients

Different commonly used predictive equations for the reaeration rate coefficient (K₂) have been evaluated using 231 data sets obtained from the literature and 576 data sets measured at different reaches of the River Kali in western Uttar Pradesh, India. The data sets include stream/channel velocity, bed slope, flow depth, cross‐sectional area and reaeration rate coefficient (K₂), obtained from the literature and generated during the field survey of River Kali, and were used to test the applicability of the predictive equations. The K₂ values computed from the predictive equations have been compared with the corresponding K₂ values measured in streams/channels. The performance of the predictive equations has been evaluated using different error estimation, namely standard error (SE), normal mean error (NME), mean multiplicative error (MME) and coefficient of determination (r²). The results show that the reaeration rate equation developed by Parkhurst and Pomeroy yielded the best agreement, with the values of SE, NME, MME and r² as 33·387, 4·62, 3·58 and 0·95, respectively, for literature data sets (case 1) and 37·567, 3·57, 2·6 and 0·95, respectively, for all the data sets (literature data sets and River Kali data sets) (case 2). Further, to minimize error estimates and improve correlation between measured and computed reaeration rate coefficients, supplementary predictive equations have been developed based on Froude number criteria and a least‐squares algorithm. The supplementary predictive equations have been verified using different error estimates and by comparing measured and computed reaeration rate coefficients for data sets not used in the development of the equations. Copyright © 2003 John Wiley & Sons, Ltd.     

Observations of the upper frequency cutoffs of the auroral kilometric radiation

Intense auroral kilometric radiation (AKR) is being frequently observed with POLRAD from the Auroral Probe (Interball-2). Observations of the abrupt upper frequency cutoffs (UFCs) in the spectra of AKR are reported. The UFCs can be observed at a frequency range from 300 to 700 kHz, corresponding to AKR generation altitudes from approximately 4800 to 2100 km, and are distributed in magnetic local time (MLT) hours similarly to the AKR events, with a maximum at 1 h MLT. The observed frequency extent of the UFCs is 12 kHz, and is often determined by the instrumental resolution (4 kHz). It is suggested that the UFC may be associated with an abrupt switching on of the generation mechanism, when the electron density becomes sufficiently low inside a plasma depletion at an altitude where the ratio of f_pe/f_ce crosses some threshold value. The steepness of the UFCs can imply a non-linear process of generation. The estimated distance of the e-folding field aligned wave amplification is between 3 and 8 km. The UFCs are sometimes, though very seldom (10%), accompanied by narrow band (less than 4 kHz) ridges of radiation observed at the cutoff frequency. They are smoothly drifting in frequncy for several minutes. The power density of radiation in the ridge can be up to 2 orders of magnitude stronger than in the accompanying wide band emission of AKR. The ridge at UFC can imply either energy concentration at the source bottom, or focusing, if specific conditions for the escape of the radiation are assumed.     

Semiannual and annual variations in the height of the ionospheric F2-peak

Ionosonde data from sixteen stations are used to study the semiannual and annual variations in the height of the ionospheric F2-peak, hmF2. The semiannual variation, which peaks shortly after equinox, has an amplitude of about 8 km at an average level of solar activity (10.7 cm flux = 140 units), both at noon and midnight. The annual variation has an amplitude of about 11 km at northern midlatitudes, peaking in early summer; and is larger at southern stations, where it peaks in late summer. Both annual and semiannual amplitudes increase with increasing solar activity by day, but not at night. The semiannual variation in hmF2 is unrelated to the semiannual variation of the peak electron density NmF2, and is not reproduced by the CTIP and TIME-GCM computational models of the quiet-day thermosphere and ionosphere. The semiannual variation in hmF2 is approximately isobaric, in that its amplitude corresponds quite well to the semiannual variation in the height of fixed pressure-levels in the thermosphere, as represented by the MSIS empirical model. The annual variation is not isobaric. The annual mean of hmF2 increases with solar 10.7 cm flux, both by night and by day, on average by about 0.45 km/flux unit, rather smaller than the corresponding increase of height of constant pressure-levels in the MSIS model. The discrepancy may be due to solar-cycle variations of thermospheric winds. Although geomagnetic activity, which affects thermospheric density and temperature and therefore hmF2 also, is greatest at the equinoxes, this seems to account for less than half the semiannual variation of hmF2. The rest may be due to a semiannual variation of tidal and wave energy transmitted to the thermosphere from lower levels in the atmosphere.     

Frequency domain interferometry mode observations of PMSE using the EISCAT VHF radar

During the summer of 1997 investigations into the nature of polar mesosphere summer echoes (PMSE) were conducted using the European incoherent scatter (EISCAT) VHF radar in Norway. The radar was operated in a frequency domain interferometry (FDI) mode over a period of two weeks to study the frequency coherence of the returned radar signals. The operating frequencies of the radar were 224.0 and 224.6 MHz. We present the first results from the experiment by discussing two 4-h intervals of data collected over two consecutive nights. During the first of the two days an enhancement of the FDI coherence, which indicates the presence of distinct scattering layers, was found to follow the lower boundary of the PMSE. Indeed, it is not unusual to observe that the coherence values are peaked around the heights corresponding to both the lower- and upper-most boundaries of the PMSE layer and sublayers. A Kelvin-Helmholtz mechanism is offered as one possible explanation for the layering structure. Additionally, our analysis using range-time-pseudocolor plots of signal-to-noise ratios, spectrograms of Doppler velocity, and estimates of the positions of individual scattering layers is shown to be consistent with the proposition that upwardly propagating gravity waves can become steepened near the mesopause.     

Modelling F2-layer seasonal trends and day-to-day variability driven by coupling with the lower atmosphere

This paper presents results from the TIME-GCM-CCM3 thermosphere–ionosphere–lower atmosphere flux-coupled model, and investigates how well the model simulates known F2-layer day/night and seasonal behaviour and patterns of day-to-day variability at seven ionosonde stations. Of the many possible contributors to F2-layer variability, the present work includes only the influence of ‘meteorological’ disturbances transmitted from lower levels in the atmosphere, solar and geomagnetic conditions being held at constant levels throughout a model year.In comparison to ionosonde data, TIME-GCM-CCM3 models the peak electron density (NmF2) quite well, except for overemphasizing the daytime summer/winter anomaly in both hemispheres and seriously underestimating night NmF2 in summer. The peak height hmF2 is satisfactorily modelled by day, except that the model does not reproduce its observed semiannual variation. Nighttime values of hmF2 are much too low, thus causing low model values of night NmF2. Comparison of the variations of NmF2 and the neutral [O/N₂] ratio supports the idea that both annual and semiannual variations of F2-layer electron density are largely caused by changes of neutral composition, which in turn are driven by the global thermospheric circulation.Finally, the paper describes and discusses the characteristics of the F2-layer response to the imposed ‘meteorological’ disturbances. The ionospheric response is evaluated as the standard deviations of five ionospheric parameters for each station within 11-day blocks of data. At any one station, the patterns of variability show some coherence between different parameters, such as peak electron density and the neutral atomic/molecular ratio. Coherence between stations is found only between the closest pairs, some 2500 km apart, which is presumably related to the scale size of the ‘meteorological’ disturbances. The F2-layer day-to-day variability appears to be related more to variations in winds than to variations of thermospheric composition.     

PUBLICATIONS OF THE INTERNATIONAL ASSOCIATION OF HYDROLOGICAL SCIENCES

Abstract

Most commonly used biochemical oxygen demand (BOD) and dissolved oxygen (DO) models have been tested for their applicability in the River Kali, which is one of the most polluted rivers in India. A total of 732 field data sets were generated during field survey from March 1999 to February 2000. The modelling of BOD and DO in the River Kali involves derivation and solution of the governing equations that describe concentration change with time and space brought on by advective, decay, settling and loading functions. However, due to continuous discharges (e.g. from wastewater treatment plants) and steady-state flow conditions in the River Kali, the dispersion effects are found to be insignificant. In the analysis, the model parameters used in BOD-DO models were optimized using the Newton-Raphson technique and the performance of different models was evaluated using correlation statistics (r ²) and error estimation, viz. standard error (SE) and mean multiplicative error (MME). The results indicate that the BOD-DO models developed after Camp (1963) yielded the best agreement with the observed values as compared with several other approaches.     

Using principal component analysis in the investigation of groundwater hydrochemistry of Upper Jezireh Basin,Syria

Principal component analysis (PCA) was applied to hydrochemical and isotopic data of 34 groundwater samples. This allowed the reduction of 20 variables to four significant PCs that explain 81.9% of the total variance; F1 (47.1%) explains the groundwater mineralization, whereas F2 (17%) shows isotopic enrichment and nitrate pollution. Based on an iso-factor scores map of F1, three water zones were delineated: Zone A (F1 < ?1), with fresh groundwater from the unconfined aquifer; Zone B (1 > F1 > ?1), with moderate mineralization from the confined–unconfined aquifer boundary; and Zone C (F1 > 1), with the most mineralized hot water from the confined aquifer. The iso-factor scores map of F2 delineates positive values representing samples from the unconfined aquifer, with freshwater and nitrate contamination associated with stable isotope enrichment, whereas negative values represent samples from the confined aquifer. The results clearly demonstrate the usefulness of PCA in groundwater hydrochemistry investigations.     

Long term ionospheric electron content variations over Delhi

Ionospheric electron content (IEC) observed at Delhi (geographic co-ordinates: 28.63°N, 77.22°E; geomagnetic co-ordinates: 19.08°N, 148.91E; dip Latitude 24.8°N), India, for the period 1975/80 and 1986/89 belonging to an ascending phase of solar activity during first halves of solar cycles 21 and 22 respectively have been used to study the diurnal, seasonal, solar and magnetic activity variations. The diurnal variation of seasonal mean of IEC on quiet days shows a secondary peak comparable to the daytime peak in equinox and winter in high solar activity. IEC_max (daytime maximum value of IEC, one per day) shows winter anomaly only during high solar activity at Delhi. Further, IEC_max shows positive correlation with F_10.7 up to about 200 flux units at equinox and 240 units both in winter and summer; for greater F_10.7 values, IEC_max is substantially constant in all the seasons. IECmax and magnetic activity (A_p) are found to be positively correlated in summer in high solar activity. Winter IEC_max shows positive correlation with A_p in low solar activity and negative correlation in high solar activity in both the solar cycles. In equinox IEC_max is independent of A_p in both solar cycles in low solar activity. A study of day-to-day variations in IEC_max shows single day and alternate day abnormalities, semi-annual and annual variations controlled by the equatorial electrojet strength, and 27-day periodicity attributable to the solar rotation.     

Low frequency waves in plasmas with spatially varying electron temperature

Low frequency electrostatic waves are studied in magnetized plasmas with an electron temperature which varies with position in a direction perpendicular to the magnetic field. For wave frequencies below the ion cyclotron frequency, the waves need not follow any definite dispersion relation. Instead a band of phase velocities is allowed, with a range of variation depending on the maximum and minimum values of the electron temperature. Simple model equations are obtained for the general case which can be solved to give the spatial variation of a harmonically time varying potential. A simple analytical model for the phenomenon is presented and the results are supported by numerical simulations carried out in a 2.5-dimensional particle-in-cell numerical simulation. We find that when the electron temperature is striated along B₀ and low frequency waves (_ci) are excited in this environment, then the intensity of these low frequency waves will be striated in a manner following the electron temperature striations. High frequency ion acoustic waves (_ci) will on the other hand have a spatially more uniform intensity distribution.     

Meteorological drought analysis using data-driven models for the Lakes District,Turkey

Abstract

Droughts may be classified as meteorological, hydrological or agricultural. When meteorological drought appears in a region, agricultural and hydrological droughts follow. In this study, the standardized precipitation index (SPI) was applied for meteorological drought analysis at nine stations located around the Lakes District, Turkey. Analyses were performed on 3-, 6-, 9- and 12-month-long data sets. The SPI drought classifications were modelled by Adaptive Neural-Based Fuzzy Inference System (ANFIS) and Fuzzy Logic, which has the advantage that, in contrast to most of the time series modelling techniques, it does not require the model structure to be known a priori. Comparison of the observed values and the modelling results shows a better agreement with SPI-12 and ANFIS models than with fuzzy logic models.     

Quiet-condition hmF2, NmF2, and B0 variations at Jicamarca and comparison with IRI-2001during solar maximum

We use the measurements of the Jicamarca digisonde to examine the variations in F2 layer peak electron density (NmF2), its height (hmF2), and the F2 layer thickness parameter (B₀) near the dip equator. The hourly ionograms during geomagnetic quiet-conditions for a 12-month period close to the maximum solar activity, April 1999–March 2000, are used to calculate the monthly averages of these parameters, for each month. The averages are compared with the International Reference Ionosphere (IRI)-2001 model values. The results show that the higher hmF2 values during daytime, associated with the upward velocity, are mainly responsible for the greater values of NmF2 and B₀; while the nighttime lower hmF2, related to the downward velocity, are responsible for the smaller NmF2 and B₀. For daytime, hmF2 and NmF2 are correlated with the solar activity in the equinoctial and summer months. The hmF2 and B₀ peaks at sunset with an associated sharp decrease in NmF2 are presented in the equinoctial and summer months, but not in the winter months. Comparison of the measured hmF2 values with the International Radio Consultative Committee (CCIR) maps used in IRI-2001 (IRI-CCIR) reveals an IRI overestimate in hmF2 during daytime. The most significant discrepancy is that the IRI-CCIR does not model the post-sunset peak in hmF2. For the NmF2 comparison, the values obtained from both the CCIR and URSI maps are generally close to the observed values. For the B₀ comparison, the highest discrepancy between the observation and the Gulyaeva option (IRI-Gulyaeva) is the location of the annual maximum for the daytime values, also the winter daytime predictions are too low. Additionally, the significant negative difference between the observation and the B₀-table option (IRI-B₀-table) provides a slightly better prediction, except for 0400–1000 LT when the model significantly overestimates. The post-sunset peak in B₀ at some months is predicted by neither the IRI-Gulyaeva nor the IRI-B₀-table options.     

4D tropospheric tomography using GPS slant wet delays

Tomographic techniques are successfully applied to obtain 4D images of the tropospheric refractivity in a local dense network of global positioning system (GPS) receivers. We show here how GPS data are processed to obtain the tropospheric slant wet delays and discuss the validity of the processing. These slant wet delays are the observables in the tomographic processing. We then discuss the inverse problem in 4D tropospheric tomography making extensive use of simulations to test the system and define the resolution and the impact of noise. Finally, we use data from the Kilauea network in Hawaii for February 1, 1997, and a local 4 × 4 × 40 voxel grid on a region of 400 km² and 15 km in height to produce the corresponding 4D wet refractivity fields, which are then validated using forecast analysis from the European Center for Medium Range Weather Forecast (ECMWF). We conclude that tomographic techniques can be used to monitor the troposphere in time and space.