<p>Thanks folks for helping.
I knew about SPK files but from Horizon's telnet inteface there is a limit of 200 bodies per request.
Making over 5000 request might be possible but not very productive.
I also looked at the  DE421 file from <a href="ftp://ssd.jpl.nasa.gov/pub/eph/planets/bsp/" rel="nofollow noreferrer">ftp://ssd.jpl.nasa.gov/pub/eph/planets/bsp/</a> but it looks like the corrections for the 8+1 planets</p>
<p>here is the output from BRIEF, a SPICE utility: <a href="https://naif.jpl.nasa.gov/naif/utilities.html" rel="nofollow noreferrer">https://naif.jpl.nasa.gov/naif/utilities.html</a></p>
<pre><code>BRIEF -- Version 4.0.0, September 8, 2010 -- Toolkit Version N0066
Summary for: de421.bsp
Bodies: MERCURY BARYCENTER (1)  SATURN BARYCENTER (6)   MERCURY (199)
    VENUS BARYCENTER (2)    URANUS BARYCENTER (7)   VENUS (299)
    EARTH BARYCENTER (3)    NEPTUNE BARYCENTER (8)  MOON (301)
    MARS BARYCENTER (4)     PLUTO BARYCENTER (9)    EARTH (399)
    JUPITER BARYCENTER (5)  SUN (10)                MARS (499)
    Start of Interval (ET)              End of Interval (ET)
    -----------------------------       -----------------------------
    1899 JUL 29 00:00:00.000            2053 OCT 09 00:00:00.000
</code></pre>
<p>I can always do a statistical estimate, any close approach being false but the overall count being close to reality.
Too bad to drop the project, I was getting nices pictures from it. The coloring is done against eccentricity, red=0</p>
<p><a href="https://i.sstatic.net/hjCsH.jpg" rel="nofollow noreferrer"><img src="https://i.sstatic.net/hjCsH.jpg" alt="enter image description here" /></a></p>

CC BY-SA 4.0
2021-04-08T13:43:52.957Z
**1.** It would be a shame to drop this project. Maybe contact Horizons, explain what you want to do, and see if they have any suggestions. They may be able to supply the SPK files in bulk (depending on what time span you want to cover), or have some other practical suggestions. **2.** Since you're investigating near encounters, I guess even the SPK files aren't adequate when two asteroids are close enough to have significant gravitational affect on each other, unless one or both are in the set of 343 asteroids included in the dynamical model.
PM 2Ring
CC BY-SA 4.0
2021-04-08T15:27:58.557Z
well, thanks PM, you can give it a ride at http://geomaitre.com/geomaitre/

You must have a graphic card that support opengl 4.5 most modern one do

Good luck
Tiblemont
CC BY-SA 4.0
2021-08-04T20:25:17.250Z
This image is quite compelling! How are the colors assigned? Have you thought of doing some animations? I see two red groups of Trojan asteroids at Sun-Jupiter L1 and L2 but are there any at L3? What made the ones at L1 and L2 red? Scott Manley has several older videos in YouTube from a previous lifetime, cf. [Asteroids In Resonance With Jupiter](https://www.youtube.com/watch?v=yt1qPCiOq-8) and also [Why did the distribution of asteroids discovered in 2010 have a radial modulation?](https://astronomy.stackexchange.com/q/25736/7982)...
uhoh
CC BY-SA 4.0
2021-08-04T20:33:52.220Z
and answers to [Orbits of Trojan Asteroids](https://astronomy.stackexchange.com/q/43322/7982) and [Code/Package for multiple planet simulation](https://astronomy.stackexchange.com/q/35056/7982) I think you are at the beginning of something really interesting and useful and can take this in several directions. Also see answers to [What is a thermal distribution of eccentricities?](https://astronomy.stackexchange.com/q/34535/7982) Maybe you can "measure the temperature" of of the asteroid belt, or plot average "temperature" as a function of radius, or see of Trojans are "hot" or "cold"?
uhoh
CC BY-SA 4.0
2021-08-04T21:11:08.227Z
See also upcoming answers (hopefully) to [Looking for that cool astronomical and spaceflight simulation website that's been discussed here](https://space.meta.stackexchange.com/q/2839/12102) and despite the rather blatant narcissism in the beginning this video tells quite an interesting story: [Why I Have a Real Planet Named After Me - (33434) Scottmanley](https://www.youtube.com/watch?v=oIQ2byD6b6o)
uhoh
CC BY-SA 4.0
2021-08-04T22:07:25.413Z
The coloring is done against eccentricity where 0 is red using HSL color  [0 1] for the hue part.
The trojan are close to Jupiter and those away from it are Hilda . see https://en.wikipedia.org/wiki/Hilda_asteroid
Tiblemont
CC BY-SA 4.0
2021-04-08T12:13:38.357Z
2021-04-08T12:31:00.387Z
2021-04-08T12:31:00.387Z
Tiblemont
<p>I finally computed positions of over 1 million asteroids for a few months from June to October 2021.
There is significantly less close approaches than I first thought, but on the average, there is a conjunction under 10 thousand kilometers almost every other day.</p>
<p>For example on 2021-08-10 20:05 asteroids (2002 GD77) and (2008 TD184) will be at 825 km according to Nasa JPL Small Planets Center and at 290 millions km from earth.</p>
<p>The figure below show  their orbits and positions relative to the sun(20x) and some major planets(1000x)</p>
<p><a href="https://i.sstatic.net/WAmln.jpg" rel="nofollow noreferrer"><img src="https://i.sstatic.net/WAmln.jpg" alt="enter image description here" /></a></p>
<p>There will be a closer approach at 745 km on 2021-08-16 17:15 between (2003 WH156) and (2014 WR26) but at 520 millions km from earth. Both conjunction have an approach speed of 5.5 km/sec.</p>
<p>My results are consistentt with JPL website and should be close to reality, but as Jon Giorgini (senior analyst at JPL, <em><strong>thanks again for your help, Jon</strong></em>) remarks, due to the uncertainty in the keplerian elements you cannot tell where <em>exactly</em> the asteroids where at a certain time.</p>
<p>Now these results bring more question than answers.</p>
<ol>
<li>Is it possible for high end amateur instruments to detect/photograph such a conjunction?</li>
<li>Is it possible to compute in advance the minimum-maximum variation in asteroid seraration?</li>
<li>Is it possible with different equipment/sofware/data to tell if that event occured or not?</li>
<li>How close should an approach be for the orbits being modified significantly?</li>
<li>If the orbits are being modified significantly, how long would it take to be detected from earth?</li>
<li>There is also a whole bunch of &quot;couples&quot; that have VERY similar orbits and travel as old companions at very close distances. Are they duplicates or realy different but bound together by some resonance?</li>
</ol>
<p>Example:</p>
<pre><code>2435 km     2021-08-10 14:27:08 (1999 FM92)  &amp; (2010 DH114)
824 km      2021-08-10 20:05:00 (2002 GD77)  &amp; (2008 TD184)
8508 km     2021-08-11 06:44:22 (2006 CK14)  &amp; (2015 CY70)
2326 km     2021-08-11 19:40:46 (2017 RD134) &amp; (2019 JH75)
852 km      2021-08-12 05:58:43 (2010 BO137) &amp; (2015 FB346)
</code></pre>
<p><div class="youtube-embed"><div>
                <iframe width="640px" height="395px" src="https://www.youtube.com/embed/VZbItaf9VAc?start=0"></iframe>
            </div></div></p>

CC BY-SA 4.0
2021-08-04T20:12:22.453Z
This is such a cool project and tool that you've created I think there are a lot of interesting things you can address. The five new questions you've enumerated are really launching points for investigations as well as new Stack Exchange questions. For item #4 there's no real cut-off for significant vs. non-significant for a $1/r^2$ force, but I am sure there's been a lot of work on *perturbation of orbital elements* vs distance and nature of close approach and that could be asked separately, and of course the masses of some asteroids are known by their perturbations on other asteroids!
uhoh
CC BY-SA 4.0
2021-08-04T20:14:36.710Z
And of course that starts to address #5 as well. You may find the following questions in Space SE interesting as well:  [Algorithmic methods or techniques to find conjunctions in high N Keplerian element ensembles?](https://space.stackexchange.com/q/29946/12102) and to a lesser extent [Algorithmic methods or techniques to find conjunctions in large ensembles of state vectors?](https://space.stackexchange.com/q/33161/12102). In artificial satellite orbital mechanics "conjunction" is a understated term for satellite near-misses or collisions.
uhoh
CC BY-SA 4.0
2021-08-05T22:02:57.157Z
*Il est fort amusant* indeed! For YouTube links you can just leave the bare url in the post and SE will automatically embed the video once you accept the change, no need to manually create a link.
uhoh
CC BY-SA 4.0
2021-08-04T18:13:18.967Z
2021-08-07T12:22:58.873Z
2021-08-07T12:22:58.873Z
Tiblemont
<p>For those interested in observing conjunctions, here is a partial list of events with approaches under 10,000 km for the coming weeks.</p>
<pre><code>km     = distance between asteroids 
km e6  = distance from earth in  million km 
km/sec = relative speed of approach 
elong  = elongation of asteroids from sun

        |    UTC  time     |asteroid A | km  | asteroid B |km e6| km/sec |elong|
        |------------------|------:----|----:|-----:------|----:|-------:|----:|
        | 2021-08-20 00:43 | 2006 CK14 |*8668 | 2015 CY70 | 480 | 0.0003 | 109 |
        | 2021-08-26 23:48 | 2000 SE89 | 9566 | 2017 XG51 | 432 | 3.4653 | 161 |
        | 2021-09-04 10:31 |    Dyer   | 3660 | 2014 WG584| 420 | 2.5170 | 143 |
        | 2021-09-09 00:09 | 2005 CO9  | 3305 | 2009 KZ1  | 387 | 6.3235 |  94 |
        | 2021-09-18 21:46 | 2008 FW83 | 6355 | 2014 EQ117| 417 | 3.7596 | 138 |
        | 2021-09-28 01:42 | 2007 TZ425| 3475 | 2008 XN66 | 363 | 3.0231 |  97 |
        | 2021-09-28 19:21 | 2001 DL93 | 2670 | 2006 AR3  | 291 | 17.150 |  13 |
        | 2021-09-29 07:02 |   Utkin   | 3575 | 2004 SN2  | 323 | 2.2074 | 154 |
        | 2021-10-07 07:57 | 2013 UH30 | 2396 | 2014 OW110| 309 | 2.5056 | 101 |
        | 2021-10-10 01:58 | 1991 LA1  | 4809 | 2011 LR7  | 424 | 3.5947 | 159 |
        | 2021-10-10 15:18 | 2002 JE35 | 3091 | 2013 PS4  | 379 | 3.7531 | 102 |
        | 2021-10-14 02:17 | 2006 BN154| 8145 | 2009 FB34 | 380 | 4.4304 |  93 |
        | 2021-10-16 16:35 | 2002 GN160| 2826 | 2010 LD118| 462 | 7.8022 | 137 | 
</code></pre>
<p>The first entry shows 2 asteroids with very similar orbits, traveling together as old chaps. You can tell by their closing speed near zero. Try your astronomy software, tracking them for a few months.</p>
<p>If you don't have any software, you can get ephemeris from JPL by typing a command like that in your web browser (no space , no line feed)</p>
<p><a href="https://ssd.jpl.nasa.gov/api/horizons.api?format=text&amp;COMMAND=%27DES=2000%20SE89%27&amp;MAKE_EPHEM=%27YES%27&amp;CENTER=%27500@10%27&amp;TABLE_TYPE=%27VEC%27&amp;VEC_TABLE=%271%27&amp;START_TIME=%272021-08-26%2023:45%27&amp;STOP_TIME=%272021-08-26%2023:50&amp;STEP_SIZE=%271%20m%27&amp;CSV_FORMAT=%27YES%27" rel="nofollow noreferrer">https://ssd.jpl.nasa.gov/api/horizons.api?format=text&amp;COMMAND=%27DES=2000%20SE89%27&amp;MAKE_EPHEM=%27YES%27&amp;CENTER=%27500@10%27&amp;TABLE_TYPE=%27VEC%27&amp;VEC_TABLE=%271%27&amp;START_TIME=%272021-08-26%2023:45%27&amp;STOP_TIME=%272021-08-26%2023:50&amp;STEP_SIZE=%271%20m%27&amp;CSV_FORMAT=%27YES%27</a></p>
<p>Note the high approach speed of 2021-09-28 19:21 |  2001 DL93 | 2670 | 2006 AR3  | 291 | 17.1502|  46 |</p>

CC BY-SA 4.0
2021-08-25T19:04:20.807Z
2021-10-16T12:08:16.787Z
2021-10-16T12:08:16.787Z
Tiblemont
<p>I am currently writing a C++ program to show asteroids in 3D, and find close approaches or collisions. I got my orbital elements from JPL <a href="https://ssd.jpl.nasa.gov/sbdb_query.cgi" rel="nofollow noreferrer">https://ssd.jpl.nasa.gov/sbdb_query.cgi</a>. So far so good, with over a million asteroids drawn in 5ms with my old Quadro GPU</p>
<p>But the problem is the variation of the elements over time. For example, CERES:</p>
<pre><code>epoch,a,e,i,om,w,ma 
2459200.5,2.766089105818,.07816842657453,10.58789954719,80.27235841368,73.72488984426,205.5454154582
2459310.5,2.765760313090,.07831877879848,10.58807660401,80.26860808947,73.73699886586,229.1146825391
</code></pre>
<p>After only 4 months, the eccentricity and all other elements have changed as seen on the second line.</p>
<p>How can I compute the change in the parameters without knowing their derivatives?</p>
<p>Or where to get those derivatives?</p>
<p>I searched all JPL small bodies site and their telnet access or email, but could not find a way to download orbital elements with derivatives or a way to compute the change.</p>
<p>I am aware that the difference might be just a pixel on the screen but it represents hundreds of thousands of kilometres. What I want to do is study the close approach between the asteroids themselves and eventually near collision. I already did that for artificial satellites using SGP4 propagator using OpenCL running on GPU video card. I can propagate 20,000 satellites (include debris) in a few milliseconds and get results exact to 1 kilometre. Compare with Celestrak Socrate. A prediction of a 1,000-kilometre approach must be possible.</p>
<p>Does anybody know how?</p>

CC BY-SA 4.0
2021-04-07T00:57:54.007Z
few quick questions: 1) if you are drawing a million asteroids will anybody know that a few are a half-pixel off? 2) wait  you downloaded orbital elements for *a million* asteroids from JPL? 3) what are the uncertainties on those osculating orbital elements? 4) if you had all their derivatives, are you confident you would know how to use them to propagate correctly? 5) could there be a reason that osculating elements don't come with derivatives?
uhoh
CC BY-SA 4.0
2021-04-07T01:53:51.300Z
Stupid answer from me: Grab every single object in the solar system at epoch `J2000` and run a huge nbody simulation over it. If you compare that to your given parameters and it is reasonably accurate, keep doing it. But remember that a butterfly flapping its wings can cause a hurricane!
slowerthanstopped
CC BY-SA 4.0
2021-04-07T01:57:37.800Z
Salut François! The equations describing the precession of nodes and the precession of apsides are quite complex. I doubt you’d want to perform them for millions of bodies! As @uhoh pointed out, nobody will notice that your asteroids are “a half-pixel off” (if that much!). Unless you want to go into the distant past or distant future, I see no use for that.
Pierre Paquette
CC BY-SA 4.0
2021-04-07T02:41:13.950Z
If the osculating orbital elements of asteroids change a lot it's mostly because of Jupiter or other big planets. Those calculations will also sometimes include effects of sunlight.
uhoh
CC BY-SA 4.0
2021-04-07T04:24:32.247Z
general references [Where can I find the positions of the planets, stars, moons, artificial satellites, etc. and visualize them?](https://astronomy.stackexchange.com/q/13488) and [Is there free, accurate and updated data on planetary body positions in the solar system?](https://space.stackexchange.com/q/13956) and for the effect of sunlight see Marsden parameters [here](https://space.stackexchange.com/a/18438) and [here](https://space.stackexchange.com/q/18474).
uhoh
CC BY-SA 4.0
2021-04-07T04:27:13.310Z
For an example of how to propagate state vectors instead of using orbital elements see [this answer](https://space.stackexchange.com/a/234090 though you would want to do a simpler propagation than this one.
uhoh
CC BY-SA 4.0
2021-04-07T18:44:18.020Z
@FrancoisBilodeau that's exactly the kind of information that you should put in the question to begin with. The more you can share about the type of solution you need, the better that answers can address your specific need. You emphasized fast screen painting in the question and said nothing about accuracy or close approach detection. If that's what you are after, great! But please emphasize that clearly by editing your question post accordingly. Many/most users won't dig down into the comments before writing an answer. *Thanks!*
uhoh
CC BY-SA 4.0
2021-04-07T18:45:40.643Z
In this case you really might want to consider switching to propagating state vectors using a gravity model of several solar system bodies, nor simply orbital elements.
uhoh
CC BY-SA 4.0
2021-04-08T05:03:25.720Z
@slowerthanstopped Hopefully, that's not necessary. From https://en.wikipedia.org/wiki/Jet_Propulsion_Laboratory_Development_Ephemeris#Construction "As of DE421, perturbations from 343 asteroids, representing about 90% of the mass of the main asteroid belt, have been included in the dynamical model.[8]"
PM 2Ring
CC BY-SA 4.0
2021-04-08T05:06:23.840Z
Do you know about the [SPK files](https://ssd.jpl.nasa.gov/?horizons_doc#spk) that Horizons can generate? It's also important to understand the [Statement of Ephemeris Limitations](https://ssd.jpl.nasa.gov/?horizons_doc#limitations)
PM 2Ring
<orbit><asteroids><orbital-elements>
Asteroids conjunction
Tiblemont
# Asteroids conjunction

(Asked by: username_0 on 2021-04-07)

I am currently writing a C++ program to show asteroids in 3D, and find close approaches or collisions. I got my orbital elements from JPL https://ssd.jpl.nasa.gov/sbdb_query.cgi. So far so good, with over a million asteroids drawn in 5ms with my old Quadro GPU
But the problem is the variation of the elements over time. For example, CERES:
epoch,a,e,i,om,w,ma 
2459200.5,2.766089105818,.07816842657453,10.58789954719,80.27235841368,73.72488984426,205.5454154582
2459310.5,2.765760313090,.07831877879848,10.58807660401,80.26860808947,73.73699886586,229.1146825391

After only 4 months, the eccentricity and all other elements have changed as seen on the second line.
How can I compute the change in the parameters without knowing their derivatives?
Or where to get those derivatives?
I searched all JPL small bodies site and their telnet access or email, but could not find a way to download orbital elements with derivatives or a way to compute the change.
I am aware that the difference might be just a pixel on the screen but it represents hundreds of thousands of kilometres. What I want to do is study the close approach between the asteroids themselves and eventually near collision. I already did that for artificial satellites using SGP4 propagator using OpenCL running on GPU video card. I can propagate 20,000 satellites (include debris) in a few milliseconds and get results exact to 1 kilometre. Compare with Celestrak Socrate. A prediction of a 1,000-kilometre approach must be possible.
Does anybody know how?

**Tags:** orbit, asteroids, orbital-elements

> **Comments:**
> **username_1**: few quick questions: 1) if you are drawing a million asteroids will anybody know that a few are a half-pixel off? 2) wait  you downloaded orbital elements for *a million* asteroids from JPL? 3) what are the uncertainties on those osculating orbital elements? 4) if you had all their derivatives, are you confident you would know how to use them to propagate correctly? 5) could there be a reason that osculating elements don't come with derivatives?  
> **username_2**: Stupid answer from me: Grab every single object in the solar system at epoch `J2000` and run a huge nbody simulation over it. If you compare that to your given parameters and it is reasonably accurate, keep doing it. But remember that a butterfly flapping its wings can cause a hurricane!  
> **username_3**: Salut François! The equations describing the precession of nodes and the precession of apsides are quite complex. I doubt you’d want to perform them for millions of bodies! As @username_1 pointed out, nobody will notice that your asteroids are “a half-pixel off” (if that much!). Unless you want to go into the distant past or distant future, I see no use for that.  
> **username_1**: If the osculating orbital elements of asteroids change a lot it's mostly because of Jupiter or other big planets. Those calculations will also sometimes include effects of sunlight.  
> **username_1**: general references [Where can I find the positions of the planets, stars, moons, artificial satellites, etc. and visualize them?](https://astronomy.stackexchange.com/q/13488) and [Is there free, accurate and updated data on planetary body positions in the solar system?](https://space.stackexchange.com/q/13956) and for the effect of sunlight see Marsden parameters [here](https://space.stackexchange.com/a/18438) and [here](https://space.stackexchange.com/q/18474).  
> **username_1**: For an example of how to propagate state vectors instead of using orbital elements see [this answer](https://space.stackexchange.com/a/234090 though you would want to do a simpler propagation than this one.  
> **username_1**: @FrancoisBilodeau that's exactly the kind of information that you should put in the question to begin with. The more you can share about the type of solution you need, the better that answers can address your specific need. You emphasized fast screen painting in the question and said nothing about accuracy or close approach detection. If that's what you are after, great! But please emphasize that clearly by editing your question post accordingly. Many/most users won't dig down into the comments before writing an answer. *Thanks!*  
> **username_1**: In this case you really might want to consider switching to propagating state vectors using a gravity model of several solar system bodies, nor simply orbital elements.  
> **username_4**: @username_2 Hopefully, that's not necessary. From https://en.wikipedia.org/wiki/Jet_Propulsion_Laboratory_Development_Ephemeris#Construction "As of DE421, perturbations from 343 asteroids, representing about 90% of the mass of the main asteroid belt, have been included in the dynamical model.[8]"  
> **username_4**: Do you know about the [SPK files](https://ssd.jpl.nasa.gov/?horizons_doc#spk) that Horizons can generate? It's also important to understand the [Statement of Ephemeris Limitations](https://ssd.jpl.nasa.gov/?horizons_doc#limitations)  

---

## 3 answers

### Answer by username_0 on 2021-04-08. Score: 3

Thanks folks for helping.
I knew about SPK files but from Horizon's telnet inteface there is a limit of 200 bodies per request.
Making over 5000 request might be possible but not very productive.
I also looked at the  DE421 file from ftp://ssd.jpl.nasa.gov/pub/eph/planets/bsp/ but it looks like the corrections for the 8+1 planets
here is the output from BRIEF, a SPICE utility: https://naif.jpl.nasa.gov/naif/utilities.html
BRIEF -- Version 4.0.0, September 8, 2010 -- Toolkit Version N0066
Summary for: de421.bsp
Bodies: MERCURY BARYCENTER (1)  SATURN BARYCENTER (6)   MERCURY (199)
    VENUS BARYCENTER (2)    URANUS BARYCENTER (7)   VENUS (299)
    EARTH BARYCENTER (3)    NEPTUNE BARYCENTER (8)  MOON (301)
    MARS BARYCENTER (4)     PLUTO BARYCENTER (9)    EARTH (399)
    JUPITER BARYCENTER (5)  SUN (10)                MARS (499)
    Start of Interval (ET)              End of Interval (ET)
    -----------------------------       -----------------------------
    1899 JUL 29 00:00:00.000            2053 OCT 09 00:00:00.000

I can always do a statistical estimate, any close approach being false but the overall count being close to reality.
Too bad to drop the project, I was getting nices pictures from it. The coloring is done against eccentricity, red=0

> **Comments:**
> **username_4**: **1.** It would be a shame to drop this project. Maybe contact Horizons, explain what you want to do, and see if they have any suggestions. They may be able to supply the SPK files in bulk (depending on what time span you want to cover), or have some other practical suggestions. **2.** Since you're investigating near encounters, I guess even the SPK files aren't adequate when two asteroids are close enough to have significant gravitational affect on each other, unless one or both are in the set of 343 asteroids included in the dynamical model.  
> **username_0**: well, thanks PM, you can give it a ride at http://geomaitre.com/geomaitre/  You must have a graphic card that support opengl 4.5 most modern one do  Good luck  
> **username_1**: This image is quite compelling! How are the colors assigned? Have you thought of doing some animations? I see two red groups of Trojan asteroids at Sun-Jupiter L1 and L2 but are there any at L3? What made the ones at L1 and L2 red? Scott Manley has several older videos in YouTube from a previous lifetime, cf. [Asteroids In Resonance With Jupiter](https://www.youtube.com/watch?v=yt1qPCiOq-8) and also [Why did the distribution of asteroids discovered in 2010 have a radial modulation?](https://astronomy.stackexchange.com/q/25736/7982)...  
> **username_1**: and answers to [Orbits of Trojan Asteroids](https://astronomy.stackexchange.com/q/43322/7982) and [Code/Package for multiple planet simulation](https://astronomy.stackexchange.com/q/35056/7982) I think you are at the beginning of something really interesting and useful and can take this in several directions. Also see answers to [What is a thermal distribution of eccentricities?](https://astronomy.stackexchange.com/q/34535/7982) Maybe you can "measure the temperature" of of the asteroid belt, or plot average "temperature" as a function of radius, or see of Trojans are "hot" or "cold"?  
> **username_1**: See also upcoming answers (hopefully) to [Looking for that cool astronomical and spaceflight simulation website that's been discussed here](https://space.meta.stackexchange.com/q/2839/12102) and despite the rather blatant narcissism in the beginning this video tells quite an interesting story: [Why I Have a Real Planet Named After Me - (33434) Scottmanley](https://www.youtube.com/watch?v=oIQ2byD6b6o)  
> **username_0**: The coloring is done against eccentricity where 0 is red using HSL color  [0 1] for the hue part. The trojan are close to Jupiter and those away from it are Hilda . see https://en.wikipedia.org/wiki/Hilda_asteroid  

---

### Answer by username_0 on 2021-08-04. Score: 3

I finally computed positions of over 1 million asteroids for a few months from June to October 2021.
There is significantly less close approaches than I first thought, but on the average, there is a conjunction under 10 thousand kilometers almost every other day.
For example on 2021-08-10 20:05 asteroids (2002 GD77) and (2008 TD184) will be at 825 km according to Nasa JPL Small Planets Center and at 290 millions km from earth.
The figure below show  their orbits and positions relative to the sun(20x) and some major planets(1000x)

There will be a closer approach at 745 km on 2021-08-16 17:15 between (2003 WH156) and (2014 WR26) but at 520 millions km from earth. Both conjunction have an approach speed of 5.5 km/sec.
My results are consistentt with JPL website and should be close to reality, but as Jon Giorgini (senior analyst at JPL, thanks again for your help, Jon) remarks, due to the uncertainty in the keplerian elements you cannot tell where exactly the asteroids where at a certain time.
Now these results bring more question than answers.

Is it possible for high end amateur instruments to detect/photograph such a conjunction?
Is it possible to compute in advance the minimum-maximum variation in asteroid seraration?
Is it possible with different equipment/sofware/data to tell if that event occured or not?
How close should an approach be for the orbits being modified significantly?
If the orbits are being modified significantly, how long would it take to be detected from earth?
There is also a whole bunch of "couples" that have VERY similar orbits and travel as old companions at very close distances. Are they duplicates or realy different but bound together by some resonance?

Example:
2435 km     2021-08-10 14:27:08 (1999 FM92)  & (2010 DH114)
824 km      2021-08-10 20:05:00 (2002 GD77)  & (2008 TD184)
8508 km     2021-08-11 06:44:22 (2006 CK14)  & (2015 CY70)
2326 km     2021-08-11 19:40:46 (2017 RD134) & (2019 JH75)
852 km      2021-08-12 05:58:43 (2010 BO137) & (2015 FB346)

> **Comments:**
> **username_1**: This is such a cool project and tool that you've created I think there are a lot of interesting things you can address. The five new questions you've enumerated are really launching points for investigations as well as new Stack Exchange questions. For item #4 there's no real cut-off for significant vs. non-significant for a $1/r^2$ force, but I am sure there's been a lot of work on *perturbation of orbital elements* vs distance and nature of close approach and that could be asked separately, and of course the masses of some asteroids are known by their perturbations on other asteroids!  
> **username_1**: And of course that starts to address #5 as well. You may find the following questions in Space SE interesting as well:  [Algorithmic methods or techniques to find conjunctions in high N Keplerian element ensembles?](https://space.stackexchange.com/q/29946/12102) and to a lesser extent [Algorithmic methods or techniques to find conjunctions in large ensembles of state vectors?](https://space.stackexchange.com/q/33161/12102). In artificial satellite orbital mechanics "conjunction" is a understated term for satellite near-misses or collisions.  
> **username_1**: *Il est fort amusant* indeed! For YouTube links you can just leave the bare url in the post and SE will automatically embed the video once you accept the change, no need to manually create a link.  

---

### Answer by username_0 on 2021-08-25. Score: 2

For those interested in observing conjunctions, here is a partial list of events with approaches under 10,000 km for the coming weeks.
km     = distance between asteroids 
km e6  = distance from earth in  million km 
km/sec = relative speed of approach 
elong  = elongation of asteroids from sun

        |    UTC  time     |asteroid A | km  | asteroid B |km e6| km/sec |elong|
        |------------------|------:----|----:|-----:------|----:|-------:|----:|
        | 2021-08-20 00:43 | 2006 CK14 |*8668 | 2015 CY70 | 480 | 0.0003 | 109 |
        | 2021-08-26 23:48 | 2000 SE89 | 9566 | 2017 XG51 | 432 | 3.4653 | 161 |
        | 2021-09-04 10:31 |    Dyer   | 3660 | 2014 WG584| 420 | 2.5170 | 143 |
        | 2021-09-09 00:09 | 2005 CO9  | 3305 | 2009 KZ1  | 387 | 6.3235 |  94 |
        | 2021-09-18 21:46 | 2008 FW83 | 6355 | 2014 EQ117| 417 | 3.7596 | 138 |
        | 2021-09-28 01:42 | 2007 TZ425| 3475 | 2008 XN66 | 363 | 3.0231 |  97 |
        | 2021-09-28 19:21 | 2001 DL93 | 2670 | 2006 AR3  | 291 | 17.150 |  13 |
        | 2021-09-29 07:02 |   Utkin   | 3575 | 2004 SN2  | 323 | 2.2074 | 154 |
        | 2021-10-07 07:57 | 2013 UH30 | 2396 | 2014 OW110| 309 | 2.5056 | 101 |
        | 2021-10-10 01:58 | 1991 LA1  | 4809 | 2011 LR7  | 424 | 3.5947 | 159 |
        | 2021-10-10 15:18 | 2002 JE35 | 3091 | 2013 PS4  | 379 | 3.7531 | 102 |
        | 2021-10-14 02:17 | 2006 BN154| 8145 | 2009 FB34 | 380 | 4.4304 |  93 |
        | 2021-10-16 16:35 | 2002 GN160| 2826 | 2010 LD118| 462 | 7.8022 | 137 | 

The first entry shows 2 asteroids with very similar orbits, traveling together as old chaps. You can tell by their closing speed near zero. Try your astronomy software, tracking them for a few months.
If you don't have any software, you can get ephemeris from JPL by typing a command like that in your web browser (no space , no line feed)
https://ssd.jpl.nasa.gov/api/horizons.api?format=text&COMMAND=%27DES=2000%20SE89%27&MAKE_EPHEM=%27YES%27&CENTER=%27500@10%27&TABLE_TYPE=%27VEC%27&VEC_TABLE=%271%27&START_TIME=%272021-08-26%2023:45%27&STOP_TIME=%272021-08-26%2023:50&STEP_SIZE=%271%20m%27&CSV_FORMAT=%27YES%27
Note the high approach speed of 2021-09-28 19:21 |  2001 DL93 | 2670 | 2006 AR3  | 291 | 17.1502|  46 |


---

