Author Archives: Garry R. Osgood

About Garry R. Osgood

Booking travel arrangements for Wordpresses on the go our specialty!

The Art of Gerrymandering – Part III

The State of Gerrymandering in these United States

1. The State of Gerrymandering in these United States

In our last post, we wrote about how to compute the Gerrymander Index of shapes, including Congressional districts. Since then we’ve fetched the U. S. Census Bureau tl_2014_us_cd114 Esri shapefile data set of the 435 Congressional Districts for the current 114th Congress, which includes the nine non-voting districts that send delegates to Congress. If you are terminally curious, download the comma-separated value text file of our results, based on Census Bureau dataset. We’re not going to discuss all 444 maps; restricting our attention to the best and the worst, the state of gerrymandering in these United States, and how the States of California and Texas are gerrymandering these days.

Of the four hundred forty four records in tl_2014_us_cd114, the most gerrymandered district is North Carolina Congressional District 12, with a Gerrymander Index of 0.0291, this based on its cartographic boundary as defined in the Census Bureau shape file.

Continue reading

The Art of Gerrymandering – Part II

gerry_00

1. Four units of rope make a square, one unit on a side, and encompasses an area of 1 square unit.

Suppose you had a piece of rope. Your aim is to encompass as much area as possible. The rope doesn’t stretch or shrink, nor break no matter how hard it’s pulled – it’s good rope. Could be made out of cytoplastic nanotubes or something.

No matter if you measure your rope in inches, miles or centimeters. Just invent a fiat measure defined as one quarter of the length of the rope, then stake out a square. Each side is one unit and the square’s perimeter is four units. The square’s area is the square of a side: 1 unit × 1 unit = 1 unit2.

That’s all well and good, but are we encompassing the most area that we can with this rope? A regular octagon suggests otherwise:

gerrymander one

2. Four units of rope also make a regular octagon, 1/2 unit to a side. In this case, the same amount of rope encloses 1.207 square units, a bit more than 1/5th additional area beyond that of a square.

For the same amount of rope, a regular octagon encompasses a little more than 1.20711 times the area encompassed by a square.

We arrived at this regular octagon by applying a modification rule to the square: we halved the sides, doubled their number and made all the interior angles the same. That gave us a polygon that encompassed more area than the square for the same perimeter, four units.

Gerrymander Three

3. Four units of rope also make a regular hexadecagon, 1/4 unit to a side, with sixteen sides equalling four units of rope. In this case, the same amount of rope encloses 1.256 square units, a wee bit better than that of an octagon.

Halving the length of the sides yet again, increasing their number by two, gives rise to a hexadecagon which encompasses an even greater area: 1.256 square units surrounded by four units of rope, divided into 16 sides of one quarter unit each.

You can see where this goes – as the number of sides increase, the regular polygon more nearly approximates a circle. So, by the miracle of calculus, we magically increase the number of sides up through a regular myriagon (10,000 sides) past the regular apeirogon (an countable infinity of sides) arriving at the limit figure: the circle, which encompasses the largest possible area for a given perimeter: about 1.27324 square units for a rope four units long.

gerrymander four

4. A circle with a radius of 2 over pi has a perimeter of four units and an area of 4 over pi, about 1.27324. Compare this to the unit square with a perimeter of the same length, but which only has an area of one unit square.

This is the essence of the Gerrymander Index (GI): any perimeter P of a given shape S that encompasses a particular area Ap, also encompasses a circle with the maximum possible area, Ac. The Gerrymander Index for that particular shape S with perimeter P is then:

[pmath size=12]GI_P = A_p/A_c[/pmath]

The Gerrymander Index is a computable property of a particular shape. It is a unitless measure arising from a ratio of areas and compares a shape’s area with that of the circle whose circumference equals the shape’s perimeter. When the shape is that circle, the Gerrymander Index is one, the ideal. Line segments, which do not have interior areas, have Gerrymander Indices of zero. All shapes that encompass some interior area, but are not circles, have Gerrymander Indices that fall somewhere between one (circle-like) and zero (line-like).

The Gerrymander Index is independent of size. The index compares the area of a shape relative to that of the circle with an equal perimeter, hence they “scale together.” We may compare Gerrymander Indices of huge, rural Congressional districts with block-sized urban districts without getting into apples-versus-oranges side debates on whether size matters.

To get a feel for this index, consider our unit square. It has a Gerrymander Index of 1.00000 (Ap) divided by 1.27324 (Ac): GIsquare = 0.785398. Most people would think of a square Congressional district as not being especially gerrymandered, so a GI of 0.7 – 0.8 can be considered “quite decent.”

The regular octagon with a perimeter of four units and area of 1.20711 square units has a Gerrymander index of 1.20711 divided by 1.27324, or 0.94806 – very nearly circular and probably too good to be seen much in the real world. The regular hexadecagon weighs in with a Gerrymander Index of 1.25684 ÷ 1.27324 = 0.98712 – a hair shy of a circle and too good to be true.

Gerrymander Five

5. A tiny bridge connects two shapes with original GerryMander indices of about 0.98, but the index for the combined shape plummets to 0.507.

Conversely, we can combine a couple of shapes, each with a pretty good Gerrymander Index, into one that doesn’t have a particularly great index. The two circular regions in Illustration 5 on their own have Gerrymander Indices of 0.98. The tiny connection bridge linking them gives rise to an overall “dumb-bell” shape with a Gerrymander Index of just 0.507.

How does the Gerrymander Index fit in with the Great National Discourse, at least insofar as Congressional Districts are concerned?

Justin Levitt, on the faculty of Loyola Law School in Los Angeles, furnishes us with a summary table of the criteria that various state level redistricting organizations follow. Thirty seven states include some sort of “compactness” guideline, but as Prof. Levitt points out, the precise meaning of “compactness” is often wanting, with definitions using language that the shape should be “regular” or that voters in a district should “live close together” or not be “far-flung.” This kind of language wanders around the concept of a numerical index, but doesn’t hit it on the head, leaving quite a bit of interpretive play. Different observers of a particular Congressional district may form different impressions of how “compact” that district is.

In contrast, the GI injects a hard number, one based solely on the geometry of a Congressional District. It enables us to discuss how much a particular district is like a circle. One is not obliged to consider the various political forces that caused a district to be shaped in a particular way. One only needs to apply a technique – taking a ratio of two areas, one that the perimeter encompasses, the other that a circle with an equivalent perimeter encompasses.

A hard number such as the Gerrymander Index allows us to consider particular thresholds. We might argue, for example, that any proposed district with a Gerrymander Index below 0.03 be disallowed as “too contorted”. Of course, that threshold is subject to debate and should be debated. We just wish to point out that at this juncture, the Gerrymander Index makes such a debate possible, as it is a concrete property of the shape.

Alternatively, we can set a threshold on the downward change in the Gerrymander Index from one re-districting effort to the next. Illustration 5 makes the point visually. Two districts with quite excellent GI’s of 0.98 are combined to produce a new district with a GI of 0.507, a downward plunge of 0.473. What if downward changes in GI were limited to a threshold of 0.2, while upward changes in any measure would be allowed? Such a policy would grandfather badly drawn districts initially, but over time, with significant GI drops disallowed, Congressional districts would all tend to compactness, with higher GI indices becoming the norm.

A word of caution is in order. The Gerrymander Index stems from the length of a shape’s perimeter. That comes from a map. To what precision is a map measured? The astute will now recognize that we teeter on the edge of the coastline paradox, attributed by Benoit Mandelbrot to mathematician Lewis Fry Richardson.

Two maps of the State of New York, each published by the U. S. government, illustrate the paradox.

The State of New York as published by the Census Bureau (tan) in mapset "tl_2014_state" Overlaid in turquoise, the State of New York as published by the National Weather Service in "mapset s_16de14". The inclusion of shoreline data causes the GI to plunge from 0.291 to 0.055.

6. The State of New York as published by the Census Bureau (tan) in mapset “tl_2014_state” Overlaid in turquoise, the State of New York as published by the National Weather Service in “mapset s_16de14”. The inclusion of shoreline data causes the GI to plunge from 0.291 to 0.055.

The Census Bureau clearly documents that their maps are for display and illustration. They deliberately simplify coastlines along large bodies of water, though political borders are carefully drawn.

For our purposes, this illustration reminds us that we cannot talk or write about Gerrymander Indices in isolation. The index is absolutely keyed to the map from which it is calculated, and in honest debate, the source of maps must always be mentioned.

As remarked in our second technical note, we use the Census Bureau map sets because those are the ones from which Congressional Districts have been published. As it so happens, the pruning of complex coastlines usually put Congressional districts in a more favorable light. For example, the 1st Congressional District of New York, currently occupying the eastern third of Long Island, would be “naturally gerrymandered” by the North and South Forks, the tiny islands between the two, and the barrier islands running along the south shore. The simplifications applied in the Census Bureau maps omits those details from the map.

So long as we are clear that we ground our Gerrymander Index on this particular map set, and agree that Census Bureau modifications serve technical purposes only, there should be no cause for “apples v. oranges” debates. Though in the polarized atmospheres that encompass much 21st Century political discourse, such an agreement could be hard to obtain in practice.

Next and Last Part: Some of our favorite Gerrymanders.

Technical Note 1: Isoperimetric Index

The Gerrymander Index is not entirely original with this author. Dr. James Case presented a similar formulation in the SIAM Journal in 2007, and his sources trace the technique back to ancient Greece, so even Pythagoreans had some notion of a gerrymandering index.

Case reports on a unitless measure arising from a ratio of areas, but for the numerator he takes the area formed by the length of a shape’s perimeter, Pand compares this with the area 4πAp, the denominator, with Ap equal to the area of the shape encompassed by P. If P happens to be circular, then the value P2 will equal 4πAp, the numerator and denominator have the same value and the ratio of the two areas becomes one.  This arises from the relation that couples a circle’s circumference (perimeter) to its area: A =πr2; P = 2πr; P2 = 4π2r2; P2 = 4π(πr2); P2 = 4πA. Thus, the ideal in Dr. Case’s “isoperimetric index” is identical to the Gerrymander Index: unity.

The usual non-circular case may be reached by holding Ap constant and pulling, pushing and twisting the perimeter P out of round so that it grows in length, encompassing the fixed area Ap less and less efficiently. It becomes more contorted and “longer,”  leaving P2 > 4πA, the isoperimetric index exceeding one.

The isoperimetric index behaves somewhat like the inverse of the Gerrymander Index, reporting divergence from the circular ideal with ever-larger numbers. This is a technical difference. Conceptually, it too is a hard number and enters into the Great Discourse the same way that the Gerrymander Index does: injecting numeric precision into a debate that suffers from fuzziness.

Technical Note 2: The Area of Arbitrary-Shaped Closed Polygons

Few shapes in the real world, Congressional Districts included, come with neat formulae that give exact areas; the world is fractal. So how does one deal with the realities of Illinois Congressional District #4?

Computational geometry gives us one method that does not constrain us too much, so long as we limit ourselves to closed polygons with sides only consisting of line segments and which do not self-intersect. That is, our polygon lays “flat” on a surface without any part of it folding over any other part. When the shape does not self-intersect, it may be stretched topologically into a circle. If reshaping entails one or more holes, like a doughnut, then the shape has intersected itself. Barring that, and with only line segments for sides, a shape may be otherwise arbitrarily convoluted.

This is the kernel operation. It produces an area fragment, Aj: , corresponding to adjacent vertices (xj, yj) and (xj+1, yj+1):

[pmath]({x_j}{y_{j+1}} – {y_j}{x_{j+1}})/2[/pmath]

The kernel operation works on pairs of adjacent vertices, j and j+1. We start with vertex zero and one, apply the kernel operation, go on to vertex one and two, apply the kernel operation, and so on, until we come to the final pair, which is the very last vertex paired with the very first. We add up the area fragments and take the absolute value of the sum. This gives us the area of the arbitrary polygon. The absolute value operation disguises the fact that walking counterclockwise around a polygon calculates a negative area. This may unsettle the casual reader. Negative areas have their uses but we’ll put such aside and just take the absolute value of the numerator.

Gerrymander Six

6. The area of an arbitrary polygon arises from a sum of “area fragments”, each composed from pairs of adjacent vertices.

Usually, polygonal datasets for congressional districts consist of thousands to tens of thousands of vertices, making this effort a bit tedious for paper-and-pencil work. That’s what computers are for.

Where do we get our data? The U. S. Census Bureau furnishes data on the shapes of 2013 Congressional districts in the form of Esri Shapefiles, which meet these criteria and are available to the public from the Census Bureau product page. With the exception of Minnesota, these shapefiles also describe the districts for the current 114th Congress. These files come in different resolutions to serve varying display purposes.  TIGER® (Topologically Integrated Geographic Encoding and Referencing)/Line shapefiles are for high resolution work; they can get quite large. Cartographic Boundary shapefiles are light-weight, low-resolution versions of the Line files: quick to down-load, easy to render, but furnish only somewhat coarse approximations of a political boundary. We use the TIGER/Line files for Gerrymander Index computations.

What do we do with our data? Get a shapefile reader, which come in a variety of shapes and sizes, and do a bit of scripting for the Gerrymander math. For this series, we use the Python scripting language to do our math and Joel Howland’s pyshp module to interface with Esri Shapefiles. This is a lightweight approach for those accustomed to scripting. It’s how we made our pictures, too.

For those who just want to load and visualize, one needs GIS software. These too come in a variety of shapes and sizes, but none, at the moment, give readings on either isoperimetric or Gerrymander indices. High-end jobs like GRASS can take add-ons written in a variety of languages, so one could, in principle, add a Gerrymander Index calculator. GRASS is a world unto itself, however, so we didn’t go that route, wishing to finish this post before the century closed. There are also online communities centered on web-based geodata. Google Maps is the best known, and lets users integrate Shapefiles. Injecting custom calculators into the mix do not seem possible at the moment, but there is always the future.

Further Reading

  1. “James Case “Flagrant Gerrymandering: Help from the Isoperimetric Theorem?”” SIAM News, Volume 4 Number 9, November 2007

  2. H. R. 1347 (114th Congress, 2015-2017) John Tanner Fairness and Independence in Redistricting Act

  3. Justin Levitt All About Redistricting “Where are the lines drawn?”

The Art of Gerrymandering – Part I

gerrymander_00

The Original Gerry Mander

The Constitution tasks the House of Representatives with setting the number of U. S. citizens that its members may represent. The Apportionment Act of 1792 fixed the House of Representatives for the Third Congress at 105 members, one Representative for 33,000 constituents. The Census of 1790, first of its kind, found the young nation numbering around 3,900,000 individuals. For purposes of computing the ratio of representatives to those represented, slaves constituted three-fifths of a free person.

112 years on, 1901, roughly midway between the Constitution’s ratification and the present day, each Representative of the 57th Congress fielded the concerns of 213,000 people and carried a six-fold increase in “representational load” over his 1792 counterpart. The House then had 357 members representing around 76 million. Had the House stayed with its 1792 ratio of one Representative to 33,000 constituents, it would have had 2,303 members in 1901, far more than what the seating in the south wing of the Capitol building could accommodate.

114 years on, the 114th Congress finds a House of 435 voting members, a number which has been fixed since the Apportionment Act of 1911. These worthies now represent about 309 million, or roughly 710,000 citizens per Representative, a four-fold increase over the 1901 representational load and a twenty-four fold increase over that of 1792. At the original ratio, the House would have almost 9,364 members, a number making for a mad house – though some think it is anyway.

Continue reading

Colony Collapse Disorder

Californians have been growing almonds for a long time. "The Almond Grove At  T. H. Selby's Residence, Fair Oaks, Cal." ca. 1870. NY Public Library

Californians have grown almonds for a long time. “Almond Grove At T. H. Selby’s Residence, Fair Oaks, Cal.” ca. 1870. NY Public Library

2006 was something of a banner year for both beekeepers and growers. A high fraction of beekeepers experienced better than 50% losses. Many growers could not rent bees at any price. Lots of almond groves in California went without pollination that year, and, as with many species of nuts, if the pollen isn’t carried over from an unrelated plant, nuts do not develop. Almonds, like many crops, (and most nuts) do not self-pollinate. Unless pollinated by an agent like a honey bee, they just don’t bear fruit. The price of almonds shot up in 2006 and has not dropped very much since.

Continue reading

The case of the Disappearing Bee

Now you see her; now you don't

Now you see her; now you don’t.

There aren’t enough bees to perform all of the pollination work in the United States. So, if you like to travel, make money doing it, and don’t mind a sting now and again, then professional beekeeping may just be for you.  And hey – it’s a seller’s market.

About 1.4 million hives get hauled to and around California every spring. They travel by plane, semi or pickup trucks. 60% of the American hives are engaged in commercial pollination. They pollinate almond trees, and other key crops:

  • Cranberries in Massachusetts and New Jersey.
  • Blueberries and peaches from Jersey to Georgia,
  • Pumpkins from Jersey to Illinois.
  • If it comes from a flowering plant, it’s pollinated by bees.

You can take your bees “On The Road” and make money without writing a word. And if you actually hate to travel, you can just send the bees.

Just make sure they all don’t all disappear.

Continue reading

Robot Rodeos Conjure Up Disasters and Pancake Contests

Extreme Hazard robot essaying an obstacle course at the Robotic Vehicle Range, Kirtland AFB, Alberquerque, NM. Sandia National Labs.

Extreme Hazard robot essaying an obstacle course at the Robotic Vehicle Range, Kirtland AFB, Albuquerque, NM. Sandia National Labs.

A sultry day was in the offing near Purnell OK, the seat of McCurtain County in the state’s southeast quadrant, just a dozen miles northwest from the triple point where Arkansas, Texas and Oklahoma meet. One hundred forty miles northeast, the National Weather Service Doppler radar station KSRX at Ft. Smith Arkansas, was monitoring a cold front approaching from the west, driven by a mass of cool dry air sweeping down from the northern plains. Typical for the late spring in the American prairie, this eastbound mass was colliding with a warm, wet air mass streaming north from the Gulf of Mexico, now roiling under a cool dry tongue at 700 mb. Buoyant but trapped under heavier cool air, supercells were forming in the humid 850 mb surface layer twenty miles west of Purnell.

Continue reading

The Celestial Shooting Gallery, Part Four: “You Have Nothing to Worry About (click) Worry About (click) Worry About (click)…”

Stability Model of an experimental distribution grid

A stability map of a simple power grid. Each point on this image represents an operating state of a simple power grid consisting of a few generators. Bluish regions constitute stable working states, red unstable and ‘salt-and-pepper’ represent chaotic behavior. One can tune a grid for stability by controlling the phasing of generators and transformers on the grid and such settings suffice for day-to-day operations. It is difficult to decide where, or by how much, abnormalities such as geomagnetic storms might push a system into red, unstable regions, or, worse, salt-and-pepper regions where the system oscillates between states. It is easy to find cases on the map where chaotic regions lie very close to stable regions, indicating that the destabilizing push need not be large at all. James Thorp, Cornell University, published in IEEE Spectrum

People paid to worry about the North American power grid regard geomagnetic storms as “high impact, low-frequency” events, spawning the inevitable acronym: HILF. Low frequency, in that a geomagnetic storm as intense as May 1921, at 5,000 nano-Teslas/minute, or the 1859 Carrington Event, best guess: 7,500 nano-Teslas/minute, might not happen in our lifetimes, the lifetimes of our children, or even our grand children. If signature traces in Arctic ice core samples are correct, these are ‘500 year events.’ When it comes to deciding where to put that preventative maintenance dollar, storm-proofing Oklahoma elementary schools against EF 5 tornadoes seems a far more practical spend than the hardening of electrical grids against a half-theoretical event that might not even happen in 500 years.

What pulls planners up short is the high impact part: the utter god-awfulness of a power grid that crashes and which then can’t boot itself up. There is a self-referential dependency: fixing a dysfunctional power grid requires it to be functional, as key aspects of the manufacturing of transformers need electricity.

Nor can one expect the cavalry to ride in anytime soon, as the vast geographic reach of geomagnetic storms means that one strong enough to take down the North American grid may very likely take down Eurasian grids as well – entire hemispheres could wind up in the toilet, and we only have two hemispheres. That and the statistical variableness to it all: the Carrington 1859 and May 1921 storms, nominally two ‘500 year events’ were, in fact, separated by only sixty-two years.

Where does the buck stop? Continue reading

Celestial Shooting Gallery, Part Three: When a CME Hits the Atmosphere

Failed GSU transformer at Salem River, NJ

A Generator Step Up (GSU) transformer failed at the Salem River Nuclear Plant during the March 1989 geomagnetic storm. The unit is depicted on the left; some of the burned 22kV primary windings are shown on the right. Though immersed in cooling oil, the windings became hot enough to melt copper, at about 2000 degrees F. John Kappenman, Metatech

Coronal Mass Ejections are mainly charged particles, protons and electrons. When a CME arrives at Earth, the charged protons and electrons come under the influence of the Earth’s own magnetic field, the magnetosphere. Charged particles spin around the lines of magnetic force that comprise the magnetosphere, which diverts most of CME harmlessly around the planet, keeping Earth’s surface tranquil.

If the ejection is large enough, however, it can distort the shape of the magnetosphere, occasionally causing magnetic flux lines to snap and reconnect. When this happens, charged particles leak in and follow the magnetosphere’s flux lines down to the Earth’s ionosphere. There, they strike oxygen and nitrogen molecules and strip them of electrons. These ionized gases glow, giving rise to the ethereal beauty of the auroras around the north and south poles. Unfortunately, these excess charged particles also produce immense electrojets.

Continue reading

Celestial Shooting Gallery, Part Two: The Physics of Geomagnetic Storms

goddard_cme_earth

On August 31, 2012 a long filament of solar material that had been hovering in the sun’s atmosphere, the corona, erupted out into space at 4:36 p.m. EDT. The coronal mass ejection, or CME, traveled at over 900 miles per second. The CME did not travel directly toward Earth, but did connect with Earth’s magnetic environment, or magnetosphere, causing aurora to appear on the night of Monday, September 3. The image above includes an image of Earth to show the size of the CME compared to the size of Earth. NASA Goddard Spaceflight Center

Thursday, May 2nd, 2013, a coronal mass ejection (CME) hurled nearly one billion tons of charged particles from the sun’s corona at an outward velocity of one million miles per hour – 270 miles per second.

In less than a half hour, 2,700 virtual Empire State Buildings, 340,000 tons apiece – give or take a few gorillas – erupted from an active region of the Sun’s surface called AR1748, a northern latitude sunspot. AR1748 had just become visible on the western limb of the Sun’s surface when it ejected this mass, so the vast bulk of it hurled outward, not toward us in Libra, but more or less toward Cancer, at right-angles to us. In practical terms, it shot wide of its mark. Still an impressive shot. The CME had been triggered by an M class solar flare, the second largest in a five step scheme (An, Bn, Cn, Mn, Xn; for n a relative magnitude). It had been the largest coronal mass ejection observed thus far in 2013.

And it was still early in the day for AR1748.

Continue reading

Celestial Shooting Gallery, Part One: The Day We Lost Quebec

Electrojets over N. America

John Kappenman reconstructed the electrojets which formed in the ionosphere late in the March 13, 1989 geomagnetic storm which compromised the Hydro-Quebec power grid in Canada. Concurrently, the eastward jet induced ground currents that severely strained the electrical distribution grid of northern continental United States, resulting in a transformer failure at the Salem Nuclear Power Plant, in New Jersey. Courtesy of Metatech

Nearly a quarter century ago, on March 13, 1989,  a geomagnetic storm led to the collapse of the Hydro-Quebec electrical grid system, which furnishes power to much of the province of Quebec, Canada. So pervasive were abnormal currents, that protective circuit breakers tripped throughout the system, bringing the entire grid to a halt in about one and a half minutes. The grid’s self-protective systems were geared toward local abnormalities happening in particular places. In contrast, ground induced currents created abnormalities everywhere. The good news was that most of the hardware protected itself. The bad news was that six million customers were without power for as long as nine hours, and where transformer damage did occur, outages continued for another week.

Further south, the United States experienced a close shave. A second surge in the March 13 storm generated similar ground induced currents in the northern United States, with large current spikes observed from the Pacific Northwest to the mid-Atlantic states, one spike destroying a large GSU transformer at the Salem Nuclear Power Plant in New Jersey. According to John Kappenman, of the Metatech Corporation “It was probably at this time that we came uncomfortably close to triggering a blackout that could have literally extended clear across the country.”

Continue reading