I. BACKGROUND
When someone who has caught the flu seeks medical attention, the examining doctor will measure their blood pressure with a sphygmomanometer and listen to their breathing to conduct a pulmonary examination. In mammals and birds infected with the avian influenza virus, the immune system is activated and chemical messengers called cytokines are released from the mast cells of the tissue to fight the virus, increasing the permeability of the vascular wall and allowing amino acids to escape from the blood vessels. A mechanism that facilitates leakage activated. However, when the number of cytokines increases, blood pressure decreases, creating a cytokine storm. In this condition, the lungs become severely edematous and breathing becomes shallow and rapid. For this reason, in humans, the examining physician or nurse monitors blood pressure and observes breathing patterns.
Several devices for humans based on radar Doppler using continuous waves reflected from the skin of the chest have been developed. Radar Doppler is not limited to medical care; it is also used in health and sleep management. Through vibrations, it can sense whether a person is awake and moving or sleeping. In humans, since the vibration of the heart is transmitted to the chest wall, the heartbeat can be obtained from vibrations of the skin. High-frequency microwaves (10−77 GHz) are reflected at the surface of the body, so the propagation of radio waves is only in the air. The amount of power attenuation is relatively small for this Doppler system, so the microwave transmission power can be very low. However, if the inputs and reflected waves (nonlinear) are simply mixed with each other, the detected signal will be nonlinear waves. To compensate for this shortcoming, a method of simultaneously interfering waves with a phase shift of λ/4 (λ: wavelength) can be used; however, the improvement is limited. Since it is not possible to obtain waveforms of heart vibrations that hidden under the breathing waves, the frequency peaks (first peak is breathing rate and second peak is heart rate) are obtained by performing the fast Fourier transform (FFT).
Since FFT is an integration of the time axis, the resolution of the time axis is reduced. Therefore, in the RR interval, which requires the accuracy of the time axis, the waveform derived from the heart motion must be extracted, IFFT performed, the original waveform drawn, and the fluctuation of the RR interval must be measured. Since the evaluation of the RR interval is the standard deviation of the interval, accuracy cannot be obtained unless the heart-derived signal is large enough. Regarding radar Doppler, there are no reports on obtaining a blood pressure waveform even if the frequency is lowered to irradiate the inside of the thoracic wall.
Measurement of esophageal impedance for chickens (monitoring of blood pressure and breathing patterns)
We previously developed a device for chickens that can monitor the cardiac output as the phase and breathing changes as the magnitude by adding a 1-MHz medium wave through electrocardiogram electrodes of the esophageal catheter [1-3].
In addition to electrocardiogram electrodes, this avian esophageal catheter incorporates an angular velocity sensor that correlates with stroke volume and a heart sound microphone that can confirm the closing sounds of the aortic and mitral valves.
Based on these physiological data, we can determine the blood flow of the cardiac output and the movement of the thoracic wall associated with respiration. Fig. 1 shows an example of measurement using an esophageal catheter.
Fig. 1 shows (from top to bottom) the invasive axillary artery pressure, magnitude (impedance), angular velocity, electrocardiogram, phase (impedance), and heart sound. The phase waveform shows the outflow from the left ventricle and blood flow into the atria. The closing sounds of the mitral valve and the aortic valve can be confirmed. Since the closing sounds of the aortic valve match the notch timing of the waveform (shape) of the phase, the phase can be determined as the blood flow originating from the heart. Regarding the peak of the axillary artery pressure (mea- sured by invasively cutting down the axillary artery), a propagation delay was observed comparison from the peak of the phase waveform. The difference (L) between heart and axillary artery is the blood flow velocity V=L/Δt. Our previous paper demonstrates that there is a correlation between the magnitude of the amplitude of the angular velocity and the magnitude of the peak of blood flow [3].
In this paper, a carrier wave of 2.4 GHz is transmitted from outside the body, and microwaves that have undergone frequency modulation and amplitude modulation effected by the heart and chest wall, and the phase (frequency) and magnitude (amplitude) are separated and demodulated by the receiving device. Even if the voice is modulated with a medium wave carrier wave on 1 MHz, the voice can be demodulated in the same way even if the voice is modulated with a microwave on 2.4 GHz. Therefore, it is estimated that even if microwaves are used, the impedance waveform will resemble that of 1 MHz of the esophageal catheter. The difference is the site of measurement, that is, the spatial difference. Due to the difference between the measurement-site at a position close to the heart or ribcage and the measurement-site at a distance outside the body, the latter is generally considered to be physically blunted in the waveform. The solution of spatial loss, that could be solved by antenna directivity and penetration of tissue. The longer the wavelength, the better the penetration, and the shorter the wavelength, the sharper the directivity, so there is a spatial gain. These two factors are relationship of tread-off. A practical product for poultry farmers to understand the health status of their birds must be contactless. Because chickens are handled in cages. With this in mind, this paper describes the method and results of how to find the trend-off conflict.
The circulation and respiratory monitoring using with microwave are methods that do not contact the individual, and are spatially separated from the heart, lungs, and air sacs. So what does the microwave waveform mean? Is the temporal timing consistent with other multimedia data? These points must be watched. Occasionally, it may be observed with delays or skeletal muscle spasms as a noise. In order to evaluate, multimedia data measured in multiple dimensions, are essential. The further the distance from the heart, the more supplementary correction data is required. Therefore, the evaluation of microwave waveforms shall be recommended with data analysis by multimedia.
Birds do not have a diaphragm, so skeletal muscles in the abdomen and chest control balloons called air sacs to direct airflow into the lungs. Avian lungs are hard, and in addition to ventilation, which draws in outside air through the trachea, airflow is also exchanged between air sacs. Therefore, airflow flows into the lungs multiple times in one ventilation. Skeletal muscles are always in a state of tension, which suggests that the sympathetic nerves, which control the excitation state of skeletal muscles, maintain a dominant state in the autonomic nerves. The bar-tailed godwits tacked by GPS, a migratory bird that flies for 14 days without rest and sleep, flying from Alaska to the off-shore of New Zealand [28]. It is extremely difficult for humans, mammals, to run a marathon 24 hours a day, but birds have a physiological structure that allows them to maintain sympathetic excitation at all times. Therefore, in respiratory monitoring, skeletal muscle spasm may appear in the air sac pressure, but this is not noise but important information controlling the flow of air between air sacs. Therefore, air sac pressure monitors must be amplified even at relatively high frequencies without noise rejection.
II. PRINCIPLE
Here, we explain the principle of monitoring the blood flow associated with the heartbeat and the breathing pattern associated with the movement of the thorax. The relationship between the free-space propagation loss and frequency can be expressed as
c: speed of light (m/s)
f: frequency (Hz)
R: distance (m)
λ: wavelength (m)
From this equation, the attenuation is 29.77 dB at 2,450 MHz and a distance of 30 cm.
The proposed device passes microwaves through the bird’s tissue and measures changes in impedance. The tissue attenuates the microwave power. Under the assumption that the depth at which power is attenuated by half is D, the relationship between frequency f, dielectric constant εr, and the dielectric dissipation factor (tan δ) is as follows.
The dielectric loss tangent is a numerical value that expresses one aspect of the ratio at which part of the energy becomes heat when an AC electric field is applied to a dielectric. If an AC electric field is applied to an ideal dielectric and only the capacitance component (capacitor) is considered, there is a 90-degree phase difference between the electric field and the current due to polarization/charge reversal.
The power consumption of the AC current is expressed as the vector product of the electric field and the current, so if the phase difference between the two is 90 degrees, the power consumption will be zero. It becomes an ideal dielectric. However, in actual materials, complete insulation is not possible. The resistance of the dielectric itself is finite, so a leakage current is generated. In addition, dielectric polarization and partial discharge create a parallel resistance component. The phase difference is thus shifted from 90 degrees. Let δ denote the shift from 90 degrees. Its tangent (tan δ), called the dielectric tangent, corresponds to the component where the electric field and current are in phase. Only this in-phase component contributes to power consumption or loss. The dielectric loss tangent is proportional to the dielectric loss [4-13]. In addition, the complex dielectric constant due to the parallel plate capacitor due to the difference in materials is described in the literature [14], the difference in the frequency characteristics of pure water and blood is described in [15], and the effects of weak microwave exposure are described in [16].
The dielectric constant of water is about 80. In living organisms, glucose in body fluids decreases the dielectric constant. Under the assumption that a body fluid has an average glucose content of 5%, the dielectric constant is 76.5 at around 2,450 MHz, as shown in Fig. 2.
Furthermore, the non-dielectric constant of the chest is lowered by the air in the lungs, connective tissue fibers, fatty tissue, etc., so the average non-dielectric constant can be considered to be 25.5 or lower. For the estimation of the attenuation of microwaves, a factor of 25.5 may be chosen to maximize attenuation.
From Fig. 3, the power half-power depth is 26.2 mm with a non-dielectric constant of 25.5, so it is 4 times this (2 to the power of 4=16) 26.2×4=104.8 mm Approximately 105 mm has 1/16 of the power. The power is attenuated by −12 dB corresponds to half-power loss (1/16=−12 dB). It is assumed here that 105 mm is approximately the width of the anterior chest of a large chicken.
For microwaves, a higher frequency leads to more energy being reflected by the surface of a substance. Less energy is reflected when the wave is incident at right angles to the surface. Whether the reflective surface of an uneven material can be regarded as flat is determined by the Rayleigh roughness criterion, expressed as
where σ is the standard deviation of the undulation amount in the first Fresnel zone and θ is the incident angle with respect to the normal direction of the surface.
If g < 1, the reflective surface is flat; if g>1, it is undulating. In chickens and pigeons, the probe is placed perpendicular to the chest wall.
Since the wavelength for 2.5 GHz waves is 12 cm and a bird’s skin is smooth, g < 1. In experiments with chickens and pigeons, the probe was intentionally set perpendicular to the chest wall, with a wavelength of λ 12 cm at a freq- uency of 2.5 GHz. Since the denominator of the above formula is large enough, the skin of birds is smooth, and g<1 is self-evident. There is no problem numerically assuming that scattered waves are relatively small at 2.5 GHz. However, the beam half width of the 10 dBi antenna is 49 degrees, so the scatter increases geometrically (it is not a pointing loss). Based on the past propagation experiments at 2.5 GHz show that the beam direction of antenna is incident perpendicularly (almost 90 degrees) and reflected or scattered by the skin, resulting in a halving of the power (i.e., a 3 dB drop).
In the equivalent circuit of avian thorax, Rbody is the change in distance mainly due to the movement of the ribcage. This change is very large. Many radar Doppler systems track only real numbers and blood flow cannot be seen logically because they calculated only magnitude (real number). Even if the change in blood flow is reflected in the real number ΔR, a small ΔR cannot be detected unless the dynamic range of the radar Doppler system is increased and the amplification circuit is not linear. The amount of change in blood flow on the real axis (Rleakage in Fig. 4) is about 1/100 of the resistance value of the thorax.
In addition, the notch caused by the aortic valve closure cannot be seen in the real axis direction change unless the resolution is high. Previous studies have thus applied the FFT power spectrum to the received waveform. The first peak is taken as the respiratory rate and the second peak is taken as the heart rate.
In birds, the heart is covered by a large sternum and thus the vibrations of the heart are only slightly transmitted to the skin. The valve closing sound is even difficult to detect using a stethoscope. Reflective microwave radar, which detects skin vibrations, thus cannot be used for birds.
Using the equivalent circuit shown in Fig. 4, it is possible to measure the change on the complex plane, where the imaginary part is the phase and the change in the real number is the magnitude, using transmission-type microwaves (continuous waves), as described later.
Specifically, when the transmit and receive antenna gains are 10 dBi, the frequency is 2,450 MHz, and the distance is 30 cm, the transmit power is −40 dBW and the received power is −67.76 dBW. The theoretical margin is 40.75 dB. Such a margin is difficult to obtain with radar Doppler.
Chickens are constantly moving and thus there is a pointing loss of the antenna toward the heart. In addition, the skeletal muscles of chickens constantly vibrate to control skeletal muscles which keep tension of air sacs [2], and those vibration causes fluctuations in the baseline of noise level. The actual measurement margin is about 30 dB, which is different from the theoretical value, but this 30 dB margin is necessary to draw the phase waveform including the notch of the aortic valve. Table 1 shows the line design of the system considering the free-space propagation loss, half-power depth, reflection/scattering, and pointing loss.
III. METHODS
The conceptual diagram and block diagram of the implemented system are shown in Fig. 5. The reflection coefficient (S11) of the antenna built into the probe is shown in Fig. 6. The implemented device is shown in Fig. 7. The experimental conditions are summarized in Table 2.
It is necessary to measure the heart and respiration with another sensor to corroborate the microwave measurements. In pigeons, an air pressure sensor and an electrocardiogram electrode were surgically sewn subcutaneously in the sagittal direction in the midline of the back under isoflurane inhalation anesthesia in advance. The implemented transmission-type microwave device for pigeons (chest antenna: horizontal polarization, abdomen antenna: vertical polarization) is shown in Fig. 7.
For chickens, an esophageal catheter can be inserted through the crop to the side of the heart. While maintaining inhalation anesthesia and spontaneous breathing, the operator guided the esophageal catheter to pass through the crop and into the lower esophagus. Isoflurane was then discontinued, sufficient air was aspirated, and the chicken was gently aroused and kept upright for the transmission-type microwave examination and data acquisition. The esophageal catheter was removed when it was noticed by the chicken. Data were then collected continuously. In addition, pressure sensors with silicon tubes were surgically inserted into the air sacs.
Respiration-induced changes in received power and heartbeat-synchronized waveforms were obtained from all target birds. Fig. 8 shows changes in the received magnitude waveform (top) and changes in the abdominal air sac pressure obtained by the pressure sensor, both of which are synchronous.
Similarly, the change in the received power obtained from the chest probe of the pigeon (Columba livia), the FFT of the signal, the inverse FFT of the frequency components below 17 Hz, and the differentiation of this waveform are shown in Fig. 9. The waveform and the electrocardiogram obtained at the same time are also shown.
A heartbeat is obtained from a chest probe. The measured signal was about 20 dB lower than the calculated value. The reason for this was clarified during the post-experiment autopsy. Homing pigeons used for racing have well-developed breast wall meat, which is not much different from chicken breast wall thickness, and their wing meat is thicker than that of chickens. Nevertheless, a margin of about 30 dB was secured.
Eight chickens were viewed using a circularly polarized probe and compared with esophageal catheter data. In Fig. 10, photograph shows the measurements taken on chest side of chickens.
The obtained magnitude waveform (black) is shown in Fig. 11. This waveform is strongly correlated with the anterior thoracic air sac pressure (yellow). In birds, the breathing pattern is inspiration, short duration, expiration, then long duration (one cycle) and both curves (black and yellow) suggest that principle. Amplitude modulation of the anterior thoracic air sac, indicated by asterisks, can be observed as turbulent airflow from other air sacs (posterior air sacs) to the precordial air sacs via the lungs. Such vibrations are due to phase differences in skeletal muscle contractures, which are described in detail in the paper [2]. Birds, unlike mammals, are efficient with more than one airflow to the lungs per breath.
The magnitude waveform has a strong correlation with the anterior thoracic air sac pressure. It is speculated that the skeletal muscles that generate the anterior thoracic air sac pressure strongly affect the movement of the chicken’s chest wall. In more detail, it can be inferred from this graph that the breathing control muscles are different between the star-marked time period and the time period of inspiration and expiration.
Fig. 12 shows the phase waveform of the microwaves (purple). There is a notch in the middle of the mountain-shaped waveform, as indicated by an arrow. This notch is synchronized with the aortic valve closing sound from the heart sound (green).
The upper part of the graph shows the angular velocity of the axial rotation of the esophagus. Previous studies [3] have shown that its amplitude correlates with blood pressure.
When the measurement point was determined at the timing of breathing and the correlation between the amplitude of the angular velocity and the amplitude of the microwave phase was obtained, a linear correlation was established at R2=0.95 (Fig. 13).
This indicates that the microwave phase reflects the blood flow in the heart and that the notch measures the closing of the aortic valve and the subsequent inflow into the left atrium.
IV. CONSIDERATIONS
Many experiments have been conducted on heartbeat/ respiratory monitors using microwaves. We have developed a monitor that has a strong correlation between the cardiac output flow and the pressure of the anterior thoracic air sac in birds using transmission-type microwaves.
We can obtain not only the cardiac output flow but also the interval between beats (RR interval). This interval reflects the condition of autonomic nerves (sympathetic and parasympathetic nerve control) and is an important parameter for confirming health conditions. For radar Doppler, it is difficult to accurately measure interval fluctuations because increasing the frequency resolution leads to frequency components at intervals of Δt that are at least 1,024, 2,048, or more times the sampling time.
Therefore, the interval of time Δt becomes at least an error factor for time domain. Furthermore, since the RR interval evaluates the standard deviation of the 100 intervals, the error in Δt is included in the result of standard deviation. Therefore, the RR interval obtained from the heart rate by FFT using reflected radar Doppler is unreliable.
Since it is not fair to discuss only our strengths unilaterally, we will also introduce the advantages of the Doppler radar system. The biggest advantage is that only one antenna is required, so the high-frequency structure is extremely small. The distance setting range is also wide, so even if the target moves, it does not matter where it is within the line gain range. However, the target must be stationary at the moment of measurement [14-18]. To increase spatial resolution, increase the frequency. In this case, the penetration rate into the tissue decreases, so the distance should be shortened and the transmission power should be increased [19-20]. Furthermore, a detection method that emphasizes phase has been proposed to suppress magnitude wave due to respiration [21-22]. The paper [27] captures the waveform correlated with the blood flow with the reflection type, blood flow of intracranial space [7], and the paper [17] also reports the paper using the transmission microwave. A comparison of radar Doppler and the proposed method is given in Table 3 [14-26].
Birds do not have a diaphragm and their lungs are hard and thus barely expand. They thus breathe by contracting balloon-like air sacs that expand in their chest and abdomen.
The avian great vessel is the right descending aorta. The cells of the left ventricle and the interventricular septum look morphologically abnormal. The muscle fibers are short and thick. The normal avian heart presents macro-histological abnormalities in macroscopic tissue sections. T these abnormalities look similar to hypertrophic myocarditis in humans.
Cardiac ejection fraction is close to 100%, so the branch of the His bund from the atrial node is running outside of left ventricular descends laterally rather than medially. The heart rate of chickens is 200−300 beats/min, Hummingbird has more than 1,000/min, which is much faster than that of humans. Birds cannot accurately transmit electrical stimulation to the left ventricle with the same conduction pathway as mammals. Therefore, the outside of the left ventricular wall is a stimulus conduction path (myocardium degenerated like a nerve fiber).
The autonomic nervous system is controlled by a battle between the sympathetic and parasympathetic nervous systems. The indices are the standard deviation of the RR interval, systolic blood pressure, and heart rate. Parasympathetic nerves in birds are weakly suppressed and many birds suffer myocardial damage and large blood vessel damage when excited. Air sacs are controlled by the movement of respiratory muscles (skeletal muscles), but the subcutaneous tissue of birds has extremely weak connections and adhesions with the fascia, making it difficult to visually detect changes associated with respiration. Therefore, it is not possible to screen infected birds based on their appearance alone.
Infection in chickens with a highly pathogenic avian influenza virus propagates in type II alveolar epithelium and causes a cytokine storm. This leads to tachycardia, low blood pressure, shallow breathing, increased breathing rate, and slower skeletal muscle contraction. A diagnosis should be confirmed by a polymerase chain reaction test. However, at poultry farms, a simple antigen-antibody reaction is used. It takes time for the influenza virus antibody titer to reach a positive level. In settings with a lot of chickens, if the cardiac output and respiratory pattern of poultry could be monitored without restraint and contact, infected birds could be screened early and isolated from the population, preventing the culling of many birds.
The hosts of avian influenza viruses are mainly waterfowl such as ducks and swans, not chickens. A low-virulence influenza virus replicates slowly in the gut rather than the lungs, with a potential prevalence of up to 30% in nature. In poultry, on the other hand, the influenza virus has been found to spread rapidly from the intestine to the lungs, and it has been reported that even low pathogenicity causes mild respiratory symptoms. Screening for early mild respiratory distress is therefore important in the poultry industry. We recommend a system that can continuously monitor cardiac output and breathing patterns without contact.
V. CONCLUSION
We developed a device that can monitor the cardiac output and anterior thoracic air sac pressure in birds using transmission-type microwaves (2.4−2.5 GHz). Compared to reflector-type radar Doppler devices, the proposed device can visualize blood flow and continuously measure breathing patterns that are correlated with the anterior thoracic air sac pressure. Data collection allows the early screening of unwell birds, which can be isolated to prevent mass infections in poultry houses and prevent pandemics.